Skip to content

方法Method

基于观测器的输出反馈与分离原理

Observer-based output feedback · Luenberger observer · Separation principle in linear control

用输出创新驱动全阶观测器,在状态与估计误差坐标中证明块三角分离,复算双积分器四极点及估计暂态的额外代价。

形式陈述 ​

考虑精确、无噪声的有限维连续时间线性模型

x′=Ax+Bu,y=Cx,

其中 A∈Rn×n、B∈Rn×m、C∈Rp×n。控制器读取输出 y,知道自己施加的输入 u,但不预设真实初态已知。

若 (A,B) 可稳定,可选 K 使 A−BK 为 Hurwitz;若 (A,C) 可检测,可选 L 使 A−LC 为 Hurwitz。采用全阶观测器和估计状态反馈

(1)x^′=Ax^+Bu+L(y−Cx^),u=−Kx^.

y−Cx^ 是实际输出与估计输出的差,称为创新。令估计误差 e=x−x^,则闭环在 (x,e) 坐标中为

(2)(x′e′)=(A−BKBK0A−LC)(xe).

因此全部闭环特征值,按代数重数计算,是 A−BK 与 A−LC 两组特征值的合并。两块都 Hurwitz 时,真实状态与估计误差均指数趋零。这就是这里的分离原理:稳定反馈与稳定误差动力学可分别设计,再由 (1) 连接。

这份结论针对准确的线性模型和稳定化目标,不声称估计误差为零、不改变暂态代价,也不自动给出含噪声或约束模型的最优控制。

直觉

观测器先用相同的动力学和已知输入预测状态,再用输出差修正预测。把真实系统与观测器相减时,共同的 Bu 消失,因此误差只由 A−LC 演化。控制器虽然使用估计值,估计误差仍像一个逐渐衰减的外部作用进入状态方程。

块三角中的右上块 BK 一般不为零:误差会影响真实状态。左下块为零,表示真实状态不反过来驱动已经相减后的误差方程。正是这个单向结构使特征多项式分成两份;它并未让实际轨道分成两个互不影响的过程。

分离原理的单向耦合与代价

图中箭头表示误差经过 BK 进入真实状态方程。下方采用本页的位置测量算例:初始估计为零时,误差虽然衰减,仍使总代价从全状态已知的二增到 9/2。

例子与边界

位置传感器下的四个闭环极点 ​

沿用双积分器

A=(0100),B=(01),C=(1 0).

执行器施加加速度,传感器只测位置。(A,B) 可控,且堆叠 C,CA 得 I2,所以位置的时间轨迹足以辨别速度。

选 LQR 页中 Q=diag(1,2)、R=1 的最优增益 K=(1 2),其闭环为

F=A−BK=(01−1−2),det⁡(sI−F)=(s+1)2.

令 L=(ℓ1,ℓ2)T,则

A−LC=(−ℓ11−ℓ20),det⁡(sI−A+LC)=s2+ℓ1s+ℓ2.

把误差极点设为 −2,−3,匹配 s2+5s+6,得到 L=(5,6)T。实际可执行的控制器为

x^1′=x^2+5(y−x^1),x^2′=u+6(y−x^1),u=−x^1−2x^2.

在 (x1,x2,e1,e2) 次序下,(2) 为

A=(0100−1−21200−5100−60),

其特征多项式为

det⁡(sI−A)=(s+1)2(s+2)(s+3).

四个极点因此确为 −1,−1,−2,−3。−1 允许出现 Jordan 块,但相应多项式乘指数项仍然衰减。

估计暂态带来的额外代价 ​

取真实初态 x(0)=(1,0)T,初始估计 x^(0)=0,故 e(0)=(1,0)T。误差方程可显式解为

e1(t)=−2e−2t+3e−3t,e2(t)=−6e−2t+6e−3t.

它满足初值,且直接求导可验证 e1′=−5e1+e2、e2′=−6e1。于是

Ke(t)=−14e−2t+15e−3t.

LQR 的价值矩阵为 P=(2112)。沿实际控制 u=−Kx^,有 u+Kx=Ke。整个闭环稳定,所以配平方恒等式中的终端项趋零,得到

J=x(0)TPx(0)+∫0∞(Ke(t))2dt=2+1964−4205+2256=2+52=92.

若一开始就令 x^(0)=x(0),误差恒零,才恢复全状态 LQR 的代价二。分离保证稳定,但初态估计误差的代价可以明确非零;这里没有在未知初态模型下声称该观测器增益具有最优性。

只测速度时无法补出恒定位置差 ​

把传感器改为 C=(0 1)。对任意 L=(ℓ1,ℓ2)T,

A−LC=(01−ℓ10−ℓ2)

都保留零特征值。位置方向 e1 不可观而且不衰减,因此系统不可检测。若真实初态为 e1、估计初态为零,则输出恒零、控制恒零,真实位置一直保持一,而观测器一直保持零。这不是把 L 选得不够大,而是整个输出轨迹不包含该位置差。

对正常的位置传感器,单纯静态输出反馈 u=−ky=−kx1 也不够:闭环多项式为 s2+k,任意实数 k 都不能使两根实部严格为负。动态观测器利用整段输出的演化获得速度信息,才把稳定状态反馈变成可实现的输出控制。

推论与应用

误差方程和坐标变换 ​

先保留一般已知输入。由真实状态和观测器方程相减,

e′=Ax+Bu−Ax^−Bu−L(Cx−Cx^)=(A−LC)e.

当 u=−Kx^=−K(x−e) 时,真实状态满足 x′=(A−BK)x+BKe,得到 (2)。这个推导同时固定了误差的符号,不能把 e=x−x^ 与相反约定混用。

也可直接核验相似变换。在 (x,x^) 坐标中,矩阵是

M=(A−BKLCA−BK−LC).

令

T=(I0I−I),T−1=T.

则 T(x,x^)T=(x,e)T,直接块乘法得到 TMT−1=A,其中 A 是 (2) 的块矩阵。因此这是真实闭环的同一组特征值,没有通过舍掉观测器状态来降低系统阶数。

为什么块耦合不破坏指数衰减 ​

由块三角行列式,特征多项式正是两个对角块的乘积。也可从轨道本身看出衰减。若

‖e(A−BK)t‖≤M1e−αt,‖e(A−LC)t‖≤M2e−βt,

则 e(t)=e(A−LC)te(0),而

x(t)=e(A−BK)tx(0)+∫0te(A−BK)(t−s)BKe(A−LC)se(0)ds.

积分范数至多为

M1M2‖BK‖‖e(0)‖te−min{α,β}t.

对任意 0<γ<min{α,β},te−(min{α,β}−γ)t 有界,故整个轨道仍有统一的 e−γt 上界。这也覆盖两组极点重合时的多项式因子。

设计条件各自负责什么 ​

可稳定性让 K 消除输入能够影响的不稳定模态;可检测性让 L 消除输出能够辨别的不稳定误差模态。不可控或不可观的稳定模态可以保留,所以不要求全部状态都可控且可观。

若要全阶观测器对所有初始误差都指数收敛,A−LC 必须 Hurwitz;不可观的不稳定或零实部模态会原样留在误差方程中。反馈稳定也同样需要 A−BK Hurwitz。这两项条件应分别检查,不可用“可控”替代“可检测”。

若实际测量是 y=Cx+v,误差方程会变成 e′=(A−LC)e−Lv;若模型或施加的输入与观测器内部不一致,还会出现额外项。此时必须分析噪声或失配的影响,不能只看无扰动的四个极点。采样实现也需要先转换到对应离散模型,再判断离散闭环稳定性。

参考资料
  • Lars Grüne, Mathematical Control Theory, Winter 2020/2021,§§4.3–4.4, Theorems 4.17、4.19, pp. 46–49:动态观测器、误差方程及输出反馈的块三角证明。原文采用不同反馈和误差符号;本文统一为 u=−Kx^、e=x−x^。
关系图谱17 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系