形式陈述
考虑实系数、有限维、连续时间常系数系统
x ′ = A x + B u , y = C x , 其中 A ∈ R n × n 、B ∈ R n × m 、C ∈ R p × n 。在任意 [ 0 , T ] 上,输入 u ∈ L 2 已知,状态是满足线性 ODE 公理库 线性常微分方程组 Linear system of ordinary differential equations 形如 x′=A(t)x+b(t) 的向量值一阶线性方程组。 的绝对连续解,输出无测量噪声。称 ( A , C ) 可观测,若任意 T > 0 上的输入、输出轨迹唯一确定初态 x 0 。
定义观测矩阵与观测 Gramian:
O = ( C C A ⋮ C A n − 1 ) , W o ( T ) = ∫ 0 T e A T t C T C e A t d t . 可观测性等价于 O 的列秩 公理库 线性映射的秩 Rank of a linear map · Matrix rank 线性映射像空间的维数,表示其保留下来的独立输出方向数。 为 n ,也等价于某个(从而每个)正 T 上 W o ( T ) 正定 公理库 正定与半正定矩阵 Positive definite matrix · Positive semidefinite matrix · PSD matrix 由二次能量严格为正或非负定义的实对称与复 Hermitian 矩阵。 。传感器不必瞬时读出整个状态,动力学带来的时间信息可以补足缺失坐标。
直觉
仅看一张位置照片,无法知道初速度;看一段准确的位置轨迹则可以。可观测性比较的是两个初态在同一已知输入下能否产生完全相同的输出。输入的贡献可以计算并扣除,因此判断初态是否可辨只依赖 A , C ,不依赖本次实际输入有多丰富。
这里 A , B , C 都已给定,所求的是初始状态。状态空间时间序列模型 公理库 状态空间时间序列模型 State-space time-series model · State space model · 状态空间模型 用隐藏状态的 Markov 演化与给定状态的观测分布分解时间依赖的序列模型。 还要描述随机状态与观测误差,未知参数的统计识别与状态的后验不确定性则需要进一步考察观测分布。
例子与边界
从不可辨方向到重建公式
由常数变易公式,扣掉已知输入贡献后得到
z ( t ) = y ( t ) − C ∫ 0 t e A ( t − s ) B u ( s ) d s = C e A t x 0 . 两初态相差 v 时,校正输出之差为 C e A t v ,且
v T W o ( T ) v = ∫ 0 T ‖ C e A t v ‖ 2 d t . 零积分由连续性等价于输出差恒为零。对恒等式在零点求导,得到 C A k v = 0 ,故 v ∈ ker O 。反向使用Cayley–Hamilton 定理 公理库 Cayley–Hamilton 定理 Cayley–Hamilton theorem 每个有限维线性算子都满足自己的特征多项式。 :前 n 项为零推出全部高次项为零,再代入指数级数。因此
对 所 有 ker W o ( T ) = ker O = { v : C e A t v = 0 对所有 t ∈ [ 0 , T ] } . 只要这个核非零,x 0 与 x 0 + v 就无法区分;核为零则初态唯一,证明了秩、正定与可观测性等价。满秩时还有具体重建式:将 z ( t ) = C e A t x 0 左乘 e A T t C T 并积分,得到
x 0 = W o ( T ) − 1 ∫ 0 T e A T t C T z ( t ) d t . 同一运动,两个传感器
令 A = ( 0 1 0 0 ) ,正常传感器 C = ( 1 0 ) 测位置。则
C e A t = ( 1 t ) , O = I 2 , W o ( T ) = ( T T 2 / 2 T 2 / 2 T 3 / 3 ) , det W o ( T ) = T 4 / 12. 校正输出 z ( t ) = a + b t 的截距和斜率分别为初始位置与速度。在 T = 1 时,若测得 z ( t ) = 2 − t ,则
∫ 0 1 ( 1 t ) z ( t ) d t = ( 3 / 2 2 / 3 ) , W o ( 1 ) − 1 = ( 4 − 6 − 6 12 ) , 相乘确得 x 0 = ( 2 , − 1 ) T 。
若改测速度,C = ( 0 1 ) ,则 C A = 0 、W o ( T ) = diag ( 0 , T ) ,不可观空间为 span ( e 1 ) 。任意初始位置差都不改变速度输出。保留执行器 B = e 2 时系统仍然可控,说明可控不推出可观;反过来取 B = e 1 , C = ( 1 0 ) ,系统可观却不可控。
重建公式要求准确、连续时间的输出轨迹。噪声或模型误差会传入重建结果,即使 Gramian 正定,接近奇异也会放大这些误差。有限采样下须按采样后的模型重新判断可观测性,再用相应的离散重建或估计方法。
推论与应用
( A , C ) 可观等价于转置输入对 ( A T , C T ) 可控 公理库 线性系统的可控性 Controllability of a linear system · Kalman 可控性秩判据 用可控矩阵与有限时域 Gramian 判断输入能否到达任意状态,并构造最小能量输入。 :二者的秩矩阵互为转置。这是输入与输出角色交换的代数对偶,不表示实际执行器和传感器可以互换。
Kalman 滤波器 公理库 Kalman 滤波器 Kalman filter · 卡尔曼滤波器 在线性高斯状态空间模型中递推计算状态预测与过滤分布均值、协方差及似然的算法。 处理带噪声状态的条件估计;可观测性提供“哪些初态差能进入数据”的结构问题。LQR 公理库 连续时间线性二次调节器 Continuous-time linear quadratic regulator · LQR · 线性二次调节器 在可稳定与代价可检测条件下,通过稳定 Riccati 解求无限时域二次代价的最优状态反馈。 使用的可检测性则允许不可观方向存在,只要求它们自身稳定;那里观察矩阵来自状态代价,未必是物理传感器 C 。
参考资料
Jonathan P. How、Emilio Frazzoli,MIT 16.30/31,2010,Topic 14: Estimators ,14–6 至 14–10,已知输入、转置对偶与观测秩条件。
Sanjay Lall,Stanford Engr210a,2001-10-30,Lecture 9 ,9–14 的输出分解。本文由该分解推导有限时域核等式与初态重建;这一有限区间论证对稳定和不稳定的 A 都成立。