原创力文档- 1、本文档共10页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
MATLAB1MATLAB1
《Matlab语言控制系统》设计
（ 2010 -- 2011 学年度 第 一 学期）
院 系： 自动控制工程系
班 级： 自动化二班
姓 名：
学 号： 0704240206
提 交 日 期： 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）