Skip to content

混合时间

Mixing time · Markov chain mixing time

Markov 链从最坏初态接近平稳分布到指定总变差精度所需的步数。

形式陈述

P 是有限状态 Markov 链的转移矩阵,π 是一个指定的平稳分布。从状态 x 出发,t 步后的分布为 Pt(x,)。定义相对于 π 的最坏初态总变差距离

d(t)=maxxPt(x,)πTV,

以及精度 ε(0,1) 下的混合时间

tmix(ε)=min{t0:d(t)ε}.

若集合为空,就约定 tmix(ε)=。不可约且非周期的有限链具有唯一平稳分布,并从任意初态收敛到它,因此这是保证混合时间有限的标准条件;不满足这些条件的链仍可代入定义,但结果可能为无穷。本文把

tmix=tmix(1/4)

作为默认约定。对标准遍历链,任何固定 ε(0,1/2) 得到的时间尺度仅相差由阈值决定的常数或对数因子,但具体数值仍依赖约定。若从初始分布 μ 出发,则相应距离是 μPtπTV;最坏初态定义保证对所有 μ 同时成立。

直觉

混合不是链在某一刻“到达”平稳分布。有限步后,Pt(x,) 往往仍与 π 不完全相等;混合时间问的是它们何时已经近到任何事件测试都难以区分,而且这个保证必须覆盖最不利的起点。

转移每执行一步都会遗忘一部分初态信息。若不同起点产生的分布逐渐靠拢,平稳分布便成为共同参照。障碍也由此清楚:多个闭沟通类会永久保留“起点在哪一类”的信息;周期链会保留时间相位;很窄的瓶颈则让概率质量跨区域移动得很慢。

例子与边界

二状态链

P=(1aab1b),0<a,b<1,

具有平稳分布 π=(b/(a+b),a/(a+b))。其非平凡特征值为 1ab,从任一初态到平稳分布的差按 |1ab|t 衰减;当 a+b 很小时,链长时间停留原状态,混合变慢。

确定交替的两状态链具有平稳分布 (1/2,1/2),却从状态 0 出发在奇偶时刻来回跳动,总变差距离不趋于零。加入以 1/2 概率原地停留的 lazy 版本可消除周期,但这产生新转移矩阵,混合步数也应按新链计算。时间平均趋向平稳分布并不能替代逐时刻边缘分布的混合。

推论与应用

总变差距离固定了“接近”的意义;不可约性非周期性排除结构性障碍。可逆链还能借助谱隙、conductance 等工具给出混合上界,而耦合法通过两条链的相遇概率直接控制 d(t)

在 MCMC 中,混合时间决定需要丢弃多少初始步骤以及有限运行带来多大偏差。它是最坏情形分布保证,不等同于某个观测量看似稳定,也不能由一条模拟轨迹的视觉平稳性替代。

参考资料
  • David A. Levin, Yuval Peres, and Elizabeth L. Wilmer, Markov Chains and Mixing Times, 2nd ed., American Mathematical Society, 2017,Ch. 4, total variation and mixing time。
  • J. R. Norris, Markov Chains, Cambridge University Press, 1997,convergence to equilibrium for discrete-time chains。