Skip to content

定理Theorem

Markov 链平稳分布收敛定理

施工中 · 内容逐步修订

Convergence to stationarity for Markov chains · Markov chain convergence theorem

不可约、正常返且非周期的可数状态 Markov 链,其逐时刻转移分布收敛到唯一平稳分布。

形式陈述 ​

设 P 是可数状态空间 S 上不可约、正常返且非周期 Markov 链的转移核,π 是其唯一平稳分布。则对任意状态 i,j∈S,

Pn(i,j)⟶π(j).

更强地,从任意初始分布 μ 出发,边际分布满足总变差距离 ‖μPn−π‖TV→0。在离散度量下,这也蕴含依分布收敛到平稳律。有限不可约链自动正常返,所以有限情形只需另加非周期性。

这条结论只描述第 n 个时刻的分布,不是沿单条样本路径求平均。后者由时间平均遍历定理控制,并不要求非周期性。

总变差距离在这里取约定

‖ν−π‖TV=12∑j∈S|ν(j)−π(j)|.

从逐点收敛到总变差收敛还用到了极限质量为 1:先选有限集 F 使 π(F)>1−ε,在 F 上逐项取极限,补集上的误差由 ν(Fc)+π(Fc) 控制。若逐点极限丢失概率质量,这一步就失败。对任意初始分布的结论则先在初始状态的有限集上求和,再用该集之外的初始概率控制尾项。

用相遇构造证明收敛 ​

耦合法把结论归结为两条链最终相遇。先考虑两条相互独立的 P 链,其联合转移为 P⊗P。不可约且非周期意味着对每个固定状态对 i,j,充分大的 n 都有 Pn(i,j)>0:在一个状态选取有限组最大公因数为 1 的返回长度,它们的非负整数和包含所有充分大的整数。具体说,模其中一个返回长度取余,这些生成元生成全部剩余类;每类选一个非负组合代表,再加该长度的倍数,就覆盖所有足够大的整数。最后拼接到达和离开的路径。因此两个坐标可以选取同一个足够大的步数,乘积链不可约。

乘积链有平稳概率 π⊗π,由不可约链的正常返判据,它几乎必然到达任意指定的 (a,a),从而两条独立链的首次相遇时间有限。现在构造另一条联合链:相遇前独立更新,相遇后从相同状态使用同一个 P 转移,始终保持相等。每个联合转移行的两个边缘仍是对应的 P 行,所以两个坐标各自都是所需的 Markov 链;相遇前的路径律也与独立乘积链一致。

取 X0=i、Y0∼π,令 τ 为首次相遇时间。每个确定时刻都有 Yn∼π,而耦合不等式给出

‖Pn(i,⋅)−π‖TV≤P(Xn≠Yn)≤P(τ>n)⟶0.

这里不需要、也没有断言随机相遇时刻的状态服从 π。最后按初始分布混合,并用前述有限集截断,就得到任意 μ 的结论。

直觉

不可约性排除永久隔离的多个状态类,正常返性让平稳概率质量能够归一化,非周期性则打散“只能在某些同余时刻回到这里”的相位。三者合在一起,初始状态留下的信息才会随步数消散,单步边际逐渐只剩平稳分布。

周期性并不阻止长期访问比例稳定,却会让每个固定时刻仍记得当前相位。时间平均把多个相位混在一起,逐时刻收敛没有这种平均作用,因此假设更强。

例子与边界

先看可直接算出速度的二状态链:

P=(3/41/41/21/2),π=(2/3,1/3).

若 qn=P(Xn=0),则 qn+1=1/2+qn/4,因而

qn−2/3=4−n(q0−2/3).

从状态 0 出发,总变差距离恰为 (1/3)4−n。这里速率来自额外计算,不能由一般定理直接读出。

二状态确定性交替链 0↔1 不可约,平稳分布为 (1/2,1/2),但从 0 出发时分布在 δ0 与 δ1 之间交替,永不趋于 π。给每个状态加入正概率自环后周期被打破,边际分布才趋稳。

整数上的简单对称随机游走不可约且常返,却是零常返,没有平稳概率分布,不能套用定理。可约链则可能有多个闭类;此时极限若存在,也可能依赖初始分布落入哪个类。无限状态情形还需核对正常返性,不能把有限链直觉直接外推。

推论与应用

混合时间量化有限链达到给定总变差精度需要多少步;本页定理只给极限存在,不自动给收敛速率。可逆链可进一步用谱隙或 conductance 控制速度,耦合法则通过两条链相遇来给出概率界。

MCMC 的 burn-in 讨论依赖逐时刻分布趋稳,而最终估计量通常是轨道平均,还同时需要遍历定理。二者分别控制初始偏差和长期平均,不能用“链看起来稳定”把其中任一项省略。

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

拖动节点调整位置。

显示关系

显示:依赖

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

使用的工具