全维状态观测器的设计
全维状态观测器例
所示。
u
x2
1 x2 s
3
1 x1 s
2
2
y
受控对象 部分
23.5
8.5
ˆ2 1 x s
3
1 s
ˆ1 x
2
ˆ y
状态观测器部分
2
图 9.24 例 9.47 全维状态观测器及受控对象状态变量图
Байду номын сангаас
(2) 给定极点对应的期望特征方程
( 10) 2 2 20 100 0
(3) 设观测器输出反馈阵 H h1
(9.203)
h2 ,则观测器系统矩阵为
T
1 h1 2h1 0 A HC 2 0 2 3 h2 2 2h2
观测器特征方程
1 3
(9.204)
I ( A HC) 2 (2h1 3) (6h1 2h2 2) 0
比较式(9.203)和(9.204)各对应项系数,有
2h1 3 20 6h1 2h2 2 100 h1 8.5, h2 23.5
全维状态观测器例
【例 9.47】 已知受控对象传递函数
Y ( s) 2 U ( s) ( s 1)( s 2)
试设计全维状态观测器,将观测器极点配置在-10,-10。 解: (1) 可观性判别 传递函数无零极点对消,故系统可控、可观测。若写出可控标准形实现,则
1 0 0 , b , C 2 0 A 2 3 1
实验 6 极点配置与全维状态观测器的设计
实验 6 极点配置与全维状态观测器的设计一、实验目的1. 加深对状态反馈作用的理解。
2. 学习和掌握状态观测器的设计方法。
二、实验原理在MATLAB 中,可以使用acker 和place 函数来进行极点配置,函数的使用方法如下:K = acker(A,B,P) A,B为系统系数矩阵,P为配置极点,K为反馈增益矩阵。
K = place(A,B,P) A,B为系统系数矩阵,P为配置极点,K为反馈增益矩阵。
[K,PREC,MESSAGE] = place(A,B,P) A,B为系统系数矩阵,P为配置极点,K为反馈增益矩阵,PREC 为特征值,MESSAGE 为配置中的出错信息。
三、实验内容1.已知系统(1)判断系统稳定性,说明原因。
(2)若不稳定,进行极点配置,期望极点:-1,-2,-3,求出状态反馈矩阵k。
(3)讨论状态反馈与输出反馈的关系,说明状态反馈为何能进行极点配置?(4)使用状态反馈进行零极点配置的前提条件是什么?1.(1)(2)代码:a=[-2 -1 1;1 0 1;-1 0 1];b=[1,1,1]';p=[-1,-2,-3]';K=acker(a,b,p)K =-1 2 4(3)讨论状态反馈与输出反馈的关系, 说明状态反馈为何能进行极点配置?在经典控制理论中,一般只考虑由系统的输出变量来构成反馈律,即输出反馈。
在现代控制理论的状态空间分析方法中,多考虑采用状态变量来构成反馈律,即状态反馈。
从状态空间模型输出方程可以看出,输出反馈可视为状态反馈的一个特例。
状态反馈可以提供更多的补偿信息,只要状态进行简单的计算再反馈,就可以获得优良的控制性能。
(4)使用状态反馈配置极点的前提是系统的状态是完全可控的。
2.已知系统设计全维状态观测器,使观测器的极点配置在12+j,12-j 。
(1)给出原系统的状态曲线。
(2)给出观测器的状态曲线并加以对比。
(观测器的初始状态可以任意选取)观察实验结果,思考以下问题:(1)说明反馈控制闭环期望极点和观测器极点的选取原则。
现代控制理论习题之状态观测设计
对应于原系统的观测器矩阵: ⎡0 ⎤ ⎡0 ⎤ P1 = V0 −1 ⎢ ⎥ = ⎢ ⎥, Po = [ p1 ⎣1 ⎦ ⎣1 ⎦
u
∑ ( A, B, C )
y
6.5
2
1 x
x1
15.3 x
x3
题 6-2 图 1
(2) 确定降维观测器的维数:m=1,n=3,则 n-m= 2。 分解输出系数矩阵 c,获得线性变换矩阵 T,对原状态空间表达式进行线性变换,使 各输出变量 y 变成各状态变量的单值函数:
f *(s) = (s + 3)(s + 4) = s2 + 7s +12 ⎡s 0⎤ ⎡−1 −1⎤ ⎡l1⎤ f (s) = sI − (A22 − LA 12) = ⎢ ⎥ −⎢ ⎥ + ⎢ ⎥[− 2 − 4] ⎣0 s⎦ ⎣−1 −1⎦ ⎢ ⎣l2⎥ ⎦ = s2 + (−4l2 − 2l1 + 2)s + (2l1 − 2l2) ⎡l ⎤ ⎡ 3.1667⎤ f *(s) = f (s) ⇒ L = ⎢ 1⎥ = ⎢ ⎥ ⎢l2⎦ ⎥ ⎣− 2.8333 ⎦ ⎣
系统能观,可设计观测器。 求希望特征多项式:
f * ( s ) = ( s + 3)( s + 4)( s + 5) = s 3 + 12 s 2 + 47 s + 60
求观测器特征多项式:
f ( s ) = sI − A + LC
5.5状态观测器设计
N =B, K =E
于是得到一特定的n 维KX 观测器。
⎧ z& = ( A − LC )z + Bu + Ly
⎨ ⎩W
=
Kz
(5-35)
称此为∑ 的一个全维KX观测器;K=I为 ∑ 的一个全维状态观测
器.
因为满足结构条件的L 不唯一,全维观测器也不唯一。全6 维
观测器设计较简单。
5.5.2 全维状态观测器设计
⎨ ⎩
W = Ez
(5-34)
称为全维观测器。
若 r < n ,M ≠ 0 ,相应观测器称为降维观测器。
对 r = n 全维观测器,参数除按通常步骤外,有特定取法:
F = A − LC ,
G=L
则 PA− FP = PA− (A − LC)P = PA− AP + LCP = LC
有 P = In 从而
y = [1 1 0 ] x
设计特征值为-3,-3和-4的全维状态观测器.
解:
⎡c⎤
⎡1 1 0⎤
rank
⎢ ⎢
cA
⎥ ⎥
=
rank
⎢ ⎢
−1
−3
−1⎥⎥ = 3 = n
8
⎢⎣cA2 ⎥⎦
⎢⎣ 0 5 0 ⎥⎦
可知系统完全观测.
⎡−1 0 1 ⎤
⎡1⎤
A
=
AT
=
⎢ ⎢
−
2
−1
0
⎥ ⎥
,
b
= cT
=
完全能控
⎡ 1 0 0⎤ ⎡0 −1 1⎤ ⎡1 0 0⎤ ⎡ 2 2 1⎤
P = ⎡⎣ A2b Ab b ⎤⎦ ⎢⎢α 2
1
利用MATLAB设计状态观测器
x%
+
⎢ ⎢
0
⎥ ⎥
u
+
⎢⎢147.3875⎥⎥
y
⎜⎝ ⎢⎣1.244 0.3965 −3.145⎥⎦ ⎢⎣544.3932⎥⎦
⎟⎠ ⎢⎣1.244⎥⎦ ⎢⎣544.3932⎥⎦
⎡ −16.855 = ⎢⎢−147.3875
⎣⎢−543.1492
1 0 0.3965
0 ⎤ ⎡ 0 ⎤ ⎡ 16.855 ⎤
统(5.9)具有任意给定的极点,这样的矩阵 L 可以应用全阶观测器的设计方法来设计。矩
阵 L 也称为是系统的降阶观测器增益矩阵。
对于降阶观测器的设计,使用 MATLAB 软件中的函数
L=(acker(Abb’,Aab’,V))’ 或
L=(place(Abb’,Aab’,V))’
可以得到观测器的增益矩阵 L 。其中的 V 是由降阶观测器的期望极点所组成的向量。
确定所需要的观测器增益矩阵。例如,对于单输入单输出系统,观测器的增益矩阵可以由函
数
L=(acker(A’,C’,V))’ 得到。其中的 V 是由期望的观测器极点所构成的向量。类似的,也可以用
L=(place(A’,C’,V))’ 来确定一般系统的观测器矩阵,但这里要求 V 不包含相同的极点。
5.3.2 降阶观测器设计
实验 5 利用 MATLAB 设计状态观测器
5.1 实验设备
同实验 1。
5.2 实验目的
1、学习观测器设计算法;
2、通过编程、上机调试,掌握基于观测器的输出反馈控制系统设计方法。
5.3 实验原理说明
5.3.1 全阶观测器设计
考虑如下的线性时不变系统
⎧x& = Ax + Bu
现代控制实验--状态反馈器和状态观测器的设计
现代控制实验--状态反馈器和状态观测器的设计-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)检验引入状态反馈后的特征值与希望极点是否一致。
全维状态观测器的设计
全维状态观测器的设计全维状态观测器(Full State Observer)是一种常用于控制系统中的重要部件,用于获取系统的全部状态信息。
它通常是通过对系统的输入输出进行观测,并通过数学模型来估算系统的状态。
全维状态观测器的设计可以通过以下步骤来完成。
第一步是系统建模。
将所要观测的系统建立数学模型,可以采用物理方程或者数学模型的方式。
常见的数学模型包括状态方程和输出方程。
状态方程描述了系统状态的时间演变规律,输出方程则描述了系统输出与状态之间的关系。
这些方程可以通过系统的运动方程,控制方程和物理特性等来建立。
第二步是选择观测器类型。
全维状态观测器有多种类型,包括基本观测器、极点配置观测器和最优观测器等。
基本观测器是使用系统的状态方程和输出方程来估算系统状态的观测器,而极点配置观测器和最优观测器则是通过最小化误差来估算系统状态,从而提高观测器的精度。
合适的观测器类型应该根据控制系统的需求来选择。
第三步是计算观测器矩阵。
观测器矩阵是观测器中用来计算系统状态的矩阵。
它可以使用系统的状态方程和输出方程来计算。
观测器矩阵需要满足一些性质,例如它需要是可观测的,并且需要保证系统状态与观测器状态的误差最小。
第五步是实现观测器。
实现观测器需要将观测器矩阵和观测器增益输入到观测器中,并对观测器的输入输出进行校验。
一旦观测器被设计并实现,它就可以用于控制系统中,并用来估算系统的全部状态信息。
总之,全维状态观测器的设计是控制系统中的重要部件,可以极大地提高控制系统的精度和稳定性。
设计一个好的全维状态观测器需要仔细分析系统模型和观测器类型,计算观测器矩阵和观测器增益,并进行实现和调试。
全维状态观测器的设计
卖验报告镌程 銭性*统理论基础 令敲可朗2016年g 刀互目 *业去飲 _______ 学号 _耐俎人 _________ 丈 检名 怂 全* 状 态 观 删 星 的 设计 评分 ______馳闻敎艸签丰 ___________一、卖軽可的1. 学习用状态观测森灰取糸统状态估计值的方出,了鮮全维状态观测森的 极点对状态的估计谖臭的彩响;2. 拿握全维状态观测森的设计方法;3. 拿提带有状态观测森的状态反馈糸统设计方出。
x = Ax + bu开环糸统彳,其中y = exa)用状态反馈配.更糸疣的诃环极点:一2土丿PJM —5 ;b) 设计全维状态观测爰.叽测乐的圾点为:一5±)2丿亍,一10; c) 研兗观测乐圾点住逍对仕计状态遍近彼仕计值的彩响;・01 0_OA =0 1 ,b =-6-11 61c = p 0 0]绘制糸统的输出阶跃响应曲线。
三.实脍环境MATLAB6.5宾唸々理图丿及歩螺利用状态及馈可以使诃环糸疣的圾盍配逍疫所希玺的住置上,其条件是必须对全部状态变量都能进行测量,仅A卖际糸统中.并不是所冇状态变量却能测量的,这就给状态反馈的卖现岌成了困难。
因此要设法利用己知的传息(綸出量y和綸入量x)>通过一个栈空重新构逡糸统状态以对状态变量进行仕计。
该楼型就称为状态观测森。
若状态观测森的阶次与务统的阶次是相同的,这擇的状态观测容就称为全维状态观测森戎全阶观测怎。
设糸疣兜全可观.則可构凌如图4・1所承的状态观測余国4-1 全维状态观測容为求出状态观測爰的反馈ke增益,与极点紀.1.类很.也可有而种方法:为比一:构逡雯换矩阵Q,使糸统支成标准能观淫.然后根据持征方程求岀匕;•资料. ^方出二:是可采用Ackermann 公 < : k c = [O 0…0 if,其中0。
为可观性矩阵。
利用对偶虑理,可使沒计问題大为简化。
管丸构凌对偶糸统'^=A T^+ C T Vhg然后可由变换法或Ackermann小式求出极点紀.1.的反馈k增益,这也可由MATLAB的place和acker函散得刊;呆后求出状态观测乐的反馈增益。
观测器设计
3
§6.1 观测器的定义
图6.1
4
§6.1 观测器的定义
观测器的形式: 对于线性定常系统 (6.1.1) x 为 n 维状态变量,u 为 p 维输入变量,y 为 q 维 A 输出变量, , B, C 分别为 n × n, n × p, q × n 实常阵。 其观测器也是一个线性系统, 状态空间描述一般可表示为
§6.1 观测器的定义
状态反馈可改变系统的性能, 从而使之有很好的优越性。 实际中由于技术或经济上的原因, 系统的全部状态不能 够直接获取. 如何来解决状态反馈性能上的优越性和物理上的不可实 现性所形成的矛盾呢? 途径之一就是通过重构系统的状态, 并且利用这个重构 状态代替系统的真实状态来实现所要求的状态反馈, 此 即状态重构问题, 也称观测器设计问题. 观测器设计问题不仅有理论意义, 而且有应用价值。
由状态反馈不改变系统的零点知, 闭环传函应为 s −1 GC ( s ) = ( s + 1)( s − 1)
23
§6.2-2 基于全维状态观测器的输出动态反馈特性
s 2 − K 2 s − K1 = ( s + 1)( s − 1) = s 2 − 1
K 2 = 0, K1 = 1, K = [1 0]
0⎤ In ⎥ ⎦
则有
16
§6.2-2 基于全维状态观测器的输出动态反馈特性
BK ⎡ A ⎤ T ⎢ ⎥T ⎣GC A − GC + BK ⎦ 0⎤⎡ A BK ⎡ In ⎤ ⎡ In =⎢ − I n I n ⎥ ⎢GC A − GC + BK ⎥ ⎢ I n ⎦⎣ ⎣ ⎦⎣
−1
0⎤ In ⎥ ⎦
n
20
§6.2-2 基于全维状态观测器的输出动态反馈特性 系统(6.2.7)直接状态反馈 u = Kx + v 的闭环传递函数为 −1 (6.2.14) Gk ( s) = C ( sI n − A − BK ) B Gc ( s ) = Gk ( s ). 结论:全维状态观测器的引入, 并不改变原系统直接状态 反馈的闭环传递函数。 s −1 例1 设开环系统的传递函数Go ( s ) = 2 , 试作全维状态 s 观测器,并构造输出动态反馈,使观测器极点 −10 ± j , 1 . GC ( s ) = 闭环传函为 s +1
现代控制理论状态反馈和状态观测器的设计实验报告
本科实验报告课程名称:现代控制理论实验项目:状态反馈和状态观测器的设计实验地点:中区机房专业班级:自动化学号:学生:指导教师:年月日现代控制理论基础一、实验目的(1)熟悉和掌握极点配置的原理。
(2)熟悉和掌握观测器设计的原理。
(3)通过实验验证理论的正确性。
(4)分析仿真结果和理论计算的结果。
二、实验要求(1)根据所给被控系统和性能指标要求设计状态反馈阵K。
(2)根据所给被控系统和性能指标要求设计状态观测器阵L。
(3)在计算机上进行分布仿真。
(4)如果结果不能满足要求,分析原因并重复上述步骤。
三、实验容(一)、状态反馈状态反馈是将系统的状态变量乘以相应的反馈系数,然后反馈到输入端与参考输入叠加形成控制作为受控系统的控制输入,采用状态反馈不但可以实现闭环系统的极点任意配置,而且也是实现解耦和构成线性最优调节器的主要手段。
1.全部极点配置给定控制系统的状态空间模型,则经常希望引入某种控制器,使得该系统的闭环极点移动到某个指定位置,因为在很多情况下系统的极点位置会决定系统的动态性能。
假设系统的状态空间表达式为(1)其中 n m C r n B n n A ⨯⨯⨯::;:;: 引入状态反馈,使进入该系统的信号为Kx r u -=(2)式中r 为系统的外部参考输入,K 为n n ⨯矩阵. 可得状态反馈闭环系统的状态空间表达式为(3)可以证明,若给定系统是完全能控的,则可以通过状态反馈实现系统的闭环极点进行任意配置。
假定单变量系统的n 个希望极点为λ1,λ2, …λn, 则可以求出期望的闭环特征方程为=)(*s f (s-λ1)(s-λ2)…(s-λn)=n n n a s a s +++- 11这是状态反馈阵K 可根据下式求得K=[])(100*1A f U c -(4)式中[]b A Ab b U n c 1-= ,)(*A f 是将系统期望的闭环特征方程式中的s 换成系统矩阵A 后的矩阵多项式。
例1已知系统的状态方程为u x x ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---=•111101101112 采用状态反馈,将系统的极点配置到-1,-2,-3,求状态反馈阵K..其实,在MATLAB的控制系统工具箱中就提供了单变量系统极点配置函数acker(),该函数的调用格式为K=acker(A,b,p)式中,p为给定的极点,K为状态反馈阵。
状态重构与状态观测器的设计_5.4_5.5
解 根据给定的受控系统, 根据给定的受控系统,求得能观测性矩阵及能控性矩阵的 秩为
C 1 1 0 rank CA = rank 1 2 1 = 3 2 CA 1 4 4
1 1 1 A2 B = rank 0 1 4 = 3 1 2 4
ˆ x = e ( A−GC )t x0 , x0 = x0 − x0 , t0 = 0
5.4 状态重构与状态观测器的设计
3. 全维状态观测器极点任意配置条件 定理5-4 可用图 所示的结构,设计全维状态观测 可用图5-6所示的结构 所示的结构, 定理
重构出系统所有的状态, 器,重构出系统所有的状态 ,并且观测器的极点可 以任意配置的充分必要条件是系统完全能观测。 以任意配置的充分必要条件是系统完全能观测。
5.4 状态重构与状态观测器的设计
4. 设计反馈矩阵 设计反馈矩阵G (1)按照极点配置的方法 )
(2)极点选取:若是选得离虚轴愈远,状态误差趋 )极点选取:若是选得离虚轴愈远, 于零的速度就愈快。 于零的速度就愈快。过于远离虚轴则状态观测器的 频带过宽,将降低状态观测器抗高频干扰的性能。 频带过宽,将降低状态观测器抗高频干扰的性能。
,上式变为: 上式变为: 上式变为
ɺ x = ( A − BK ) x +BKx +Bv
(4)同时 观测器的状态误差方程 )同时,
ɺ x = ( A − GC )x
(5) 上两式= 0 ɺ
(6)特征方程 )
λ I − A + BK
0
BK x B x + 0 v A- GC
全维、降维观测器
本文通过具体的例子阐明如何在 MATLAB系统中进行全维状态观测器和降维状态观测器的设计。
MATLAB 为状态空间设计提供了很多有用的函数,方便了矩阵方程的求解,其中的MATLAB里面提供的库函数对全维状态观测器和降维状态观测器的设计也显得非常地方便。
现通过例子说明如何用 MATLAB 设计状态观测器。
为了评价 MATLAB 所设计的状态观测器的性能,本文通过在 SIMULINK 环境下来仿真一个三阶状态观测器,来说明用 MATLAB 设计状态观测器的准确性。
1、全维观测器的设计已知三阶系统的状态空间方程为:首先检验系统的是否完全能观A=[2 -1 1;0 -1 2;1 0 -2];C=[2 1 0];N=[C;C*A;C*A*A]rank(N)rank(N)ans = 3 ,说明系统是完全能观的。
下面选择观测器需要配置的期望极点为:s1 =-12 s2,3 =-3±0.88i 由此求出观测器增益矩阵G:A=[2 -1 1;0 -1 2;1 0 -2];C=[2 1 0];P =[-12;-3+0.88*i;-3-0.88*i];G = acker(A',C',P);求得G = [11.6527 -6.3054 1.0619]可得全维观测器的方程为:下面可依据上式构建simulink图,据此观察观测器的跟踪能力:跟踪效果图如下:X1X2X3从图中可以看出状态观测器的状态X1,X2,X3 能够完全估计原三阶系统的状态。
如果原三阶系统和其状态观测器的初始状态不同,状态观测器的状态X1,X2,X3 不能同原三阶系统状态完全一致,但能很快跟踪原三阶系统状态(主要取决于状态观测器的响应速度,即状态观测器的极点配置。
)2、降维观测器设计从上面的全维观测器输出方程可以看出,此系统输出就等于第一个状态,因此该状态可以有输出求得,即变换矩阵P为单位阵,而降维观测器的阶次为2。
降维观测器的期望特征根选为-3±0.88i据此求观测器增益a22=[-1 1;0 -1];a12=[-2 -2];pe=[-3+1i*2*7^(1/2)/3;-3-1i*2*7^(1/2)/3];lt=acker(a22',a12',pe);l=lt'求得,得到引入中间变量得降维观测器的状态方程为下面可依据上式构建simulink图,据此观察观测器的跟踪能力X2X3从降维观测器仿真的输出图形中可以看出,系统状态变量X2、X3的跟踪能力比较理想。
现代控制理论第五章讲义1
对于q维输出系统,有q个输出变量可直接由 传感器测得,若选取该q个输出作为状态变 量,它们便无需由观测器作出估计,观测器 只需估计(n-q)个状态变量,称为降维观 测器。它是(n-q)维子系统,结构简单, 工程上易于实现。为此,需要由受控对象动 态方程导出(n-q)维子系统动态方程,建 立降维观测器的观测模型。
g1 8.5 3 2 g1 20 2 4 g1 2 g 2 100 g 2 32
状态观测器为
g1 G g2
ˆ ˆ x [ A GC ]x bu Gy ˆ ˆ y Cx
5.5 状态观测器的设计
四、降维观测器
第六节 状态观测器实现状态反馈
在前面几节中,我们讲述了利用状态观测器 解决受控系统的维数重构问题从而使得状态 反馈系统得以实现,本节主要讨论利用观测 器进行状态估值反馈的系统与状态直接反馈 的系统之间的区别。
5.6 利用状态观测器实现状态反馈 一、系统结构与状态空间表达式
在一个带有全 维状态观测器 的状态反馈系 统中,设能控 能观的受控系 统∑0=(A、B、 C)为
* g 0 a 0 a0 * g1 a1 a1 g a* a ˆn n 1 n 1
5.5 状态观测器的设计
例、已知系统 1 1 0 x x 1u 0 2 y 2 0x 试设计一个状态观测器 ,使其极点为- , 10。 10
1
sI A HC
1
1
B 0
C 0
sI A BK B
1
0
C sI A BK B
状态观测器设计
状态观测器设计利用状态反馈实现闭环系统的极点配置,需要利用系统的全部状态变量。
然而系统的状态变量并不都是能够易于用物理方法量测出来的,有些根本就无法量测;甚至一些中间变量根本就没有常规的物理意义。
此种情况下要在工程上实现状态反馈,就需要对系统的状态进行估计,即构造状态观测器。
状态观测器,是一个在物理上可以实现的动态系统,它利用待观测系统的可以量测得到的输入和输出信息来估计待观测系统的状态变量,以便用该组状态变量的估计值来代替待观测系统的真实状态变量进行状态反馈设计,实现闭环系统极点的再配置。
1. 全维状态观测器当对象的所有状态均不可直接量测时,若要进行状态反馈设计,就需对全部状态变量进行观测。
这时构造的状态观测器,其阶次与对象的阶次相同,被称为全维状态观测器。
考虑如下n阶单输出线性定常离散系统(1)其中,A为n×n维系统矩阵,B为n×r输入矩阵,C为n×1维输出矩阵。
系统结构图如图1所示。
图1 全维状态观测器构造一个与受控系统具有相同参数的动态系统(2)当系统(1)与(2)的初始状态完全一致时,则两个系统未来任意时刻的状态也应完全相同。
但在实际实现时,不可能保证二者初始状态完全相同。
为此,应引入两个系统状态误差反馈信号构成状态误差闭环系统,通过极点配置使误差系统的状态渐趋于零。
由于原受控系统状态不可直接量测,故用二个系统的输出误差信号代替。
引入了输出误差的状态观测器状态方程为(3)其中,H为状态观测器的输出误差反馈系数矩阵,有如下形式定义状态估计误差为,用式(7.65)与(7.67)相减可得(4)即(5)通过式(5)可以看出,若选择合适的输出误差反馈矩阵H 使得状态估计误差系统(5)的所有极点均位于z平面单位圆内,则误差可在有限拍内趋于零,即状态估计值在有限拍内可以跟踪上真实状态,且极点越靠近原点状态估计误差趋于零的速度越快,反之越慢。
可见,能否逼近x(k)以及逼近速度是由H阵决定的。
全维状态观测器的设计
实 验 报 告课程 线性系统理论基础 实验日期 2016年 6月 6 日 专业班级 姓名 学号 同组人 实验名称全维状态观测器的设计评分批阅教师签字一、实验目的1. 学习用状态观测器获取系统状态估计值的方法,了解全维状态观测器的极点对状态的估计误差的影响;2. 掌握全维状态观测器的设计方法;3. 掌握带有状态观测器的状态反馈系统设计方法。
二、实验内容开环系统⎩⎨⎧=+=cxy bu Ax x,其中[]0100001,0,10061161A b c ⎡⎤⎡⎤⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦a) 用状态反馈配置系统的闭环极点:5,322-±-j ; b) 设计全维状态观测器,观测器的极点为:10,325-±-j ; c) 研究观测器极点位置对估计状态逼近被估计值的影响; d) 求系统的传递函数(带观测器及不带观测器时);绘制系统的输出阶跃响应曲线。
三、实验环境 MATLAB6.5四、实验原理(或程序框图)及步骤利用状态反馈可以使闭环系统的极点配置在所希望的位置上,其条件是必须对全部状态变量都能进行测量,但在实际系统中,并不是所有状态变量都能测量的,这就给状态反馈的实现造成了困难。
因此要设法利用已知的信息(输出量y 和输入量x),通过一个模型重新构造系统状态以对状态变量进行估计。
该模型就称为状态观测器。
若状态观测器的阶次与系统的阶次是相同的,这样的状态观测器就称为全维状态观测器或全阶观测器。
设系统完全可观,则可构造如图4-1所示的状态观测器图4-1 全维状态观测器为求出状态观测器的反馈ke 增益,与极点配置类似,也可有两种方法: 方法一:构造变换矩阵Q ,使系统变成标准能观型,然后根据特征方程求出k e ;方法二:是可采用Ackermann 公式:[]Toe Q A k 1000)(1-Φ=,其中O Q 为可观性矩阵。
利用对偶原理,可使设计问题大为简化。
首先构造对偶系统⎩⎨⎧=+=ξηξξTT T b v c A 然后可由变换法或Ackermann 公式求出极点配置的反馈k 增益,这也可由MATLAB 的place 和acker 函数得到;最后求出状态观测器的反馈增益。
离散控制系统中的状态观测器设计
离散控制系统中的状态观测器设计离散控制系统是指系统的输入和输出是离散的,并且在时间上以离散的方式进行测量和控制。
状态观测器是离散控制系统中重要的组成部分,用于估计系统的状态变量,从而实现对系统的控制。
本文将介绍离散控制系统中状态观测器的设计方法及其应用。
一、状态观测器的概念和作用状态观测器是一种用于估计系统状态的装置或算法。
在离散控制系统中,通过观测系统的输出值和输入值,结合系统的数学模型,状态观测器能够推断出系统的状态变量,从而实现对系统的监测和控制。
状态观测器在离散控制系统中具有重要的作用。
首先,通过对系统状态的估计,可以实现对系统的运行状态的实时监测,减少故障的发生。
其次,状态观测器可以提供系统未知状态变量的估计值,从而实现对系统的控制。
因此,状态观测器在离散控制系统中具有广泛的应用。
二、状态观测器的设计方法状态观测器的设计方法可以分为两类:基于传统观测器设计方法和基于最优观测器设计方法。
1. 基于传统观测器设计方法基于传统观测器设计方法的核心思想是通过系统的输出值来估计系统的状态变量。
最常用的传统观测器设计方法有:(1)全阶观测器设计:全阶观测器是指观测器的状态向量与系统的状态向量具有相同的维数。
全阶观测器可以通过系统的输出值和输入值来准确地估计系统的状态变量。
(2)低阶观测器设计:低阶观测器是指观测器的状态向量比系统的状态向量的维数低。
低阶观测器设计方法通过将系统的状态变量投影到一个低维的观测空间中来实现对系统状态的估计。
2. 基于最优观测器设计方法基于最优观测器设计方法的核心思想是通过优化问题来设计状态观测器,使得估计误差最小。
最优观测器能够最大程度地准确估计系统的状态变量。
最常用的最优观测器设计方法是卡尔曼滤波器。
卡尔曼滤波器能够通过系统的输出值和输入值来估计系统的状态变量,并且可以自适应地调整观测器的参数,以最小化估计误差。
三、状态观测器的应用状态观测器在离散控制系统中有广泛的应用。
全维状态观测器及其设计方法
《现代控制理论》课程课件节选
6.5.1 全维状态观测器及其设计方 法
1
全维状态观测器及其设计方法(1/1)
6.5.1 全维状态观测器及其设计方法
下面分别介绍 ➢ 开环状态观测器 ➢ 渐近状态观测器
2
开环状态观测器(1/6)
1. 开环状态观测器
设线性定常连续系统的状态空间模型为(A,B,C),即为
xˆ Axˆ Bu
yˆ
Cxˆ
其中 xˆ 为被控系统状态变量x(t)的估计值.
4
开环状态观测器(3/6)
该状态估计系统称为开环状态观测器, ➢ 简记为 ˆ (A, B,C),
➢ 其结构如下图所示.
u
+
B
x' ∫ x C y
+ A
+ B
xˆ
xˆ
∫
+
开环状态观测器
A
yˆ
C
xˆ
图6-8 开环状态观测器的结构图
➢ 为了和下面讨论的状态观测器区分开来,通常把该观测器 称为开环状态观测器.
8
5
开环状态观测器(4/6)
比较系统(A,B,C)和 ˆ (A, B,C)的状态变量,有
x&(t) xˆ&(t) A x(t) xˆ(t)
全维观测器的原理
全维观测器的原理
全维观测器(全局观测器)是一种系统控制理论中常用的观测器设计方法,它用于估计一个系统的未测量状态变量。
全维观测器的原理基于系统状态方程和输出方程,通过在系统中引入一个观测器来估计系统的状态变量。
观测器的结构与实际系统的结构相似,由一个状态方程、一个输出方程和一个观测误差方程组成。
全维观测器的状态方程是由实际系统的状态方程推导得到的,但是观测器的参数是通过一定的设计方法确定的。
观测器的输出方程与实际系统的输出方程相同,即观测器可以输出与实际系统完全相同的测量结果。
观测器的原理是基于对实际系统的估计误差进行反馈修正的原理。
通过观测器的输出与实际系统的输出之间的误差,可以得到对实际系统状态变量估计的误差。
利用这个误差,可以通过一定的修正算法来更新观测器的参数,使得观测器的估计结果逐渐接近实际系统的状态变量。
全维观测器的设计方法有很多,常见的方法包括最小二乘法、Kalmman滤波器等。
其中,最小二乘法是利用观测器的输出与实际系统的输出之间的误差最小化来确定观测器的参数。
Kalmman滤波器则是一种利用贝叶斯定理来估计系统状态的方法,它通过观测器的输出与实际系统的输出之间的协方差矩阵来确定观测器的参数。
总之,全维观测器通过利用系统的状态方程、输出方程和观测误差方程来估计系统的状态变量。
利用观测器的输出与实际系统的输出之间的误差,可以对观测器的参数进行修正,从而逐渐接近实际系统的状态变量。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
实 验 报 告
课程 线性系统理论基础 实验日期 2016年 6月 6 日 专业班级 姓名 学号 同组人 实
验
名
称
全
维
状
态
观
测
器
的
设
计
评分
批阅教师签字 一、实验目的
1、 学习用状态观测器获取系统状态估计值的方法,了解全维状态观测器的极点对状态的估计误差的影响;
2、 掌握全维状态观测器的设计方法;
3、 掌握带有状态观测器的状态反馈系统设计方法。
二、实验内容
开环系统⎩
⎨⎧=+=cx y bu Ax x
&,其中
[]0100001,0,10061161A b c ⎡⎤⎡⎤⎢⎥⎢⎥===⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦
a) 用状态反馈配置系统的闭环极点:5,322-±-j ; b) 设计全维状态观测器,观测器的极点为:10,325-±-j ; c) 研究观测器极点位置对估计状态逼近被估计值的影响; d) 求系统的传递函数(带观测器及不带观测器时);
绘制系统的输出阶跃响应曲线。
三、实验环境 MATLAB6、5
四、实验原理(或程序框图)及步骤
利用状态反馈可以使闭环系统的极点配置在所希望的位置上,其条件就是必须对全部状态变量都能进行测量,但在实际系统中,并不就是所有状态变量都能测量的,这就给状态反馈的实现造成了困难。
因此要设法利用已知的信息(输出量y 与输入量x),通过一个模型重新构造系统状态以对状态变量进行估计。
该模型就称为状态观测器。
若状态观测器的阶次与系统的阶次就是相同的,这样的状态观测器就称为全维状态观测器或全阶观测器。
设系统完全可观,则可构造如图4-1所示的状态观测器
图4-1 全维状态观测器
为求出状态观测器的反馈ke 增益,与极点配置类似,也可有两种方法: 方法一:构造变换矩阵Q,使系统变成标准能观型,然后根据特征方程求出k e ;
方
法
二
:
就
是
可
采用Ackermann 公式
:
[]T
o e Q A k 1000)(1
Λ-Φ=,其中O Q 为可观性矩阵。
利用对偶原理,可使设计问题大为简化。
首先构造对偶系统
⎩⎨⎧=+=ξ
ηξξ
T
T T b v c A & 然后可由变换法或Ackermann 公式求出极点配置的反馈k 增益,这也可
由MATLAB的place与acker函数得到;最后求出状态观测器的反馈增益。
五、程序源代码、实验数据、结果分析
(a)源程序:
A=[0 1 0;0 0 1;-6 -11 6];
B=[0;0;1];
C=[1 0 0];D=0;
P1=[-2+2*sqrt(3)*i;-2-2*sqrt(3)*i;-5];
K1=place(A,B,P1)
sysnew=ss(A-B*K1,B,C,D)
运行结果:
K1 =
74、0000 25、0000 15、0000
a =
x1 x2 x3
x1 0 1 0
x2 0 0 1
x3 -80 -36 -9
b = u1
x1 0
x2 0
x3 1
c = x1 x2 x3
y1 1 0
d = u1
y1 0
(b)源程序:
A=[0 1 0;0 0 1;-6 -11 6];
B=[0;0;1];
C=[1 0 0];D=0;
P2=[-5+2*sqrt(3)*i;-5-2*sqrt(3)*i;-10];
K2=acker(A',C',P2);L=K2'
Anew=A-L*C
运行结果:
L =
26
282
1770
Anew =
-26 1 0
-282 0 1
-1776 -11 6
(c)研究观测器极点位置对估计状态逼近被估计值的影响:
观测器极点距离虚轴越近,估计状态逼近被估计值得速度越快。
(d)不带观测器:
源程序:
A=[0 1 0;0 0 1;-6 -11 6];
B=[0;0;1];
C=[1 0 0];D=0;
P1=[-2+2*sqrt(3)*i;-2-2*sqrt(3)*i;-5];
K1=place(A,B,P1)
sysnew=ss(A-B*K1,B,C,D);
[num,den]=ss2tf(A-B*K1,B,C,D);
Gb=tf(num,den)
step(Gb)
grid on;
title('不带观测器的系统的阶跃响应曲线');
运行结果:
K1 =
74、0000 25、0000 15、0000
Transfer function:
7、105e-015 s^2 + 1、208e-013 s + 1
--------------------------------------------
s^3 + 9 s^2 + 36 s + 80
带观测器:
源程序:
A=[0 1 0;0 0 1;-6 -11 6];
B=[0;0;1];
C=[1 0 0];D=0;
P1=[-2+2*sqrt(3)*i;-2-2*sqrt(3)*i;-5];
K1=place(A,B,P1);
sysnew=ss(A-B*K1,B,C,D);
P2=[-5+2*sqrt(3)*i;-5-2*sqrt(3)*i;-10];
K2=acker(A',C',P2);L=K2';
An=[A -B*K1;L*C A-B*K1-L*C]
Bn=[B;B]
Cn=[C 0 0 0]
Dn=0;
[num,den]=ss2tf(An,Bn,Cn,Dn);
Go=tf(num,den)
step(Go)
grid on;
title('带观测器的系统的阶跃响应曲线');
运行结果:
An =
1、0e+003 *
0 0、0010 0 0 0 0
0 0 0、0010 0 0 0
-0、0060 -0、0110 0、0060 -0、0740 -0、0250 -0、0150
0、0260 0 0 -0、0260 0、0010 0
0、2820 0 0 -0、2820 0 0、0010
1、7700 0 0 -1、8500 -0、0360 -0、0090 Bn =
1
1
Cn =
1 0 0 0 0 0
Transfer function:
-1、137e-013 s^4 + s^3 + 20 s^2 + 137 s + 370
-------------------------------------------------------------------------------
s^6 + 29 s^5 + 353 s^4 + 2403 s^3 + 9862 s^2 + 2、428e004 s + 2、96e004
Time (sec)
A m p l i t u d e
结果分析:
σ%=10、8% tp=1、15s ts=1、63s 原系统方框图
Gain2
10
3
2.5
2
1.5
1
0.5
-0.5
原系统阶跃响应加观测器的方框图:
Scope1:
Scope2:
Scope3:。