《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,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):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);