MATLAB2MATLAB2.doc

《Matlab语言控制系统》设计 （ 2010 -- 2011 学年度 第 一 学期） 院 系： 自动控制工程系 班 级： 自动化二班 姓 名： 学 号： 0704240205 提 交 日 期： 2010 年 12 月 22 日 实验一 求状态空间表达式 实验题目： 一、控制系统微分方程为： y⑷+10y⑶+35y⑵+50y⑴+24y=u⑶+7u⑵+24u⑴+u y⑷+3y⑶+2y⑴=-u⑴ 求其状态空间表达式。 实验原理： 介绍运用MATLAB对微分方程组的仿真方法。首先给出一阶、二阶的解,并利用SIMULINK建立仿真模型,给出仿真结果图形;然后给出一个具体的应用实例。仿真算法具有模型设计简单、结果直观等特点有很强的使用价值。 实验程序运行如下： To get started, select MATLAB Help or Demos from the Help menu. The element type name must be terminated by the matching end-tag /name. Could not parse the file: c:\matlab7\toolbox\ccslink\ccslink\info.xml  format compact f=[1 7 24 24]; g=[1 10 35 50 24]; printsys(f,g,s) num/den = s^3 + 7 s^2 + 24 s + 24 --------------------------------- s^4 + 10 s^3 + 35 s^2 + 50 s + 24 disp() z=roots(f) z = -2.7306 + 2.8531i -2.7306 - 2.8531i -1.5388 p=roots(g) p = -4.0000 -3.0000 -2.0000 -1.0000 k=f(min(find(f(:)~=0)))/g(1) k = 1 [z1,p1,k1]=tf2zp(f,g) z1 = -2.7306 + 2.8531i -2.7306 - 2.8531i -1.5388 p1 = -4.0000 -3.0000 -2.0000 -1.0000 k1 = 1 disp() [r,p,h]=residue(f,g) r = 4.0000 -6.0000 2.0000 1.0000 p = -4.0000 -3.0000 -2.0000 -1.0000 h = [] disp() [A,B,C,D]=tf2ss(f,g) A = -10 -35 -50 -24 1 0 0 0 0 1 0 0 0 0 1 0 B = 1 0 0 0 C = 1 7 24 24 D = 0 printsys(A,B,C,D) a =  x1 x2 x3 x4 x1 -10.00000 -35.00000 -50.00000 -24.00000 x2 1.00000 0 0 0 x3 0 1.00000 0 0 x4 0 0 1.00000 0 b = u1 x1 1.00000 x2 0 x3 0 x4