- 1、本文档共15页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一. 状态空间模型和模型变换
1.1 已知SISO系统的传递函数为
(1)给出系统的一个状态空间实现,写出状态空间模型的完整表达式。
(2)求出系统特征值,比较系统的特征值和极点是否一致,为什么?所得的状态空间实现是否是最小实现?
解:(1)
G(s)==
将其输入到MATLAB 工作空间:
代码
num=[1,2];
den=[1,12,45,58,24];
TFG=tf(num,den);
SG=ss(TFG)
运行结果
a =
x1 x2 x3 x4
x1 -12 -5.625 -1.813 -0.75
x2 8 0 0 0
x3 0 4 0 0
x4 0 0 1 0
b =
u1
x1 0.25
x2 0
x3 0
x4 0
c =
x1 x2 x3 x4
y1 0 0 0.125 0.25
d =
u1
y1 0
Continuous-time model.
根据运行结果,得到系统的一个状态空间实现,状态空间模型的完整表达式为:
= 1 \* GB3 \* MERGEFORMAT ①在MATLAB 工作空间中,求系统特征值的代码为
A=[-12,-5.625,-1.813,-0.75;8,0,0,0;0,4,0,0;0,0,1,0];
V=eig(A)
运行结果
V =
-6.0019
-3.9964
-1.0335
-0.9681
故系统的特征值为
由SISO系统的传递函数知系统的极点为p1=p2=-1,p3=-4,p4=-6。
由上可知系统的特征值和极点是一致的。因为系统矩阵的特征值方程和系统传递函数的特征方程是等价的,特征值也是特征方程D(s)=0的根。在经典控制理论中,系统传递函数的特征方程D(s)=0解出的根具有负实部,系统就是稳定的;在现代控制理论中,特征值具有负实部系统也是稳定的。故系统的特征值和极点是一致的。
= 2 \* GB3 \* MERGEFORMAT ②在MATLAB 工作空间中,验证所得的状态空间实现是否是最小实现的代码为:
A=[-12,-5.625,-1.813,-0.75;8,0,0,0;0,4,0,0;0,0,1,0];
B=[0.25;0;0;0];
C=[0,0,0.125,0.25];
D=[0];
G=ss(A,B,C,D);
Gm=minreal(G)
运行结果
a =
x1 x2 x3 x4
x1 -12 -5.625 -1.813 -0.75
x2 8 0 0 0
x3 0 4 0 0
x4 0 0 1 0
b =
u1
x1 0.25
x2 0
x3 0
x4 0
c =
x1 x2 x3 x4
y1 0 0 0.125 0.25
d =
u1
y1 0
Continuous-time model.
该系统的传递函数没有零极点对消,并且由运行结果可知所得的状态空间实现是最小实现。
1.2 已知SISO系统的状态空间表达式为
(1)计算系统在等价变换下的状态空间方程,其中。比较变换前后的系统特征值是否一致,为什么?
(2)求系统的传递函数。
解:(1)
= 1 \* GB3 \* MERGEFORMAT ①计算系统在等价变换下的状态空间方程,其中。
将其输入到MATLAB 工作空间:
代码
A=[0 1 0;3 0 2;-12 -7 -6];
B=[2 1 7];
C=[1 1 1];
G=ss(A,B,C,0);
P=[2 1 1;1 0 1;3 2 -1];
G1=ss2ss(G,P)
运行结果
a =
x1 x2 x3
x1 -5
文档评论(0)