Skip to content

文档日期

2026-09-02 · 让稀有事件变成典型事件 ​

0. 导读与元结构 (Orientation & Meta-structure) ​

起点问题。 弱大数律告诉我们:若 X1,X2,… 独立同分布且 μ=EX1 存在,则经验均值

X―n:=1n∑i=1nXi

偏离 μ 的概率趋于零。但它不回答衰减多快。指数型 Markov 不等式容易给出

P(X―n≥x)≤exp⁡{−nI(x)},

真正困难的是:这个指数是否只是一个方便上界,还是概率的真实对数尺度?今日的答案是 Cramér theorem(Cramér 定理):在适当指数矩条件下,它确实是真实速率;匹配下界的发动机不是再做一次 Markov 估计,而是 change of measure(变化测度)。

对象与量词分四层:

  1. 固定一个实值概率分布 ν,令 Xi∼iidν;
  2. 对每个样本规模 n,考察 X―n 的分布;
  3. 对每个闭集 F⊂R 与开集 G⊂R,分别比较概率的上、下指数;
  4. 令 n→∞,以 n−1log⁡P(X―n∈⋅) 读取速度 n 的衰减率。

完整量词链是

∀F⊂R closed,lim supn→∞1nlog⁡P(X―n∈F)≤−infx∈FI(x),∀G⊂R open,lim infn→∞1nlog⁡P(X―n∈G)≥−infx∈GI(x).

先修:Markov 不等式、独立随机变量、弱大数律、凸函数与导数。主导 proof engines:指数 Markov 给局部上界;有限覆盖与 exponential tightness(指数紧性)把它升级为闭集上界;exponential tilting(指数倾斜)使目标均值在新测度下典型;弱大数律再给匹配下界。去向:Gärtner–Ellis theorem、Sanov theorem、relative entropy(相对熵)、统计物理中的自由能对偶,以及 rare-event importance sampling(稀有事件重要抽样)。

本期是三日 anchor micro-arc 的 Day 1。它诚实深化 2026-08-09:当日已把 Chernoff–Legendre 优化用于随机图割,只需固定候选的指数上界;今日删去候选计数与 union bound,补齐 Cramér 原理的集合表述和变化测度下界。Day 2 将把“调均值所支付的代价”改写成相对熵变分问题。

1. 核心概念、定理与谱系 (Concept, Theorem & Lineage) ​

1.1 对数矩母函数与凸对偶 ​

设 X 为实值随机变量。定义 log moment-generating function(对数矩母函数)

Λ(θ):=log⁡EeθX∈(−∞,+∞],θ∈R.

今日采用一个清楚而略强的版本:

(C)∃η>0,Eeη|X|<∞.

它保证 Λ(θ)<∞ 对所有 |θ|≤η 成立,特别地 μ=EX 有限。定义 Legendre–Fenchel transform(Legendre–Fenchel 变换)

I(x):=Λ∗(x)=supθ∈R{θx−Λ(θ)}.

逐项检查基本性质:

  • I 是仿射函数族 x↦θx−Λ(θ) 的逐点上确界,故凸且下半连续;
  • 取 θ=0 得 I(x)≥0;
  • Jensen 不等式给出 Λ(θ)≥θμ,故 I(μ)=0;
  • 若 x>μ,则负的 θ 不会改善上确界,所以
    I(x)=supθ≥0{θx−Λ(θ)};
    x<μ 时相应只需 θ≤0。

当 Λ 在某个 θx 可微且

Λ′(θx)=x,

凸对偶的切线条件给出

I(x)=θxx−Λ(θx).

这不是纯形式优化:θx 将定义一个均值恰为 x 的新概率分布。

1.2 Cramér 定理 ​

定理(Cramér large-deviation principle,一维版)。 设 X1,X2,… 为实值 i.i.d. 随机变量,并满足条件 (C)。令 X―n=n−1∑i=1nXi,令 I=Λ∗。则 {X―n} 以速度 n 满足大偏差原理:

  • 对每个闭集 F⊂R,
    lim supn→∞1nlog⁡P(X―n∈F)≤−infx∈FI(x);
  • 对每个开集 G⊂R,
    lim infn→∞1nlog⁡P(X―n∈G)≥−infx∈GI(x).

约定 log⁡0=−∞、inf∅=+∞。在一维 Cramér 定理中 I 是 good rate function(好率函数):每个次水平集 {x:I(x)≤c} 都紧。

若 x>μ 且 I 在 x 附近有限连续,则取 F=[x,∞) 与 G=(x,∞) 可得

limn→∞1nlog⁡P(X―n≥x)=−I(x).

这才是“P(X―n≥x)≈e−nI(x)”的精确定义:等式只在 n−1log 尺度上成立,次指数因子未被描述。

1.3 上界:Markov 之后还缺什么 ​

固定 θ≥0 与 x∈R。由独立性,

P(X―n≥x)=P(eθnX―n≥eθnx)≤e−θnxEeθ∑i=1nXi=exp⁡{−n(θx−Λ(θ))}.

对 θ≥0 取上确界便得到上尾 Chernoff bound。下尾类似。这一步只控制半直线;闭集上界还要处理两个问题。

  1. 局部化。 对每个 x∈F 与每个 a<I(x),可选某个 θx 使
    θxx−Λ(θx)>a.
    由连续性,该不等式在 x 的某个邻域仍成立,给出该邻域的指数上界。
  2. 从局部到整体。 有界部分 F∩[−R,R] 是紧集,可由有限多个上述邻域覆盖,有限和不改变主指数。尾部则由条件 (C) 和 Markov 不等式控制:
    P(|X―n|>R)≤exp⁡{−n(ηR−Λ(η))}+exp⁡{−n(ηR−Λ(−η))}.
    令 R→∞,得到 exponential tightness,封住无界尾部。

因此“优化 Markov”给的是局部代价;闭集上界还依赖紧致覆盖与指数紧性。若忽略后两步,便尚未证明完整 LDP。

1.4 下界:把目标均值变成新世界里的典型值 ​

先处理 regular point(正则点):固定 x,假设存在 θx 使 Λ 在该点可微且 Λ′(θx)=x。定义 tilted law(倾斜分布)νθx:

dνθxdν(y)=exp⁡{θxy−Λ(θx)}.

归一化成立,因为

∫eθxy−Λ(θx)ν(dy)=1.

在新分布下,随机变量的均值是

EθxX=E[XeθxX]EeθxX=Λ′(θx)=x.

令 Pθx 表示 n 个倾斜样本的乘积分布。对任意 δ>0,在事件

An,δ:={|X―n−x|<δ}

上,似然比满足

dPdPθx=exp⁡{−nθxX―n+nΛ(θx)}.

若 θx≥0,则 X―n<x+δ 给出

P(An,δ)=Eθx[dPdPθx1An,δ]≥exp⁡{−n[θx(x+δ)−Λ(θx)]}Pθx(An,δ)=exp⁡{−n[I(x)+θxδ]}Pθx(An,δ).

θx<0 时把误差写成 |θx|δ 即可。关键发生在最后一项:在倾斜分布下,x 已是均值,因此弱大数律给出

Pθx(An,δ)⟶1.

所以

lim infn→∞1nlog⁡P(An,δ)≥−I(x)−|θx|δ.

令 δ↓0,得到 −I(x) 的匹配下界。一般开集 G 中选取一个 x∈G,再取足够小的区间 (x−δ,x+δ)⊂G;最后对 x 取优化。非正则边界点由凸性、逼近与下半连续性补齐。

证明的核心可压成

原测度下稀有→exponential tilt新测度下典型→likelihood costI(x).

上界的 I(x) 来自优化;下界说明同一个优化参数真的制造出该偏差,因而率函数不是估计伪影。

1.5 例:Gaussian 分布的二次代价 ​

若 X∼N(μ,σ2) 且 σ>0,则

Λ(θ)=μθ+σ2θ22.

因此

I(x)=supθ∈R{θ(x−μ)−σ2θ22}.

唯一极大点是

θx=x−μσ2,

故

I(x)=(x−μ)22σ2.

倾斜后的分布仍为 Gaussian,只是均值从 μ 平移到

μ+σ2θx=x,

方差保持 σ2。这给出一个完全可见的模型:指数倾斜不是抽象技巧,而是在最小修改分布形状的同时把目标偏差变成中心。

1.6 谱系:从正态尺度之外到现代 LDP ​

Harald Cramér 1938 年的法文论文 Sur un nouveau théorème-limite de la théorie des probabilités 研究独立同分布随机变量之和在比中央极限定理的 n 尺度更远处如何衰减。论文明确要求生成函数在零点邻域存在,并使用后来称为 Esscher / Cramér transform 的变化测度。其结果是现代 Cramér 定理的历史起点;现代闭集—开集 LDP 表述则把“尾概率的一个精确渐近”升级为适用于一般集合、可通过连续映射传递的框架。

历史谱系不能压成“旧定理已经过时”:原始问题、指数倾斜和主指数都存续下来;后来变化的是假设的整理、集合语言、凸分析表述及从 i.i.d. 和推广到经验测度、Markov 过程与路径空间的范围。

Graduate / modern extension. Cramér 定理是 contraction principle(伸缩原理)的基本输入:若 X―n 以率函数 I 满足 LDP,且 f:R→R 连续,则 f(X―n) 的率函数为

J(y)=inf{I(x):f(x)=y}.

这说明大偏差率函数像一种“指数尺度上的能量”,在投影或粗粒化时沿纤维取最小值。Day 2 将把同一结构提升到经验分布,并由 relative entropy 解释这种最小代价。

2. 固定模板习题 (Template Exercise) ​

题目 ​

设随机变量 X 满足

a≤X≤ba.s.,EX=0,

其中 a<b。证明 Hoeffding lemma(Hoeffding 引理)

∀λ∈R,log⁡EeλX≤λ2(b−a)28.

进而设 X1,…,Xn 相互独立,ai≤Xi≤bi a.s.,证明对每个 t>0,

P(∑i=1n(Xi−EXi)≥t)≤exp⁡(−2t2∑i=1n(bi−ai)2).

模板识别 ​

看到“独立、有界、控制和的尾概率”,应识别 bounded-mgf → exponential Markov → optimize(有界矩母函数—指数 Markov—参数优化) 模板。真正需要证明的不是尾界本身,而是每个中心化变量的 log-mgf 有统一二次包络;独立性随后把这些包络相加。

分级提示 ​

  1. 由 x↦eλx 的凸性,对每个 x∈[a,b] 用端点割线控制 eλx,把一般 X 约化为一个只取 a,b 的 Bernoulli 型变量。
  2. 令 p=−a/(b−a),并把得到的上界写成某个函数 g(u),其中 u=λ(b−a);核验 g(0)=g′(0)=0。
  3. 证明 g″(u)≤1/4;这是倾斜 Bernoulli 分布的方差至多 1/4。两次积分得到 g(u)≤u2/8。
  4. 对中心化的 Xi 应用引理,使用独立性相乘 mgf,再对 λ>0 优化一个形如 −λt+λ2L/8 的二次式,其中 L=∑i(bi−ai)2。

3. 大家之言 (Mathematician in Their Own Words) ​

人物与语境。 Terence Tao 在其第一人称课程博客 254A, Notes 0: A review of probability theory(2010-01-01)回顾 Markov、Chebyshev、矩方法与概率尾界。下列短引已在原文页面逐字核对。

核验短引(原文):

“Chebyshev’s inequality gives a simple but still useful large deviation inequality”

中文翻译:

“Chebyshev 不等式给出一个简单但仍然有用的大偏差不等式。”

评论。 Tao 的原文是在搭建概率工具箱,不是陈述 Cramér 定理。短引有价值之处恰在于它把 Chebyshev 放在正确层级:若只知方差,二阶矩法已经能证明偏差概率变小;但对独立有界和,它通常只给多项式上界,无法读取真实的指数衰减率。今日的 log-mgf 与指数倾斜并不否定 Chebyshev,而是提高问题的分辨率:从“偏差罕见”推进到“每单位样本支付多少对数代价”。

这也给出一种方法论上的克制:先问目标需要方差尺度、中心极限尺度,还是远尾指数尺度,再选择工具。Chebyshev、中央极限定理与大偏差原理并非互相替代;它们分别回答粗尾界、中心波动与远尾指数三个不同问题。短引只说明这一工具层级,不能替代 Section 1 的变化测度证明。

4. 反例与边界 (Counterexample & Boundary) ​

4.1 有限均值不保证速度 n 的正大偏差指数 ​

令 X 服从 Pareto 分布

P(X>t)=t−3,t≥1.

逐项核验:

  1. X≥1,且
    EX=∫0∞P(X>t)dt=1+∫1∞t−3dt=32<∞.
    因而弱大数律仍适用。
  2. 对每个 θ>0,
    EeθX=∞,
    因为指数增长压过密度的多项式衰减。正向 Cramér condition 失效。
  3. 固定 x>μ=3/2。若 X1>nx,则其余变量非负,所以 X―n>x。因此
    P(X―n>x)≥P(X1>nx)=(nx)−3.
  4. 另一方面该概率至多 1,故
    −3log⁡(nx)n≤1nlog⁡P(X―n>x)≤0.
    由夹逼,
    limn→∞1nlog⁡P(X―n>x)=0.

因此偏差概率确实趋于零,但不是 e−cn 级;在速度 n 上,正率函数退化为零。主导机制是 single big jump(单个大跳),而不是所有样本共同受到指数倾斜。

4.2 这暴露的不是大数律边界,而是 Cramér 机制边界 ​

该例同时满足:

  • X―n→μ in probability;
  • P(X―n>x)→0;
  • 但没有正的速度 n 指数。

所以“收敛”与“指数收敛”不可混同。指数矩假设的作用正是排除过重的尾部,使整体偏移比单个巨值更经济;没有它时,Legendre 优化所需的正参数根本不可用。

4.3 为什么定理写闭集上界、开集下界 ​

即使指数矩良好,也不能对每个任意 Borel 集 A 都承诺极限

limn→∞1nlog⁡P(X―n∈A)

存在。格点分布和随 n 可达性变化的不规则集合会让概率在零与非零之间振荡。闭集上界允许从外部封闭,开集下界允许在内部放入小邻域;二者正是稳定于极限的拓扑形式。只有当

infA∘I=infA―I

时,上下界才自动合成一个集合级指数极限。

5. 今日一问与数学地图 (Question & Map) ​

今日一问。 设 X1,X2,… 为 i.i.d. 标准 Gaussian 随机变量,已知

X―n的率函数为I(x)=x22.

令

Yn=(X―n)2.

不用直接计算 Yn 的密度,尝试仅从 Cramér 定理推出 Yn 的大偏差原理及其率函数 J:R→[0,∞]。你需要回答:

  1. 对给定 y>0,连续映射 f(x)=x2 有两个原像,率函数为何应在它们之间取最小?
  2. y=0 处应怎样处理?
  3. y<0 为什么必须令 J(y)=+∞?
  4. 这个例子怎样展示 contraction principle 中“沿纤维取最小代价”的含义?

方向提示。 先写

J(y)=inf{x22:x2=y},

但先不要急着化简;逐项处理空原像、单原像与双原像。

数学地图。

Markov 不等式、独立和、弱大数律↓log-mgf Λ 与凸对偶 I=Λ∗↙↘指数上界与指数紧性指数倾斜与匹配下界↓↓Cramér 大偏差原理↙↘contraction principle、Gärtner--EllisSanov、相对熵、重要抽样↓统计物理自由能、随机图变分问题、稀有事件算法

附录 A:习题完整解答 (Full Solution) ​

展开完整解答

模板:bounded-mgf → exponential Markov → optimize(有界矩母函数—指数 Markov—参数优化)。

先证 Hoeffding 引理。由 EX=0 且 a≤X≤b,必有 a≤0≤b。若 a=b=0 结论平凡;以下设 a<b。

对每个 x∈[a,b],凸函数 u↦eλu 位于连接端点值的割线之下:

eλx≤b−xb−aeλa+x−ab−aeλb.

取期望并用 EX=0,得

EeλX≤bb−aeλa+−ab−aeλb.

令

p:=−ab−a∈[0,1],ν:=λ(b−a).

则 a=−p(b−a)、b=(1−p)(b−a),上式右侧等于

(1−p)e−pν+pe(1−p)ν=e−pν(1−p+peν).

定义

g(ν):=−pν+log⁡(1−p+peν).

直接计算 g(0)=g′(0)=0。再令

qν:=peν1−p+peν∈[0,1].

则

g″(ν)=qν(1−qν)≤14.

对 ν≥0,两次积分得

g(ν)=∫0ν(ν−s)g″(s)ds≤14∫0ν(ν−s)ds=ν28.

对 ν<0,同一 Taylor 积分余项或将积分方向反转,仍得 g(ν)≤ν2/8。所以

log⁡EeλX≤g(λ(b−a))≤λ2(b−a)28.

Hoeffding 引理得证。

现在令

Yi:=Xi−EXi.

平移不改变区间长度,所以 EYi=0,且 Yi 落在长度 bi−ai 的某个区间内。对任意 λ>0,由指数 Markov 不等式、独立性与刚证的引理,

P(∑i=1nYi≥t)≤e−λtEexp⁡(λ∑i=1nYi)=e−λt∏i=1nEeλYi≤exp⁡(−λt+λ28∑i=1n(bi−ai)2).

记

L:=∑i=1n(bi−ai)2.

若 L=0,每个 Xi 都几乎处处为常数,左侧概率为 0。若 L>0,对二次式

−λt+Lλ28

求最小值,唯一最优参数为

λ∗=4tL>0.

代入得

−λ∗t+Lλ∗28=−2t2L.

因此

P(∑i=1n(Xi−EXi)≥t)≤exp⁡(−2t2∑i=1n(bi−ai)2).

这套模板的可复用点是:先把单变量 log-mgf 控制成二次函数;独立性把二次预算相加;最后只做一维凸优化。

资料与来源 (Sources) ​

  1. David Gamarnik, “Lecture 2: Large Deviations for i.i.d. Random Variables”, MIT 6.265/15.070J Advanced Stochastic Processes, Fall 2013, dated 2013-09-09. 包含指数 Markov、log-mgf 与 Legendre 变换的推导。
    https://ocw.mit.edu/courses/15-070j-advanced-stochastic-processes-fall-2013/ed35842f2e902421e4a6cdd4d9eca6fe_MIT15_070JF13_Lec2.pdf

  2. David Gamarnik, “Lecture 3: Cramér’s Theorem”, MIT 6.265/15.070J Advanced Stochastic Processes, Fall 2013, dated 2013-09-11. 包含闭集上界、开集下界、率函数性质与指数变化测度证明。
    https://ocw.mit.edu/courses/15-070j-advanced-stochastic-processes-fall-2013/f536bc29826255afee543b03b605bf4d_MIT15_070JF13_Lec3.pdf

  3. Harald Cramér, “Sur un nouveau théorème-limite de la théorie des probabilités”, Colloque consacré à la théorie des probabilités, Actualités scientifiques et industrielles 736, Hermann, Paris, 1938, pp. 2–23; Hugo Touchette English translation, 2018. 这是历史原始论文的可核验法英对照版本,记录指数矩条件、变化分布与大偏差尺度。
    https://appliedmaths.sun.ac.za/~htouchette/archive/docs/cramer1.pdf

  4. Terence Tao, “254A, Notes 0: A review of probability theory”, What’s new, 2010-01-01. 本期 Section 3 的第一人称课程笔记与核验短引来源。
    https://terrytao.wordpress.com/2010/01/01/254a-notes-0-a-review-of-probability-theory/