MATLAB练习题和答案.docx

《MATLAB练习题和答案.docx》由会员分享，可在线阅读，更多相关《MATLAB练习题和答案.docx（31页珍藏版）》请在优知文库上搜索。

1、限制系统仿真试验MatIab部分试验结果目录试验一MATLAB基本操作1试验二Matlab编程5试验三Matlab底层图形限制6试验四限制系统古典分析12试验五限制系统现代分析15试验六PID限制器的设计19试验七系统状态空间设计23试验九直流双闭环调速系统仿真25试验一MATLAB基本操作1用MATLAB可以识别的格式输入下面两个矩阵1233一-l4z43678-A二2357233554+2/1357B=26+7/534232391895431894一再求出它们的乘积矩阵C,并将C矩阵的右下角2X3子矩阵赋给D矩阵。赋值完成后，调用相应的叮嘱查看MATLAB工作空间的占用状况。A=1Z2,3

2、Z3;2Z3,5Z7;1Z3,5Z7;3,2,3Z9;1Z8,9Z4;B=l+4i,4,3,6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;l,8,9,5,4,3;C=A*B;D=C(4:5,4:6);whosNameSizeBytesClassAttributesA5x4160doubleB4x6384doublecomplexC5x6480doublecomplexD2x396doublecomplex2选择合适的步距绘制出下面的图形sin(lO,其中(TJ)t=-l:0.1:1;y=sin(l.t);plt(t,y)3对下面给出的各个矩阵求取矩阵的行列式、秩、特征

3、多项式、范数、特征根、特征向量和逆矩阵。-7.53.500-57658334.1071087A=B=09103-1.568109003.719.3579101234-3-3-2456785-518C=,D91011121185-7131415165-1-3-1A=7.5,3.5f0,0;8,33,4.1,0;0,9z103z-1.5;0,0,3.7,19.3;B=5,7,6z5;7x10,8,7;6,8,10,9;5,7,9,10;C=l:4;5:8;9:12;13:lrtf6;D=3z-3,-2r45,-5zl,8llz8,5z-75z-lz-3,-l;det（八）;det(B)；det(C

4、);det(D);rank（八）;rank(B);rank(C);rank(D);a=oly（八）;b=oly(B);c=poly(C);d=poly(D);norm（八）;norm(B);norm(C);norm(D);v,d=eig(A,nobalance,);v,d=eig(B,nobalance,);vrd=eig(C,nobalance,);v,d=eig(D,nobalance,);m=inv（八）;n=inv(B);p=inv(C);q=inv(D);4求解下面的线性代数方程，并验证得出的解真正满足原方程。B=4;7;-l;0;X=AB;C=A*X;(b)A=1,3,2,13;7

5、,2,1,-2;9,15,3,-2;-2,-2,11,5;B=9z0;6z4;llz7;-2,-l;X=AB;C=A*X;5.(1)初始化-10*10矩阵,其元素均为1ones(10,10);(2)初始化-10*10矩阵，其元素均为0zeros(10,10);(3)初始化10*10对角矩阵v=l:10;diag(v);(4)输入A=75;256；315,B=l11;222;333,执行下列叮嘱，理解其含义A(2,3)表示取A矩阵第2行、第3列的元素；A(:,2)表示取A矩阵的第2列全部元素；A(3,:)表示取A矩阵第3行的全部元素；A(:,1:2:3)表示取A矩阵第1、3列的全部元素；A(:,

6、3).*B(:,2)表示A矩阵第3列的元素点乘B矩阵第2列的元素A(:,3)*B(2,:)表示A矩阵第3列的元素乘以B矩阵第2行A*B矩阵AB相乘A.*B矩阵A点乘矩阵BA2矩阵A的平方A2矩阵表示求矩阵A的每一个元素的平方值B/A表示方程AX=B的解XB./A表示矩阵B的每一个元素点除矩阵A的元素6在同一坐标系中绘制余弦曲线y:cos(t-0.25)和正弦曲线y=sin(L0.5),t0,2,用不同颜色，不同线的类型予以表示，留意坐标轴的比例限制。t=0:0.01:2*pi;yl=cos(t-0.25);plot(tzylz,r-)holdony2=sin(t-0.5);plot(t,y2z

7、k,)试验二Matlab编程1分别用for和while循环结构编写程序，求出63X=2,=1+2+22+23+262+263f=0并考虑一种避开循环的简洁方法来进行求和。(a) j=l;n=0;sm=l;forn=n+l:63fori=l:nj=j*2;endsum=sum+j;j=l；endsum(b) j=l;n=l;sum=l;whilen=64i=l;whilein+lj=j*2;i=i+l;endn=n+l;sum=sum+j;j=l；endSum(c) i=0:63;k=sum(2.i);2计算l+2+n2000时的最大n值s=0;m=0;while(sDy=(x)=hlDx,xZ

8、)h,XVD存放于文件ff.m中，令D=3,h=l求出，1.5),40.5),5).D=3;h=l;x=-2*D:l/2:2*D;y=-h*(x=-D)&(xD);plot(x,y);gridonfl=y(find(x=-1.5)f2=y(find(x=0.5)f3=y(find(x=5)试验三Matlab底层图形限制1在MATLAB叮嘱行中编程得到y=sin(t)和yl=cos函数,plot(t，y);figure(10);plot(t,yl);t=-pi:0.05:pi;y=sin(t);yl=cos;lot(t,y)figure(IO);lot(t,yl)2在MATLAB叮嘱行中键入h=

9、get(O),查看根屏幕的属性，h此时为根屏幕句柄的符号表示，0为根屏幕对应的标号。h=get(O)h=BeingDeleted:offBusyAction:queueButtonDownFcn:CallbackObject:UChildren:2xldoubleClipping:,on,CommandWindowSize:8927JCreateFcn:CurrentFigure:1DeleteFcn:Diary:,offDiaryFile:diary,Echo:offFixedWidthFontName:CourierNewFormat:,short,FonnatSpacing:looseH

10、andIeVisibility:on,HitTest:,on,Interruptible:,on1.anguage:zh-cn.gbk,MonitorPositions:1I1440900More:offParent:PointerLocation:1048463PointerWindow:0RecursionLimit:500ScreenDepth:32ScreenPixelsPerInch:96ScreenSize:111440900Selected:offSeIectionHighlight:,onShowHiddenHandles:offTag:Type:rootUIContextMe

11、nu:Units:pixelsUserData:Visible:,on,3hl=get(l);h2=get(10)51,10分别为两图形窗口对应标号，其中1为MatIab自动支配，标号10已在AgUre(Io)中指定。查看hl和h2属性，留意CurrentAxes和CurrenObject属性。h1=get(1)hl=Alphamap:1x64doubleBeingDelcted:,offBusyAction:queueButtonDownFcn:Children:170.0012Clipping:,onCIoseRequestFcn:,closereqColor:0.80000.80000.

12、8000Colormap:64x3doubleCreateFcn:CurrentAxes:170,0012CunentCharacter:,CurrentObject:CurrcntPoint:00DeleteFcn:DockControls:onFileName:FixcdColors:10x3doubleHandleVisibiIity:on,HitTest:onIntegerHandIe:,on,Interruptible:onInvcrtHardcopy:on,KeyPressFcn:KeyReleaseFcn:MenuBar:figureMinColormap:64Name:Next

13、PIot:,add,NumberTitle:on,PapcrOrientation:,portraitPaperPosition:0.63456.345220.304615.2284PaperPositionMode:,manualPaperSize:20.984029.6774PaperType:A4,PaperUnits:centimetersParent:0Pointer:arrow,PointerShapcCData:16x16doublePointerShapeHotSpot:11Position:440378560420Renderer:painters*RendererMode:auto,Resize:,on,ResizeFcn:Selected:,offSelectionHighlight:on