路径: 本科数学/ 进阶数学/ 矩阵分析与谱图论/随机过程与鞅论:马尔可夫链、布朗运动与鞅
30

随机过程与鞅论:马尔可夫链、布朗运动与鞅

Stochastic Processes & Martingales · 随机世界怎么演化

【章首引子】明天会不会下雨、用户会不会点广告、股价怎么走——这些都是"随时间演化的随机量",叫随机过程。本章抓三个最有用的:马尔可夫链(下一步只看现在)、泊松过程(随机事件计数)、布朗运动(连续随机游走),再加一个贯穿始终的工具——鞅(公平赌博)。强化学习、Diffusion、MCMC、随机梯度下降都在这套语言里。

卡必背公式 / 记忆口诀
马氏链:P 转移矩阵,π=πP 平稳
泊松:N(t)~Poisson(λt),独立增量
布朗运动:W(t)~N(0,t),独立平稳增量
鞅:E(M_{n+1}|Fₙ)=Mₙ
口诀:马氏只看现在,泊松数事件,布朗瞎走,鞅是公平赌

① 是什么:四类随机过程

啥随机过程家族

① 随机过程定义:一族随机变量 {X(t) : t∈T},t 是时间。离散时间 X₀,X₁,X₂,…;连续时间 X(t)。平稳过程:统计性质不随时间平移改变;独立增量:不重叠时间区间上的增量相互独立。

② 马尔可夫链:离散状态 + 离散时间,核心是马尔可夫性——下一步只依赖当前状态,与历史无关:P(X_{n+1}=j | Xₙ=i, X_{n−1},…)=P(X_{n+1}=j|Xₙ=i)=p_{ij}。排成转移矩阵 P=(p_{ij}),n 步转移是 Pⁿ。

③ 平稳分布:若 πP=π(即 π 是 P 的左特征向量、特征值 1)且 π 非负归一,则 π 是平稳分布。不可约(任意两状态互通)+ 非周期 ⇒ 遍历:无论初始分布如何,Pⁿ→π。吸收态:一旦进入就出不来(p_{ii}=1)。

④ 泊松过程:N(t)=[0,t] 内事件数。N(t)~Poisson(λt),即 P(N(t)=k)=e^{−λt}(λt)ᵏ/k!。独立增量、无记忆(等待时间指数分布 Exp(λ))。

⑤ 布朗运动 W(t):连续轨道、W(0)=0、独立平稳增量、W(t)−W(s)~N(0,t−s)。它是连续时间的"随机游走极限",是伊藤积分与随机微分方程(SDE)的基石。

图 30.1:布朗运动三条样本路径(都是 W(t),但涨落剧烈不同)
t W 路径1 路径2 路径3(虚线)

② 怎么想到的

思解题心法

1

状态随时间跳、只看现在 → 马尔可夫链:写转移矩阵 P,求平稳分布 π=πP。

2

数"事件发生了几次" → 泊松过程:N(t)~Poisson(λt),间隔时间独立指数分布。

3

连续时间随机游走 → 布朗运动:方差随时间线性增长 Var W(t)=t。

4

"无偏估计/公平游戏" → 鞅:未来期望等于当下,E(M_{n+1}|Fₙ)=Mₙ。

证核心定理

定理:① 马氏链遍历定理:不可约 + 非周期 + 有限状态 ⇒ 存在唯一平稳分布 π,且 Pⁿ(i,j)→π(j)。② 平稳分布方程:πⱼ=Σᵢ πᵢ p_{ij}(细致平衡 πᵢp_{ij}=πⱼp_{ji} 是充分条件)。③ 泊松过程:P(N(t)=k)=e^{−λt}(λt)ᵏ/k!;等待时间 T₁~Exp(λ),P(T₁>t)=e^{−λt}(无记忆)。④ 布朗运动:W(t)~N(0,t);W(t)−W(s)~N(0,t−s) 与过去独立。⑤ 鞅定义:E|Mₙ|<∞ 且 E(M_{n+1}|Fₙ)=Mₙ;可选抽样定理:对有界停时 τ,E(M_τ)=E(M₀)——公平游戏在任何"聪明"的停时策略下仍公平。

推导思路(平稳分布):πP=π 表示"按 π 随机抽一个状态,再走一步,分布还是 π"。若链遍历,长期频率收敛到 π,所以 π 也是极限分布。细致平衡 πᵢp_{ij}=πⱼp_{ji} 直接相加即 πP=π,是找 π 的常用捷径。

推导思路(鞅直觉):鞅的"鞅停住"性质:赌徒用任何策略(见好就收、追跌加仓)都不能把公平游戏变成正期望——E(M_τ)=E(M₀)。这是 MCMC、赌博系统、期权定价公平性的数学基础。

伊藤积分直觉:布朗运动路径处处连续但处处不可导,普通积分 ∫f(W)dW 定义不了(dW 抖动太厉害)。伊藤积分用"右端点"采样把它严格化,得到 dW²=dt 的神奇规则,是 SDE dX=f dt+g dW 的根基。

③ 例题

例1:两状态马氏链的平稳分布(晴雨模型)
【审题】P=[[0.7,0.3],[0.4,0.6]](晴→晴0.7,晴→雨0.3;雨→晴0.4,雨→雨0.6)。
思路:解 πP=π,π₁+π₂=1。
逐步:π₁=0.7π₁+0.4π₂ ⇒ 0.3π₁=0.4π₂ ⇒ π₁=(4/3)π₂。又 π₁+π₂=1 ⇒ (7/3)π₂=1 ⇒ π₂=3/7≈0.429,π₁=4/7≈0.571。长期:57% 晴天、43% 雨天。
例2:泊松过程 λ=2/小时,2 小时内到 5 个事件的概率
【审题】N(2)~Poisson(4)。
思路:套 Poisson 公式。
逐步:P(N(2)=5)=e^{−4}·4⁵/5! = e^{−4}·1024/120 ≈ 0.0183×8.533 ≈ 0.156。
例3:布朗运动 W(1) 与 W(2) 的分布
【审题】W(t)~N(0,t),独立增量。
思路:直接用定义。
逐步:W(1)~N(0,1);W(2)~N(0,2),标准差 √2≈1.414。增量 W(2)−W(1)~N(0,1) 且与 W(1) 独立。直觉:时间越长,随机游走扩散得越宽(方差∝t)。
随机过程核心公式 马氏链:Pⁿ = n 步转移;平稳 π=πP;不可约+非周期 ⇒ Pⁿ→π
泊松:N(t)~Poisson(λt),P(N(t)=k)=e^{−λt}(λt)ᵏ/k!;等待时间 Exp(λ)
布朗运动:W(t)~N(0,t),独立平稳增量,W(0)=0,连续路径
鞅:E(M_{n+1}|Fₙ)=Mₙ;可选抽样 E(M_τ)=E(M₀)(有界停时)

④ AI 落点

算法
AI 落点 1:MCMC 采样

Metropolis-Hastings 构造一条平稳分布为目标后验 π 的马氏链,长时间跑样本即来自 π——贝叶斯推断的核心引擎。

算法
AI 落点 2:Diffusion 模型

前向加噪是布朗运动/高斯扩散过程 Xₜ=√(t)ε;反向学 score 去噪。整条链就是随机过程理论的工程化。

稳定性
AI 落点 3:强化学习 MDP

马尔可夫决策过程 = 马氏链 + 动作 + 奖励;贝尔曼方程就是在马氏链上做期望。Q-Learning 是随机逼近。

优化
AI 落点 4:SGD 的鞅噪声

随机梯度 = 真梯度 + 鞅差噪声;分析 SGD 收敛要用到鞅的大数定律与 Azuma 不等式。

⑤ 延展

展知识衔接地图

往本科走:概率论的期望、条件期望、大数定律是本章的语言。

往算法走:马氏链 → MCMC/PageRank;布朗运动 → SDE/扩散模型;鞅 → 随机优化与金融数学。

思维陷阱

以为平稳分布一定存在。可约链(有多个闭类)或周期链不一定有唯一平稳分布;必须不可约+非周期才有遍历极限。

把独立增量当平稳增量。独立增量说的是"不重叠区间增量独立",平稳增量说的是"增量分布只依赖时间差"——布朗运动两者都有,一般过程未必。

忽略停时有界性。可选抽样定理要求停时 τ 有界(或满足矩条件);无限停时下 E(M_τ)=E(M₀) 可能失效(赌徒输光原理)。

练习

【基础】泊松过程 λ=3/分钟,1 分钟内恰好 3 个事件的概率?

查看解答P(N(1)=3)=e^{−3}·3³/3! = e^{−3}·27/6 ≈ 0.224。
【自评反馈】对 → 下一题。

【进阶】W(0.5) 的分布?W(0.5)−W(0.2) 的方差?

查看解答W(0.5)~N(0,0.5),标准差 √0.5≈0.707;增量 W(0.5)−W(0.2)~N(0,0.3),方差 0.3。
【自评反馈】对 → 下一题。

【挑战】说明为什么 E(W(t)|W(s))=W(s)(t>s),即 W 是鞅。

查看解答W(t)=W(s)+(W(t)−W(s)),取条件期望 E(W(t)|W(s))=W(s)+E(W(t)−W(s)|W(s))=W(s)+0=W(s),因为增量独立且均值 0。
记
小结卡

① 马氏链:P 转移、Pⁿ n 步、π=πP 平稳;不可约非周期则遍历收敛。

② 泊松过程 N(t)~Poisson(λt),独立无记忆;布朗运动 W(t)~N(0,t),连续随机游走。

③ 鞅 E(M_{n+1}|Fₙ)=Mₙ,可选抽样 E(M_τ)=E(M₀);是公平游戏与随机优化的理论语言。