Equivalence of Continuous-Time Markov Chains and Linear Dynamical Systems

技术摘要：连续时间马尔可夫链与线性动力系统的等价性

问题陈述
本短篇笔记探讨了连续时间线性动力系统与连续时间马尔可夫链（CTMC）之间的结构关系。虽然在离散时间设定下，已知一个 $d$ 状态马尔可夫链的动力学受一个维度至多为 $d-1$ 的线性动力系统支配（反之亦然），作者研究了这种等价性是否在连续时间领域同样成立。具体而言，本文旨在建立以下等价关系：一个 $d$ 状态 CTMC 的演化可以简化为一个维度为 $d-\ell$ 的线性动力系统（其中 $\ell$ 是零特征值的重数），并且反之，一个维度为 $d-1$ 的线性动力系统可以被嵌入到一个 $d$ 状态的 CTMC 中。

方法论
本文采用线性代数技术，特别是谱分解和矩阵指数化，来推导这种等价性。该方法论依赖于以下核心步骤：

1. 生成矩阵的谱分析： 作者首先建立了一个关于随机生成矩阵（$A$）的引理。证明了 $A$ 的 0 始终是其特征值，且其代数重数与几何重数相等（均为 $\ell$），而所有其他特征值的实部均严格为负。这是通过将 $A$ 与离散随机矩阵 $M = I + aA$ 联系起来并应用 Perreron-Frobenius 定理得出的。
2. 子空间分解： 将状态空间 $\mathbb{C}^d$ 分解为两个正交子空间：$V_0$（由特征值 0 对应的特征向量所张成）和 $W$（由所有非零特征值的广义特征向量所张成）。
3. 矩阵变换： 利用变换矩阵 $R = [P \mid Q]$（其中 $P$ 和 $Q$ 的列分别构成 $V_0$ 和 $W$ 的基），将生成矩阵 $A$ 转换为块对角形式。这使得矩阵指数 $\exp(At)$ 可以被分离为一个常数项（与稳态相关）和一个动态项（与瞬态行为相关）。
4. 嵌入构造： 对于逆向过程，本文构造了一个特定的嵌入矩阵 $Q$ 和标量 $\rho$，将一个低维生成矩阵 $C$ 映射到一个高维随机生成矩阵 $A$，从而确保生成的系统满足 CTMC 的性质（非负非对角项且列和为零）。

主要贡献与结果
本文提出了两个主要定理来形式化这种等价性：

• 定理 1（归约）： 给定一个分布 $v$ 和一个特征值 0 的重数 $\ell \ge 1$ 的随机生成矩阵 $A \in \mathbb{Q}^{d \times d}$，可以计算出一个归约后的生成矩阵 $B \in \mathbb{Q}^{(d-\ell) \times (d-\ell)}$、一个嵌入矩阵 $Q$、一个稳态分布 $s$ 以及一个向量 $u$，使得对于所有 $t \ge 0$：
  $$\exp(At)v = s + Q \exp(Bt)u$$
  至关重要的是，$Q$ 的选择独立于初始分布 $v$。项 $Q \exp(Bt)u$ 在 $t \to \infty$ 时收敛于 0，从而仅留下稳态 $s$。

• 定理 2（嵌入）： 给定一个向量 $u \in \mathbb{Q}^{d-1}$ 和一个生成矩阵 $C \in \mathbb{Q}^{(d-1) \times (d-1)}$，可以计算参数 $\rho, \eta$，一个分布 $v$，以及一个随机生成矩阵 $A \in \mathbb{Q}^{d \times d}$，使得对于所有 $t \ge 0$：
  $$\exp(At)v = s + \eta e^{-\rho t} \cdot Q \exp(Ct)u$$
  此处，$s$ 和 $Q$ 是固定的选择（$s = \mathbf{1}_d/d$ 且 $Q$ 为特定的块矩阵）。标量 $\rho$ 被选择得足够大，以确保构造出的生成矩阵的非对角项保持非负。

意义与主张
作者明确指出，本笔记的目的是“记录”已知于离散时间线性动力系统（特别参考了文献 [1] 和 [2] 中的工作）的核心结果在连续时间设定下同样成立。本文并未提出新的应用或实验验证。相反，其意义在于提供了一个严谨的数学证明，证明了这种结构上的等价性：即一个 $d$ 状态连续时间马尔可夫链的动力学在根本上受一个维度至多为 $d-1$ 的线性动力系统支配，并且任何此类低维系统都可以被嵌入到马尔可夫链中。其证明过程镜像了所引用离散时间文献中的逻辑，并针对生成矩阵和矩阵指数的特性进行了适配。