- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第十四次课MATLAB.ppt
* * 内容: 三 极点配置 基于状态反馈的极点配置法就是通过状态反馈将系统的闭环极点配置到期望的极点位置上,从而使系统特性满足要求。 1. 极点配置原理 假设原系统的状态空间模型为: 若系统是完全可控的,则可引入状态反馈调节器,且: 这时,闭环系统的状态空间模型为: 系统方框图如下: 2. 极点配置的MATLAB函数 在MATLAB控制工具箱中,直接用于系统极点配置的函数有acker()和place()。调用格式为: K=acker(A,B,P) 用于单输入单输出系统 其中:A,B为系统矩阵,P为期望极点向量,K为反馈增益向量。 K=place(A,B,P) (K,prec,message)=place(A,B,P) place()用于单输入或多输入系统。Prec为实际极点偏离期望极点位置的误差;message是当系统某一非零极点偏离期望位置大于10%时给出的警告信息。 3. 极点配置步骤: (1)获得系统闭环的状态空间方程; (2)根据系统性能要求,确定系统期望极点分布P; (3)利用MATLAB极点配置设计函数求取系统反馈增益K; (4)检验系统性能。 极点配置实例分析 例1:已知控制系统的系数矩阵为: 闭环系统的极点为s=-1,-2,-3,对其进行极点配置。 解:用acker()函数对系统进行极点配置,给出程序ExamAcker1.m %极点配置 A=[-2.0 -2.5 -0.5;1 0 0;0 1 0]; B=[1;0;0]; P=[-1 -2 -3]; %极点配置 K=acker(A,B,P) %闭环系统矩阵 Ac=A-B*K %验证闭环系统牲特征值 eig(Ac) 程序运行结果如下: K = 4.0000 8.5000 5.5000 Ac = -6 -11 -6 1 0 0 0 1 0 ans = -3.0000 -2.0000 -1.0000 由运行结果可知,配置结果与题目要求相符,配置过程正确。 所以状态反馈控制器为: K=[4 8.5 5.5] 例2:已知控制系统系数矩阵为: 闭环系统的极点为s=-1±5i,-10,对其进行极点配置。 解:用place()函数对系统进行极点配置,给出程序ExamPlace.m %极点配置 A=[-0.1 5 0.1;-5 -0.1 5;0 0 -10];B=[0;0;10]; P=[-1-5*i,-1+5*i,-10]; %极点配置 K=place(A,B,P); %显示极点配置信息 [K,Prec,Mes]=place(A,B,P) 程序运行结果如下: K = -0.1404 0.3754 0.1800 Prec = 15 Mes = 由运行结果可知,配置过程中没有出错和警告信息。 所以状态反馈控制器为: K=[-0.1404 0.3754 0.18]. 例3: 已知控制系统的传递函数为 ,试判别系统的 可控性,并设计反馈控制器,使闭环系统极点为-2,-1±i。 解:首先判别系统的可控性,然后对可控系统进行极点配置。程序 ExamPlace2.m如 下: %判别系统的可控性 num=10;den=conv(conv([1,0],[1 1]),[1,2]); [a,b,c,d]=tf2ss(num,den);n=3; CAM=ctrb(a,b); if det(CAM)~=0; rcam=rank(CAM); if rcam==n disp(系统是可控的) p=[-2,-1+i,-1-i]; K=place(a,b,p); %极点配置 disp(系统反馈矩阵K,实际极点偏离期望极点的误差Prec,警告信息Mes各为:) [K,Prec,Mes]=place(a,b,p) sys=ss(a-b*K,b,c,d); poles=pole(sys); step(sys/dcgain(sys),10); elseif rcamn dis
文档评论(0)