形式陈述
固定状态维数 n ≥ 1 。给定实常矩阵 A ∈ R n × n 、B ∈ R n × m ,考虑连续时间线性方程组 公理库 线性常微分方程组 Linear system of ordinary differential equations 形如 x′=A(t)x+b(t) 的向量值一阶线性方程组。
x ′ = A x + B u , x ( 0 ) = x 0 . 在有限区间 [ 0 , T ] 上取 u ∈ L 2 ( [ 0 , T ] ; R m ) ,不限制输入幅值与变化率。状态取绝对连续函数,方程几乎处处成立;唯一解由常数变易公式给出。称 ( A , B ) 可控,若对每个 T > 0 和任意 x 0 , x T ∈ R n ,都存在输入使 x ( T ) = x T 。
定义可控矩阵与有限时域可控 Gramian:
C = [ B A B ⋯ A n − 1 B ] , W c ( T ) = ∫ 0 T e A t B B T e A T t d t . 以下条件等价:( A , B ) 可控;C 的秩 公理库 线性映射的秩 Rank of a linear map · Matrix rank 线性映射像空间的维数,表示其保留下来的独立输出方向数。 等于 n ;对某个 T > 0 ,W c ( T ) 正定 公理库 正定与半正定矩阵 Positive definite matrix · Positive semidefinite matrix · PSD matrix 由二次能量严格为正或非负定义的实对称与复 Hermitian 矩阵。 ;对每个 T > 0 ,W c ( T ) 正定。有限时域结论不要求 A 稳定。
直觉
B 给出输入当下能推动的方向,A B , A 2 B , … 则记录动力学怎样把这些作用传递到其他坐标。执行器少于状态维数不意味着不可控:一个改变速度的力,经过时间积累也能改变位置。真正的问题是这些传播方向能否张成整个状态空间。
图片加载失败 输入可达方向与输出不可辨方向 正常执行器 B = e 2 经动力学传播能张成二维空间,而退化执行器 B = e 1 只能产生位置方向的增量。左图蓝色方向 e A t e 2 = ( t , 1 ) T 随时间改变,两个不同时间的方向已线性无关;红色箭头表示始终不变的 e 1 ,为避免与坐标轴重叠而稍向下平移,其竖直分量仍为零。右图只测速度,同一水平虚线上的状态具有相同输出,表达的是可观测性 公理库 线性系统的可观测性 Observability of a linear system · Kalman 可观测性秩判据 用观测 Gramian 重建初态,再分离不可观子空间,证明可检测性与稳定输出注入的等价关系。 中的初态不可辨,而非输入方向的缺失。
例子与边界
为什么秩与 Gramian 给出同一个答案
对任意 v ∈ R n ,
v T W c ( T ) v = ∫ 0 T ‖ B T e A T t v ‖ 2 d t . 因此 W c 半正定。若上式为零,由被积函数连续可知 B T e A T t v 恒为零;在 t = 0 逐阶求导,得到 B T ( A T ) k v = 0 ,特别是 0 ≤ k < n 时成立。反过来,若这前 n 项为零,Cayley–Hamilton 定理 公理库 Cayley–Hamilton 定理 Cayley–Hamilton theorem 每个有限维线性算子都满足自己的特征多项式。 将更高次幂化为前 n 次幂的线性组合,故所有项均为零,矩阵指数的级数也为零。于是
ker W c ( T ) = ker C T . 这对每个正 T 都成立,从而证明了秩条件与两个正定条件等价。
再把目标与自由演化的差记为 d = x T − e A T x 0 ,并令
L T u = ∫ 0 T e A ( T − s ) B u ( s ) d s . 若 W c 可逆,取
u ∗ ( s ) = B T e A T ( T − s ) W c ( T ) − 1 d , 直接代入得到 L T u ∗ = d ,所以目标可达。若 W c 奇异,取非零 v ∈ ker W c ,所有输入都满足 v T L T u = 0 ,目标增量 d = v 就不可能实现。这完成全部等价关系的证明,也说明从零出发的全空间可达性在本模型中足够判断任意初末状态的可控性。
最小输入能量与双积分器
上述 u ∗ 还在所有到达同一终点的输入中唯一地(按几乎处处相等)最小化 ∫ 0 T ‖ u ‖ 2 d t 。事实上,任意可行输入可写成 u = u ∗ + w ,其中 L T w = 0 ,而
∫ 0 T u ∗ T w d s = d T W c − 1 L T w = 0. 勾股分解因此给出
∫ 0 T ‖ u ‖ 2 d s = d T W c − 1 d + ∫ 0 T ‖ w ‖ 2 d s . 取位置、速度状态 x = ( x 1 , x 2 ) T ,令
A = ( 0 1 0 0 ) , B = ( 0 1 ) , e A t = ( 1 t 0 1 ) . 此时 C = ( 0 1 1 0 ) 满秩,且
W c ( T ) = ( T 3 / 3 T 2 / 2 T 2 / 2 T ) , det W c ( T ) = T 4 / 12 > 0. 要求在一单位时间内从 ( 0 , 0 ) 到 ( 1 , 0 ) ,有
W c ( 1 ) − 1 = ( 12 − 6 − 6 4 ) , u ∗ ( s ) = 6 − 12 s . 积分得到速度 x 2 ( s ) = 6 s − 6 s 2 、位置 x 1 ( s ) = 3 s 2 − 2 s 3 ,端点确为 ( 1 , 0 ) ,能量为 ∫ 0 1 ( 6 − 12 s ) 2 d s = 12 。控制先加速、后减速,最终位置改变而速度归零。
若改成 B = e 1 ,则 A B = 0 ,可达增量只有 span ( e 1 ) ,W c ( T ) = diag ( T , 0 ) ;速度不受输入影响。这里所有结论均在无约束实输入下成立,加入饱和或非负输入约束后,满秩不能保证任意时限内到达任意目标。
推论与应用
可达子空间与不可控块
令 R = im C 。Cayley–Hamilton 定理给出 A R ⊆ R ,且 im B ⊆ R 。选取 R 的一组基,再补成全空间基,以这些列组成可逆矩阵 T 。在坐标 x = T z 下,
(1) A ~ = T − 1 A T = ( A c A 12 0 A u ) , B ~ = T − 1 B = ( B c 0 ) . 左下零块来自 R 的不变性,输入的下块为零来自 im B ⊆ R 。所有 A ~ j B ~ 的上块为 A c j B c 、下块为零;它们张成全部前 r = dim R 个坐标。因此 ( A c , B c ) 可控,更高次幂再由 A c 的 Cayley–Hamilton 关系约化。
任意负反馈 公理库 线性状态反馈与稳定化 Linear state feedback · 线性反馈稳定化 用全状态反馈改变闭环谱,证明单输入可控系统的极点配置,并区分可控与可稳定。 u = − K ~ z ,写 K ~ = ( K c K u ) ,都会给出
A ~ − B ~ K ~ = ( A c − B c K c A 12 − B c K u 0 A u ) . 所以不可控块 A u 的特征值完全不随反馈改变。系统可稳定当且仅当 A u 为 Hurwitz,即其全部特征值实部严格为负;下面补足“可控块总能稳定”的充分性。
多输入可控块的一个稳定增益
若 r > 0 ,选 α > ‖ A c ‖ 2 ,令 H = A c + α I 。− H 为 Hurwitz。平移和取负号不改变输入的 Krylov 张成空间,故 ( − H , B c ) 仍可控。于是
W = ∫ 0 ∞ e − H t B c B c T e − H T t d t 收敛且正定:收敛来自指数衰减,正定可由任意正时长的有限 Gramian 正定得到。对被积矩阵求导并积分,有
H W + W H T = B c B c T . 取 K c = B c T W − 1 、F = A c − B c K c ,则
F W + W F T = − 2 α W − B c B c T ≺ 0. 左右乘 W − 1 ,得到 F T W − 1 + W − 1 F ≺ 0 。对任意复特征向量 F v = λ v ,左侧的二次型为 2 Re λ v ∗ W − 1 v ,故 Re λ < 0 。这直接证明可控块可以稳定,且不要求单输入。
若 A u 已稳定,令 K u = 0 即使两个对角块都稳定,原坐标增益为 K = K ~ T − 1 。若 r = 0 ,没有可控块,结论退化为 A 自身必须稳定;若 r = n ,不可控块为空,其稳定条件自动满足。
例如 A = diag ( 0 , − 2 ) 、B = e 1 只有第一方向可控,但第二方向自行指数衰减。取 K = ( 1 0 ) 即得闭环 diag ( − 1 , − 2 ) 。把第二个对角元改成零或正数,就无法通过这个执行器实现指数稳定。半正定代价的 LQR 例子 公理库 连续时间线性二次调节器 Continuous-time linear quadratic regulator · Continuous-time LQR · 连续时间LQR 从零终端代价的有限时域 Riccati 解构造稳定 CARE 极限,在半正定状态代价下证明最优性与唯一性,并复算双积分器及双曲正切轨道。 继续使用前一种系统,明确计算有限时域到无限时域的价值矩阵。
与可观测性的转置对应
上述有限维、实常系数、任意正有限时域的可控性判据,经 ( A , B ) ↦ ( A T , B T ) ,恰与可观测性判据 公理库 线性系统的可观测性 Observability of a linear system · Kalman 可观测性秩判据 用观测 Gramian 重建初态,再分离不可观子空间,证明可检测性与稳定输出注入的等价关系。 对应。转置系统把 B T 当输出矩阵,其观测矩阵是 C T ,观测 Gramian 也等于本页的 W c ( T ) 。所以 ( A , B ) 可控当且仅当该转置系统可观。这个等价交换的是给定矩阵的输入与输出角色,并不把原系统的实际执行器当成传感器。
到达能量与稳定调节
满秩回答“能否到达”,d T W c − 1 d 则回答“到达要花多少输入能量”。Gramian 很小的特征值对应昂贵方向,因此理论可控不等于数值上容易控制。数值求解宜使用线性方程求解器计算 W c − 1 d ,不必显式求逆。
线性状态反馈 公理库 线性状态反馈与稳定化 Linear state feedback · 线性反馈稳定化 用全状态反馈改变闭环谱,证明单输入可控系统的极点配置,并区分可控与可稳定。 进一步让输入随当前状态变化,以改变未来轨道;线性二次调节器 公理库 连续时间线性二次调节器 Continuous-time linear quadratic regulator · Continuous-time LQR · 连续时间LQR 从零终端代价的有限时域 Riccati 解构造稳定 CARE 极限,在半正定状态代价下证明最优性与唯一性,并复算双积分器及双曲正切轨道。 还同时惩罚状态偏差与输入。本页求固定终点的有限时域最小能量输入,LQR 则求无限时域的最优调节反馈;各自的最优解由相应时域、约束和代价共同决定。
参考资料