第 15 章 · 遍历理论与大偏差

15.2 大偏差原理与后续研究

Large Deviations and Further Study
中心极限定理描述均值附近的波动,大偏差理论研究固定距离的稀有偏离。本节准确陈述开集下界、闭集上界与 Cramér 定理,从 Chernoff 方法推导上界,再用指数倾斜说明下界机制。伯努利例子完整计算相对熵速率,并比较有限样本概率与指数界;最后明确本书覆盖的研究生核心与仍需专门学习的方向。

15.2.1 大数、中心极限与大偏差的尺度

设 \(X_i\) iid,样本均值为 \(\overline X_n\)、总体均值为 \(\mu\)。大数定律说明 \(\overline X_n\to\mu\);有限正方差的中心极限定理刻画 \(\sqrt n(\overline X_n-\mu)\) 的分布。大偏差理论则关心固定距离的偏离,例如 \(P(\overline X_n\ge a)\),其中 \(a>\mu\) 不随 \(n\) 靠近均值。

CLT 在固定标准化阈值处给出极限;上述事件的标准化阈值随 \(\sqrt n\) 增大,不能仅由 CLT 推出准确的远尾指数。大偏差原理研究 \(n^{-1}\ln P(\overline X_n\in A)\) 的极限范围。本节采用自然对数,速率的单位与第 14.1 节使用的比特不同。

大偏差原理(large deviation principle,LDP)

称实随机变量序列 \(Z_n\) 满足速度 \(n\)、速率函数 \(I:\mathbb R\to[0,\infty]\) 的 LDP,如果 \(I\) 下半连续,且对每个闭集 \(F\) 与开集 \(G\),

\[\limsup_{n\to\infty}\frac1n\ln P(Z_n\in F)\le-\inf_{x\in F}I(x),\qquad\liminf_{n\to\infty}\frac1n\ln P(Z_n\in G)\ge-\inf_{x\in G}I(x).\]
(15.2-1)

约定 \(\ln0=-\infty\)、\(\inf\varnothing=+\infty\)。若每个有限 \(c\) 对应的集合 \(\{I\le c\}\) 都紧,则称为好速率函数(good rate function)。一般速度也可以是趋于无穷的 \(a_n\),本节集中讨论独立样本均值的速度 \(n\)。

对一般 Borel 集 \(A\),用其内部与闭包可得上下界。只有二者上的速率下确界相同,才可据此断言对数概率极限等于该共同值。即使 \(P(Z_n\in A)=\exp[-nI_*+o(n)]\),也只控制指数尺度;概率比值未必趋于 \(1\)。

15.2.2 Cramér 定理与速率函数

Cramér 定理 · 实值 iid 均值

设存在 \(\delta>0\),使 \(E e^{\theta X_1}<\infty\) 对 \(|\theta|<\delta\) 成立。定义对数矩母函数及其 Legendre–Fenchel 变换:

\[\Lambda(\theta)=\ln Ee^{\theta X_1},\qquad I(x)=\sup_{\theta\in\mathbb R}\{\theta x-\Lambda(\theta)\}.\]
(15.2-2)

矩母函数发散处约定 \(\Lambda=+\infty\),相应参数不贡献上确界。则 \(\overline X_n\) 满足速度 \(n\)、好速率函数 \(I\) 的 LDP。零点邻域的有限性已经保证全部有限阶矩存在,但不要求矩母函数在整条实轴有限。

\(I\) 是连续仿射函数族的上确界,因而凸且下半连续。取 \(\theta=0\) 得 \(I\ge0\);Jensen 不等式给出 \(\Lambda(\theta)\ge\theta\mu\),故 \(I(\mu)=0\)。在假设的邻域内任选 \(d>0\),有

\[I(x)\ge d|x|-\max\{\Lambda(d),\Lambda(-d)\}.\]
(15.2-3)

因此其有限水平集有界,再结合下半连续得到紧性。凸性还使 \(I\) 在均值右侧不减、左侧不增。下面完整推导闭集上界与可由指数倾斜实现的点处的下界;对矩母函数有效域边界等一般情形的下界延拓,作为 Cramér 定理的技术部分引用文献 [1–2],不把内部点的证明当作所有边界的证明。

15.2.3 从 Chernoff 界到闭集上界

尾部与局部概率的上界

对 \(\theta\ge0\),独立性和指数 Markov 不等式给出

\[P(\overline X_n\ge a)\le\exp[-n(\theta a-\Lambda(\theta))].\]
(15.2-4)

对所有使 \(\Lambda\) 有限的非负参数取下确界。若 \(a\ge\mu\),非正参数的 \(\theta a-\Lambda(\theta)\le\theta(a-\mu)\le0\),故所得上确界就是 \(I(a)\)。左尾用非正参数处理。更一般地,任取 \(x\)、半径 \(r>0\) 和可用参数 \(\theta\),当 \(|\overline X_n-x|<r\) 时有 \(\theta\overline X_n\ge\theta x-|\theta|r\),所以

\[P(|\overline X_n-x|<r)\le\exp[-n(\theta x-\Lambda(\theta)-|\theta|r)].\]
(15.2-5)
紧覆盖与指数紧性

令 \(K\) 为紧集,取任意实数 \(L<\inf_KI\)。对每个 \(x\in K\),上确界定义保证能选一个参数 \(\theta_x\),使 \(\theta_xx-\Lambda(\theta_x)>L\),再选小邻域保持指数大于 \(L\)。紧性给出有限个邻域覆盖,因此式 (15.2-5) 和并集界给出 \(\limsup n^{-1}\ln P(\overline X_n\in K)\le-L\)。

对一般闭集 \(F\),还要控制无界部分。固定式 (15.2-3) 的 \(d\),两侧 Chernoff 界给出

\[P(|\overline X_n|>R)\le2\exp\{-n[dR-\max(\Lambda(d),\Lambda(-d))]\}.\]
(15.2-6)

令 \(R\) 足够大,可以让右侧指数衰减率超过任意指定的 \(L\),这称为指数紧性。将 \(F\) 分成紧集 \(F\cap[-R,R]\) 与外部部分,两个上界相加仍有对数上极限不超过 \(-L\)。最后令 \(L\uparrow\inf_FI\),就证明了式 (15.2-1) 的闭集上界;若下确界为无穷,可让 \(L\) 任意增大。空集情形直接成立。

15.2.4 指数倾斜给出局部下界

对位于 \(\{\Lambda<\infty\}\) 内部的参数 \(\theta\),以 \(P\) 表示单个 \(X\) 的分布,定义倾斜分布(exponential tilting)

\[\frac{dP_\theta}{dP}(x)=e^{\theta x-\Lambda(\theta)}.\]
(15.2-7)

其积分为 \(1\)。参数附近的指数可积上界允许对矩母函数求导,得到 \(E_\theta X=\Lambda'(\theta)\)、\(\operatorname{Var}_\theta(X)=\Lambda''(\theta)\)。这给出了选择新分布的方法:让原本稀有的均值 \(x\) 成为新分布的平均值。

下界推导 · 固定一个可实现的均值

令 \(x=\Lambda'(\theta)\),在 \(P_\theta\) 下独立生成 \(n\) 项。由大数定律,对 \(A_{n,r}=\{|\overline X_n-x|<r\}\) 有 \(P_\theta^{\otimes n}(A_{n,r})\to1\)。把测度改回原分布,得到

\[\begin{aligned}P^{\otimes n}(A_{n,r})&=E_\theta\left[e^{-n\theta\overline X_n+n\Lambda(\theta)}\mathbf1_{A_{n,r}}\right]\\&\ge e^{-n[\theta x-\Lambda(\theta)+|\theta|r]}P_\theta^{\otimes n}(A_{n,r}).\end{aligned}\]
(15.2-8)

凸性使 \(\theta\) 在式 (15.2-2) 中实现 \(I(x)=\theta x-\Lambda(\theta)\)。取对数、除以 \(n\),第二个因子的贡献趋零,因此下极限至少为 \(-I(x)-|\theta|r\)。对包含 \(x\) 的开集,选任意足够小的邻域并令 \(r\downarrow0\),即得下界 \(-I(x)\)。对这些 \(x\) 再优化便得到相应的开集下界。有效域边界或不能表示为 \(\Lambda'(\theta)\) 的点,需要前述额外延拓,不能省略。

同一个变换也解释重要性采样为何能用于稀有事件:在新分布下让事件常见,再乘回似然比。大偏差推导关注指数尺度;模拟算法还必须检查二阶矩与计算成本。

15.2.5 伯努利均值:相对熵作为速率

例 1 · 完整求出速率函数

设 \(X_i\sim\operatorname{Bernoulli}(p)\) iid,\(0<p<1\)。有 \(\Lambda(\theta)=\ln(1-p+pe^\theta)\)。对 \(0<x<1\),解 \(x=\Lambda'(\theta)=pe^\theta/(1-p+pe^\theta)\),得到

\[\theta_x=\ln\frac{x(1-p)}{p(1-x)},\qquad I(x)=x\ln\frac{x}{p}+(1-x)\ln\frac{1-x}{1-p}.\]
(15.2-9)

在端点按连续延拓定义 \(I(0)=-\ln(1-p)\)、\(I(1)=-\ln p\);在 \([0,1]\) 外为无穷。端点事件的概率分别为 \((1-p)^n,p^n\),与该速率精确一致。所有内部点都可由有限参数倾斜实现,所以本例没有遗漏内部下界。倾斜后的单项恰为 \(\operatorname{Bernoulli}(x)\)。

这个速率是使用自然对数的 \(D(\operatorname{Bernoulli}(x)\Vert\operatorname{Bernoulli}(p))\),等于第14.1节比特单位相对熵乘以 \(\ln2\)。求导可得 \(I'(p)=0\)、\(I''(p)=1/[p(1-p)]\),所以在均值附近 \(I(p+t)=t^2/[2p(1-p)]+o(t^2)\)。这解释了局部二次形状与 CLT 方差的联系,但单靠这个 Taylor 展开不能证明所有随 \(n\) 变化的中等偏差结论。

伯努利参数零点三的速率与局部二次式 I(x)0100.31x实线:相对熵速率;虚线:均值附近的二次式最小值在 x = p = 0.3;远离均值时须使用完整速率函数
图 15.2-1 两条曲线按精确公式计算,横轴为样本比例。虚线 \((x-0.3)^2/(2\cdot0.3\cdot0.7)\) 只表达局部二阶展开。
例 2 · 固定远尾的数值比较

取 \(p=0.3\)、\(a=0.6\)。有 \(I(0.6)\approx0.192042\)、倾斜参数 \(\theta_a=\ln3.5\)。Chernoff 上界给出 \(P(\overline X_{100}\ge0.6)\le e^{-100I(0.6)}\approx4.5680\times10^{-9}\)。直接二项求和 \(\sum_{k=60}^{100}\binom{100}{k}0.3^k0.7^{100-k}\approx5.12995\times10^{-10}\)。上界有效,却不是该有限样本概率的精确值。

由于 \(I\) 在 \([p,1]\) 递增且在 \(a\in(p,1)\) 连续,闭尾集的上界与开尾集的下界在 \(a\) 处吻合,因此 \(\lim_{n\to\infty}n^{-1}\ln P(\overline X_n\ge a)=-I(a)\)。这一结论允许有限样本中仍有显著的前置因子差异。

指数倾斜把成功概率零点三变为零点六原分布 Bernoulli(0.3)倾斜分布 Bernoulli(0.6)0.70.30.40.60101θ = ln 3.5,所有柱高按概率绘制
图 15.2-2 倾斜后,样本比例趋近 \(0.6\) 成为大数定律描述的典型事件;将概率换回原分布时,指数似然比给出其稀有程度。

15.2.6 习题与后续学习范围

习题 15.2-1 · 正态均值的速率函数

若 \(X_i\sim N(\mu,\sigma^2)\) iid,\(\sigma>0\),求 \(I(x)\) 和右尾固定阈值 \(a>\mu\) 的对数概率极限。

解答。\(\Lambda(\theta)=\mu\theta+\sigma^2\theta^2/2\)。对 \(\theta(x-\mu)-\sigma^2\theta^2/2\) 配方或求导,最大值在 \(\theta=(x-\mu)/\sigma^2\),故 \(I(x)=(x-\mu)^2/(2\sigma^2)\)。速率在 \(\mu\) 右侧递增且连续,因此开闭尾集上下界吻合,极限为 \(-(a-\mu)^2/(2\sigma^2)\)。正态模型中的二次速率对所有 \(x\) 成立;一般分布只在均值附近具有二次展开。

习题 15.2-2 · 矩母函数不必处处有限

若 \(X_i\) iid 且为速率 \(1\) 的指数变量,求全部实数上的 \(I(x)\),并求 \(a=2\) 的右尾速率。

解答。\(\Lambda(\theta)=-\ln(1-\theta)\) 当 \(\theta<1\),其他处为无穷。对 \(x>0\),优化 \(\theta x+\ln(1-\theta)\) 得 \(\theta=1-1/x\),故 \(I(x)=x-1-\ln x\)。若 \(x=0\),令 \(\theta\to-\infty\),对数项趋于无穷;\(x<0\) 时线性项也趋于正无穷。因此 \(I(x)=\infty\) 于 \(x\le0\)。\(I(2)=1-\ln2\approx0.306853\),相应对数右尾极限为其相反数。矩母函数在零点邻域有限,已足以满足定理假设。

习题 15.2-3 · 为什么必须区分开集与闭集

对伯努利均值,取一个位于 \((0,1)\) 的无理数 \(a\) 和集合 \(A=\{a\}\)。比较真实概率与 \(e^{-nI(a)}\),说明是否违反 LDP。

解答。每个样本均值都是有理数 \(k/n\),所以 \(P(\overline X_n\in A)=0\),对数概率为负无穷;而 \(I(a)\) 有限。这不违反 LDP:闭集上界只要求负无穷不超过 \(-I(a)\),成立;\(A\) 的内部为空,其下界为 \(-\infty\),没有给出有限下界。对任意集合都写一个相等的对数极限会遗漏这种格点限制。

习题 15.2-4 · 删除独立性后的失败

令 \(Z\sim\operatorname{Bernoulli}(0.3)\),每个 \(X_i=Z\)。尽管边缘矩母函数与例 1 相同,计算 \(n^{-1}\ln P(\overline X_n\ge0.6)\) 的极限。

解答。样本均值恒等于 \(Z\),尾概率恒为 \(0.3\),所以极限为 \(\lim_n(\ln0.3)/n=0\),并非 \(-0.192042\)。指数 Markov 不等式本身仍然成立,但联合矩母函数不再是单项矩母函数的 \(n\) 次方。因而 iid 条件不能只用边缘同分布代替。

本书的覆盖范围与继续学习

前面的本科基础、测度积分、条件期望、鞅、马尔可夫过程、布朗运动和随机分析入门,构成研究生概率核心的学习路线;它们并不穷尽概率论的所有研究方向。本节也没有展开一般大偏差理论中的全部边界处理、函数空间紧性或无限维原理。

表 15.2-1 后续专题及所需基础
后续方向需要继续学习的内容本书可用的基础
随机场与空间模型Gibbs 测度、渗流、空间相关与相变条件分布、Markov 性、遍历性
Lévy 过程与跳跃分析无限可分分布、Lévy–Khintchine 公式、跳跃补偿特征函数、泊松过程、鞅
更深入的随机分析一般半鞅积分、测度变换、局部时布朗运动、伊藤积分与随机微分方程
随机偏微分方程无穷维噪声、函数空间估计与方程解的正则性测度、Hilbert 空间与随机分析;另需偏微分方程
函数空间大偏差Sanov、Schilder、收缩原理与变分方法弱收敛、紧性、速率函数与指数倾斜

关键术语

  • 大偏差原理:用开集下界与闭集上界刻画指数尺度概率。
  • 速度:对数概率的归一化尺度,本节为 n。
  • 速率函数:量化各位置的指数衰减代价。
  • 好速率函数:有限水平集都紧的速率函数。
  • Legendre–Fenchel 变换:对线性函数减对数矩母函数取上确界。
  • 指数紧性:无界区域的概率可获得任意强度的指数控制。
  • 指数倾斜:按指数权重改变概率测度。
  • 对数等价:控制概率的主要指数,不保证普通比值趋于一。

参考文献

  1. Durrett, R. (2019). Probability: Theory and Examples, 5th ed., §2.7。作者公开教材:独立和的大偏差及指数变换。
  2. MIT 6.265 / 15.070J (2013). Lecture 3, Large Deviations Theory. Cramér’s Theorem。官方讲义:开闭集表述、速率性质和换测度下界;讲义展开的证明采用额外的有效域与支撑假设。
  3. University of Oxford, C8.6 Limit Theorems and Large Deviations in Probability (2021–22)。官方课程纲要:函数空间弱收敛、指数紧性及 Sanov、Schilder、收缩原理等后续专题。