Skip to content

定理Theorem

连续时间链的 Kolmogorov 方程

Kolmogorov forward and backward equations

由短时间生成矩阵和半群性质推出前向与后向方程,明确行向量约定和无限状态边界。

形式陈述 ​

局部故障速率是每小时一个数,而我们需要的是运行 t 小时后的故障概率。Kolmogorov 方程把局部速率积累成整个转移矩阵。

设有限状态连续时间 Markov 链生成矩阵为 Q,P(t)=(Pij(t)),采用行概率向量约定,所有导数逐个矩阵元素理解;t=0 处取右导数。转移半群满足

P(0)=I,P(s+t)=P(s)P(t),P′(t)=P(t)Q=QP(t).

P′(t)=P(t)Q 称为前向方程,P′(t)=QP(t) 称为后向方程,逐列应用常系数线性 ODE的初值唯一性,解为 P(t)=etQ。若初始分布为行向量 p(0),则 p(t)=p(0)P(t),所以 p′(t)=p(t)Q。

有限状态是本页完整陈述的范围;可数状态下需要检查无穷求和与微分交换、非爆炸及方程解的唯一性,不能只把有限矩阵改成无穷矩阵便自动得到所有结论。

直觉

前向方程固定起点,观察最后一小段:时刻 t 在状态 k 的概率乘以从 k 流入终点 j 的速率,再对 k 求和。后向方程则先看起点的第一小段跳向哪里,之后再运行剩余时间。它们分别沿路径的尾部和头部分解同一概率。

证明也反映这一顺序。由 P(h)=I+hQ+o(h) 和半群性质,

P(t+h)−P(t)h=P(t)P(h)−Ih→P(t)Q.

交换 P(t) 与 P(h) 的时间分解次序,便得后向式。对 t>0,再用 P(t)−P(t−h)=P(t−h)(P(h)−I) 得到相同的左导数;有限随机矩阵的有界性和 P(h)→I 保证这里 P(t−h)→P(t)。这是导数的极限计算,不是任意矩阵之间都可交换;P(t) 与 Q 可交换来自同一个齐次半群。

例子与边界

对工作/维修链,令故障速率 α>0、修复速率 β>0,从工作态出发。设 d(t)=P(X(t)=D)。流入维修态的概率速率为 α(1−d(t)),流出为 βd(t),所以

d′(t)=α−(α+β)d(t),d(0)=0.

乘积分因子 e(α+β)t 后积分,得到

d(t)=αα+β(1−e−(α+β)t).

它在 t=0 的导数为 α,符合局部故障速率;长期趋于 α/(α+β),符合停留时间占用比。完整矩阵为

P(t)=1α+β(β+αe−ktα(1−e−kt)β(1−e−kt)α+βe−kt),k=α+β.

每行和为一、非对角元非负,都是解完后应检查的概率约束。

短时间近似 I+hQ 未必在大步长下合法:若 h>1/α,工作态对角元 1−αh 会为负。它是局部 Euler 近似,不是精确转移矩阵。精确的 etQ 则保持非负性与行和。

非齐次链 Q(t) 的前向方程为 ∂tP(s,t)=P(s,t)Q(t)。若不同时刻的生成矩阵不交换,不能写成 exp⁡(∫stQ(u)du);时间顺序会影响结果。

推论与应用

有限链的不变分布可由 πQ=0、∑iπi=1 求得,它保证 πetQ=π。这只是在求稳态;从任意初态是否收敛仍需链的结构条件。数值计算若想保留概率非负性并给出截断误差,可采用均匀化的 Poisson 加权级数。

参考资料
关系图谱16 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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