MATLAB报告王志文.docx

《MATLAB报告王志文.docx》由会员分享，可在线阅读，更多相关《MATLAB报告王志文.docx（18页珍藏版）》请在课桌文档上搜索。

1、请创建如图所示的结构数组姓名编号指标江明顺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.00yy=-0.44701.97803.11(X)5.25005.020

2、04.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.9043xx=0:0.0):l;yy=polyval(P.xx);plot(xx.yy,-b,x0,y0.Y.Markcrsizc,20);Iegend(拟合曲线,惊始数据Y1.OCatiOn.SouthEasCxlabel(x)试验二一.试睑目的及试验环境1 .试验目的驾驭矩阵,多项式，导函数的计算在MAT1.AB上如何应用。2 .试验环境MAT1.AB二.试验内容习题1：设矩阵A和B满足关系式AB=A+2B已知A=R

3、2R求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.01.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.formalralpI=conv(1t0,2.conv(1,4,1,1

4、);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).

5、(1-69)I错误：网括号或方括号不对称或异样。p=conv(I-4,conv(15).(1-69)P=Columns1through5I-5-17129-180px=poly2str(p.x)PX=x4-5x3-I7x2+129x-180x=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)+5xx=-3:pi

6、/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实现离散数据可

7、视化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(%),仿照运行，分析衣现形式不同的缘

8、由。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

9、=x2,);(ext(6*pb4.-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);plX(x,y3);titlc(y3=xA24cos(24xH；yl=x.A2;subpl(2.2.1):bar:IiHayl=x八2条形图,)：subplo(2.2.2);stairs(x.yl):HtkeyI=小2阶梯

10、图,)：SUbPloo2在一5r.*(x0)z=Cosxcosye4Je中X的21个值匀称分布在卜551范围,y的31个位匀称分布在0,10卜要求运用subplot(2.1.Iillsubplot(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-);subplx(2,l,2);COlMOUr3(x.yz5O);补充：符号计算试骐【习题1】试求习JH2求

11、y11.(2*-j)2kJ【习J3用solve函数解线性方程也-x1+2.v2=22.v1.r2+x3=34.r15x2+Ix3=0xl+x2+5工=-5【习题5】解方程组a=1;b=2;c=3;d=4：M1=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,1

12、1),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):plotcx(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.

13、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).*(x0);plot(x,y)X=IinSP

14、aCe(-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);上机补充题：MiBsymsx;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)

15、,lxdfd(2=diff(f,t,2)dfdt2=0.6t10,0)dfdxdt=difftdiff(f.x).t)dfdxdt=I0,0)I-sin(x),0第3symsk;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=I234M2=ahc.d:eval(M2)ans=I234M3=Sym(Ia.bed):eval(M3)ans=I234第5题=(-1.2.0;A=-120A=-1.2,0;2.1.1;4.5.7;.,5;B=2:3:0:-5):f=AB13/111

16、(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四.窝试数据及运行结果“HiQXxOX目1!r*11三Qv.ZF11QU.O117U.1CaU.OA：.1.：v.0：个个iAu.*三米三:米三T0-五.总结心得体会：通过这次忒验，我熟识了MAI1.AB软件应用环境，理解了MAT1.AB的基本原理，驾驭了基本的MAT1.AB函数操作及其运用，并懂得/常用窗口的功能和运用方法。在试验得过程中学会如何调试和验证结果