生物医学网络科学Biomedical Network Science
网络可控性与稳定性08 / 17

可观测性与观测节点

阅读
约 3 分钟

在控制一个系统之前,先要知道它处在状态空间的什么位置。细胞内成千上万种代谢物、蛋白不可能同时测量,临床上也只能采集少数指标。可观测性问:从少数可测变量(输出)的时间序列,能否唯一推断系统的全部内部状态(等价地,推断初始状态 𝒙(0)\mathbf x(0))?

可观测性:系统 𝒙̇=A𝒙+B𝒖\dot{\mathbf x}=A\mathbf x+B\mathbf u,𝒚=C𝒙\mathbf y=C\mathbf x,若由一段时间内的 𝒚(t)\mathbf y(t) 与 𝒖(t)\mathbf u(t) 能唯一确定 𝒙(0)\mathbf x(0),则称 (A,C)(A,C) 可观测。线性判据为可观测性矩阵满秩: 𝒪=(CCA⋮CAN−1),rank⁡𝒪=N.\mathcal O=\begin{pmatrix}C\\CA\\\vdots\\CA^{N-1}\end{pmatrix},\qquad \operatorname{rank}\mathcal O=N .

对偶原理(Kalman):(A,C)(A,C) 可观测 ⇔\iff (A⊤,C⊤)(A^\top,C^\top) 可控。在网络语言中,网络的可观测性等于把每条边反向后网络的可控性。

因此对线性系统,最少观测节点数 = 反向网络的最少驱动节点数,可直接用最大匹配求得;可观测性 Gramian Wo=∫0∞eA⊤τC⊤CeAτdτW_o=\int_0^\infty e^{A^\top\tau}C^\top Ce^{A\tau}d\tau 则刻画状态估计的难易(观测噪声被放大的程度)。

非线性系统:推断图与根强连通分支

对非线性系统,对偶性不再成立,而精确判据(基于李导数的秩条件)只能用于很小的系统,穷举观测组合需检验约 2N2^N 种可能。Liu、Slotine 与 Barabási(2013,PNAS)提出图论方法:

  1. 推断图:若变量 xjx_j 出现在 xix_i 的微分方程中(∂fi/∂xj≢0\partial f_i/\partial x_j\not\equiv0),画一条边 xi→xjx_i\to x_j,表示测量 xix_i 能获得关于 xjx_j 的信息(推断图是系统影响图的反向图)。
  2. 强连通分支分解:把推断图分解为极大强连通分支(SCC);同一 SCC 内的节点彼此“可推断”。
  3. 选择观测节点:没有入边的 SCC 称为根 SCC,每个根 SCC 至少放一个观测节点。

观测节点的必要条件:对任意(线性或非线性)动力学,监测每个根 SCC 中至少一个节点是可观测的必要条件;根 SCC 的个数是最少观测节点数的严格下界。任何状态观测器若不监测这些节点都会失败。

这个条件对线性系统一般不充分:结构对称性(扩张)也会造成不可观测,需要更多观测节点(代码示例:200 节点随机网络中根 SCC 有 44 个,而线性对偶给出的最少观测节点是 68 个;只放 44 个时可观测性矩阵的秩为 168/200)。但对生化反应等大型非线性系统,状态变量之间的对称性极少,Liu 等(2013)发现根 SCC 给出的观测集合往往既必要又充分。

目标可观测性与生物标志物

很多时候只需推断少数目标变量(例如与疾病相关、但难以直接测量的代谢物)。目标可观测性要求:观测节点到目标节点在推断图中存在有向路径(测量 xsx_s 能推断 xtx_t 的前提是 xs→⋯→xtx_s\to\cdots\to x_t)。据此可以为每个目标找出最合适的可测变量——这为生物标志物设计提供了系统的方法:标志物不必是疾病相关分子本身,而是在网络中能“看到”它的分子。

参数可辨识性与网络重建。把未知参数(如反应速率常数)看作导数为零的额外状态变量,参数辨识就成了扩展系统的可观测性问题;若参数编码了连接是否存在,问题就推广为网络重建。这把可观测性与基因调控网络推断与共表达分析中的网络推断联系在一起。

推断图依赖已知的动力学方程结构(谁出现在谁的方程里)。用统计相关构建的网络(共表达、功能连接)没有方向和机制含义,不能直接当作推断图使用;在这类网络上计算“观测节点”只是启发式的。