茆书复习笔记 04|特征函数、大数定律与中心极限定理
导读
大数定律和中心极限定理,最容易混淆的不是结论,而是条件:哪个要求独立同分布?哪个只要两两不相关?哪个连方差都不需要?
这一篇先整理特征函数这个“工具”,再把几个极限定理的条件放进表格里对比,最后给出用特征函数证明辛钦大数定律和林德伯格-莱维中心极限定理的主线。
1. 特征函数
定义(P216)
\[\varphi(t)=E(e^{itX})=\begin{cases}\sum_k e^{itx_k}P(X=x_k),&\text{离散}\\[4pt]\int_{-\infty}^{\infty}e^{itx}p(x)\,dx,&\text{连续}\end{cases}\]记住欧拉公式 \(e^{i\theta}=\cos\theta+i\sin\theta\)。因为 \(\vert e^{itX}\vert =1\),所以任何随机变量的特征函数都存在(矩母函数就不一定)。
性质
- \(\vert \varphi(t)\vert \le\varphi(0)=1\),\(\varphi(-t)=\overline{\varphi(t)}\)
- 若 \(Y=aX+b\),则 \(\varphi_Y(t)=e^{ibt}\varphi_X(at)\)
- 若 \(X\) 与 \(Y\) 独立,则 \(\varphi_{X+Y}(t)=\varphi_X(t)\varphi_Y(t)\)
- 若 \(E(X^k)\) 存在,则 \(\varphi^{(k)}(0)=i^kE(X^k)\)
-
逆转公式:若 \(x_1<x_2\) 是 \(F\) 的连续点,则
\[F(x_2)-F(x_1)=\lim_{T\to\infty}\frac1{2\pi}\int_{-T}^{T}\frac{e^{-itx_1}-e^{-itx_2}}{it}\varphi(t)\,dt\] -
唯一性定理:特征函数与分布函数一一对应。若 \(\int\vert \varphi(t)\vert \,dt<\infty\),则 \(X\) 为连续型,且
\[p(x)=\frac1{2\pi}\int_{-\infty}^{\infty}e^{-itx}\varphi(t)\,dt\]
狄利克雷积分(P223,证明逆转公式时要用)
\[D(a)=\frac1\pi\int_0^{\infty}\frac{\sin at}{t}\,dt=\begin{cases}\tfrac12,&a>0\\0,&a=0\\-\tfrac12,&a<0\end{cases}\]常用分布的特征函数(P220,指导书 P227)
| 分布 | \(\varphi(t)\) |
|---|---|
| 单点 \(P(X=a)=1\) | \(e^{ita}\) |
| 0-1 分布 \(b(1,p)\) | \(pe^{it}+q\) |
| 二项 \(b(n,p)\) | \((pe^{it}+q)^n\) |
| 泊松 \(P(\lambda)\) | \(e^{\lambda(e^{it}-1)}\) |
| 几何 \(Ge(p)\) | \(\dfrac{pe^{it}}{1-qe^{it}}\) |
| 负二项 \(Nb(r,p)\) | \(\left(\dfrac{pe^{it}}{1-qe^{it}}\right)^r\) |
| 正态 \(N(\mu,\sigma^2)\) | \(\exp\left\{i\mu t-\frac{\sigma^2t^2}2\right\}\) |
| 标准正态 \(N(0,1)\) | \(e^{-t^2/2}\) |
| 柯西 \(Cau(\mu,\lambda)\) | \(\exp\{i\mu t-\lambda\lvert t\rvert\}\) |
| 均匀 \(U(a,b)\) | \(\dfrac{e^{itb}-e^{ita}}{it(b-a)}\) |
| 均匀 \(U(-a,a)\) | \(\dfrac{\sin at}{at}\) |
| 指数 \(Exp(\lambda)\) | \(\left(1-\frac{it}{\lambda}\right)^{-1}\) |
| 伽玛 \(Ga(\alpha,\lambda)\) | \(\left(1-\frac{it}{\lambda}\right)^{-\alpha}\) |
| 卡方 \(\chi^2(n)\) | \((1-2it)^{-n/2}\) |
| 贝塔 \(Be(a,b)\) | \(\dfrac{\Gamma(a+b)}{\Gamma(a)}\sum_{j=0}^{\infty}\dfrac{\Gamma(a+j)\,(it)^j}{\Gamma(a+b+j)\,\Gamma(j+1)}\) |
推导示例:泊松与指数
泊松:\(\varphi(t)=\sum_{k=0}^{\infty}e^{itk}\frac{\lambda^k}{k!}e^{-\lambda}=e^{-\lambda}\sum_{k=0}^{\infty}\frac{(\lambda e^{it})^k}{k!}=e^{-\lambda}e^{\lambda e^{it}}\)。
指数:\(\varphi(t)=\int_0^{\infty}e^{itx}\lambda e^{-\lambda x}\,dx=\frac{\lambda}{\lambda-it}=\left(1-\frac{it}\lambda\right)^{-1}\)。
记忆技巧:伽玛是 \(\alpha\) 个“指数”相加,所以是指数的 \(\alpha\) 次方;卡方是 \(Ga(\frac n2,\frac12)\),代入即得。
2. 收敛性之间的关系
定理 4.1.1(P209):若 \(X_n\xrightarrow{P}a\),\(Y_n\xrightarrow{P}b\),则
- \[X_n\pm Y_n\xrightarrow{P}a\pm b\]
- \[X_nY_n\xrightarrow{P}ab\]
- \(X_n/Y_n\xrightarrow{P}a/b\)(\(b\ne0\))
更一般地,若 \(g\) 在 \(a\) 处连续,则 \(g(X_n)\xrightarrow{P}g(a)\)。
定理 4.1.2(P212):\(X_n\xrightarrow{P}X\ \Rightarrow\ X_n\xrightarrow{L}X\)(反之一般不成立)。
定理 4.1.3:若 \(c\) 为常数,则 \(X_n\xrightarrow{P}c\iff X_n\xrightarrow{L}c\)。
连续性定理:\(F_n\to F\)(弱收敛)\(\iff\varphi_n(t)\to\varphi(t)\) 对每个 \(t\) 成立。这是用特征函数证明极限定理的关键工具。
3. 大数定律
统一形式:
\[\lim_{n\to\infty}P\left(\left\vert \frac1n\sum_{i=1}^nX_i-\frac1n\sum_{i=1}^nE(X_i)\right\vert <\varepsilon\right)=1.\]| 名称 | 页码 | 条件 | 掌握要求 |
|---|---|---|---|
| 伯努利大数定律 | P230 | \(S_n\sim b(n,p)\),结论为 \(\frac{S_n}n\xrightarrow{P}p\) | 会证 |
| 切比雪夫大数定律 | P233 | \(\{X_i\}\) 两两不相关,方差存在且有共同上界 | |
| 马尔可夫大数定律 | P234 | \(\frac1{n^2}\mathrm{Var}\left(\sum_{i=1}^nX_i\right)\to0\)(马尔可夫条件) | |
| 辛钦大数定律 | P235 | \(\{X_i\}\) 独立同分布,\(E(X_i)\) 存在(不要求方差存在) | 会证 |
直观理解:伯努利大数定律说明“频率稳定于概率”,辛钦大数定律说明“样本均值稳定于总体均值”,这是矩估计和蒙特卡罗方法的依据。
证明:伯努利(用切比雪夫不等式)
\(E\left(\frac{S_n}{n}\right)=p\),\(\mathrm{Var}\left(\frac{S_n}{n}\right)=\frac{p(1-p)}{n}\),所以
\[P\left(\left\vert \frac{S_n}n-p\right\vert \ge\varepsilon\right)\le\frac{p(1-p)}{n\varepsilon^2}\le\frac{1}{4n\varepsilon^2}\to0.\]切比雪夫、马尔可夫大数定律的证明完全相同,只是方差的上界换一下。
证明:辛钦(用特征函数)
设 \(E(X_i)=a\),\(\varphi(t)\) 为 \(X_i\) 的特征函数。因为 \(\varphi'(0)=ia\),所以 \(\varphi(t)=1+iat+o(t)\)。
\(Y_n=\frac1n\sum X_i\) 的特征函数为
\(e^{iat}\) 是单点分布 \(P(X=a)=1\) 的特征函数,所以 \(Y_n\xrightarrow{L}a\),再由定理 4.1.3 得 \(Y_n\xrightarrow{P}a\)。
习题 4.1.18:\(\{X_n\}\) 独立同分布,\(E(X_n)=0\),\(\mathrm{Var}(X_n)=\sigma^2\),则 \(\frac1n\sum_{k=1}^nX_k^2\xrightarrow{P}\sigma^2\)。(对 \(X_k^2\) 用辛钦大数定律,\(E(X_k^2)=\sigma^2\)。)
4. 中心极限定理
| 名称 | 页码 | 条件 | 结论 |
|---|---|---|---|
| 林德伯格-莱维(会证) | P240 | \(\{X_n\}\) 独立同分布,\(E(X_i)=\mu\),\(\mathrm{Var}(X_i)=\sigma^2>0\) 存在 | \(\dfrac{\sum X_i-n\mu}{\sigma\sqrt n}\xrightarrow{L}N(0,1)\) |
| 棣莫弗-拉普拉斯 | P242 | \(S_n\sim b(n,p)\)(伯努利) | \(\dfrac{S_n-np}{\sqrt{npq}}\xrightarrow{L}N(0,1)\) |
证明:林德伯格-莱维(特征函数法)
令 \(Y_i=\frac{X_i-\mu}{\sigma}\),其特征函数 \(\varphi(t)=1-\frac{t^2}{2}+o(t^2)\)(因为 \(E Y_i=0\),\(E Y_i^2=1\))。
\(Z_n=\frac{1}{\sqrt n}\sum Y_i\) 的特征函数为
这正是 \(N(0,1)\) 的特征函数,由连续性定理即得结论。
应用:二项分布的正态近似(P243)
\[P(a\le S_n\le b)\approx\Phi\left(\frac{b+0.5-np}{\sqrt{npq}}\right)-\Phi\left(\frac{a-0.5-np}{\sqrt{npq}}\right)\](加减 0.5 是连续性修正。)
直观理解:不管单个 \(X_i\) 长什么样,只要独立同分布且方差有限,大量叠加之后就会“磨成”钟形。这也是正态分布在自然界中无处不在的原因。
本系列目录
- 01|概念简答索引
- 02|常用分布
- 03|多维分布、协方差与条件期望
- 04|特征函数、大数定律与中心极限定理(本篇)
- 05|抽样分布、次序统计量与充分统计量
- 06|参数估计
- 07|区间估计、假设检验与习题结论
← 上一篇:多维分布、协方差与条件期望 · 下一篇:抽样分布、次序统计量与充分统计量 →
本系列是我复习茆诗松、程依明、濮晓龙《概率论与数理统计教程》时整理的笔记,共 7 篇。文中 P××× 为茆书页码,“指导书”为配套的学习指导与习题解答。如有错漏,欢迎指出。