利用MATLAB设计状态观测器—现代控制理论实验报告.docx

利用MATLAB设计状态观测器现代控制理论实验报告实验六利用MATLA般计状态观测器*******1121*****实验目的：1、学习观测器设计算法：2、通过编程、上机调试，掌握基于观测器的输出反馈控制系统设计方法实验原理：1、全阶观测器模型：x=Ax+Bu4-L(y—Cr)=(A—LC)x+Bu+Ly由极点配置和观测器设计问题的对偶关系，也可以应用MATLAB中极点配置的函数来确定所需要的观测器增益矩阵例如，对于单输入单输出系统，观测器的增益矩阵可以由函数L^ackcrCA*C,¥))'得到其中的V是由期望的观测器极点所构成的向量类似的，也可以用L=(place(A',C',V))'来确定一般系统的观测器矩阵，但这里要求V不包含相同的极点2、降阶观测器模型：•人AAvv=Aw+Bv+Fuxh=Ly基于降阶观测器的输出反馈控制器是：由二(】一户勺)访十住一户+u=-Khw-(Ka+KhL)y对于降阶观测器的设计，使用MATLAB软件中的函数L=(acker(Abb\Aab‘9V))'或L-(place(Abb\Aab',V))，可以得到观测器的增益矩阵Lo其中的V是由降阶观测器的期望极点所组成的向量。

实验要求1 .在运行以上例程序的基础上，考虑图6.3所示的调节器系统，试针对被控对象设计基于全阶观测器和降阶观测器的输出反馈控制器设极点配置部分希望的闭环极点是42二一2士；2x/3,希望的观测器极点是(a)对于全阶观测器，4二一8和人=一8；(b)对于降阶观测器，//=-8o比较系统对下列指定初始条件的响应：(a)对于全阶观测器：x,(0)=l,x2(0)=0,el(0)=l,e2(0)=0(b)对于降阶观测器：芭(0)=1,/(0)=0,巧(0)=1进一步比较两个系统的带宽U»图6.3调节器系统2 .假设SISO受控系统的开环传递函数为G(s)(1)若根据系统的性能指标要求，希望配置的系统极点为1 3, 2,32 j2L(s+2)得,L=[3.5;9.0]J=[-2+j*2*sqrt (3)2-j*2*sqrt (3)];求受控系统的状态反馈矩阵2)设计观测器反馈系数矩阵L,使全维状态观测器的希望极点均为-3.实验结果一、设计基于全阶观测器和降阶观测器的输出反馈控制器1、全阶观测器:1)计算全阶观测器的增益矩阵由图6.3所示的调节器系统1/s*执行以下的M-文件:a=[01;0-2];b=[0;1];c=[40]；v=[-8-8];l=(acker(a',c',v))'result:3.50009.0000计算得出，全阶观测器的增益矩阵相应的全阶观测器是：-????????????x(ALC)xBuLy=[999999]x+[??]u+「？？?]y-????-??????2)计算全阶观测器的状态反馈矩阵K及其状态响应A=[01;0-2];B=[0;1];K=place(A,B,J);sys=ss(A+B*K,eye(2),eye(2),eye(2));t=0:0.01:20;x=initial(sys,[1;0],t);x1=[10]*x';x2=[01]*x';subplot(2,1,1);plot(t,x1);gridxlabel('t(sec)');ylabel('xl');subplot(2,1,2);plot(t,x2);gridxlabel('t(sec)');ylabel('x2');可以得到：K=[-16-2]得到状态反馈后的状态变量的状态响应曲线:t(sec)t(sec)于是有，全阶观测器的输出反馈控制器为：?=(A-LC+BK)x+Ly=[—??????]x+[????]y-????-????2、降阶观测器1)计算降阶观测器的增益矩阵L执行以下的M-文件：Aaa=[0];Aab=[1];Aba=[0];Abb=[-2];Ba=[0]；Bb=[1]；v=[-8];l=(acker(Abb',Aab',v))'Ahat=Abb-l*AabBhat=Ahat*l+Aba-l*AaaFhat=Bb-l*Baresult:l=6Ahat=-8Bhat=-48Fhat=1计算得出，降阶观测器的增益矩阵L=[6]相应的降阶观测器是：wAwByFU=[-8]w+[-48]y+[1]u2)计算全阶观测器的状态反馈矩阵K及其状态响应Ahat=[-8];Bhat=[-48];Fhat=[1];K=[-2];sys=ss(Ahat-Fhat*K,eye(1),eye(1),eye(1));t=0:0.01:20;w=initial(sys,[-1],t);plot(t,w);xlabel('t(sec)');ylabel('w');得到状态响应曲线如下图所示：CJi।i।।i।।।- 0.1=J=- 021□3--- 0.4■■-豆-0.5■--- 0,0---□7- □8--- 0.9--hilltill■O246a10112141S1020t(seeJi二、开环传递函数1)求受控系统的状态反馈矩阵执行以下应用函数acker编制的M-文件：A=[010;001;000];B=[0；0；1]；J=[-3-2+j*22-j*2];K=acker(A,B,J)得到K=[-24-20-7]2)设计观测器反馈系数矩阵La=[010;001;000];b=[0;0;1];c=[100];v=[-3-3-3];l=(acker(a',c',v))'l=92727??计算得出，全维观测器的增益矩阵L=[????????。