- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
系统的能控性与能观性分析及状态反馈极点配置
实 验 报 告
课程 自动控制原理 实验日期 12 月26 日
专业班级 姓名 学号
实验名称 系统的能控性与能观性分析及状态反馈极点配置 评分
批阅教师签字
一、实验目的
加深理解能观测性、能控性、稳定性、最小实现等观念,掌握状态反馈极点配置方法,掌握如何使用MATLAB进行以下分析和实现。
1、系统的能观测性、能控性分析;
2、系统的最小实现;
3、进行状态反馈系统的极点配置;
研究不同配置对系统动态特性的影响。mineral;
(b)已知连续系统的传递函数模型,,当a 分别取-1,0,1时,判别系统的能控性与能观测性;
(c)已知系统矩阵为,,,判别系统的能控性与能观测性;
(d)求系统的最小实现。
2.实验内容
原系统如图所示。图中,X1和X2是可以测量的状态变量。
图试设计状态反馈矩阵使系统加入状态反馈后其动态性能指标满足给定的要求 (1) 已知:K=10,T=1秒,要求加入状态反馈后系统的动态性能指标为:
σ%≤20%,ts≤1秒。
(2) 已知:K=1,T=0.05秒,要求加入状态反馈后系统的动态性能指标为:
σ%≤5%,ts≤0.5秒。
状态反馈后的系统,如图所示:
图状态反馈后分别观测状态反馈前后两个系统的阶跃响应曲线,并检验系统的动态性能
指标是否满足设计要求。 (1-1)
其中A为n×n维状态矩阵;B为n×m维输入矩阵;C为p×n维输出矩阵;D为p×m维传递矩阵,一般情况下为0。
系统的传递函数阵和状态空间表达式之间的关系如式(1-2)所示:
(1-2)
式(1-2)中,表示传递函数阵的分子阵,其维数是p×m;表示传递函数阵的分母多项式,按s降幂排列的后,各项系数用向量表示。
系统的能控性、能观测性分析是多变量系统设计的基础,包括能控性、能观测性的定义和判别。
系统状态能控性定义的核心是:对于线性连续定常系统(1-1),若存在一个分段连续的输入函数u(t),在有限的时间(t1-t0) (1-3)
系统状态能观测性的定义:对于线性连续定常系统(1-1),如果对t0时刻存在ta,t0ta,根据[t0,ta]上的y(t)的测量值,能够唯一地确定系统在t0时刻的任意初始状态x0,则称系统在t0时刻是状态完全能观测的,或简称系统在[t0,ta]区间上能观测。
状态能观测性判别方法也分为2种:一般判别和直接判别法,后者是针对系统的系数阵A是对角标准形或约当标准形的系统,状态能观性判别时不用计算,应用公式直接判断,是一种直接简易法;前者状态能观测性分为一般判别是应用最广泛的一种判别法。
状态能观测性判别式为:
(1-4)
系统的传递函数阵和状态空间表达式之间的有(1-2)式所示关系。已知系统的传递函数阵表述,求其满足(1-2)式所示关系的状态空间表达式,称为实现。实现的方式不唯一,实现也不唯一。其中,当状态矩阵A具有最小阶次的实现称为最小实现,此时实现具有最简形式。
2、状态反馈极点配置
一个受控系统只要其状态是完全能控的,则闭环系统的极点可以任意配置。极点配置有两种方法:①采用变换矩阵T,将状态方程转换成可控标准型,然后将期望的特征方程和加入状态反馈增益矩阵K后的特征方程比较,令对应项的系数相等,从而决定状态反馈增益矩阵K;②基于Carlay-Hamilton理论,它指出矩阵满足自身的特征方程,改变矩阵特征多项式的值,可以推出增益矩阵K这种方法推出增益矩阵K的方程式叫Ackermann公式。 num=[1 -1];den=[1 10 27 18];[a,b,c,d]=tf2ss(num,den);
Qc=ctrb(a,b)
Qc =
1 -10 73
0 1 -10
0 0 1
rank(Qc)
ans =
3
Qo=obsv(a,c)
Qo =
0 1 -1
1 -1 0
-11 -27 -18
rank(Qo)
ans =
3
num=[1 0];den=[1 10 27 18];[a,b,c,d]=tf2ss(num,den);
Qc=ctr
文档评论(0)