形式陈述
局部故障速率是每小时一个数,而我们需要的是运行 小时后的故障概率。Kolmogorov 方程把局部速率积累成整个转移矩阵。
设有限状态连续时间 Markov 链公理库连续时间 Markov 链Continuous-time Markov chain · CTMC以指数停留时间和跳转概率构造连续时间链,并将局部速率与非爆炸条件分开。生成矩阵为 ,,采用行概率向量约定,所有导数逐个矩阵元素理解; 处取右导数。转移半群满足
称为前向方程, 称为后向方程,逐列应用常系数线性 ODE公理库线性常微分方程组Linear system of ordinary differential equations形如 x′=A(t)x+b(t) 的向量值一阶线性方程组。的初值唯一性,解为 。若初始分布为行向量 ,则 ,所以 。
有限状态是本页完整陈述的范围;可数状态下需要检查无穷求和与微分交换、非爆炸及方程解的唯一性,不能只把有限矩阵改成无穷矩阵便自动得到所有结论。
直觉
前向方程固定起点,观察最后一小段:时刻 在状态 的概率乘以从 流入终点 的速率,再对 求和。后向方程则先看起点的第一小段跳向哪里,之后再运行剩余时间。它们分别沿路径的尾部和头部分解同一概率。
证明也反映这一顺序。由 和半群性质,
交换 与 的时间分解次序,便得后向式。对 ,再用 得到相同的左导数;有限随机矩阵的有界性和 保证这里 。这是导数公理库导数Derivative函数增量比在步长趋零时的极限。的极限计算,不是任意矩阵之间都可交换; 与 可交换来自同一个齐次半群。
例子与边界
对工作/维修链,令故障速率 、修复速率 ,从工作态出发。设 。流入维修态的概率速率为 ,流出为 ,所以
乘积分因子 后积分,得到
它在 的导数为 ,符合局部故障速率;长期趋于 ,符合停留时间占用比。完整矩阵为
每行和为一、非对角元非负,都是解完后应检查的概率约束。
短时间近似 未必在大步长下合法:若 ,工作态对角元 会为负。它是局部 Euler 近似,不是精确转移矩阵。精确的 则保持非负性与行和。
非齐次链 的前向方程为 。若不同时刻的生成矩阵不交换,不能写成 ;时间顺序会影响结果。
推论与应用
有限链的不变分布可由 、 求得,它保证 。这只是在求稳态;从任意初态是否收敛仍需链的结构条件。数值计算若想保留概率非负性并给出截断误差,可采用均匀化公理库连续时间链的均匀化Uniformization · Randomization method for CTMC · 均匀化用共同 Poisson 时钟和自环表示有界速率连续时间链,并以 Poisson 尾控制数值截断误差。的 Poisson 加权级数。
参考资料