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

《2022年9《Matlab及应用》上机实践报告.docx》由会员分享，可在线阅读，更多相关《2022年9《Matlab及应用》上机实践报告.docx（18页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精品学习资源Matlab 及应用上机实践报告系 别：专业班级：学 号：姓 名：一、简易运算题： 第 1-4 小题每道题 5 分，第 5 小题 10 分，共 30 分1、运算，当和时的值；方法一：编写一个函数文件 humps.m ，即：functiony=humpsxy=x.3+x-0.98.2/x+1.35.3-5*x+1./x；然后调用，在命令窗口中输入： y1=humps2 y1 =-4.4723 y2=humps4 y2 =42.8096 y=fevalhumps,2,4 y =-4.442342.8077方法二： f=inliney=x.3+x-0.98.2/x+1.35.3-5*x+

2、1./x ； f2ans =-4.4723 f4ans =42.8096 y=fevalf,2,4 y =-4.442342.80772、运算； cos60*pi/180-9-sqrt21/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%求特点多项式的

3、系数p1 =1.0000 -22.0000 -37.0000 122.0000 rootsp1%求特点根ans =23.3603-3.06451.7042欢迎下载精品学习资源所以，原对称实矩阵的特点根为5、设 为符号变量， 试进行如下运算1）；2）； 3）；4）； 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% fx-gx ans =x4-3*x2-7-x3-5*x% expandfx*gx ans =x7+4*x6+6*x5+12*x4+6*x3+12*x2+5*x+8% expandfx/gxa

4、ns = 1/x3+4*x2+5*x+8*x4+1/x3+4*x2+5*x+8*x2+1/x3+4*x2+5*x+8% factorfx ans =x2+x+1*x2-x+1%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,3p1 =%3阶公式拟合后的多项式系数-0.06670.3700-0.78830.7350 x1=linspace0.8

5、,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 x=446

6、,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=

7、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

8、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=sy

9、msin2*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=quadx./x.4+4,1,10 y =0.2718方法二：第一建立函数文件 myfun.m，即functiony=myfunx y=x./x.4+4；然后调用，在命令窗口中输入： y=quadmyfun,1,10 y =0.2718方法三： f=inlinex./x.4+4； y=quadf,1,10

10、 y =0.2718方法四： f=xx./x.4+4； y=quadf,1,10 y =0.27187、运算二重积分；方法一： S=dblquadx-y.2+sinx+y.2,pi,2*pi,0,pi S =118.5787方法二：S=dblquadinlinex-欢迎下载精品学习资源y.2+sinx+y.2,pi,2*pi,0,pi S =118.5787方法三：第一建立函数文件 jifen2.m，即functionz=jifen2x,yz=x-y.2+sinx+y.2；然后调用，在命令窗口中输入： S=dblquadjifen2,pi,2*pi,0,pi S =118.5787四、编程画图

11、题： 每道题 6 分，共 12 分1、设， 把区间分为 125 点，画出以为横坐标，为纵坐标的曲线； x=linspace0,2*pi,125； y=cosx.*0.5+3*sinx./1+x.2；欢迎下载精品学习资源 ploty； titlemy first plot； xlabelx； ylabely； 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

12、.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, 画出三维

13、图形及其在三个平面上的投影；（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=ode45rossler,0, 500,12,2,9 ；subplot2,2,1 ploty:,1,y:,2 xlabelindependent Variable x ylabeldependent Variable y subplot2,2,2ploty: ,1,y:,3xlabelindependent Variable x y

14、labeldependent Variable z subplot2,2,3ploty: ,2,y:,3欢迎下载精品学习资源xlabelindependent Variable y ylabeldependent Variable z subplot2,2,4 plot3y:,1,y:,2,y:,3xlabelVariable xylabelVariable y zlabelVariable zfunction 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；2a=0.3欢迎下载精品学习资源欢迎下载