状态观测器的设计

利用MATLAB 设计状态观测器本节将介绍用MATLAB 设计状态观测器的若干例子。

我们将举例说明全维状态观测器和最小阶状态观测器设计的MATLAB 方法。

------------------------------------------------[例1] 考虑一个调节器系统的设计。

给定线性定常系统为Cxy Bu Ax x =+=& 式中]01[,10,06.2010=⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡=C B A且闭环极点为)2,1(==i s i μ，其中4.28.1,4.28.121j j −−=+−=μμ期望用观测-状态反馈控制，而不是用真实的状态反馈控制。

观测器的期望特征值为821−==μμ试采用MATLAB 确定出相应的状态反馈增益矩阵K 和观测器增益矩阵e K 。

[解]对于题中给定的系统，可利用如下MATLAB Program 1来确定状态反馈增益矩阵K和观测器增益K。

矩阵eMATLAB Program 1% Pole placement and design of observer ------% ***** Design of a control system using pole-placement% technique and state observer. Solve pole-placement% problem *****% ***** Enter matrices A,B,C,and D *****A=[0 1;20.6 0];B=[0;1]C=[1 0];D=[0];% ***** Check the rank of the controllability matrix Q *****Q=[B A*B];Rank(Q)ans=2% ***** Since the rank of the controllability matrix Q is 2, % arbitrary pole placement is possible *****% ***** Enter the desired characteristic polynomial by% defining the following matrix J and computingpoly(J) *****J=[-1.8+2.4*i 0;0 -1.8-2.4*i];Poly(J)ans=1.000 3.6000 9.0000% ***** Enter characteristic polynomial Phi *****Phi=polyvalm(poly(J),A);% ***** State feedback gain matrix K can be given by ***** K=[0 1]*inv(Q)*PhiK=29.6000 3.6000% ***** The following program determines the observer matrix Ke *****% ***** Enter the observability matrix RT and check its rank *****RT=[C’ A’*C’];rank(RT)ans=2% ***** Since the rank of the observability matrix is 2, design of% the observer is possible *****% **** Enter the desired characteristic polynomial by defining % the following matrix J0 and entering statement poly(JO) *****JO=[-8 0;0 -8];Poly(JO)ans=1 16 64% ***** Enter characteristic polynomial Ph ***** Ph=polyvalm(ply(JO),A);% ***** The observer gain matrix Ke is obtained from ***** Ke=Ph*(inv(RT’))*[0;1]Ke=16.000084.60000求出的状态反馈增益矩阵K 为[]6.36.29=K观测器增益矩阵e K 为⎥⎦⎤⎢⎣⎡=6.8416e K 该观测-状态反馈控制系统是4阶的，其特征方程为0=+−+−C K A sI BKA sI e通过将期望的闭环极点和期望的观测器极点代入上式，可得5764.3746.1306.19)8)(4.28.1)(4.28.1(2342++++=+++−+=+−+−s s s s s j s j s C K A sI BK A sI e这个结果很容易通过MATLAB得到，如MATLAB Program 2所示(MATLAB Program 2是K已MATLAB Program1的继续。

• 考察
T & u x = Ax+B z = A z +CTυ & 的对偶系统 n = BT z y =Cx
(4)Biblioteka 且定义 v=r-Kz • 则
T z =(A −CTk)z +CTr & n = BT z
(5)
• 注意到 AT −CTk =(A−kTC)T
(A−kTC)T 和 A−kTC • 而
(4).带观测器的状态反馈系统 带观测器的状态反馈系统. 带观测器的状态反馈系统 • 在状态反馈中 不采样原系统的状态进行反 在状态反馈中,不采样原系统的状态进行反 馈而采用状态观测器估计的状态进行反馈, 馈而采用状态观测器估计的状态进行反馈 其结构图如下图所示. 其结构图如下图所示
• 状态估计器
2.极点配置条件 极点配置条件 • 若被控系统 Σ0(A, B) 是状态完全能控的,那么 是状态完全能控的 那么 反馈系统的极点必是可以任意配置的,或者 反馈系统的极点必是可以任意配置的 或者 说,能使闭环系统极点任意配置的条件是被 能使闭环系统极点任意配置的条件是被 控系统完全可控. 控系统完全可控
极点配置(仅讨论单输入 单输出系统) 二.极点配置 仅讨论单输入 单输出系统 极点配置 仅讨论单输入/单输出系统 1.什么是极点配置 什么是极点配置. 什么是极点配置 A− • 如果 Σk[(A−Bk), B,C]的全部 个)极点可以通过 的全部(n个 极点可以通过 选择状态反馈矩阵K的各元素而移至 的各元素而移至S平面 选择状态反馈矩阵 的各元素而移至 平面 任意指定的位置,称该系统是极点可任意配 任意指定的位置 称该系统是极点可任意配 置的。 置的。
• 由(1)和(2)得 和 得

实验报告课程践性系统理抡基碣实验日期2016年6月6日 专业册级学号同组人 实验名称全维狀盗现測器的投廿坪分 批阅敎斷签字一、实鲨目的1. 学习用状态观测器获取系统状态IS it li ft 方法，了解全细状态观测器的 根点对状态的估it 误差的靈响；2. 拿捋全维状态观測器的设it 方法;3. 拿捋带有狀态规測器的状态反饋系躱设廿方法。

二、实验容a ） 用状态反馈配置系説的用环根点：一2±丿・2朽,一5；b ） 设廿全细狀态规測器，规测器的极点为：一5±/2巧,一10；0研究规测器tUE 置对估计状态逼近被估it 值的影响； d ）求系统的传递函数（带规测器及不帑观测器时）;绘制系躱的输出阶跃响应曲线。

三. 实螫环境MATLAB6.5・01 o'oA =0 0 1,b = 0-6-1161c = p 0 0]开环系统x = Ax + bu y = ex0.实验眾理（或程序框图）及步驟利用状态反饋可以便闭环系就的根点配置在所希里的位置上，以条件是必须对全部状态变量邵能iSliH量，但在实际系统中,并不是所有狀盗变量那能测量的，这就给狀态反锁的实现it成了困难。

SUtl设法利用已*11的信息（输岀量y和输人量x）,通过一个模塑車新构适系貌状态以对状态变量进fiIfiito该模塑就枕为狀态规測器。

若状态观测器的阶次与系貌的阶次是曲同的，这样的状态规汹器就祢力全维狀态观测器或全阶观測器。

设系説完全可观，剧可fiiia图4"所示的狀态规测器图4-1全维状态範測器为求岀状态焜测器的反備ke熾益，与板点配置类做，也可有两种方法:方法一：构造变换矩阵Q，便系统变成标旌能规里，然后根齬特征方程求HI ke；方法二：是可采用Ackermann公式：& =①(A)Q：［o 0・・・0 if, 其中0。

力可规性矩眸。

利用对偶原理，可便设廿冋题大为简化。

首先构造对偶系貌然后可由变换沫或Ackermann公式求岀様点配置的反馈k增益，逹也可由MATLAB的place和acker函数得到;晟后求出狀态规测器的反饋增益。

现代控制实验--状态反馈器和状态观测器的设计-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN状态反馈器和状态观测器的设计一、实验设备PC 计算机，MATLAB 软件，控制理论实验台，示波器二、实验目的（1）学习闭环系统极点配置定理及算法，学习全维状态观测器设计法；（2）掌握用极点配置的方法（3）掌握状态观测器设计方法（4）学会使用MATLAB工具进行初步的控制系统设计三、实验原理及相关知识（1）设系统的模型如式所示若系统可控，则必可用状态反馈的方法进行极点配置来改变系统性能。

引入状态反馈后系统模型如下式所示：（2）所给系统可观，则系统存在状态观测器四、实验内容（1）某系统状态方程如下10100134326x x u •⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦[]100y x =理想闭环系统的极点为[]123---.（1）采用 Ackermann 公式计算法进行闭环系统极点配置；代码：A=[0 1 0;0 0 1;-4 -3 -2];B=[1; 3; -6];P=[-1 -2 -3];K=acker(A,B,P)Ac=A-B*Keig(Ac)（2）采用调用 place 函数法进行闭环系统极点配置；代码：A=[0 1 0;0 0 1;-4 -3 -2];B=[1;3;-6];eig(A)'P=[-1 -2 -3];K=place(A,B,P)eig(A-B*K)'（3）设计全维状态观测器，要求状态观测器的极点为[]---123代码：a=[0 1 0;0 0 1;-4 -3 -2];b=[1;3;-6];c=[1 0 0];p=[-1 -2 -3];a1=a';b1=c';c1=b';K=acker(a1,b1,p);h=(K)'ahc=a-h*c（2）已知系统状态方程为：10100134326x x u •⎡⎤⎡⎤⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥----⎣⎦⎣⎦[]100y x =（1）求状态反馈增益阵K ，使反馈后闭环特征值为[-1 -2 -3]；代码：A=[0 1 0;0 0 1;4 -3 -2];b=[1;3;-6];p=[-1 -2 -3];k=acker(A,b,p)A-b*keig(A-b*k)（2）检验引入状态反馈后的特征值与希望极点是否一致。