形式陈述
考虑精确、无噪声的有限维连续时间线性模型
x ′ = A x + B u , y = C x , 其中 A ∈ R n × n 、B ∈ R n × m 、C ∈ R p × n 。控制器读取输出 y ,知道自己施加的输入 u ,但不预设真实初态已知。
若 ( A , B ) 可稳定 公理库 线性状态反馈与稳定化 Linear state feedback · 线性反馈稳定化 用全状态反馈改变闭环谱,证明单输入可控系统的极点配置,并区分可控与可稳定。 ,可选 K 使 A − B K 为 Hurwitz;若 ( A , C ) 可检测 公理库 线性系统的可观测性 Observability of a linear system · Kalman 可观测性秩判据 用观测 Gramian 重建初态,再分离不可观子空间,证明可检测性与稳定输出注入的等价关系。 ,可选 L 使 A − L C 为 Hurwitz。采用全阶观测器和估计状态反馈
(1) x ^ ′ = A x ^ + B u + L ( y − C x ^ ) , u = − K x ^ . y − C x ^ 是实际输出与估计输出的差,称为创新。令估计误差 e = x − x ^ ,则闭环在 ( x , e ) 坐标中为
(2) ( x ′ e ′ ) = ( A − B K B K 0 A − L C ) ( x e ) . 因此全部闭环特征值,按代数重数计算,是 A − B K 与 A − L C 两组特征值的合并。两块都 Hurwitz 时,真实状态与估计误差均指数趋零。这就是这里的分离原理 :稳定反馈与稳定误差动力学可分别设计,再由 (1) 连接。
这份结论针对准确的线性模型和稳定化目标,不声称估计误差为零、不改变暂态代价,也不自动给出含噪声或约束模型的最优控制。
直觉
观测器先用相同的动力学和已知输入预测状态,再用输出差修正预测。把真实系统与观测器相减时,共同的 B u 消失,因此误差只由 A − L C 演化。控制器虽然使用估计值,估计误差仍像一个逐渐衰减的外部作用进入状态方程。
块三角中的右上块 B K 一般不为零:误差会影响真实状态。左下块为零,表示真实状态不反过来驱动已经相减后的误差方程。正是这个单向结构使特征多项式分成两份;它并未让实际轨道分成两个互不影响的过程。
图片加载失败 分离原理的单向耦合与代价 图中箭头表示误差经过 B K 进入真实状态方程。下方采用本页的位置测量算例:初始估计为零时,误差虽然衰减,仍使总代价从全状态已知的二增到 9 / 2 。
例子与边界
位置传感器下的四个闭环极点
沿用双积分器
A = ( 0 1 0 0 ) , B = ( 0 1 ) , C = ( 1 0 ) . 执行器施加加速度,传感器只测位置。( A , B ) 可控,且堆叠 C , C A 得 I 2 ,所以位置的时间轨迹足以辨别速度。
选 LQR 页 公理库 连续时间线性二次调节器 Continuous-time linear quadratic regulator · LQR · 线性二次调节器 从零终端代价的有限时域 Riccati 解构造稳定 CARE 极限,在半正定状态代价下证明最优性与唯一性,并复算双积分器及双曲正切轨道。 中 Q = diag ( 1 , 2 ) 、R = 1 的最优增益 K = ( 1 2 ) ,其闭环为
F = A − B K = ( 0 1 − 1 − 2 ) , det ( s I − F ) = ( s + 1 ) 2 . 令 L = ( ℓ 1 , ℓ 2 ) T ,则
A − L C = ( − ℓ 1 1 − ℓ 2 0 ) , det ( s I − A + L C ) = s 2 + ℓ 1 s + ℓ 2 . 把误差极点设为 − 2 , − 3 ,匹配 s 2 + 5 s + 6 ,得到 L = ( 5 , 6 ) T 。实际可执行的控制器为
x ^ 1 ′ = x ^ 2 + 5 ( y − x ^ 1 ) , x ^ 2 ′ = u + 6 ( y − x ^ 1 ) , u = − x ^ 1 − 2 x ^ 2 . 在 ( x 1 , x 2 , e 1 , e 2 ) 次序下,(2) 为
A = ( 0 1 0 0 − 1 − 2 1 2 0 0 − 5 1 0 0 − 6 0 ) , 其特征多项式 公理库 特征多项式 Characteristic polynomial 由 det(λI−T) 定义并编码线性算子特征值的多项式。 为
det ( s I − 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 。误差方程可显式解为
e 1 ( t ) = − 2 e − 2 t + 3 e − 3 t , e 2 ( t ) = − 6 e − 2 t + 6 e − 3 t . 它满足初值,且直接求导可验证 e 1 ′ = − 5 e 1 + e 2 、e 2 ′ = − 6 e 1 。于是
K e ( t ) = − 14 e − 2 t + 15 e − 3 t . LQR 的价值矩阵为 P = ( 2 1 1 2 ) 。沿实际控制 u = − K x ^ ,有 u + K x = K e 。整个闭环稳定,所以配平方恒等式中的终端项趋零,得到
J = x ( 0 ) T P x ( 0 ) + ∫ 0 ∞ ( K e ( t ) ) 2 d t = 2 + 196 4 − 420 5 + 225 6 = 2 + 5 2 = 9 2 . 若一开始就令 x ^ ( 0 ) = x ( 0 ) ,误差恒零,才恢复全状态 LQR 的代价二。分离保证稳定,但初态估计误差的代价可以明确非零;这里没有在未知初态模型下声称该观测器增益具有最优性。
只测速度时无法补出恒定位置差
把传感器改为 C = ( 0 1 ) 。对任意 L = ( ℓ 1 , ℓ 2 ) T ,
A − L C = ( 0 1 − ℓ 1 0 − ℓ 2 ) 都保留零特征值。位置方向 e 1 不可观而且不衰减,因此系统不可检测。若真实初态为 e 1 、估计初态为零,则输出恒零、控制恒零,真实位置一直保持一,而观测器一直保持零。这不是把 L 选得不够大,而是整个输出轨迹不包含该位置差。
对正常的位置传感器,单纯静态输出反馈 u = − k y = − k x 1 也不够:闭环多项式为 s 2 + k ,任意实数 k 都不能使两根实部严格为负。动态观测器利用整段输出的演化获得速度信息,才把稳定状态反馈变成可实现的输出控制。
推论与应用
误差方程和坐标变换
先保留一般已知输入。由真实状态和观测器方程相减,
e ′ = A x + B u − A x ^ − B u − L ( C x − C x ^ ) = ( A − L C ) e . 当 u = − K x ^ = − K ( x − e ) 时,真实状态满足 x ′ = ( A − B K ) x + B K e ,得到 (2)。这个推导同时固定了误差的符号,不能把 e = x − x ^ 与相反约定混用。
也可直接核验相似变换。在 ( x , x ^ ) 坐标中,矩阵是
M = ( A − B K L C A − B K − L C ) . 令
T = ( I 0 I − I ) , T − 1 = T . 则 T ( x , x ^ ) T = ( x , e ) T ,直接块乘法得到 T M T − 1 = A ,其中 A 是 (2) 的块矩阵。因此这是真实闭环的同一组特征值,没有通过舍掉观测器状态来降低系统阶数。
为什么块耦合不破坏指数衰减
由块三角行列式,特征多项式正是两个对角块的乘积。也可从轨道本身看出衰减。若
‖ e ( A − B K ) t ‖ ≤ M 1 e − α t , ‖ e ( A − L C ) t ‖ ≤ M 2 e − β t , 则 e ( t ) = e ( A − L C ) t e ( 0 ) ,而
x ( t ) = e ( A − B K ) t x ( 0 ) + ∫ 0 t e ( A − B K ) ( t − s ) B K e ( A − L C ) s e ( 0 ) d s . 积分范数至多为
M 1 M 2 ‖ B K ‖ ‖ e ( 0 ) ‖ t e − min { α , β } t . 对任意 0 < γ < min { α , β } ,t e − ( min { α , β } − γ ) t 有界,故整个轨道仍有统一的 e − γ t 上界。这也覆盖两组极点重合时的多项式因子。
设计条件各自负责什么
可稳定性让 K 消除输入能够影响的不稳定模态;可检测性让 L 消除输出能够辨别的不稳定误差模态。不可控或不可观的稳定模态可以保留,所以不要求全部状态都可控且可观。
若要全阶观测器对所有初始误差都指数收敛,A − L C 必须 Hurwitz;不可观的不稳定或零实部模态会原样留在误差方程中。反馈稳定也同样需要 A − B K Hurwitz。这两项条件应分别检查,不可用“可控”替代“可检测”。
若实际测量是 y = C x + v ,误差方程会变成 e ′ = ( A − L C ) e − L v ;若模型或施加的输入与观测器内部不一致,还会出现额外项。此时必须分析噪声或失配的影响,不能只看无扰动的四个极点。采样实现也需要先转换到对应离散模型,再判断离散闭环稳定性。
参考资料