《MATLAB语言及应用上机实习报告.docx》由会员分享,可在线阅读,更多相关《MATLAB语言及应用上机实习报告.docx(50页珍藏版)》请在优知文库上搜索。
1、MAT1.AB软件与基础数学试验指导老师:报告人:学号:班级:报全部代码均关纳MAT1.AB2016a编写.和之前版本不同.它首次攵持实时编辑,以下内容在可能在吩喟行窗口、场蜉器或实时防.饵热中埸写.试验1MAT1.AB基本特性与基本运算例1.1求12+2(7-4)13:的分术运算结果.(12+2*(7-4)V32ans=2例-2计算5!并把运算结果赋给变Iityy=fac(orial(5)Y=120例1-2计算2开平方2-(1/2)ans=1.4142例1-4计算2开平方并肮值给变量区(不显示)会看X的肮依状况x=2(12):xX=1.4142Sin(Ial+|)例1-5设=-24”,方=7
2、5计籁Jlan(I+/”)的值,a=-24180zb=7518():sin(abs(八)+abs(b)(tan(abs(a+b)ans=1.7951例1-6设三角形三边长为“=4,8=3,c=2.求此三角形的面枳.a=4:b=3x=2;p=(a+b+c)/2:s=(p*(p-a)”(pb).(p)(l2)2.9047=I23456B=-I2()-113例1-7设I0192II.计算八+8.八Al*.A.A=13146;1.OJl;B=|-1.2.O:1J3:2JJ);+Bans=043569312A*Bans=779131921norm(八)ans=9.5758A(-I)ans=0.83330
3、.33330.5000-0.33330.3333-I.(XX)O0.8333-0.33330.5000例1-8显示上例中矩阵A的第2行第3列元素,并对其进行修改.A(2,3)ans=6A(2,3)=I:A(2,3)ans=_sinx例IT分别而出函数)=x%sx和“-X在区间卜6,6)上的图形.x=(-6)*pi:e.ei:6*pi;y(x.A2).cos();zsin(x).;plot(,y);figure,plot(x,z);functionmianfunction()disp(nm2str(yl(1,2,3”)dis(num2str(y2)functiony=yl(arbzny-(a+b
4、)An;endfunctiony三y2(arbzny=(a-b)An;endend运行显示:11uinfuchon27-1f(x)=+1-6例1-13设(-v-03)+001(.v-0,9)+0.04,试亘出的0.2)上的曲线段。%加坐标网格=llnsace(,2);y=l./(x-0.3).A2*0.01)+l./(x-0.9).A2+0.e4)-6;plot(,y);gridon;例如:对于例胚HI3中所定义的曲),求其零点c.X=SOlvefl/(.O,3)2O.OIl(.O.92O.O4)-6=O,)X=例如:求一元函数量小值(fminbnd吩咐)由于fminbnd只能求一个最小值,通
5、过视察函数可知微小值点为8和+8x=fninbnd(y.O.inf)X=Infy=(X)1/(x-,);l-solve(s-5,x)1,root(4.z.1)root(o.z2)root(4.z3)JOO求无理数川2的近似值.(提示:e-2.7182818284.,In2=0.6931471806.)x=i:0.(X!:2:y=I:trapz(x.y)ans=0.69311.sin(x+)-sin.rIim例117求极限aah.SymSXh;Iimil(Sin(X+h)-sin(x)Vh.h.0)ans=COS(X)ff2f2f例1.18:设)=xy+SinG).求dx6、.2原办symsXy
6、n;f=sym(xAn*y+sinans-COS()Mdiff(diff(fxx)zxaB-n2y(n-l)diff(difff,y),ya11s-Sinbr)diff(diff(fxx),yans-11*,例“:求信-1.*,匐点心二时:(x+)z)dz.ASymSxyZ%声明符号更受,留近变量间必需用空将分开symsXy2:;intrx)/Ians2inn(xysynsu11Mf.kJ.i11f)ans=+=+1.+二T(ans=3/2、)%求纨敦33?3symska;f=(kW3”:symsum(f,k,O.inf)ans=(3*a)2泰勒绽开(taylor)SymSXfy=l(ltay
7、lor(f.x.ordcr.8)ans=xa7+x6-x4+x3X+i+x+xa2)求fx对自变量x(默认)在X=O点(我认傣勒就开树6项退认)SymSx;f=(x)l/(l+x+xA2);taylor(f.x)ans=-x4+x3-X+I4ifx对白变量x(默认)在x=l点奏勒归开式讲8项taylor(f.x,l,o11fer,.8)ans=(2*(x-1)a2)9-x/3-(X-1/3/9+(x-1/4/27-(X-IP681+(x-I)781+2/3方程求根(solve)fx=sym(a*xA2+b*x+c);%虚立符号葩我方包fx=()的符号解SymSabcx;y=a*x2+b*x+c
8、=0:solvey)ans=-A(IM2*a)(b-(b2-4*a*c)A(l/2y(2*a)求方程fx=O关于变Itb的杼号解solvcy.b)ans=-(a*xA2+c)x微分方程(组)求解ans=x22C2求方程y=I+y,y(0)=l,y(O=O的特Wsymsy(t);dsolve(D2y=I+Dyy=I.Dy(O)=O*)ans=exp()-(x,=xy术方V姐)=2的通解.默认白变量为tsyms,(t)x();z=dslve(Dx=x+y,.Dy=2x);Z.Xans=C6*exp2*l)-(C5i*exp.t)V2z.yans=C54cxXt)+C6*cxp(2*t)试验2MAT
9、1.AB绘制二维、三维图形例2-1在子图形窗口中画出2上正弦、余弦曲线。XIinspacc(0,2pi);y!=in(x);y2=:subplot(2.1.1):plot(x,yl);subplot(2,1,2):plot(x,y2);例2-2画出2加上正弦、余弦曲城并对城鞭加粗、点型加大,重新定冏坐标系以及加注相关说明和注拜.=linsace(,2*pi);yl-sin();y2=cos();subplot(2,l,l);plot(,yl,1.ineWidth,l);title(y-sinx,);latwl(,);ylatcl(y);subplot(2,l,2);plot(,y2,oMarkerSize,j3);axis(-l,8,-2,2);title(,ycosx,);xlatM?lCx);ylabl(y);例2-3分别在两个图形新口画出地充正