形式陈述
对实连续时间常系数系统 x ′ = A x + B u ,假设当前全状态 x 可准确获得。固定负反馈记号
u = − K x , x ′ = ( A − B K ) x . K ∈ R m × n 为常增益,闭环矩阵是 A − B K 。若采用另一约定 u = F x ,则 F = − K ,闭环相应写作 A + B F 。称 ( A , B ) 可稳定,若存在 K 使 A − B K 为 Hurwitz 矩阵,即全部特征值实部严格为负,从而原点全局指数稳定 公理库 平衡点稳定性 Stability of an equilibrium 初值受到小扰动时轨道保持接近或最终回到平衡点的性质。 。
单输入时写 B = b ∈ R n 。若 ( A , b ) 可控 公理库 线性系统的可控性 Controllability of a linear system · Kalman 可控性秩判据 用可控矩阵与有限时域 Gramian 判断输入能否到达任意状态,并构造最小能量输入。 ,则任意给定的实首一 n 次多项式 q 都可实现为 A − b K 的特征多项式。因此可以配置任意按代数重数计数、关于复共轭封闭的 n 个极点。这里给出单输入结论的完整证明。
直觉
开环输入是一条预先选好的时间函数;状态反馈则不断根据偏离原点的状态修正作用。在双积分器中,位置反馈产生回复力,速度反馈产生阻尼,两者共同决定闭环特征多项式。
把极点移入左半平面保证衰减,却未决定控制消耗有多大、测量噪声会被放大多少,也不表示越快越好。全状态可用还是一个实现前提:只测一个输出时,不能直接把未测坐标当作已知值代入反馈律。
例子与边界
单输入极点配置为什么能逐项解出
令 p ( s ) = det ( s I − A ) ,adj ( M ) 表示矩阵 M 的伴随矩阵。按列的多线性展开 det ( M + b K ) :凡选入 b K 的两列或更多列,这些列都平行于 b ,对应项为零;只选一列的项由余子式展开合为 K adj ( M ) b 。因此
det ( s I − A + b K ) = p ( s ) + K adj ( s I − A ) b . 右边附加项是次数至多 n − 1 的多项式,因此定义一个从 n 维行向量空间到 n 维多项式空间的线性映射
K ⟼ K adj ( s I − A ) b . 只需证明它单射。取任一从属于向量范数的矩阵范数;当 | s | > ‖ A ‖ 时,‖ A / s ‖ < 1 ,矩阵几何级数(Neumann 级数)收敛,给出 ( s I − A ) − 1 = s − 1 ∑ j ≥ 0 ( A / s ) j 。若上述多项式为零,在这个区域除以 p ( s ) ,利用伴随矩阵的求逆公式得到
0 = K ( s I − A ) − 1 b = ∑ j ≥ 0 K A j b s j + 1 . 幂级数系数唯一性给出 K b = K A b = ⋯ = K A n − 1 b = 0 。可控矩阵 [ b A b ⋯ A n − 1 b ] 可逆,故 K = 0 。等维线性映射单射即满射,所以可令附加项等于 q − p ,匹配其 n 个系数,得到所需实增益。由此,配置特征多项式 公理库 特征多项式 Characteristic polynomial 由 det(λI−T) 定义并编码线性算子特征值的多项式。 的系数就成为一个可逆的线性问题。
两个极点与一个无法移动的方向
对 A = ( 0 1 0 0 ) 、B = ( 0 1 ) ,写 K = ( k 1 k 2 ) ,则
A − B K = ( 0 1 − k 1 − k 2 ) , det ( s I − A + B K ) = s 2 + k 2 s + k 1 . 目标极点 − 1 , − 2 对应 ( s + 1 ) ( s + 2 ) = s 2 + 3 s + 2 ,所以 K = ( 2 3 ) ,控制为 u = − 2 x 1 − 3 x 2 。
一般地,若非零复行向量 w ∗ 满足 w ∗ A = λ w ∗ 、w ∗ B = 0 ,则对任意 K ,
w ∗ ( A − B K ) = λ w ∗ . 这个模态无法被反馈移动。将双积分器执行器改成 B = e 1 ,任意闭环都为 ( − k 1 1 − k 2 0 0 ) ,必保留零特征值,因而不可稳定。
不可控系统也可能可稳定。例如 A = diag ( 1 , − 2 ) , B = e 1 的第二方向不可控,但已自行衰减;取 K = ( 2 0 ) 得闭环 diag ( − 1 , − 2 ) 。一般多输入系统同样有判据:可稳定当且仅当所有不可控模态严格位于左半平面;其完整结论见参考资料,以上配置证明限定于单输入情形。
推论与应用
可控性允许配置所有模态,可稳定性只要求消除不稳定模态。两者的区别在欠驱动系统中特别重要:已稳定且不受输入影响的方向不妨碍调节到零,却妨碍到达任意终点。
线性二次调节器 公理库 连续时间线性二次调节器 Continuous-time linear quadratic regulator · LQR · 线性二次调节器 在可稳定与代价可检测条件下,通过稳定 Riccati 解求无限时域二次代价的最优状态反馈。 在稳定反馈之间引入状态和输入的代价权衡。上述 K = ( 2 3 ) 将在同一个双积分器上与最优增益比较,得到“两个反馈都稳定,但代价不同”的精确算例。实现该反馈需要位置和速度;若只测速度,初始位置不可观,必须补充信息才能获得所需全状态。
参考资料