MATLAB报告王志文.docx

请创建如图所示的结构数组姓名编号指标江明顺071023身高：176.体用：82于越忠060134身高：168,体术：74邓拓050839身高：182,体玳：77三.方案设计根据书上的例子，熟识MAT1.AB的运用方法，解决问题四.浦试数据及运行结果习题1.»Xo=O:0.1:1;»yO=(-.447,1,978,3.11.5.25.5.02.4.66.4.01.4.58.3.45.5.35.9.22;»xx=00.1(XX)0.2(XX)0.30000.40000.5(XX)0.60000.70000.80000.90()01.00>>yy=-0.44701.97803.11(X)5.25005.02004.66004.01(X)4.58003.45(X)5.35(X)9.22(X)»n=3;»P=POlyfh(Xo.y,n)P=56.6915-87.117440.0070-0.9043»xx=0:0.0):l;»yy=polyval(P.xx);»plot(xx.yy,"-b',x0,y0.'Y.'Markcrsizc',20);»Iegend('拟合曲线,惊始数据Y1.OCatiOn'.SouthEasC»xlabel("x')试验二一.试睑目的及试验环境1 .试验目的驾驭矩阵,多项式，导函数的计算在MAT1.AB上如何应用。2 .试验环境MAT1.AB二.试验内容习题1：设矩阵A和B满足关系式AB=A+2B已知A=R2R求BI 1OII 23习题2:术($2+4K$+1)的W余“多项式.习题3：将表达式(x-4)(x+5)(x*6x+9)展开为多项式形式，并求其对应的一元n次方程的根。习题4设/()=Vv3+2.vj-x+12+iJx+5+5,v+2求导函数，并画出导函数在-3,3区间内的图像。三.方案设计习题1.»=(4.2,3J.1.0>1.2,3);»B=inv(A-2*cyc(sizc（八）)*A;»BB=3.0000-8.0000600002.0(KX)-9.0000-6.0000-2.(XXX)12.(XMX)9.(MXM)习题2.»formalral»pI=conv(1t0,2.conv(1,4,1,1);»p2=(IOll;»q,r=deconv(pl,2);»eg=，商多项式为'©='余多项式为：»disp(cq.poly2str(q.'s,),disp(cr.poly2str(r.,s,)disp(cq,poly2str(q,s,),disp(cr,poly2str(r,)l)IError:Unbalancedorunexpectedparenthesisorbracket.»disp(cq.poly2str(q,s,),disp(cr.poly2str(r,'s")商多项式为s+5余多项式为5s2+4s+3习题3.»p=conv(1-4.conv(I5).(1-69)p=conv(l-4.conv(l5).(1-69)I错误：网括号或方括号不对称或异样。»p=conv(I-4,conv(15).(1-69)P=Columns1through5I-5-17129-180»px=poly2str(p.'x')PX=x4-5x3-I7x2+129x-180»x=roots(p)-5433习题4.»symsx;»f(x)=sqrt(xA3+2*xA2-x+l2)+(x+5)A(l/6)+5*x+2;»f(x)=sqrt(x3+2*x2-x+12)+(x+5).(I6)+5*x+2:»p=diff(f(x),x)I(6*(x+5)(56)+(3*x2+4*x-X(2*(x3+2*x2-x+I2)(l2)+5»xx=-3:pi/2OO:3;»yx=subs(p.x.xx):»plot(xx,yx);四.测试数据及运行结果习期4五.总结心得体会：通过这一周的试脸，我感受到了MAT1.AB在数学运用画图方面的强大功能，让我感受颇深同时也对MAT1.AB这门课程产生了深厚的学习爱好。试验三一.试睑目的及试验环境1 .试验目的体验数据可视化，加深对矩阵的运算，还有求极限。2 .试验环境MAT1.AB二.试验内容1.仿照运行，体会数据可视化方法。%用SIem实现离散数据可视化n=0:l2:y=1.abs(n-6):stcm(n,y)I.已知”=(M,12.m-6l运行下面程序,体会离散数据可视化方法.%Hlplot实现离散数据可视化n=0:12;%产生一组自变量数据y=17abs(n-6):%计算相应点的函数值plot(n,y.,r*'.'MarkcrSzc',20)%用红花标出数据点gridon%画坐标方格说明：PIo«和MCm指令均可以实现离散数据的可视化，但通常PkH更常用于连续函数中特别点的标记；而stem广泛运行与数字信号处理中离放点的图示。用户在运行上面例程时会发觉在吩时窗口出现警告：Warning:Dividebyze!即警告程序中出现非零数除以O的指令.MAT1.AB对于这种状况并不中止程序只是给该项减值为inf以做标记.2.下面时用图形表示连续调制波形F=Sin")sin(%),仿照运行，分析衣现形式不同的缘由。V=(0.5+4叱COSX3,设I+，在x=02n区间取点,绘制函数曲线X=linspacc(0,2*pi.101);y=(0.5+3*sin(x),(l+x.*x).*cos(x);PlO(X,y)4.已知y1=x2.>2=cos(2x).y3=yly2.完成下列操作：(I)在同一坐标系下用不同的颜色和线型绘制三条曲规:(2)分别用条形图、阶梯图、杆图和填充图绘制三条曲战。(I)x=linspacc(-24pi.24,pi.l()();yl=x.2;y2=cos(2x);y3=yl产y2:plot(x,y1.T>-,.x.y2rf.x.y3.,y-'):tcxt(4,)6,McftanowyI=x2,);(ext(6*pb'4.-1downaowy2=cos(2;tx)r);x=linspace(-2tpi.2fpi.lOO);yl=x2;y2=cos(24x):y3=yl.4y2;subplo<(1.3.l);plo<(x.yl):titlc('yl=xA2'):SUbPIO<(1,3,2);plo(x,y2);titlc('y2-cos(2*xy);SUbPkH(133);pl<X(x,y3);titlc('y3=xA24cos(24xH；yl=x.A2;subpl<(2.2.1):bar<x.yl>:IiHa'yl=x八2条形图,)：subplo<(2.2.2);stairs(x.yl):HtkeyI=小2阶梯图,)：SUbPlo<(2.2,3);sem(xtyl);iUe('yl=W2杆图subp!o(2.2.4):fll(x,y1.Y);3.已知x+6y=Ie2'X0)x+F77V>o2在一5<r<=5区间绘制函数曲雄.x=S0.01:5y=(x+sqrt(pi)/exp(2>.*(x<=0)+l/2*log(x+sqrl(l+x.A2).*(x>0)z=Cosxcosye4Je中X的21个值匀称分布在卜551范围,y的31个位匀称分布在0,10卜要求运用subplot(2.1.I>illsubplot(2,1,2>将产生的曲面图和劈高图画在同一个也口上.X=IinSPaee(-5.5.21);y=linspacc(0.10,31);x,y=mcshgrid(x.y);z=cos(x)?cos(y).*exp(-sqrt(x,A2+y.A2)M)；SUbPkM(2.1.1):surf(x,yj-);subpl<x(2,l,2);COlMOUr3(x.yz5O);补充：符号计算试骐【习题1】试求习JH2求y11.(2*-j)2kJ【习J3用solve函数解线性方程也-x1+2.v2=22.v1÷.r2+x3=34.r1÷5x2+Ix3=0xl+x2+5工=-5【习题5】解方程组»a=1;b=2;c=3;d=4：»M1=<,bjc,c»M2三'a1bc,d'val(M2)»M3=«ymfa.b;c,dD;vaM3)【习题4】指出下面的M1.M2,M3分别是什么.并上机发证.【习题6】解循分方程26.世常0分方用及It傅个初蛤条f!为翳ecoj(2jr)-F%(°)=°F(O)=I三求解读效分方程(今解.=.方案设计第IjS»n=0:l2;»y=l.abs(n-6);»ploi(n,y,'r*',MarkerSize',20)»gridon»n=0:12;»y=l.abs(n-6);»stcm(n.y)第2题»11=(0=ll)ll*pi;»yI=sin(tl).*sin(9tl);>>t2=(0:100)"00*pi;»y2=sin(t2).*sin(9*t2);>>subploi(2.2.1).ploi(tl,yl.'r')axis(0.PiJljD.UUe('子图(I),):>>subploi(2,2,2),plo(i2,y2,'r,),axis(0,pi,-l,lIhiiilefJ1IS(2)');>>subplo(2,2,3),lo(tI,yJ,y1,"r,),axis(),pi,-1,11),EleC子图(3)>>subpkx(224),pkx(t2,y2)axis(0,pi,-1,1DJiUeC子图(4),)第3题»x=linspacc(0.2*pi.IO1);»y=(O.5+3*siXx)(I+x.*x).*cos(x):>>plo<(x.y)第4题»x=1inspace(-2*pi,2*i,I(X);»yl=x.a2;»y2=cos(2*x):»y3=yI.*y2;»pot(x.y1.blX.y2Hy3,'y./):»text(4.16.,Meftarrowyl=x2');>>tcx(6*pi4,-1.Mownarrowy2=cos(2*x)r);>>tcxt(-1.5*pi.-2.25*pi*piupaowy3=y1*y21);»x=linspace(-2*pi,2*i,I(M);»yl=x.a2;»y2=cos(2*x):>>y3=yl.*y2;>>subplot(1.3.l);>>lot(x.yl);>>titlc('y=x2');>>subplot(1.3.2);>>lot(x,y2);>>titlc('y2=cos(x*2)');>>subplot(1.3.3);>>lot(x,y3);>>titlc(,y3=x2*cos(2*x);>>subplot(2,2.l);>>bar(x,yl);>>title(,y1=x2条形图')：>>subplot(2,2.2);>>stairs(x,yl);>>title('yl=x2阶梯图)：subplot(2.2.3);stem(x,yl);titlc('yl=x2杆图为subplot(2.2.4);fiU(x,yl,'r');title('yl=x2填充图工»X=-5:0.01:5:y=(x+sqr(i)exp(2).*(x<=0)+12*log(x+sqrt(I+x.2).*(x>0);plot(x,y)»»X=IinSPaCe(-5.5,21)：y=IinsPaCe(0,10,31);x.y=ncshgrid(x.y):z=cos(x).*cos(y).*exp(-sqrt(x.A2+y.A2)/4)：SUbPlOl(2,1,1);surf(x.>z);subpiot(2.1,2);contour3(x,y50);»上机补充题：MiB»symsx;SymSk;f=(l-lx)(k*x);1.f=Iimit(CXJnf);»1.fl=vpa(subs(1.f,k,syn(,-1,),48)1.fl=第2题»symsatx：»f=(a.t3：t*cos(x).log(x)：>>df=diff(Odf=I0.0)-t*sin(x),lx»dfd(2=diff(f,t,2)dfdt2=0.6t10,0)»dfdxdt=difftdiff(f.x).t)dfdxdt=I0,0)I-sin(x),0第3»symsk;f=l/(2*k-l)A2,(-l)Ak/k):a=symun(f,k,Un11Ipi28.-log(2)第4题»a=i;b=2;c=3;d=4;Ml=|a,b：c,dMl=I234»M2='ahc.d:»eval(M2)ans=I234»M3=Sym(Ia.bed'):eval(M3)ans=I234第5题>>=(-1.2.0;A=-120»A=-1.2,0;2.1.1;4.5.7;.,5;B=2:3:0:-5):f=AB13/111(X)77-123/77第6题>>y=dsolve(,D2y=cos(2*x)-yVDy(0)=(),'y(0)=l')y=cos(2*x)-cos(t)*(cos(2*x)-1)»y=simplc(y)>Insym.sinplcat41y=I-4<!sin(V2)A2*sin(x)A2四.窝试数据及运行结果“HiQXxO®X目1!r*11三Qv.ZF11QU.O117U.1CaU.OA：.×1.：v.0：个个iAu.*三米三:米三T0-五.总结心得体会：通过这次忒验，我熟识了MAI1.AB软件应用环境，理解了MAT1.AB的基本原理，驾驭了基本的MAT1.AB函数操作及其运用，并懂得/常用窗口的功能和运用方法。在试验得过程中学会如何调试和验证结果