资源描述：
《现代控制实验报告》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、实验1利用MATLAB进行传递函数和状态空间模型间的转换例1.1代码：A=[010;001;-5-25-5];B=[0;25;-120];C=[100];D=[0];[num,den]=ss2tf(A,B,C,D)结果：num=00.000025.00005.0000den=1.00005.000025.00005.0000例1.2代码：A=[01;-25-4];B=[11;01];C=[10;01];D=[00;00];[num1,den1]=ss2tf(A,B,C,D,1)[num2,den2]=ss2tf(A,B,C,D,2)结果：num1=01.00004.000000-25
2、.0000den1=1.00004.000025.0000num2=01.00005.000001.0000-25.0000den2=1.00004.000025.0000例1.3代码：num=[001010];den=[16510];[A,B,C,D]=tf2ss(num,den)结果：A=-6-5-10100010B=100C=01010D=0实验要求代码：num=[0012;0153]den=[1234];[A,B,C,D]=tf2ss(num,den)结果：A=-2-3-4100010B=100C=012153D=00实验2利用MATLAB对状态空间模型进行分析例2.1代码：
3、A=[01;-10-5];B=[0;0];D=B;C=[10;01];x0=[2;1];[y,x,t]=initial(A,B,C,D,x0);plot(t,x(:,1),t,x(:,2))gridtitle('ResponsetoInitialCondition')xlabel('Time(sec)')ylabel('x1,x2')text(0.55,1.15,'x1')text(0.4,-2.9,'x2')图形：例2.2代码：A=[-1-1;6.50];B=[11;10];C=[10;01];D=[00;00];step(A,B,C,D)图形：例2.3代码：A=[0-2;1-3]
4、;B=[2;0];C=[10];D=[0];sys=ss(A,B,C,D);x0=[1;1];t=[0:0.01:20];u=cos(t);[y,T,x]=lsim(sys,u,t,x0);subplot(2,1,1),plot(T,x(:,1))xlabel('Time(sec)'),ylabel('X_1')subplot(2,1,2),plot(T,x(:,2))xlabel('Time(sec)'),ylabel('X_2')图形：实验要求代码：A=[01;-10];B=[0;0];D=B;C=[10;01];x0=[0;1];[y,x,t]=initial(A,B,C,D,
5、x0);plot(t,x(:,1),t,x(:,2))gridtitle('ResponsetoInitialCondition')xlabel('Time(sec)')ylabel('x1,x2')text(0.55,1.15,'x1')text(0.4,-2.9,'x2')图形：实验3利用MATLAB导出连续状态空间模型的离散化模型例3.1代码：A=[01;-25-4];B=[0;1];[G,H]=c2d(A,B,0.05)结果：G=0.97090.0448-1.12120.7915H=0.00120.0448实验要求T=1时的代码：A=[-10;0-2];B=[0;1];[G,
6、H]=c2d(A,B,1)结果：G=0.3679000.1353H=00.4323T=0.5时的代码：A=[-10;0-2];B=[0;1];[G,H]=c2d(A,B,0.5)结果：G=0.6065000.3679H=00.3161实验4利用MATLAB求解极点配置问题例4.1代码：A=[010;001;-1-5-6];B=[0;0;1];J=[-2+j*4-2-j*4-10];K=acker(A,B,J)结果：K=199558代码：A=[010;001;-1-5-6];B=[0;0;1];J=[-2+j*4-2-j*4-10];K=place(A,B,J);sys=ss(A-B*
7、K,[0;0;0],eye(3),0);t=0:0.01:4;x=initial(sys,[1;0;0],t);x1=[100]*x';x2=[010]*x';x3=[001]*x';subplot(3,1,1);plot(t,x1),gridtitle('ResponsetoInitialCondition')ylabel('x1')subplot(3,1,2);plot(t,x2),gridylabel('x2')subplot(3,1,3);plot
