Skip to content

定理Theorem

线性系统的可观测性

Observability of a linear system · Kalman 可观测性秩判据

由已知输入和无噪声输出轨迹判断初始状态能否唯一辨别,并用 Gramian 重建初态。

形式陈述 ​

考虑实系数、有限维、连续时间常系数系统

x′=Ax+Bu,y=Cx,

其中 A∈Rn×n、B∈Rn×m、C∈Rp×n。在任意 [0,T] 上,输入 u∈L2 已知,状态是满足线性 ODE的绝对连续解,输出无测量噪声。称 (A,C) 可观测,若任意 T>0 上的输入、输出轨迹唯一确定初态 x0。

定义观测矩阵与观测 Gramian:

O=(CCA⋮CAn−1),Wo(T)=∫0TeATtCTCeAtdt.

可观测性等价于 O 的列秩为 n,也等价于某个(从而每个)正 T 上 Wo(T) 正定。传感器不必瞬时读出整个状态,动力学带来的时间信息可以补足缺失坐标。

直觉

仅看一张位置照片,无法知道初速度;看一段准确的位置轨迹则可以。可观测性比较的是两个初态在同一已知输入下能否产生完全相同的输出。输入的贡献可以计算并扣除,因此判断初态是否可辨只依赖 A,C,不依赖本次实际输入有多丰富。

这里 A,B,C 都已给定,所求的是初始状态。状态空间时间序列模型还要描述随机状态与观测误差,未知参数的统计识别与状态的后验不确定性则需要进一步考察观测分布。

例子与边界

从不可辨方向到重建公式 ​

由常数变易公式,扣掉已知输入贡献后得到

z(t)=y(t)−C∫0teA(t−s)Bu(s)ds=CeAtx0.

两初态相差 v 时,校正输出之差为 CeAtv,且

vTWo(T)v=∫0T‖CeAtv‖2dt.

零积分由连续性等价于输出差恒为零。对恒等式在零点求导,得到 CAkv=0,故 v∈ker⁡O。反向使用Cayley–Hamilton 定理:前 n 项为零推出全部高次项为零,再代入指数级数。因此

ker⁡Wo(T)=ker⁡O={v:CeAtv=0 对所有 t∈[0,T]}.

只要这个核非零,x0 与 x0+v 就无法区分;核为零则初态唯一,证明了秩、正定与可观测性等价。满秩时还有具体重建式:将 z(t)=CeAtx0 左乘 eATtCT 并积分,得到

x0=Wo(T)−1∫0TeATtCTz(t)dt.

同一运动,两个传感器 ​

令 A=(0100),正常传感器 C=(1 0) 测位置。则

CeAt=(1 t),O=I2,Wo(T)=(TT2/2T2/2T3/3),det⁡Wo(T)=T4/12.

校正输出 z(t)=a+bt 的截距和斜率分别为初始位置与速度。在 T=1 时,若测得 z(t)=2−t,则

∫01(1t)z(t)dt=(3/22/3),Wo(1)−1=(4−6−612),

相乘确得 x0=(2,−1)T。

若改测速度,C=(0 1),则 CA=0、Wo(T)=diag(0,T),不可观空间为 span(e1)。任意初始位置差都不改变速度输出。保留执行器 B=e2 时系统仍然可控,说明可控不推出可观;反过来取 B=e1,C=(1 0),系统可观却不可控。

重建公式要求准确、连续时间的输出轨迹。噪声或模型误差会传入重建结果,即使 Gramian 正定,接近奇异也会放大这些误差。有限采样下须按采样后的模型重新判断可观测性,再用相应的离散重建或估计方法。

推论与应用

(A,C) 可观等价于转置输入对 (AT,CT) 可控:二者的秩矩阵互为转置。这是输入与输出角色交换的代数对偶,不表示实际执行器和传感器可以互换。

Kalman 滤波器处理带噪声状态的条件估计;可观测性提供“哪些初态差能进入数据”的结构问题。LQR使用的可检测性则允许不可观方向存在,只要求它们自身稳定;那里观察矩阵来自状态代价,未必是物理传感器 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 都成立。
关系图谱21 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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