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可编辑资料 - - - 欢迎下载精品名师归纳总结可编辑资料 - - - 欢迎下载
