Skip to content

模型Model

Markov 链

Markov chain

未来条件分布在给定当前状态后与更早历史无关的随机过程。

形式陈述 ​

本页以离散时间、有限或可数状态空间为默认范围;有限链是可数链的特例。一般可测状态空间版本使用概率核与条件期望,不能把矩阵公式无条件照搬。设 S 为非空有限或可数集,配备全体子集构成的 σ-代数。取值于 S 的随机过程 (Xn)n≥0 称为 Markov 链,若对一切 n 与一切使条件事件具有正概率的 i0,…,in−1,i,j∈S,条件概率满足

P(Xn+1=j∣X0=i0,…,Xn−1=in−1,Xn=i)=P(Xn+1=j∣Xn=i).

若存在固定的行随机数组 P=(pij)i,j∈S,使每当 P(Xn=i)>0 时都有 pij=P(Xn+1=j∣Xn=i),便称链为时间齐次。这里 pij≥0 且 ∑jpij=1;一个给定初始分布下从未访问的状态,其转移行不能由这些条件概率唯一恢复,仍须作为模型的一部分指定。

每一行给出下一步的概率分布,因而 K(i,A)=∑j∈Apij 是离散空间上的概率核。有限 S 时 P 是通常的矩阵;可数无限时,下面的乘法记号表示非负可数求和,不是有限维矩阵运算。初始分布 μ 经 n 步演化为 μPn,多步转移满足 Chapman–Kolmogorov 方程

Pm+n=PmPn.

初始分布与转移矩阵还确定每条有限路径的概率:

P(X0=i0,…,Xn=in)=μ(i0)∏k=0n−1pikik+1.

这条乘积公式把一步条件分布组装成整条路径分布;反过来,对中间状态求和便得到矩阵幂中的多步转移。

在一般可测状态空间 (S,S) 上,时间齐次版本改写为

P(Xn+1∈A∣X0,…,Xn)=K(Xn,A),A∈S,

其中 K 是概率核,左侧是指标变量关于历史信息的条件期望;对每个固定 A,等式按几乎处处理解。以下转移数组、状态可达性和返回时间术语均按有限或可数版本理解。

范围与术语约定 ​

本库统一使用:irreducible 译为不可约,aperiodic 译为非周期,recurrent/transient 分别译为常返/暂留,positive recurrent/null recurrent 分别译为正常返/零常返。在可数状态空间中,不可约表示任意两状态可经某个正概率多步路径互达;状态周期是正概率返回步数集合的最大公因数,周期为 1 即非周期;常返状态若首次正返回时间期望有限,则为正常返。这些性质回答不同问题,其中正常返蕴含常返;不能把它们合并成没有量词的“遍历性”。

一般状态空间也使用不可约、Harris 常返和正常返等词,但定义需要参考测度、可达集合与返回集合等额外数据。除非页面明确声明一般状态空间版本,本库相邻 Markov 链条目的这些术语都采用上述有限或可数状态空间含义。

直觉

Markov 性要求当前状态已经概括预测未来所需的全部历史信息。过去没有消失;与未来有关的部分被压缩进状态变量。因此“无记忆”说的是给定当前状态后的条件分布,不表示相邻状态独立,也不保证状态停留时间服从无记忆分布。

Markov 性依赖建模方式。同一个现象若状态取得太粗,遗漏信息仍会影响未来,Markov 性便失败;把必要记忆并入状态后,往往可以恢复。只保留正概率转移的支撑,可以把它画成无外部输入的状态机;这个对应只记录允许的路径,丢掉边上的概率后便无法计算路径概率。

初始分布与转移核共同决定整条路径分布。有限状态时,矩阵编码一步规律,矩阵乘法汇总中间状态,矩阵幂给出多步转移。可数链使用相应的无穷求和,一般状态空间则用概率核的积分。

例子与边界

整数轴上的简单随机游走:pi,i+1=pi,i−1=1/2,是时间齐次 Markov 链。两步转移一算便知:从 i 出发以 1/4 到达 i±2、以 1/2 回到 i,这正是 P2 的第 i 行,与 Chapman–Kolmogorov 方程一致。

Markov 性依赖状态的选取。只记录天气"晴/雨"时,若真实机制是"连续晴天越久,明天转雨概率越大",那么这个二状态过程不满足 Markov 性——预测明天需要的不止是今天;把状态扩充为(今日天气,已连续天数)后 Markov 性即可恢复。一般的教训是:对 Markov 链的状态做归并(取链的函数)通常会破坏 Markov 性,因为压缩可能把必要的记忆挤出状态之外。

还应把 Markov 性与独立性分开:随机游走中 Xn+1 与 Xn 高度相关,Markov 性只约束"给定现在后过去失效",在已经满足 Markov 性的链中,若各转移行相同,下一步便独立于整个过去;但相邻变量独立本身不足以推出 Markov 性。例如取独立公平比特 A,B,令 X0=A,X1=B,X2=A⊕B,三者两两独立,给定完整的前两步却能确定第三步。另外,定义中的条件式只对正概率历史提出要求;齐次性是额外假设,转移规律随时间变化的非齐次链同样存在,只是许多结论需要在齐次假设下陈述。

若把两个状态依次记为 0,1,并取转移矩阵

P=(0.90.10.40.6)

则两步转移概率不是把矩阵元素逐项平方,而是计算

P2=(0.850.150.600.40).

例如从状态 0 出发,两步后仍在 0 的概率 0.85 来自两条互斥路径:连续留在 0,概率 0.92;或先到 1 再返回,概率 0.1×0.4。矩阵乘法正是在自动汇总所有中间状态,因此 Pn 才能编码 n 步转移。

从状态 0 出发,两步后返回 0 的概率来自 0→0→0 与 0→1→0 两条互斥路径之和,而不是对转移矩阵逐项平方。
推论与应用

条件概率给出 Markov 性的精确定义,矩阵与谱分析则控制长期行为。状态空间先由沟通类与不可约性分解,随后用常返、暂留与周期性判断质量是否逃逸以及边缘分布是否振荡;平稳分布刻画不变概率质量,细致平衡则提供一种更强的逐边对称结构。这些概念各自回答不同问题,不应在 Markov 性定义里合并成笼统的“遍历”。

极限问题应沿不同出口继续。时间平均遍历定理研究一条长轨道上的经验平均;平稳分布收敛定理研究第 n 步边际分布是否趋稳;混合时间再量化有限步距离。三者的对象、假设和量词不同,不能用“链最终稳定”一句话互相替代。

把步数换成连续时间后,连续时间链还需指定指数停留速率,并检查非爆炸性。Kolmogorov 方程由生成矩阵计算转移概率;速率有共同有限上界时,均匀化再把它表示成 Poisson 时钟驱动、允许自环的离散链。

若转移时刻受到删失,转移概率可以先从各状态风险集估计局部风险,再按时间顺序组合。Aalen–Johansen 估计给出这一乘积矩阵,并区分 Markov 条件转移概率与更宽情形下的边缘状态占据概率。

布尔立方体上逐坐标随机翻转给出一个可显式对角化的例子。布尔噪声算子把一次转移写成条件平均,并以奇偶字符为特征函数,特征值为 ρ|S|;连续施加两次转移使相关参数相乘。它展示了谱结构怎样区分不同阶数的信息衰减。

参考资料
  • J. R. Norris, Markov Chains, Cambridge University Press, 1997,Ch. 1, transition matrices and the Markov property。
  • Rick Durrett, Probability: Theory and Examples, 5th ed., Cambridge University Press, 2019,Ch. 5, discrete-time Markov chains。
关系图谱27 个相邻概念 · 3 类关系

拖动节点调整位置。

显示关系

显示:依赖

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

上位 / 更一般

下位 / 直接特例

暂未标注直接特例。

类型化关系