2022年9《Matlab及应用》上机实践报告.docx

精品学习资源Matlab 及应用上机实践报告系 别：专业班级：学 号：姓 名：一、简易运算题： 第 1-4 小题每道题 5 分，第 5 小题 10 分，共 30 分>1、运算，当和时的值；方法一：编写一个函数文件 humps.m ，即：functiony=humpsx>y=x.3+x-0.98>.2/x+1.35>.3-5*x+1./x>；然后调用，在命令窗口中输入：>> y1=humps2> y1 =-4.4723>> y2=humps4> y2 =42.8096<或者把 y1 和 y2 放在一起，通过以下方式求解）>> y=feval'humps',2,4> y =-4.442342.8077方法二：>> f=inline'y=x.3+x-0.98>.2/x+1.35>.3-5*x+1./x>'> ；>> f2>ans =-4.4723>> f4>ans =42.8096<或者把 y1 和 y2 放在一起，通过以下方式求解）>> y=fevalf,2,4> y =-4.442342.80772、运算；>> cos60*pi/180>-9-sqrt2>>1/3>欢迎下载精品学习资源ans =-1.46493、设，求，；>> A=1 4 8；-3 6 -5；2 -7 -12；>> B=5 4 -2；6 -2 -8；-1 -9 7；>> C=A*BC =21-76222621-77-20130-32>> D=A.*B D =516-16-18-1240-263-84所以，；4、设对称实矩阵， 求其特点根；>> a=2 4 9；4 2 4；9 4 18；>> p1=polya>%求特点多项式的系数p1 =1.0000 -22.0000 -37.0000 122.0000>> rootsp1>%求特点根ans =23.3603-3.06451.7042欢迎下载精品学习资源所以，原对称实矩阵的特点根为5、设 为符号变量， 试进行如下运算<1）；<2）； <3）；<4）；<5） 对进行因式分解>> syms x fx gx>> fx=x4+x2+1；>> gx=x3+4*x2+5*x+8；>> fx+gx ans =x4+5*x2+9+x3+5*x%<1）的结果>> fx-gx ans =x4-3*x2-7-x3-5*x%<2）的结果>> expandfx*gx> ans =x7+4*x6+6*x5+12*x4+6*x3+12*x2+5*x+8%<3）的结果>> expandfx/gx>ans = 1/x3+4*x2+5*x+8>*x4+1/x3+4*x2+5*x+8>*x2+1/x3+4*x2+5*x+8>%<4）的结果>> factorfx> ans =x2+x+1>*x2-x+1>%<5）的结果二、数据拟合与插值： 每道题 6 分，共 12 分>1、用三阶公式运算在和处的导数值，的值由下表欢迎下载精品学习资源给出：1.01.11.21.31.40.250.22680.20660.18900.1736>> x=1.0,1.1,1.2,1.3,1.4；>> y=0.25,0.2268,0.2066,0.1890,0.1736；>> p1=polyfitx,y,3>p1 =%3阶公式拟合后的多项式系数-0.06670.3700-0.78830.7350>> x1=linspace0.8,1.6,100>； %在0.8， 1.6区间取 100 个点>> y1=polyvalp1,x1>； %运算拟合多项式的值>> plotx,y,'*b',x1,y1,'r'>； %绘制原始数据和拟合曲线>> legend'原始数据 ','3阶多项式 '> ；>> k=polyderp1>%对拟合后的 3 阶多项式进行求导k =-0.20000.7400-0.7883欢迎下载精品学习资源>> s1=polyvalk,1.0>%对求导后的多项式运算在的值s1 =-0.2483>> s2=polyvalk,1.2>%对求导后的多项式运算在的值s2 =-0.18832、在某处测得海洋不同深处的水温如下： 深度44671495014221634<m）水温7.044.283.402.542.13利用分段线性插值函数，求在深度为500m、1000m、1500m 处的水温；>> x=446,714,950,1422,1634；>> y=7.04,4.28,3.40,2.54,2.13；>> y1=interp1x,y,500,'linear'>%运算在 500 处的插值结果y1 =6.4839>> y2=interp1x,y,1000,'linear'>%运算在 1000 处的插值结果y2 =3.3089>> y3=interp1x,y,1500,'linear'>%运算在 1500 处的插值结果y3 =2.3892三、解方程与微积分： 每道题 4 分，共 28 分> 1、求下面线性方程组的解；方法一：欢迎下载精品学习资源>> A=1 1 1 -1；2 1 -3 -14；-3 2 1 -5；7 -4 -3 2；>> B=5；-1 ；3 ；-2；>> X=AB X =1.0000-2.00005.0000-1.0000所以，原方程组的解为方法二：>> A=1 1 1 -1；2 1 -3 -14；-3 2 1 -5；7 -4 -3 2；>> B=5；-1 ；3 ；-2；>> invA> ans =0.33330.0505-0.2121-0.01010.33330.4444-1.6667-0.88890.3333-0.43431.42420.68690.00000.0606-0.4545-0.2121>> X=invA>*B X =1.0000-2.00005.0000-1.0000所以，原方程组的解为2、 求代数方程的根；>> y=3 4 7 2 9 12；>> rootsy> ans =欢迎下载精品学习资源-0.8612 + 1.4377i-0.8612 - 1.4377i0.6737 + 1.0159i0.6737 - 1.0159i-0.9583所以，原代数方程的 5 个根分别为：3、 设方程的根为， 求它们对应的关于的线性代数方程；>> r=-3 -5 -8 -9；>> y=polyr> y =125223831108所以它们对应的关于的线性代数方程为：；4、已知，求；方法一：>> syms x y>> y=sin2*x>*logx> y =sin2*x>*logx>>> diffy> ans =2*cos2*x>*logx>+sin2*x>/x方法二：>> y=sym'sin2*x>*logx>'>；>> diffy> ans =2*cos2*x>*logx>+sin2*x>/x所以5、求>> syms x fx欢迎下载精品学习资源>> fx=tanx>-sinx>>/1-cos2*x>> fx =tanx>-sinx>>/1-cos2*x>>>> f=limitfx> f =06、运算一重积分方法一：>> y=quad'x./x.4+4>',1,10> y =0.2718方法二：第一建立函数文件 myfun.m，即functiony=myfunx> y=x./x.4+4>；然后调用，在命令窗口中输入：>> y=quad'myfun',1,10> y =0.2718方法三：>> f=inline'x./x.4+4>'>；>> y=quadf,1,10> y =0.2718方法四：>> f=x>x./x.4+4>；>> y=quadf,1,10> y =0.27187、运算二重积分；方法一：>> S=dblquad'x-y>.2+sinx+y>>.2',pi,2*pi,0,pi> S =118.5787方法二：>>S=dblquadinline'x-欢迎下载精品学习资源y>.2+sinx+y>>.2'>,pi,2*pi,0,pi> S =118.5787方法三：第一建立函数文件 jifen2.m，即functionz=jifen2x,y>z=x-y>.2+sinx+y>>.2；然后调用，在命令窗口中输入：>> S=dblquad'jifen2',pi,2*pi,0,pi> S =118.5787四、编程画图题： 每道题 6 分，共 12 分>1、设， 把区间分为 125 点，画出以为横坐标，为纵坐标的曲线；>> x=linspace0,2*pi,125>；>> y=cosx>.*0.5+3*sinx>./1+x.2>>；欢迎下载精品学习资源>> ploty>；>> title'my first plot'>；>> xlabel'x'>；>> ylabel'y'>；>> grid on；2、绘制函数在，上的曲面图；>> syms x y>> z=1/2*pi>*exp-x2+y2>+2>；>> ezsurfx,y,z,-3,3,-3,3>；欢迎下载精品学习资源五、综合编程题： 每道题 9 分，共 18 分>1、编写一个主程序文件，文件名称为sum，求的值；建立主程序文件 sum.m ，即clcclear s=0 ；sign=-1；for i=1:100sign=-sign；term=sign./i；s=s+term；end out=s；运行以后，输出结果如下：欢迎下载精品学习资源>> out out =0.68822 、编写一个函数文件，实现求解， 要求函数名为 fac， 并在命令窗口中对该函数进行调用，求的值；编写一个函数文件 fac.m，即：functiony=facn> s=1 ；for i=1:ns=s*i；endy=s ；运行以后，输出结果如下：>> y=fac10> y =36288004、 rossler 模型是（1） 取 a=0.32,b=3,c=6, 画出三维图形及其在三个平面上的投影；（2） 取 b=2,c=4 转变 a 将得到如下图形：a=0.3极限环 >； a=0.35周期 2>；a=0.375 周期 4> ； a=0.386四带混沌吸引子 >； a=0.3909周期 6>； a=0.398.> ；欢迎下载精品学习资源t,y=ode45'rossler',0, 500,12,2,9> ；subplot2,2,1> ploty:,1>,y:,2>> xlabel'independent Variable x'> ylabel'dependent Variable y'> subplot2,2,2>ploty: ,1>,y:,3>>xlabel'independent Variable x'> ylabel'dependent Variable z'> subplot2,2,3>ploty: ,2>,y:,3>>欢迎下载精品学习资源xlabel'independent Variable y'> ylabel'dependent Variable z'> subplot2,2,4> plot3y:,1>,y:,2>,y:,3>>xlabel'Variable x'>ylabel'Variable y'> zlabel'Variable z'>function dy=rosslert,y> ；a=0.32；b=3；c=6 ；dy=zeros3,1> ；dy1>=-y2>-y3> ；dy2>=y1>+a*y2> ；dy3>=b-c*y3>+y1>*y3>；2>a=0.3欢迎下载精品学习资源欢迎下载