matlab上机题.doc

第一次上机： 1 (1)利用matlab画一张正弦波形，x=0:0.01:10*pi;y=sin(x);plot(x,y) (2)画出的曲线。t=0:pi/50:10*pi;plot3(sin(t),cos(t),t) 2利用subplot（n，m，k）函数，将画图区分成4个区域，分别作图：1画一正弦波形；2画一余弦波形；3画y=x的直线；4画的曲线。 x=-2*pi:0.01:2*pi; subplot(2,2,1);plot(x,sin(x)); subplot(2,2,2);plot(x,cos(x)); subplot(2,2,3);plot(x,x); subplot(2,2,4);plot(x,x.^2); 3画出的曲线。ezplot('x^2*sin(x+y^2)+y^2*exp(x+y)+5*cos(x^2+y)'); 4已知 A=[1 2 3;4 5 6;7 8 9],求矩阵A的秩。A=[1 2 3;4 5 6;7 8 9];k=rank(A) 5利用diff函数，求解的导函数。 y='x^2+3*x';dy=diff(y,'x') 6 利用dsolve函数，求解微分方程。y=dsolve('Dy+2*x*y=x*exp(-x)') 第二次上机： 1 已知系统传递函数为，求系统状态空间的表达式。 num=[0 5 1 1];den=[1 6 11 6];[A,B,C,D]=tf2ss(num,den) 2 已知状态空间表达式：A=[0 2; -2 -3] B=[1 0 ;1 1] C=[1 0;1 1] D=0 ,求系统的传递函数。 A=[0 2; -2 -3]; B=[1 0 ;1 1];C=[1 0; 1 1];D=zeros(2,2); [num1,den1]=ss2tf(A,B,C,D,1);[num2,den2]=ss2tf(A,B,C,D,2) 3已知两个子系统的传函分别是：，，（1）求两个子系统串联后系统传函。num1=1;den1=[1,4,3];num2=1;den2=[1,5,0];[num,den]=series(num1,den1,num2,den2) （2）求两个子系统并联后系统传函。num1=1;den1=[1,4,3];num2=1;den2=[1,5,0]; [num,den]=parallel(num1,den1,num2,den2) 4 已知系统的状态空间表达式为：A=[0 1;1 3] B=[0;1] C=[1 1] D=zeros[1 1]。求系统的传递函数。A=[0 1;1 3];B=[0;1];C=[1 1];D=zeros(1,1);[num1,den1]=ss2tf(A,B,C,D,1) 5求系统的系统特征值、自然频率和阻尼系数。 format;num=[5 6 1];den=[1 2 3];[wn,zeta]=damp(den),roots(den) 6 num1=[10];den1=[1 1]; num2=[1];den2=[2 0.5]; num3=[100];den3=[1]; num4=[0.1];den4=[1]; [na,da]=series(num1,den1,num2,den2); [nb,db]=feedback(na,da,num4,den4,-1); [nc,dc]=series(num3,den3,nb,db); [num,den]=cloop(nc,dc,-1); printsys(num,den) 已知系统的框图如上图示，求系统的传递函数 7 num=[1];den=conv([1,0],[2,0.5]);T=1; [numd1,dend1]=c2dm(num,den,T); printsys(numd1,dend1,'z') 已知系统如图所示，求系统在T=1s，开环脉冲传函G(s)。 8利用matlab生成一个3输入3输出的4阶的稳定续系统. [A,B,C,D]=rmodel(4,3,3) 第三次上机： 1已知某系统A=[0 1;-1 -1]，利用eye函数判断系统稳定性。 2已知某系统的闭环传函为 ，判断系统的稳定性，并求出系统的闭环极点。 num=[3,6,5,4,8,2];den=[6,7,1,4,3,9,1];r=roots(den),pzmap(num,den) 3已知某系统的开环传函为 ，求系统在单位反馈下的阶跃响应曲线和超调量。 num0=30; den0=[1 5 6];[num,den]=cloop(num0,den0);t=0:0.1:10;[y,x,t]=step(num,den,t); plot(t,y);M=((max(y