MATLAB语言及应用上机实习报告.docx
MAT1.AB软件与基础数学试验指导老师:报告人：学号：班级：报全部代码均关纳MAT1.AB2016a编写.和之前版本不同.它首次攵持实时编辑,以下内容在可能在吩喟行窗口、场蜉器或实时防.饵热中埸写.试验1MAT1.AB基本特性与基本运算例1.1求12+2×(7-4)1÷3:的分术运算结果.>>(12+2*(7-4)V32ans=2例-2计算5!并把运算结果赋给变Iity»y=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”,方=75'计籁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=I23'456B=-I2()-113例1-7设I0192II.计算八+8.八Al*.A'.»A=13146;1.OJl;»B=|-1.2.O:1J3:2JJ);>>+Bans=043569312»A*Bans=779131921»norm（八）ans=9.5758»A(-I)ans=0.83330.33330.5000-0.33330.3333-I.(XX)O0.8333-0.33330.5000例1-8显示上例中矩阵A的第2行第3列元素,并对其进行修改.>>A(2,3)ans=6>>A(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<1,2r3)>)functiony=yl(arbzn>y-(a+b)An;endfunctiony三y2(arbzn>y=(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)2÷O.OIl(.O.9>2÷O.O4)-6=O,)X=例如：求一元函数量小值(fminbnd吩咐)由于fminbnd只能求一个最小值，通过视察函数可知微小值点为8和+8»x=fninbnd(y.O.inf)X=Inf»y=(X)1/(x-<).3)a2+O.()I)+1/(x).9)2-K).O46:»ix.y=fminbnd(y,0j11f)X=InfY=-6»x,yJ=fminbnd(y.-inf,O)%运行该吩咐时会报错.因为这里的X为NAN,这可能是因为fmindnd或DPlmg泣函数的代码有bug导致,f(.tu例如：求例题1-13中所定义以)的Oj上的定枳分JUJ.»symsx;»int(I(xO.3)2÷O.OI)÷1(x-O,9)2+O,O4).6.O.2)ans=134p+atan(35641981496012/6755()630868559)-12xyxljj(xey+z?)dxdydz例1-14求：虫积分中21及三正枳分。,3似忸。同。»SyrnSXy;»int(int(x4y,xJ.2),y,0J)ans=3/4»symsXyz;»inl(inl(int(%xp(y)+z八2.x0.l).y.0.1)z.l)ans=CXP(IN21/6例1-15已知'-5/'+6r+5,设该曲线在区间o上所用曲边梯形面积为s,试求当S分别为5,IO时的X的值.symsXt;y=$ym(tA3-5*tA26xt+5);s=int(y,t>,×);×l-solve(s-5,x)×1«，root(4.z.1)'root(o.z2)root(4.z3)JOO<(ol%4)wherel=12?*20z-20x2solve(sl,×)x2«'root(1.z.1)'root(1.z.2)root(1.z,3)rool(4.z.4)/whereOI=Zj邛.12z2*20z-40例1-16利用MAT1.AB吩咐求解无理数的近似值.(1)用函数零点吩咐(f«ro)求无理数6的近似侪：»f=(x)x-cxp(l);»fzero(f.2)ans=2.7183用定积分计算吩咐(trapz,quad.quadl>求无理数川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:设')=x"y+SinG').求dx6、.2原办symsXyn;f=sym(xAn*y+sin<y);diff(frx)ansdiff(fzy>ans-COS()Mdiff(diff(fxx)zx>aB-n2y(n-l)diff(difff,y),y>a11s-Sinbr)diff(diff(fxx),y>ans-11*,例“：求信-1.'*',匐点心二时:'(x+)z)dz.ASymSxyZ%声明符号更受，留近变量间必需用空将分开symsXy2:;int<x-y/(l+xA2>rx)</Ians2inn(xy<l，乂A2),y,0,t)ans-32(i2*nint(int(x*y/(l*xA2),y,0,xA(l/2),x,0,l)1 .£ans-28int(int(int(x+y+zrz,l-x-y),yrO,l-x),xr0,1)ans-124级数求和synsum)%求级23k(ans=inf即)»symsk;»f=(k)l/k:»symsum(f,kjjnf)ans=Inf%求汲数击品k"ME)»symsk;»f=(k)i/(k*(k+l)：>>synsu11Mf.kJ.i11f)ans=+=+1.+二T(ans=3/2、)%求纨敦33?3»symska;»f=(kW3”：»symsum(f,k,O.inf)ans=(3*a)2泰勒绽开(taylor)SymSXfy=l(l»taylor(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=0:»solve<y)ans=-<b+(bA2-4*a*c>A(IM2*a)(b-(b2-4*a*c)A(l/2>y(2*a)求方程fx=O关于变Itb的杼号解»solvc<y.b)ans=-(a*xA2+c)x微分方程(组)求解<ds«lvc)求方衽户=5的通解,戏认白变量为I»synsy(t);»dsolvc(dif(y)=5)ans二Cl+5*1求方程y'=x的通解,指定自变也为X»SymSy(x);»dslve(diff(y.x)x>ans=x22C2求方程y'=I+y,y(0)=l,y'(O>=O的特W»symsy(t);»dsolve(D2y=I+Dyy<O>=I.'Dy(O)=O*)ans=exp()-(x,=xy术方V姐)'=2'的通解.默认白变量为t>>syms>,(t)x();»z=dslve('Dx=x+y,.Dy=2x');»Z.Xans=C6*exp<2*l)-(C5i*exp<.t)V2»z.yans=C54cxXt)+C6*cxp(2*t)试验2MAT1.AB绘制二维、三维图形例2-1在子图形窗口中画出°2上正弦、余弦曲线。XIinspacc(0,2pi);y!=in(x);y2=<os(x>: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(,);ylat>cl(y);subplot(2,l,2);plot(×,y2,o'MarkerSize,j3);axis(-l,8,-2,2);title(,ycosx,);xlatM?lCx');ylabl(y);例2-3分别在两个图形新口画出地充正
