04-第四章 大数定律与中心极限定理
依赖于
被以下题目直接调用
正文部分
§4.1 随机变量序列的两种收敛性
- 依概率收敛 设 {Xn} 为一随机变量序列,X 为一随机变量.如果对任意的 ε>0,有
n→∞limP{∣Xn−X∣<ε}=1,
则称 {Xn} 依概率收敛于 X,记作
XnPX.
- 依概率收敛与服从大数定律间的关系 设 {Xn} 为一随机变量序列,记
Yn=n1i=1∑nXi,E(Yn)=n1i=1∑nE(Xi).
则 {Xn} 服从大数定律等价于
Yn−E(Yn)P0.
- 依概率收敛的四则运算 如果
XnPa,YnPb,
则有
1.
Xn±YnPa±b;
-
XnYnPab;
-
Xn/YnPa/b(b=0).
- 按分布收敛,弱收敛 设 {Fn(x)} 是随机变量序列 {Xn} 的分布函数列,F(x) 为 X 的分布函数.若对 F(x) 的任一连续点 x,都有
n→∞limFn(x)=F(x),
则称 {Fn(x)} 弱收敛于 F(x),记作
Fn(x)WF(x).
也称 {Xn} 按分布收敛于 X,记作
XnLX.
- 依概率收敛与按分布收敛间的关系
-
XnPX⟹XnLX;
-
XnPc⟺XnLc(其中 c 为常数).
习题与解答 4.1
如果 XnPX,且 XnPY.试证:P(X=Y)=1.
解
对任意的 ε>0,有
{∣X−Y∣≥ε}⊂({∣X−Xn∣≥ε/2}∪{∣Xn−Y∣≥ε/2}),
故当 n→∞ 时,有
0≤P{∣X−Y∣≥ε}≤P{∣X−Xn∣≥ε/2}+P{∣Xn−Y∣≥ε/2}→0.
即对任意的 ε>0,有
P{∣X−Y∣≥ε}=0.
于是有
P(X=Y)=P(k=1⋃∞{∣X−Y∣≥k1})≤k=1∑∞P{∣X−Y∣≥k1}=0,
从而
P(X=Y)=1
成立,结论得证.
如果 XnPX,YnPY.试证:
- Xn+YnPX+Y;
- XnYnPXY.
解
- 因为
{∣Xn+Yn−X−Y∣≥ε}⊂({∣Xn−X∣≥ε/2}∪{∣Yn−Y∣≥ε/2}),
故当 n→∞ 时,有
0≤P{∣Xn+Yn−X−Y∣≥ε}≤P{∣Xn−X∣≥ε/2}+P{∣Yn−Y∣≥ε/2}→0,
即
Xn+YnPX+Y.
进一步由 YnPY 可得 −YnP−Y,所以又有
Xn−YnPX−Y.
- 先证明
Xn2PX2.
对任意的 ε>0,δ>0,取 M 足够大(譬如 ε/M≤1),使有
P{∣X∣>(M−1)/2}<δ
成立;对取定的 M,存在 N,当 n>N 时,有
P{∣Xn−X∣≥1}<δ,P{∣Xn−X∣≥ε/M}<δ.
这时有
P{∣Xn+X∣>M}≤P{∣Xn−X∣+∣2X∣>M}=P({∣Xn−X∣+∣2X∣>M}∩{∣Xn−X∣<1})+P({∣Xn−X∣+∣2X∣>M}∩{∣Xn−X∣≥1})≤P{∣2X∣>M−1}+P{∣Xn−X∣≥1}<2δ.
从而有
P{∣Xn2−X2∣≥ε}=P{∣Xn−X∣∣Xn+X∣≥ε}=P({∣Xn−X∣∣Xn+X∣≥ε}∩{∣Xn+X∣≤M})+P({∣Xn−X∣∣Xn+X∣≥ε}∩{∣Xn+X∣>M})≤P{∣Xn−X∣≥ε/M}+P{∣Xn+X∣>M}<3δ.
由 ε,δ 的任意性知
Xn2PX2.
同理可证
Yn2PY2.
由上面(1)得
2XnYn=(Xn+Yn)2−Xn2−Yn2P(X+Y)2−X2−Y2=2XY,
即
XnYnPXY.
如果 XnPX,g(x) 是直线上的连续函数,试证:g(Xn)Pg(X).
解
若 g(x) 是 m 次多项式函数,即
g(x)=i=0∑maixi,
则由上一题知有
g(Xn)Pg(X).
下证一般情况.对任意的 ε>0,δ>0,取 M 充分大,使有
P{∣X∣>M}<δ,
又选取 N1 充分大,使当 n≥N1 时,有
P{∣Xn−X∣>1}<δ,
于是有
P{∣Xn∣>M+1}≤P({∣X∣>M}∪{∣Xn−X∣>1})<2δ.
对取定的 M,因为 g(x) 是连续函数,所以可以用多项式去逼近 g(x),并且在任意有界区间上还可以是一致的,因而存在 m 次多项式 gm(x),使得当
−M−1≤x≤M+1
时,有
∣g(x)−gm(x)∣<ε/3.
对取定的 m 次多项式 gm(x),因为
gm(Xn)Pgm(X),
所以存在 N2,使当 n≥N2 时,有
P{∣gm(Xn)−gm(X)∣≥ε/3}<δ.
又因为
P{∣g(Xn)−g(X)∣≥ε}=P({∣g(Xn)−g(X)∣≥ε}∩[{∣X∣>M}∪{∣Xn∣>M+1}])+P({∣g(Xn)−g(X)∣≥ε}∩{∣X∣≤M}∩{∣Xn∣≤M+1})=I1+I2,
当 n≥max{N1,N2} 时,有
I1≤P{∣X∣>M}+P{∣Xn∣>M+1}<3δ,
又因为
{∣g(Xn)−g(X)∣≥ε}⊂{∣g(Xn)−gm(Xn)∣≥ε/3}∪{∣gm(Xn)−gm(X)∣≥ε/3}∪{∣gm(X)−g(X)∣≥ε/3},
且
P({∣g(Xn)−gm(Xn)∣≥ε/3}∩{∣X∣≤M}∩{∣Xn∣≤M+1})=0,
P({∣gm(X)−g(X)∣≥ε/3}∩{∣X∣≤M}∩{∣Xn∣≤M+1})=0,
所以
I2≤P({∣gm(Xn)−gm(X)∣≥ε/3}∩{∣X∣≤M}∩{∣Xn∣≤M+1})≤P{∣gm(Xn)−gm(X)∣≥ε/3}<δ.
从而有
P{∣g(Xn)−g(X)∣≥ε}=I1+I2<4δ.
由 ε,δ 的任意性即知
g(Xn)Pg(X),
结论得证.
如果 XnPa,则对任意常数 c,有
cXnPca.
解
记
g(x)=cx,
则 g(x) 是连续函数,由上一题即可得
cXn=g(Xn)Pg(a)=ca.
试证:XnPX 的充要条件为:当 n→∞ 时,有
E(1+∣Xn−X∣∣Xn−X∣)→0.
解
先证充分性.令
f(x)=1+xx,x>0,
则
f′(x)=(1+x)21>0,x>0,
故 f(x) 是 x 的严格单调增函数,因而对任意的 ε>0,有
{∣Xn−X∣>ε}⊂{1+∣Xn−X∣∣Xn−X∣>1+εε}.
于是对任意的 ε>0,当 n→∞ 时,有(参见 \S 2.3 第 12 题)
P{∣Xn−X∣>ε}≤P{1+∣Xn−X∣∣Xn−X∣>1+εε}≤ε1+εE(1+∣Xn−X∣∣Xn−X∣)→0.
充分性得证.
下证必要性.对任意的 ε>0,令
Aε={∣Xn−X∣>ε}.
因为 XnPX,故存在充分大的 N,使得当 n≥N 时,有
P(Aε)<ε,
于是有
E(1+∣Xn−X∣∣Xn−X∣)=E(1+∣Xn−X∣∣Xn−X∣IAε)+E(1+∣Xn−X∣∣Xn−X∣IAε)≤P(Aε)+ε<2ε.
由 ε 的任意性知,当 n→∞ 时,有
E(1+∣Xn−X∣∣Xn−X∣)→0.
结论得证.
设 D(x) 为退化分布:
D(x)={0,1,x<0,x≥0.
试问下列分布函数列的极限函数是否仍是分布函数?(其中 n=1,2,…)
- {D(x+n)};
- {D(x+1/n)};
- {D(x−1/n)}.
解
- 因为此时的极限函数为
G(x)=1,−∞<x<∞,
不满足分布函数的基本性质:
x→−∞limF(x)=0,
所以不是分布函数.
- 因为此时的极限函数为
G(x)={0,1,x<0,x≥0,
所以是分布函数.
- 因为此时的极限函数为
G(x)={0,1,x≤0,x>0,
不满足分布函数的右连续性,所以不是分布函数.
设分布函数列 {Fn(x)} 弱收敛于连续的分布函数 F(x),试证:{Fn(x)} 在 (−∞,∞) 上一致收敛于分布函数 F(x).
解
对任意的 ε>0,取 M 充分大,使有
当 x≥M 时,有 1−F(x)<ε;当 x≤−M 时,有 F(x)<ε.
对上述取定的 M,因为 F(x) 在闭区间 [−M,M] 上一致连续,故可取它的 k 个分点
−M=x1<x2<⋯<xk−1<xk=M,
使有
F(xi+1)−F(xi)<ε.
再令
x0=−∞,xk+1=+∞,
则对这些分点有
F(xi+1)−F(xi)<ε.
这时存在 N,使得当 n>N 时,有
∣Fn(xi)−F(xi)∣<ε.
对任意的 x∈(−∞,∞),必存在某个 i,使得
x∈(xi,xi+1],
由此知,当 n>N 时,有
F(xi)−ε<Fn(xi)≤Fn(x)≤Fn(xi+1)<F(xi+1)+ε.
由上式立刻得
Fn(x)−F(x)<F(xi+1)−F(x)+ε≤F(xi+1)−F(xi)+ε<2ε,
同时
Fn(x)−F(x)>F(xi)−F(x)−ε≥F(xi)−F(xi+1)−ε>−2ε.
即有
∣Fn(x)−F(x)∣<2ε,
结论得证.
如果 XnLX,且数列 an→a,bn→b.试证:anXn+bnLaX+b.
解
先证明
aXnLaX.
因为 a=0 时为显然,所以不妨设 a>0(a<0 时的修改为显然).因为
FaXn(x)=FXn(x/a),
且当 x 是 FaX(⋅) 的连续点时,则 x/a 是 FX(⋅) 的连续点,于是有
n→∞limFaXn(x)=n→∞limFXn(x/a)=FX(x/a)=FaX(x),
此即
aXnLaX.
由本节第 10 题知
Xn(an−a)P0.
再由本节第 2 题(1)知
(an−a)Xn+bnPb.
于是由前述结论及本节第 9 题知
anXn+bn=aXn+(an−a)Xn+bnLaX+b,
结论得证.
如果 XnLX,YnPa,试证:Xn+YnLX+a.
解
对任意的 ε>0,首先考虑 Xn+Yn 的分布函数
P(Xn+Yn≤z)=P(Xn+Yn≤z, ∣Yn−a∣≤ε)+P(Xn+Yn≤z, ∣Yn−a∣>ε)≤P(Xn≤z−a+ε, ∣Yn−a∣≤ε)+P(∣Yn−a∣>ε)≤P(Xn≤z−a+ε)+P(∣Yn−a∣>ε).
因此
n→∞limsupP(Xn+Yn≤z)≤n→∞limP(Xn≤z−a+ε)+n→∞limP(∣Yn−a∣>ε)=F(z−a+ε),
其中 F(⋅) 为 X 的分布函数.
类似有
P(Xn+Yn≤z)≥P(Xn+Yn≤z, ∣Yn−a∣≤ε)≥P(Xn≤z−a−ε, ∣Yn−a∣≤ε)=P(Xn≤z−a−ε)−P(Xn≤z−a−ε, ∣Yn−a∣>ε)≥P(Xn≤z−a−ε)−P(∣Yn−a∣>ε).
因此
n→∞liminfP(Xn+Yn≤z)≥n→∞limP(Xn≤z−a−ε)−n→∞limP(∣Yn−a∣>ε)≥F(z−a−ε).
由上述两个关系式,再考虑到 ε 的任意性,即可得
n→∞limP(Xn+Yn≤z)=F(z−a).
这就意味着
Xn+YnLX+a.
证毕.
如果 XnLX,YnP0.试证:XnYnP0.
解
记 Xn 与 X 的分布函数分别为 Fn(x) 和 F(x).对任给的 ε>0,取足够大的 a>0 和 b>0,使 −a,b 是 F(x) 的连续点,且
1−F(b)<ε,F(−a)<ε.
因为 Fn(x)WF(x),故存在 N1,使当 n≥N1 时,有
1−Fn(b)<2ε,Fn(−a)<2ε.
令
M=max{a,b},
因为 YnP0,故存在 N2,使当 n≥N2 时,有
P{∣Yn∣>ε/M}<ε.
而
P{∣XnYn∣>ε}=P({∣XnYn∣>ε}∩{−a<Xn≤b}∩{∣Yn∣≤ε/M})+P({∣XnYn∣>ε}∩[{Xn≤−a}∪{Xn>b}∪{∣Yn∣>ε/M}])=I1+I2.
由 M 的定义即知
{∣XnYn∣>ε}∩{−a<Xn≤b}∩{∣Yn∣≤ε/M}=∅,
所以有 I1=0.而对于 I2,当 n≥max{N1,N2} 时,有
P({∣XnYn∣>ε}∩{−a<Xn≤b})≤P({−a<Xn≤b})=P({Xn≤−a}∪{Xn>b})=Fn(−a)+[1−Fn(b)]<4ε,
且
P({∣XnYn∣>ε}∩{∣Yn∣>ε/M})≤P{∣Yn∣>ε/M}<ε,
因而
P{∣XnYn∣>ε}=I2<5ε.
由 ε 的任意性知
XnYnP0,
结论得证.
如果 XnLX,YnPa,且 Yn=0,常数 a=0.试证:Xn/YnLX/a.
解
先证
1/YnP1/a.
不妨设 a>0.对任意的 0<ε<a,当 ∣Yn−a∣<ε 时,有
∣Yna∣=∣a2+a(Yn−a)∣≥a2−aε.
因而
{YnaYn−a≥ε}⊂{a2−aεYn−a≥ε}.
于是当 n→∞ 时,有
0≤P{Yn1−a1≥ε}=P{YnaYn−a≥ε}=P({YnaYn−a≥ε}∩{∣Yn−a∣<ε})+P({YnaYn−a≥ε}∩{∣Yn−a∣≥ε})≤P{a2−aεYn−a≥ε}+P(∣Yn−a∣≥ε)→0.
所以
1/YnP1/a.
于是由本节第 10 题,有
Xn(Yn1−a1)P0.
又由本节第 8 题的证明知
Xn/aLX/a,
因而由本节第 9 题知
YnXn=Xn(Yn1−a1)+aXnLaX,
结论得证.
设随机变量 Xn 服从柯西分布,其密度函数为
pn(x)=π(1+n2x2)n,−∞<x<∞.
试证:XnP0.
解
对任意的 ε>0,当 n→∞ 时,有
P(∣Xn∣≤ε)=∫−εεπ(1+n2x2)ndx=∫−nεnεπ(1+t2)1dt→1.
即
XnP0,
结论得证.
设随机变量序列 {Xn} 独立同分布,其密度函数为
p(x)={1/β,0,0<x<β,其他,
其中常数 β>0,令
Yn=max{X1,X2,…,Xn},
试证:YnPβ.
解
因为当 x<0 时,有
P(Yn≤x)=0;
当 x≥β 时,有
P(Yn≤x)=1;
而当 0≤x<β 时,有
P(Yn≤x)=i=1∏nP(Xi≤x)=i=1∏n∫0xβ1dx=(βx)n.
所以,对任意的 ε>0(ε<β),当 n→∞ 时,有
P(∣Yn−β∣≥ε)=P(Yn≤β−ε)=(ββ−ε)n→0,
所以有
YnPβ,
结论得证.
设随机变量序列 {Xn} 独立同分布,其密度函数为
p(x)={e−(x−α),0,x≥α,x<α.
令
Yn=min{X1,X2,…,Xn},
试证:YnPα.
解
因为 Xi 的分布函数为
F(x)={1−e−(x−α),0,x≥α,x<α.
所以当 x≥α 时,有
P(Yn≥x)=i=1∏nP(Xi≥x)=[1−F(x)]n=e−n(x−α).
对任意的 ε>0,当 n→∞ 时,有
P(∣Yn−α∣≥ε)=P(Yn−α≥ε)=e−nε→0,
即
YnPα,
结论得证.
设随机变量序列 {Xn} 独立同分布,且 Xi∼U(0,1).令
Yn=(i=1∏nXi)1/n,
试证明:
YnPc,
其中 c 为常数,并求出 c.
解
因为
lnYn=n1i=1∑nlnXi,
且
E(lnXi)=∫01lnxdx=[xlnx−x]01=−1,
E(lnXi)2=∫01(lnx)2dx=x[(lnx)2]01−2∫01lnxdx=2,
从而
Var(lnXi)=1.
所以由切比雪夫不等式得,对任意 ε>0,有
P(n1i=1∑nlnXi−(−1)≥ε)≤nε21→0.
即
lnYnP−1.
再由本节第 3 题知
Yn=(i=1∏nXi)1/nPe−1,
即
c=e−1.
设分布函数列 {Fn(x)} 弱收敛于分布函数 F(x),且 Fn(x) 和 F(x) 都是连续、严格单调函数,又设 ξ 服从 (0,1) 上的均匀分布,试证:Fn−1(ξ)PF−1(ξ).
解
对任意的 ε>0,δ>0,存在充分大的 M,使有
F(M)−F(−M)>1−δ.
对取定的 M,可选取正整数 k 和 N,使有
1/k<ε,N/k>M.
对取定的 N,存在
h>0,h<δ/[2(2N+1)].
对取定的 h,因为 Fn(x)→F(x) 关于 x 是一致的(见本节第 7 题),因而存在 N1,使当 n≥N1 时,任对 x∈(−∞,∞),有
∣Fn(x)−F(x)∣<h.
因此有
P(∣Fn−1(ξ)−F−1(ξ)∣≤ε)≥P(∣Fn−1(ξ)−F−1(ξ)∣≤1/k)≥i=−N∑NP({Fn−1(ξ)−ki≤2k1}∩{F−1(ξ)−ki≤2k1})=i=−N∑NP({Fn(ki−2k1)≤ξ≤Fn(ki+2k1)}∩{F(ki−2k1)≤ξ≤F(ki+2k1)})≥i=−N∑NP(F(ki−2k1)+h≤ξ≤F(ki+2k1)−h)=i=−N∑N[F(ki+2k1)−F(ki−2k1)−2h]=F(kN+2k1)−F(−kN−2k1)−2(2N+1)h≥1−2δ.
由 ε,δ 的任意性知
Fn−1(ξ)PF−1(ξ).
结论得证.
设随机变量序列 {Xn} 独立同分布,数学期望、方差均存在,且 E(Xn)=μ.试证:
n(n+1)2k=1∑nkXkPμ.
解
已知
E(Xn)=μ,
记
Var(Xn)=σ2,
令
Yn=n(n+1)2k=1∑nkXk,
则
E(Yn)=n(n+1)2k=1∑nkμ=μ,
Var(Yn)=n2(n+1)24k=1∑nk2σ2=(n+1)24σ2k=1∑nn2k2≤n+14σ2.
对任意的 ε>0,由切比雪夫不等式得
P(∣Yn−μ∣≥ε)≤ε2Var(Yn)≤ε21⋅n+14σ2→0,n→∞.
即
YnPμ,
结论得证.
设随机变量序列 {Xn} 独立同分布,数学期望、方差均存在,且
E(Xn)=0,Var(Xn)=σ2.
试证:
n1k=1∑nXk2Pσ2.
解
这时 {Xn2} 仍独立同分布,且
E(Xn2)=Var(Xn)=σ2<∞,
由辛钦大数定律(见 \S 4.3)知结论成立.
设随机变量序列 {Xn} 独立同分布,且 Var(Xn)=σ2 存在,令
Xn=n1i=1∑nXi,Sn2=n1i=1∑n(Xi−Xn)2.
试证:Sn2Pσ2.
解
不妨设 E(Xn)=0,否则令
Xn′=Xn−E(Xn),
并以 {Xn′} 代替 {Xn},这时 σ2,Sn2 均保持不变.
易知
Sn2=n1k=1∑n(Xk−Xn)2=n1k=1∑nXk2−(Xn)2.
由本节第 18 题知
n1k=1∑nXk2Pσ2.
又因为
Xn=n1i=1∑nXiPn1i=1∑nE(Xi)=0,
所以由本节第 2 题(2)知
(Xn)2P0.
最后由本节第 2 题(1)知
Sn2=n1k=1∑nXk2−(Xn)2Pσ2,
结论得证.
将 n 个编号为 1 至 n 的球放入 n 个编号为 1 至 n 的盒子中,每个盒子只能放一个球,记
Xi={1,0,编号为 i 的球放入编号为 i 的盒子,反之.
Sn=i=1∑nXi,
试证明:
nSn−E(Sn)P0.
解
这是一个配对模型,由习题 3.4 第 25 题知:
E(Sn)=1,Var(Sn)=1.
所以由切比雪夫不等式,对任意 ε>0,有
P(nSn−E(Sn)≥ε)≤n2ε21→0.
即
nSn−E(Sn)P0.
§4.2 特征函数
- 特征函数的定义 设 X 是一个随机变量,称
φ(t)=E(eitX)
为 X 的特征函数,其表达式如下
φ(t)=E(eitX)=⎩⎨⎧∑ieitxiP(X=xi),∫−∞∞eitxp(x)dx,在离散场合,在连续场合,−∞<t<∞.
由于
∣eitx∣=cos2tx+sin2tx=1,
所以随机变量 X 的特征函数 φ(t) 总是存在的.
利用欧拉公式
eiθ=cosθ+isinθ
可把不少问题中的复变函数问题转化为实变函数问题进行处理.
上述由密度函数 p(x) 求其特征函数的公式常称傅里叶变换.
-
特征函数的性质
-
∣φ(t)∣≤φ(0)=1;
-
φ(−t)=φ(t),
其中 φ(t) 表示 φ(t) 的共轭;
-
若 Y=aX+b,其中 a,b 是常数,则
φY(t)=eibtφX(at);
- 若 X 与 Y 是相互独立的随机变量,则
φX+Y(t)=φX(t)φY(t);
- 若 E(Xℓ) 存在,则 φ(t) 可 ℓ 次求导,且对 1≤k≤ℓ,有
φ(k)(0)=ikE(Xk);
- 一致连续性:特征函数 φ(t) 在 (−∞,∞) 上一致连续;
- 非负定性:特征函数 φ(t) 是非负定的,即对任意正整数 n,及 n 个实数 t1,t2,…,tn 和 n 个复数 z1,z2,…,zn,有
k=1∑nj=1∑nφ(tk−tj)zkzj≥0;
- 逆转公式 设 F(x) 和 φ(t) 分别为随机变量 X 的分布函数和特征函数,则对 F(x) 的任意两个连续点 x1<x2,有
F(x2)−F(x1)=T→∞lim2π1∫−TTite−itx1−e−itx2φ(t)dt;
- 唯一性定理 随机变量的分布函数由其特征函数唯一决定;
- 若连续随机变量 X 的密度函数为 p(x),特征函数为 φ(t).如果
∫−∞∞∣φ(t)∣dt<∞,
则
p(x)=2π1∫−∞∞e−itxφ(t)dt.
这个公式又称傅里叶逆变换;
- 分布函数序列 {Fn(x)} 弱收敛于分布函数 F(x) 的充要条件是 {Fn(x)} 的特征函数序列 {φn(t)} 收敛于 F(x) 的特征函数 φ(t).
- 常用分布的特征函数表
分布单点分布0-1 分布二项分布 b(n,p)泊松分布 P(λ)几何分布 Ge(p)分布列 pk 或分布密度 p(x)P(X=a)=1pk=pkq1−k, q=1−p, k=0,1pk=(kn)pkqn−k, k=0,1,…,npk=k!λke−λ, k=0,1,…pk=pqk−1, k=1,2,…特征函数 φ(t)eitapeit+q(peit+q)neλ(eit−1)1−qeitp
续表
| 分布 | 分布列 pk 或分布密度 p(x) | 特征函数 φ(t) |
|---|
| 负二项分布 Nb(r,p) | pk=(r−1k−1)prqk−r,k=r,r+1,⋯ | (1−qeitp)r |
| 均匀分布 U(a,b) | p(x)=b−a1,a<x<b | it(b−a)eibt−eiat |
| 均匀分布 U(−a,a) | p(x)=2a1,−a<x<a | atsinat |
| 正态分布 N(μ,σ2) | p(x)=2πσ1exp{−2σ2(x−μ)2} | exp{iμt−2σ2t2} |
| 标准正态分布 N(0,1) | p(x)=2π1exp{−2x2} | e−t2/2 |
| 指数分布 Exp(λ) | p(x)=λe−λx,x≥0 | (1−λit)−1 |
| 伽马分布 Ga(α,λ) | p(x)=Γ(α)λαxα−1e−λx,x≥0 | (1−λit)−α |
| χ2(n) 分布 | p(x)=Γ(n/2)2n/2xn/2−1e−x/2,x≥0 | (1−2it)−n/2 |
| 贝塔分布 Be(a,b) | p(x)=Γ(a)Γ(b)Γ(a+b)xa−1(1−x)b−1,0<x<1 | Γ(a)Γ(a+b)k=0∑∞k!Γ(a+b+k)(it)kΓ(a+k) |
| 柯西分布 Cau(0,1) | p(x)=π(1+x2)1,−∞<x<∞ | e−∣t∣ |
习题与解答 4.2
设离散随机变量 X 的分布列如下,试求 X 的特征函数。
XP00.410.320.230.1
解
φX(t)=0.4+0.3eit+0.2ei2t+0.1ei3t.
设离散随机变量 X 服从几何分布
P(X=k)=(1−p)k−1p,k=1,2,⋯.
试求 X 的特征函数,并以此求 E(X) 和 \Var(X)。
解
记 q=1−p,则
φ(t)=E(eitX)=k=1∑∞eitkqk−1p=peitk=1∑∞(eitq)k−1=1−qeitpeit.
它的前二阶导数为
φ′(t)=(1−qeit)2ipeit,
φ′′(t)=(1−qeit)4−peit(1−qeit)2−2peit(1−qeit)qeit,
由此可算得几何分布的期望和方差为
E(X)=i1φ′(0)=(1−q)2p=p1,
E(X2)=i21φ′′(0)=(1−q)4p(1−q)2+2pq(1−q)=p21+q,
\Var(X)=E(X2)−[E(X)]2=p21+q−(p1)2=p2q.
设离散随机变量 X 服从帕斯卡分布
P(X=k)=(r−1k−1)pr(1−p)k−r,k=r,r+1,⋯.
试求 X 的特征函数。
解
设 X1,X2,⋯,Xr 是相互独立同分布的随机变量,且都服从参数为 p 的几何分布 Ge(p),则由上一题知 Xj 的特征函数为
φXj(t)=1−qeitpeit,
其中 q=1−p。又因为
X=X1+X2+⋯+Xr,
所以 X 的特征函数为
φX(t)=j=1∏rφXj(t)=(1−qeitpeit)r.
求下列分布函数的特征函数,并由特征函数求其数学期望和方差:
1.
F1(x)=2a∫−∞xe−a∣t∣dt(a>0);
-
F2(x)=πa∫−∞xt2+a21dt(a>0).
解
**(1)**因为此分布的密度函数为
p1(x)=2ae−a∣x∣,−∞<x<∞.
所以此分布的特征函数为
φ1(t)=2a∫−∞0eitxeaxdx+2a∫0∞eitxe−axdx=2a∫−∞0(costx+isintx)eaxdx+2a∫0∞(costx+isintx)e−axdx=a∫0∞costx⋅e−axdx=a2+t2a2.
又因为
φ1′(t)=−(a2+t2)22ta2,φ1′(0)=0,
φ1′′(t)=(a2+t2)32a2(3t2−a2),φ1′′(0)=−a22,
所以
E(X)=i1φ1′(0)=0,\Var(X)=E(X2)=i21φ1′′(0)=a22.
**(2)**因为此分布的密度函数为
p2(x)=πa⋅x2+a21,−∞<x<∞.
所以此分布的特征函数为
φ2(t)=πa∫−∞∞x2+a2eitxdx=π2a∫0∞x2+a2costxdx,
又因为当 t>0 时,有(见菲赫金哥尔茨《微积分学教程》第二卷第三分册或查积分表)
∫0∞x2+a2costxdx=2aπe−at,
所以当 t>0 时,有
φ2(t)=π2a⋅2aπe−at=e−at.
而当 t<0 时,有
φ2(t)=φ2(−t)=e−a∣t∣,
所以
φ2(t)=e−a∣t∣.
又因为 φ2(t) 在 t=0 处不可导,故此分布(柯西分布)的数学期望不存在。
设随机变量 X∼N(μ,σ2),试用特征函数的方法求 X 的 3 阶及 4 阶中心矩。
解
因为正态分布 N(μ,σ2) 的特征函数为
φ(t)=eiμt−σ2t2/2,
所以
φ′(0)=iμ,E(X)=iφ′(0)=μ,
φ′′(0)=−μ2−σ2,E(X2)=i2φ′′(0)=μ2+σ2,
φ′′′(0)=−iμ3−3iμσ2,E(X3)=i3φ′′′(0)=μ3+3μσ2,
φ(4)(0)=μ4+6μ2σ2+3σ4,E(X4)=i4φ(4)(0)=μ4+6μ2σ2+3σ4.
由此得 X 的 3 阶及 4 阶中心矩为
E(X−E(X))3=E(X3)−3E(X2)μ+3E(X)μ2−μ3=0,
E(X−E(X))4=E(X4)−4E(X3)μ+6E(X2)μ2−4E(X)μ3+μ4=3σ4.
试用特征函数的方法证明二项分布的可加性:若随机变量 X∼b(n,p),Y∼b(m,p),且 X 与 Y 独立,则 X+Y∼b(n+m,p)。
解
证 记 q=1−p,因为
φX(t)=(peit+q)n,φY(t)=(peit+q)m,
所以由 X 与 Y 的独立性得
φX+Y(t)=φX(t)φY(t)=(peit+q)n+m,
这正是二项分布 b(n+m,p) 的特征函数,由唯一性定理知
X+Y∼b(n+m,p).
试用特征函数的方法证明泊松分布的可加性:若随机变量 X∼P(λ1),Y∼P(λ2),且 X 与 Y 独立,则 X+Y∼P(λ1+λ2)。
解
证 因为
φX(t)=eλ1(eit−1),φY(t)=eλ2(eit−1),
所以由 X 与 Y 的独立性得
φX+Y(t)=φX(t)φY(t)=e(λ1+λ2)(eit−1),
这正是泊松分布 P(λ1+λ2) 的特征函数,由唯一性定理知
X+Y∼P(λ1+λ2).
试用特征函数的方法证明伽马分布的可加性:若随机变量 X∼Ga(α1,λ),Y∼Ga(α2,λ),且 X 与 Y 独立,则 X+Y∼Ga(α1+α2,λ)。
解
证 因为
φX(t)=(1−λit)−α1,φY(t)=(1−λit)−α2,
所以由 X 与 Y 的独立性得
φX+Y(t)=φX(t)φY(t)=(1−λit)−(α1+α2),
这正是伽马分布 Ga(α1+α2,λ) 的特征函数,由唯一性定理知
X+Y∼Ga(α1+α2,λ).
试用特征函数的方法证明 χ2 分布的可加性:若随机变量 X∼χ2(n),Y∼χ2(m),且 X 与 Y 独立,则 X+Y∼χ2(n+m)。
解
证 因为
φX(t)=(1−2it)−n/2,φY(t)=(1−2it)−m/2,
所以由 X 与 Y 的独立性得
φX+Y(t)=φX(t)φY(t)=(1−2it)−(n+m)/2,
这正是 χ2 分布 χ2(n+m) 的特征函数,由唯一性定理知
X+Y∼χ2(n+m).
设随机变量 Xi 独立同分布,且 Xi∼Exp(λ),i=1,2,⋯,n。试用特征函数的方法证明:
Yn=i=1∑nXi∼Ga(n,λ).
解
证 因为
φXi(t)=(1−λit)−1,
所以由诸 Xi 的相互独立性得 Yn 的特征函数为
φYn(t)=(1−λit)−n,
这正是伽马分布 Ga(n,λ) 的特征函数,由唯一性定理知
Yn∼Ga(n,λ).
设连续随机变量 X 服从柯西分布,密度函数如下:
p(x)=π1⋅λ2+(x−μ)2λ,−∞<x<∞,
其中参数 λ>0,−∞<μ<∞,常记为 X∼Cau(λ,μ)。
- 试证 X 的特征函数为 exp{iμt−λ∣t∣},且利用此结果证明柯西分布的可加性;
- 当 μ=0,λ=1 时,记 Y=X,试证 φX+Y(t)=φX(t)⋅φY(t),但是 X 与 Y 不独立;
- 若 X1,X2,⋯,Xn 相互独立,且服从同一柯西分布,试证:
n1(X1+X2+⋯+Xn)
与 X1 同分布。
解
**(1)**因为 Y=X−μ 的密度函数为
p(y)=π1⋅λ2+y2λ,−∞<y<∞,
所以由本节第 4 题(2)知 Y 的特征函数为
φY(t)=exp{−λ∣t∣}.
由此得 X=Y+μ 的特征函数
φX(t)=φY+μ(t)=exp{iμt}φY(t)=exp{iμt−λ∣t∣}.
下面利用此式证明柯西分布的可加性。设 Xi (i=1,2) 服从参数为 μi,λi 的柯西分布,其密度函数为
pi(x)=π1⋅λi2+(x−μi)2λi,−∞<x<∞,i=1,2.
若 X1 与 X2 相互独立,则
φX1+X2(t)=φX1(t)φX2(t)=exp{i(μ1+μ2)t−(λ1+λ2)∣t∣},
这正是参数为 μ1+μ2,λ1+λ2 的柯西分布的特征函数。所以由唯一性定理知,X1+X2 服从参数为 μ1+μ2,λ1+λ2 的柯西分布。
**(2)**当 μ=0,λ=1 时有
φX(t)=exp{−∣t∣},φY(t)=exp{−∣t∣},
所以
φX+Y(t)=φ2X(t)=φX(2t)=exp{−2∣t∣}=exp{−∣t∣}exp{−∣t∣}=φX(t)φY(t).
由于 Y=X,当然 X 与 Y 不独立。此题说明,由 φX+Y(t)=φX(t)φY(t) 不能推得 X 与 Y 独立。
**(3)**设 Xi 都服从参数为 μ,λ 的柯西分布,则特征函数为
φ(t)=exp{iμt−λ∣t∣}.
由相互独立性得:
n1i=1∑nXi
的特征函数为
[φ(t/n)]n=exp{iμt−λ∣t∣},
即
n1i=1∑nXi
与 X1 具有相同的特征函数,由唯一性定理知它们具有相同的分布。
设连续随机变量 X 的密度函数为 p(x),试证:p(x) 关于原点对称的充要条件是它的特征函数是实的偶函数。
解
证 记 X 的特征函数为 φX(t)。
先证充分性。若 φX(t) 是实的偶函数,则
φX(−t)=φX(t),
又因为 φX(−t)=φ−X(t),这表明 X 与 −X 有相同的特征函数,从而 X 与 −X 有相同的密度函数,而 −X 的密度函数为 p(−x),所以得
p(x)=p(−x),
即 p(x) 关于原点是对称的。
再证必要性。若 p(x)=p(−x),则 X 与 −X 有相同的密度函数,所以 X 与 −X 有相同的特征函数。由于 −X 的特征函数为 φX(−t),所以
φX(t)=φX(−t)=φX(t),
故 φX(t) 是实的偶函数。
设随机变量 X1,X2,⋯,Xn 独立同分布,且都服从 N(μ,σ2) 分布,试求
X=n1i=1∑nXi
的分布。
解
因为 Xj 的特征函数为
φj(t)=eiμt−σ2t2/2,
所以由诸 Xj 的相互独立性得 X 的特征函数为
φX(t)=[φj(t/n)]n=eiμt−σ2t2/(2n),
这正是正态分布 N(μ,σ2/n) 的特征函数,所以由唯一性定理知
X=n1i=1∑nXi∼N(μ,nσ2).
利用特征函数方法证明如下的泊松定理:设有一列二项分布 {b(k,n,pn)},若
n→∞limnpn=λ,
则
n→∞limb(k,n,pn)=k!λke−λ,k=0,1,2,⋯.
解
证 二项分布 b(k,n,pn) 的特征函数为
φn(t)=(pneit+qn)n,其中 qn=1−pn.
因为
φn(t)=[1+nnpn(eit−1)]n,
又当 n→∞ 时,npn→λ,则
n→∞limφn(t)=exp{λ(eit−1)}=φ(t).
而 φ(t) 正是泊松分布的特征函数,故得证。
设随机变量 X∼Ga(α,λ),证明:当 α→∞ 时,随机变量
αλX−α
按分布收敛于标准正态变量。
解
证 令
Yα=αλX−α,
则由 X 的特征函数
φX(t)=(1−λit)−α
可得
φYα(t)=e−itα(1−αit)−α.
两边取对数,并将 ln(1−it/α) 展开为级数形式,可得
lnφYα(t)=−itα−αln(1−αit)=−itα−α(−αit+2αt2−3ααit3+o(ααt3))=−2t2+3αit3+α⋅o(ααt3)→−2t2(α→∞).
所以
φYα(t)→e−t2/2,
而 e−t2/2 正是 N(0,1) 的特征函数。由特征函数的唯一性定理及判断弱收敛的方法知结论成立。
§4.3 大数定律
\textbf{1. 随机变量序列 {Xn} 服从大数定律} 设 {Xn} 为随机变量序列,若对任意的 ε>0,有
n→∞limP(n1i=1∑nXi−n1i=1∑nE(Xi)<ε)=1,
则称 {Xn} 服从大数定律。
2. 伯努利大数定律 设 Sn 为 n 重伯努利试验中事件 A 发生的次数,p 为每次试验中 A 出现的概率,则对任意的 ε>0,有
n→∞limP(nSn−p<ε)=1,
这是第一个大数定律,它表明:事件发生的频率依概率收敛于该事件的概率,这就是“频率稳定于概率”的含义,也是“用频率去估计概率”的依据。
3. 切比雪夫大数定律 设 {Xn} 为一列两两不相关的随机变量序列,若每个 Xi 的方差存在,且有共同的上界,即 \Var(Xi)≤c, i=1,2,⋯,则 {Xn} 服从大数定律。
4. 马尔可夫大数定律 对随机变量序列 {Xn},若有
n21\Var(i=1∑nXi)→0(n→∞),
则 {Xn} 服从大数定律。上式称为马尔可夫条件。
5. 辛钦大数定律 设 {Xn} 为一独立同分布的随机变量序列,若 Xi 的数学期望存在,则 {Xn} 服从大数定律。
习题与解答 4.3
设 {Xk} 为独立随机变量序列,且
P(Xk=±lnk)=21,k=1,2,⋯.
证明 {Xk} 服从大数定律。
解
因为 X1,X2,⋯ 相互独立,且 E(Xk)=0, \Var(Xk)=lnk,所以
\Var(k=1∑nXk)=k=1∑nlnk≤n×lnn,
由此可得马尔可夫条件
n21\Var(k=1∑nXk)≤nlnn→0(n→∞).
由马尔可夫大数定律知 {Xn} 服从大数定律。
设 {Xk} 为独立随机变量序列,且
P(Xk=±2k)=22k+11,P(Xk=0)=1−22k1,k=1,2,⋯,
证明 {Xk} 服从大数定律。
解
因为 X1,X2,⋯ 相互独立,且 E(Xk)=0, \Var(Xk)=1,由此可得马尔可夫条件
n21\Var(k=1∑nXk)=n1→0(n→∞).
由马尔可夫大数定律知 {Xn} 服从大数定律。
设 {Xn} 为独立随机变量序列,且 P(X1=0)=1,
P(Xn=±n)=n1,P(Xn=0)=1−n2,n=2,3,⋯,
证明 {Xn} 服从大数定律。
解
因为 X1,X2,⋯ 相互独立,且
E(Xn)=0,\Var(Xn)=2,n=2,3,⋯,
故可得马尔可夫条件
n21\Var(i=1∑nXi)=n22(n−1)→0(n→∞).
由马尔可夫大数定律知 {Xn} 服从大数定律。
在伯努利试验中,事件 A 出现的概率为 p,令
Xn={1,0,若在第 n 次及第 n+1 次试验中 A 都出现,其他.
证明 {Xn} 服从大数定律。
解
{Xn} 为同分布随机变量,其共同分布为
XnP01−p21p2
且
E(Xn)=E(Xn2)=p2,
从而
\Var(Xn)=p2(1−p2)≤1.
又当 ∣i−j∣≥2 时,Xi 与 Xj 独立,所以
n21\Var(i=1∑nXi)=n21[i=1∑n\Var(Xi)+2i=1∑n−1\Cov(Xi,Xi+1)].
又因为
∣\Cov(Xi,Xi+1)∣≤\Var(Xi)\Var(Xi+1)=p2(1−p2),
于是有
n21\Var(i=1∑nXi)≤n21[np2(1−p2)+2(n−1)p2(1−p2)]→0(n→∞),
即马尔可夫条件成立,故 {Xn} 服从大数定律。
设 {Xn} 为独立的随机变量序列,且
P(Xn=1)=pn,P(Xn=0)=1−pn,n=1,2,⋯.
证明 {Xn} 服从大数定律。
解
因为
E(Xn)=pn,\Var(Xn)=pn(1−pn)≤41,
所以由 X1,X2,⋯ 的独立性可得
n21\Var(k=1∑nXk)≤4n1→0(n→∞).
由马尔可夫大数定律知 {Xn} 服从大数定律。
设 {Xn} 为独立同分布的随机变量序列,其共同的分布函数为
F(x)=21+π1arctanax,−∞<x<∞.
试问:辛钦大数定律对此随机变量序列是否适用?
解
此为柯西分布的分布函数,而柯西分布的数学期望不存在。因为辛钦大数定律要求数学期望存在,所以辛钦大数定律对此随机变量序列不适用。
设 {Xn} 为独立同分布的随机变量序列,其共同分布为
P(Xn=k22k)=2k1,k=1,2,⋯.
试问 {Xn} 是否服从大数定律?
解
因为
E(Xn)=k=1∑∞k22k⋅2k1=k=1∑∞k21=6π2<∞,
即 E(Xn) 存在,所以由辛钦大数定律知 {Xn} 服从大数定律。
设 {Xn} 为独立同分布的随机变量序列,其共同分布为
P(Xn=k)=k2⋅lg2kc,k=2,3,⋯,
其中
c=(k=2∑∞k2⋅lg2k1)−1,
试问 {Xn} 是否服从大数定律?
解
因为
E(Xn)=k=2∑∞k⋅k2lg2kc=ck=2∑∞klg2k1,
由柯西积分判别法知上述级数收敛,故 E(Xn) 存在,所以由辛钦大数定律知 {Xn} 服从大数定律。
设 {Xn} 为独立同分布的随机变量序列,其中 Xn 服从参数为 n 的泊松分布,试问 {Xn} 是否服从大数定律?
解
因为
n21\Var(i=1∑nXi)≤n2nn→0(n→∞).
所以由马尔可夫大数定律知 {Xn} 服从大数定律。
设 {Xn} 为独立的随机变量序列,证明:若诸 Xn 的方差 σn2 一致有界,即存在常数 c,使得
σn2≤c,n=1,2,⋯,
则 {Xn} 服从大数定律。
解
证 因为
n21\Var(i=1∑nXi)=n21i=1∑nσi2≤nc→0(n→∞).
所以由马尔可夫大数定律知 {Xn} 服从大数定律。
(泊松大数定律)设 Sn 为 n 次独立试验中事件 A 出现的次数,而事件 A 在第 i 次试验中出现的概率为 pi, i=1,2,⋯,n,⋯,则对任意的 ε>0,有
n→∞limP(nSn−n1i=1∑npi<ε)=1.
解
证 记
Xi={1,0,第 i 次试验时 A 出现,反之.
则
Sn=i=1∑nXi,
且
E(Sn)=i=1∑npi,\Var(Sn)=i=1∑npi(1−pi)≤4n.
所以由切比雪夫不等式,对任意的 ε>0,有
P(nSn−n1i=1∑npi≥ε)≤ε2\Var(Sn/n)≤4nε21→0,
即
n→∞limP(nSn−n1i=1∑npi<ε)=1.
(伯恩斯坦(Bernstein)大数定律)设 {Xn} 是方差一致有界的随机变量序列,且当 ∣k−l∣→∞ 时,一致地有 \Cov(Xk,Xl)→0,证明 {Xn} 服从大数定律。
解
证 记 \Var(Xn)=σn2≤σ2, \Corr(Xk,Xl)=ρ(k,l)。对任意 ε>0,存在 M>0,当 ∣k−l∣>M 时,有
∣ρ(k,l)∣<ε,
所以
n21\Var(i=1∑nXi)=n21i=1∑nσi2+k,l=1∑nσkσlρ(k,l)≤n21nσ2+σ2k=1∑n∣k−l∣≤M∑∣ρ(k,l)∣+∣k−l∣>M∑∣ρ(k,l)∣≤n21[nσ2+σ2k=1∑n(2M+1+2nε)]=nσ2+n2M+1σ2+2σ2ε.
由 ε 的任意性知
n21\Var(i=1∑nXi)→0(n→∞).
所以由马尔可夫大数定律知 {Xn} 服从大数定律。
(格涅坚科(Gnedenko)大数定律)设 {Xn} 是随机变量序列,若记
Yn=n1i=1∑nXi,an=n1i=1∑nE(Xi),
则 {Xn} 服从大数定律的充要条件是
n→∞limE[1+(Yn−an)2(Yn−an)2]=0.
解
证 先证充分性。对任意 ε>0,注意到 t>0 时,f(t)=t2/(1+t2) 是增函数,故当 ∣y−an∣≥ε 时,有
1+∣y−an∣2∣y−an∣2≥1+ε2ε2,ε21+ε2⋅1+∣y−an∣2∣y−an∣2≥1.
因此有
P(∣Yn−an∣≥ε)=∫∣y−an∣≥εdFYn(y)≤ε21+ε2∫∣y−an∣≥ε1+(y−an)2(y−an)2dFYn(y)≤ε21+ε2E[1+(Yn−an)2(Yn−an)2].
所以当
n→∞limE[1+(Yn−an)2(Yn−an)2]=0
时,有
n→∞limP(∣Yn−an∣≥ε)=0,
故 {Xn} 服从大数定律。
再证必要性。设 {Xn} 服从大数定律,即
n→∞limP(∣Yn−an∣≥ε)=0,
则对任意 ε>0,存在 N,当 n>N 时,有
P(∣Yn−an∣≥ε)<ε.
因为函数 f(t)=t2/(1+t2) 是增函数及 0<f(t)<1,得
0≤E[1+(Yn−an)2(Yn−an)2]≤1+ε2ε2P(∣Yn−an∣<ε)+1×P(∣Yn−an∣≥ε)≤ε2+ε.
由于 ε 的任意性,所以
n→∞limE[1+(Yn−an)2(Yn−an)2]=0.
设 {Xn} 为独立同分布的随机变量序列,方差存在。又设 ∑n=1∞an 为绝对收敛级数。令
Yn=i=1∑nXi.
证明 {anYn} 服从大数定律。
解
证 不妨设 E(Xn)=0,否则令 Xn′=Xn−E(Xn),并讨论 {Xn′} 即可。由 E(Xn)=0 知
\Var(Xn)=E(Xn2)=σ2.
又因为 ∑n=1∞an 为绝对收敛级数,可记
c=n=1∑∞∣an∣<∞.
因为
i=1∑naiYi=i=1∑nai(j=1∑iXj)=j=1∑nXj(i=j∑nai),
故有
n21\Var(i=1∑naiYi)=n21E[j=1∑nXj(i=j∑nai)]2=n2σ2j=1∑n(i=j∑nai)2≤nc2σ2→0(n→∞).
所以由马尔可夫大数定律知 {anYn} 服从大数定律。
设 {Xn} 为独立同分布的随机变量序列,方差存在,令
Yn=i=1∑nXi.
又设 {an} 为一列常数,如果存在常数 c>0,使得对一切 n 有
∣nan∣≤c,
证明 {anYn} 服从大数定律。
解
证 与上一题一样,不妨设 E(Xn)=0,则 \Var(Xn)=E(Xn2)=σ2。对任意的 n≥k,有
i=k∑n∣ai∣≤nc(n−k+1),
因而
k=1∑n(i=k∑n∣ai∣)2≤n2c2k=1∑n(n−k+1)2=n2c2⋅6n(n+1)(2n+1).
仿上一题的证明有
n21\Var(i=1∑naiYi)≤n2σ2k=1∑n(i=k∑n∣ai∣)2≤6c2σ2⋅n3(n+1)(2n+1)→0(n→∞).
所以由马尔可夫大数定律知 {anYn} 服从大数定律。
设 {Xn} 为独立同分布的随机变量序列,其方差有限,且 Xn 不恒为常数。如果
Sn=i=1∑nXi,
试证:随机变量序列 {Sn} 不服从大数定律。
解
证 记
Tn=i=1∑nSi=[k=1∑n(n+1−k)Xk],
令 \Var(Xn)=σ2,则
\Var(Tn)=k=1∑n(n+1−k)2σ2=σ26n(n+1)(2n+1)=σ262n3+3n2+n.
由此得
3σ2n3<\Var(Tn)≤σ2n3.
倘若 {Sn} 服从大数定律,则对任意的 ε>0,有
n→∞limP(∣Tn−E(Tn)∣≥nε)=0.
于是,当 n 充分大时,有
P(∣Tn−E(Tn)∣≥nε)<σ2ε2.
记
A={∣Tn−E(Tn)∣<nε},
则
P(A)<σ2ε2,
由此得
\Var(Tn)=E[Tn−E(Tn)]2=P(A)E{[Tn−E(Tn)]2∣A}+P(A)E{[Tn−E(Tn)]2∣A}≤P(A)⋅n2ε2+P(A)⋅σ2n3≤n2ε2+n3ε2≤2n3ε2.
由 ε 的任意性,不妨取 2ε2<σ2/3,则当 n 充分大时,有
\Var(Tn)<3σ2n3,
这与前面推出的
\Var(Tn)>3σ2n3
相矛盾,所以 {Sn} 不服从大数定律。
分别用随机投点法和平均值法计算下列定积分:
J1=∫01e−1ex−1dx,J2=∫−11exdx.
解
**(1)**先计算 J1。
随机投点法:先用计算机产生 (0,1) 上均匀分布的 2n 个随机数(譬如 n=104),构成 n 个数对 (xi,yi), i=1,2,⋯,n。记
f1(x)=e−1ex−1.
以 n1 记满足不等式 yi≤f1(xi) 的数对个数,则
J1≈nn1.
平均值法:先用计算机产生 n 个 (0,1) 上均匀分布的随机数 xi, i=1,2,⋯,n,然后对每个 xi 计算 f1(xi),最后得 J1 的估计值为
J1≈n1i=1∑nf1(xi).
**(2)**对于第二个积分 J2,先将其化成 [0,1] 区间上的积分。令
f2(y)=e−e−11(e−1+2y−e−1),y∈[0,1].
则 0≤f2(y)≤1。此时有
J2=∫−11exdx=S0∫01f2(y)dy+2e−1,
其中
S0=2(e−e−1).
随机投点法:先用计算机产生 (0,1) 上均匀分布的 2n 个随机数(譬如 n=104),构成 n 个数对 (xi,yi), i=1,2,⋯,n。以 n2 记满足不等式 yi≤f2(xi) 的数对个数,则
J2≈S0⋅nn2+2e−1.
平均值法:先用计算机产生 n 个 (0,1) 上均匀分布的随机数 xi, i=1,2,⋯,n,然后对每个 xi 计算 f2(xi),最后得 J2 的估计值为
J2≈S0⋅n1i=1∑nf2(xi)+2e−1.
§4.4 中心极限定理
1. 中心极限定理 研究独立随机变量和的极限分布在什么条件下为正态分布的问题。
2. 林德伯格-莱维中心极限定理 设 {Xn} 是独立同分布的随机变量序列,且 E(Xi)=μ, \Var(Xi)=σ2>0。记
Yn∗=σnX1+X2+⋯+Xn−nμ,
则对任意实数 y,有
n→∞limP(Yn∗≤y)=Φ(y)=2π1∫−∞ye−t2/2dt.
3. 二项分布的正态近似
**(1)**棣莫弗-拉普拉斯中心极限定理 设 n 重伯努利试验中,事件 A 在每次试验中出现的概率为 p (0<p<1),记 Sn 为 n 次试验中事件 A 出现的次数,且记
Yn∗=np(1−p)Sn−np.
则对任意实数 y,有
n→∞limP(Yn∗≤y)=Φ(y)=2π1∫−∞ye−t2/2dt.
**(2)**近似中的修正 因为二项分布是离散分布,而正态分布是连续分布,所以用正态分布作为二项分布的近似计算中,作些修正可以提高精度。若 k1<k2 均为整数,一般先作如下修正后再用正态近似
P(k1≤Sn≤k2)=P(k1−0.5<Sn<k2+0.5)
=Φ(np(1−p)k2+0.5−np)−Φ(np(1−p)k1−0.5−np).
4. 三类近似计算问题 若记 β=Φ(y),则由棣莫弗-拉普拉斯中心极限定理给出的近似式
P(Yn∗≤y)≈Φ(y)=β
可用来解决三类计算问题:
- 已知 n,y 求 β(求概率);
- 已知 n,β 求 y(求分位数);
- 已知 y,β 求 n(求样本量)。
5. 独立不同分布下的中心极限定理 设 {Xn} 为独立随机变量序列,且 E(Xi)=μi, \Var(Xi)=σi2, i=1,2,⋯,记
Yn=i=1∑nXi,Bn=\Var(Yn)=σ12+σ22+⋯+σn2.
**(1)**林德伯格条件 若诸 Xi 为连续随机变量,其密度函数为 pi(x),对任意的 τ>0,称
n→∞limτ2Bn21i=1∑n∫∣x−μi∣>τBn(x−μi)2pi(x)dx=0
为林德伯格条件。
**(2)**林德伯格中心极限定理 设独立随机变量序列 {Xn} 满足林德伯格条件,则对任意的 x,有
n→∞limP(Bn1i=1∑n(Xi−μi)≤x)=2π1∫−∞xe−t2/2dt.
**(3)**假如独立随机变量序列 {Xn} 具有同分布和方差有限的条件,则必定满足林德伯格条件,即林德伯格-莱维中心极限定理是林德伯格中心极限定理的特例。
**(4)**李雅普诺夫中心极限定理 设 {Xn} 为独立随机变量序列,若存在 δ>0,满足
n→∞limBn2+δ1i=1∑nE(∣Xi−μi∣2+δ)=0,
则对任意的 x,有
n→∞limP[Bn1i=1∑n(Xi−μi)≤x]=2π1∫−∞xe−t2/2dt.
习题与解答 4.4
某保险公司多年的统计资料表明,在索赔户中被盗索赔户占 20%,以 X 表示在随意抽查的 100 个索赔户中因被盗向保险公司索赔的户数。
- 写出 X 的分布列;
- 求被盗索赔户不少于 14 户且不多于 30 户的概率的近似值。
解
(1)X 服从 n=100,p=0.2 的二项分布 b(100,0.2),即
P(X=k)=(k100)0.2k0.8100−k,k=0,1,2,⋯,100.
**(2)**利用棣莫弗-拉普拉斯中心极限定理,有
P(14≤X≤30)=P(13.5<X<30.5)≈Φ(100×0.2×0.830.5−100×0.2)−Φ(100×0.2×0.813.5−100×0.2)=Φ(2.625)−Φ(−1.625)=Φ(2.625)−1+Φ(1.625)=0.99565−1+0.948=0.9437.
这表明:被盗索赔户在 14 与 30 户之间的概率近似为 0.9437。
某电子计算机主机有 100 个终端,每个终端有 80% 的时间被使用。若各个终端是否被使用是相互独立的,试求至少有 15 个终端空闲的概率。
解
记 X 为 100 个终端中被使用的终端个数,则
X∼b(100,0.8).
利用棣莫弗-拉普拉斯中心极限定理,所求概率为
P(X≤85)=P(X<85.5)≈Φ(100×0.8×0.285.5−100×0.8)=Φ(1.375)=0.9155.
这表明至少有 15 个终端空闲的概率近似为 0.9155。
有一批建筑房屋用的木柱,其中 80% 的长度不小于 3 m,现从这批木柱中随机地取出 100 根,问其中至少有 30 根短于 3 m 的概率是多少?
解
记 X 为 100 根木柱中长度不小于 3 m 的根数,则
X∼b(100,0.8).
利用棣莫弗-拉普拉斯中心极限定理,所求概率为
P(X≤70)=P(X<70.5)≈Φ(100×0.8×0.270.5−100×0.8)=1−Φ(2.375)=0.0088.
掷一颗骰子 100 次,记第 i 次掷出的点数为 Xi, i=1,2,⋯,100,点数之平均为
X=1001i=1∑100Xi,
试求概率 P(3≤X≤4)。
解
由题意可得
E(Xi)=3.5,\Var(Xi)=1235,E(X)=3.5,\Var(X)=2407.
利用林德伯格-莱维中心极限定理,可得
P(3≤X≤4)≈Φ(7/2404−3.5)−Φ(7/2403−3.5)=2Φ(2.9277)−1=0.9966.
这表明:掷 100 次骰子点数之平均在 3 到 4 之间的概率近似为 0.9966,很接近于 1。
连续地掷一颗骰子 80 次,求点数之和超过 300 的概率。
解
记 Xi 为第 i 次投掷时出现的点数,i=1,2,⋯,80,则 P(Xi=k)=1/6, k=1,2,3,4,5,6,且
E(Xi)=27,\Var(Xi)=1235.
由林德伯格-莱维中心极限定理,所求概率为
P(i=1∑80Xi>300)≈1−Φ(80×35/12300−80×7/2)=1−Φ(1.309)=0.096.
有 10 个灯泡,设每个灯泡的寿命服从指数分布,其平均寿命为 60 天。每次用一个灯泡,当使用的灯泡坏了以后立即换上一个新的,求这些灯泡总共可使用 450 天以上的概率。
解
记 Xi 为第 i 个灯泡的寿命(单位:天),i=1,2,⋯,10,则 Xi∼Exp(1/60),且
E(Xi)=60,\Var(Xi)=602.
由林德伯格-莱维中心极限定理,所求概率为
P(i=1∑10Xi>450)≈1−Φ(10×602450−60×10)=Φ(0.79)=0.7852.
设 X1,X2,⋯,X48 为独立同分布的随机变量,共同分布为 U(0,5)。其算术平均为
X=481i=1∑48Xi,
试求概率 P(2≤X≤3)。
解
由均匀分布 U(0,5) 可算得
E(Xi)=2.5,\Var(Xi)=1225,E(X)=2.5,\Var(X)=12×4825.
利用林德伯格-莱维中心极限定理,可得
P(2≤X≤3)≈Φ(5/243−2.5)−Φ(5/242−2.5)=2Φ(2.4)−1=0.9836.
这表明:来自均匀分布 U(0,5) 的 48 个随机数的平均在 2 到 3 之间的概率近似为 0.9836,较接近于 1。
某汽车销售点每天出售的汽车数服从参数为 λ=2 的泊松分布。若一年 365 天都经营汽车销售,且每天出售的汽车数是相互独立的,求一年中售出 700 辆以上汽车的概率。
解
记 Xi 为第 i 天出售的汽车辆数,则
Y=X1+X2+⋯+X365
为一年的总销量。由
E(Xi)=\Var(Xi)=2,
知
E(Y)=\Var(Y)=365×2=730.
利用林德伯格-莱维中心极限定理,可得
P(Y>700)=1−P(Y≤700)≈1−Φ(730700−730)=1−Φ(−1.11)=0.8665.
这表明:该销售点一年售出 700 辆以上汽车的概率近似为 0.8665。
某餐厅每天接待 400 名顾客,设每位顾客的消费额(单位:元)服从 (20,100) 上的均匀分布,且顾客的消费额是相互独立的。试求:
- 该餐厅每天的平均营业额;
- 该餐厅每天的营业额在平均营业额 ±760 元内的概率。
解
记 Xi 为第 i 位顾客的消费额,则
Xi∼U(20,100),
所以
E(Xi)=60,\Var(Xi)=31600.
而该餐厅每天的营业额为
Y=i=1∑400Xi.
**(1)**该餐厅每天的平均营业额为
E(Y)=i=1∑400E(Xi)=400×60=24000(元).
**(2)**利用林德伯格-莱维中心极限定理,可得
P(−760<Y−24000<760)≈2Φ(400×1600/3760)−1
=2Φ(1.645)−1=0.90.
这表明:该餐厅每天营业额在 23240 到 24760 元之间的概率近似为 0.90。
一仪器同时收到 50 个信号,其中第 i 个信号的长度为 Ui, i=1,2,⋯,50。设 Ui 是相互独立的,且都服从 (0,10) 内的均匀分布,试求
P(i=1∑50Ui>300).
解
因为
E(Ui)=5,\Var(Ui)=325.
所以
E(i=1∑50Ui)=250,\Var(i=1∑50Ui)=31250.
利用林德伯格-莱维中心极限定理,可得
P(i=1∑50Ui>300)≈1−Φ(1250/3300−250)=1−Φ(2.45)
=1−0.9929=0.0071.
这表明:50 个信号长度之和超过 300 的概率近似为 0.0071。
计算机在进行加法运算时对每个加数取整数(取最为接近它的整数)。设所有的取整误差是相互独立的,且它们都服从 (−0.5,0.5) 上的均匀分布。
- 若将 1500 个数相加,求误差总和的绝对值超过 15 的概率;
- 最多几个数加在一起可使得误差总和的绝对值小于 10 的概率不小于 90%。
解
记 Xi 为第 i 个加数的取整误差,则
Xi∼U(−0.5,0.5),E(Xi)=0,\Var(Xi)=121.
**(1)**由
E(i=1∑1500Xi)=0,\Var(i=1∑1500Xi)=121500=125,
得所求概率为
P(i=1∑1500Xi>15)≈2−2Φ(12515)=2−2Φ(1.34)=0.1802.
**(2)**由题意可列出概率不等式
P(i=1∑nXi<10)≥0.90,
利用林德伯格-莱维中心极限定理,可改写为
2Φ(n/1210)−1≥0.90,或Φ(n203)≥0.95.
查表得
n203≥1.645,或n≤1.645203.
由此得 n≤443.45。这表明:至多 443 个数相加,才能使它们的误差总和的绝对值小于 10 的概率不小于 90%。
设各零件的质量都是随机变量,它们相互独立,且服从相同的分布,其数学期望为 0.5 kg,标准差为 0.1 kg,问 5000 只零件的总质量超过 2510 kg 的概率是多少?
解
记 Xi 为第 i 只零件的质量,由
E(Xi)=0.5,σ(Xi)=0.1,
得
E(i=1∑5000Xi)=2500,\Var(i=1∑5000Xi)=50.
利用林德伯格-莱维中心极限定理,所求概率为
P(i=1∑5000Xi>2510)≈1−Φ(502510−2500)=1−Φ(1.414)=0.0787.
这表明:5000 只零件的总质量超过 2510 kg 的概率近似为 0.0787。
某种产品由 20 个相同部件连接而成,每个部件的长度是均值为 2 mm、标准差为 0.02 mm 的随机变量。假如这 20 个部件的长度相互独立同分布,且规定产品总长为 (40±0.2) mm 时为合格品,求该产品的不合格品率。
解
记 Xi 为第 i 个部件的长度,则
Y=X1+X2+⋯+X20
为总长度,且
E(Xi)=2,σ(Xi)=0.02.
可用林德伯格-莱维中心极限定理近似得合格品率
P(39.8≤Y≤40.2)≈Φ(0.0080.2)−Φ(−0.0080.2)
=2Φ(2.236)−1=0.9746.
所以不合格品率为 0.0254。
一个保险公司有 10000 个汽车投保人,每个投保人平均索赔 280 元,标准差为 800 元。求总索赔额超过 2700000 元的概率。
解
记 Xi 为第 i 个投保人的索赔额,i=1,2,⋯,10000,则
E(Xi)=280,\Var(Xi)=8002.
由林德伯格-莱维中心极限定理,所求概率为
P(i=1∑10000Xi>2700000)≈1−Φ(10000×80022700000−10000×280)
=1−Φ(−1.25)=Φ(1.25)=0.8944.
有两个班级同时上一门课,甲班有 25 人,乙班有 64 人。该门课程期末考试平均成绩为 78 分,标准差为 14 分。试问:甲班的平均成绩超过 80 分的概率大,还是乙班的平均成绩超过 80 分的概率大?
解
记 Xi 为甲班第 i 个学生的成绩,i=1,2,⋯,25;Yj 为乙班第 j 个学生的成绩,j=1,2,⋯,64。因
E(Xi)=E(Yj)=78,\Var(Xi)=\Var(Yj)=142,
所以由林德伯格-莱维中心极限定理,甲班平均成绩超过 80 分的概率为
P(251i=1∑25Xi>80)≈1−Φ(25×14225×80−25×78)
=1−Φ(75)=0.2374.
同理可计算乙班平均成绩超过 80 分的概率为
P(641j=1∑64Yj>80)≈1−Φ(64×14264×80−64×78)
=1−Φ(78)=0.1261.
所以甲班的平均成绩超过 80 分的概率大。
进行独立重复试验,每次试验中事件 A 发生的概率为 0.25。试问能以 95% 的把握保证 1000 次试验中事件 A 发生的频率与概率差多少?此时 A 发生的次数在什么范围内?
解
记 Yn 为 1000 次试验中事件 A 发生的次数,则
Yn∼b(1000,0.25),
且
E(Yn)=250,\Var(Yn)=187.5.
设事件 A 发生的频率 Yn/1000 与概率 0.25 的差为 k,根据题意,可得不等式
P(1000Yn−0.25≤k)≥0.95,
或
P(∣Yn−250∣≤1000k)≥0.95.
利用棣莫弗-拉普拉斯中心极限定理和修正项可得
P(187.5Yn−250≤187.51000k+0.5)≥0.95,
即
2Φ(187.51000k+0.5)−1≥0.95.
由此得
Φ(187.51000k+0.5)≥0.975.
查表得
187.51000k+0.5≥1.96.
从中解得 k≥0.027,这表明能以 95% 的把握保证在 1000 次试验中事件 A 发生的频率与概率相差不大于 0.027。或者说,在 1000 次试验中事件 A 发生的次数在 (250±27) 次间,即在 223 次到 277 次间。
设某生产线上组装每件产品的时间服从指数分布,平均需要 10 min,且各件产品的组装时间是相互独立的。
- 试求组装 100 件产品需要 15 h 至 20 h 的概率;
- 保证有 95% 的可能性,问 16 h 内最多可以组装多少件产品?
解
记 Xi 为组装第 i 件产品的时间(单位:min),则
Xi∼Exp(λ),E(Xi)=1/λ=10,\Var(Xi)=1/λ2=100.
**(1)**根据题意所求概率如下,再用林德伯格-莱维中心极限定理可得
P(15×60≤i=1∑100Xi≤20×60)≈Φ(100×1001200−100×10)−Φ(100×100900−100×10)
=Φ(2)−Φ(−1)=Φ(2)+Φ(1)−1=0.8185.
**(2)**设 16 h 内最多可以组装 k 件产品。则根据题意可列出概率不等式
P(i=1∑kXi≤16×60)≥0.95,
再用林德伯格-莱维中心极限定理可得
Φ(100k960−10k)≥0.95.
由此查表得
10k960−10k≥1.645,
从中解得 k=81。
某种彩票的奖金额 X 由摇奖决定,其分布列为
X/万元P50.2100.2200.2300.1400.1500.11000.1
若一年中要开出 300 个奖,问需要多少奖金总额,才有 95% 的把握能够发放奖金。
解
记 Xi 为第 i 次摇奖的奖金额,则可得
E(Xi)=29,\Var(Xi)=764.
设奖金总额为 k(万元),根据题意可列出不等式
P(i=1∑300Xi≤k)≥0.95,
再用林德伯格-莱维中心极限定理可得
Φ(300×764k−300×29)≥0.95.
由此查表得
229200k−8700≥1.645,
从中解得
k≥9487.54.
取 k=9488(万元)即可。这表明:该种彩票一年开出 300 个奖需要准备 9488 万元,才有 95% 的把握够发放奖金。
一家有 500 间客房的大旅馆的每间客房装有一台 2 kW(千瓦)的空调机。若开房率为 80%,需多少千瓦的电力才有 99% 的可能性保证有足够的电力供使用空调机?
解
记
Xi={1,0,第 i 间客房开房,第 i 间客房未开房,
则 Xi∼b(1,0.8),由此得
Y=X1+X2+⋯+X500∼b(500,0.8).
设共有 k kW 的电力可供使用,根据题意可列不等式
P(2Y≤k)=P(Y≤k/2)≥0.99,
再用棣莫弗-拉普拉斯中心极限定理和修正项可得
Φ(500×0.8×0.2k/2+0.5−500×0.8)≥0.99,
或
Φ(85k−799)≥0.99,
由此查表得
85k−799≥2.33,
从中解得 k≥840.68,取 k=841 kW 即可。这表明:该旅馆每天需要 841 kW 电力,才能以 99% 的把握保证空调机用电。
设某元件是某电气设备的一个关键部件,当该元件失效后立即换上一个新的元件。假定该元件的平均寿命为 100 小时,标准差为 30 小时,试问:应该有多少备件,才能有 0.95 以上的概率,保证这个系统能连续运行 2000 小时以上?
解
记 Xi 为第 i 个元件的寿命,i=1,2,…,n,则
E(Xi)=100,Var(Xi)=302.
根据题意可列不等式
P(i=1∑nXi>2000)≥0.95,
再由林德伯格-莱维中心极限定理可得
Φ(n×302100n−2000)≥0.95.
由此查表得
n×302100n−2000≥1.645,
从中解得 n≥22.33,所以取 n=23,即应有 23 个此种元件,可有 0.95 以上的概率保证这个系统能连续运行 2000 小时以上。
独立重复地对某物体的长度 a 进行 n 次测量,设各次测量结果 Xi 服从正态分布 N(a,0.22)。记 X 为 n 次测量结果的算术平均值,为保证有 95% 的把握使平均值与实际值 a 的差异小于 0.1,问至少需要测量多少次?
解
因为
E(X)=a,Var(X)=n0.04,
所以根据题意可列不等式
P(∣X−a∣<0.1)≥0.95,
再用林德伯格-莱维中心极限定理可得
2Φ(0.04/n0.1)−1≥0.95,
或
Φ(2n)≥0.975.
由此查表得 n/2≥1.96,从中解得 n≥15.3664,取 n=16 即可以 95% 的把握使平均值与实际值 a 的差异小于 0.1。
某工厂每月生产 10000 台液晶投影机,但它的液晶片车间生产液晶片合格品率为 90%,为了以 99.7% 的可能性保证出厂的液晶投影机都能装上合格的液晶片,试问该液晶片车间每月至少应该生产多少片液晶片?
解
设每月至少应该生产 n 片液晶片,其中合格品数记为 X,则有
X∼b(n,0.90).
下求 n,使概率不等式
P(X≥10000)≥0.997,或P(X<10000)≤0.003
成立。利用二项分布的正态近似,可得
Φ(0.9×0.1×n10000+0.5−0.9n)≤0.003,
查表可得
0.3n10000.5−0.9n≤−2.75,
由此解得 n≥11209,即每月至少应该生产 11209 片液晶片。
某产品的合格品率为 99%,问包装箱中应该装多少个此种产品,才能有 95% 的可能性使每箱中至少有 100 个合格产品。
解
设包装箱中装有 n 个产品,其中合格品数记为 X,则有
X∼b(n,0.99).
下求 n,使
P(X≥100)≥0.95,或P(X<100)≤0.05
成立。利用二项分布的正态近似,可得
Φ(0.99×0.01×n100+0.5−0.99n)≤0.05,
查表可得
0.0995n100.5−0.99n≤−1.645,
由此解得 n>103.19,即每箱装有 104 个产品,能有 95% 的可能性使每箱中至少有 100 个合格产品。
为确定某城市成年男子中吸烟者的比例 p,任意调查 n 个成年男子,记其中的吸烟人数为 m,问 n 至少为多大才能保证 m/n 与 p 的差异小于 0.01 的概率大于 95%。
解
因为
m∼b(n,p),
所以
E(m/n)=p,Var(m/n)=np(1−p).
根据题意有
0.95<P(nm−p<0.01)≈2Φ(p(1−p)0.01n)−1,
由此得
0.975<Φ(p(1−p)0.01n).
查表得
p(1−p)0.01n≥1.96,或n≥1962×p(1−p).
因为 p(1−p)≤1/4,所以当
n≥41962=9604
时,必可满足要求,因此至少抽 9604 个成年男子,可使其吸烟频率 m/n 与实际成年人中吸烟率 p 的误差 ∣m/n−p∣ 小于 0.01 的概率大于 95%。
设 X∼Ga(n,1),试问 n 应该多大,才能满足
P(nX−1>0.1)<0.01.
解
因为
E(X)=n,Var(X)=n,
所以由中心极限定理得
P(nX−1≤0.1)≐2Φ(0.1n)−1≥0.99,
即
Φ(0.1n)≥0.995.
查标准正态分布函数值表得
0.1n≥2.575,
所以得 n≥663.06,取 n=664 即可满足要求。
设 {Xn} 为一独立同分布的随机变量序列,已知 E(Xik)=αk, k=1,2,3,4。试证明:当 n 充分大时,
Yn=n1i=1∑nXi2
近似服从正态分布,并指出此正态分布的参数。
解
因为 {Xn} 为独立同分布的随机变量序列,所以 {Xn2} 也是独立同分布的随机变量序列。根据林德伯格-莱维中心极限定理,
Yn=n1i=1∑nXi2
近似服从正态分布,其参数为
E(Yn)=n1i=1∑nE(Xi2)=α2,
Var(Yn)=n21i=1∑nVar(Xi2)=nα4−α22.
用概率论的方法证明:
n→∞lim(1+n+2!n2+⋯+n!nn)e−n=21.
解
设 {Xn} 为独立同分布的随机变量序列,其共同分布为参数 λ=1 的泊松分布 P(λ)。故
E(Xn)=Var(Xn)=1.
又由泊松分布的可加性知,
Yn=i=1∑nXi
服从参数 λ=n 的泊松分布。由林德伯格-莱维中心极限定理知
P(nYn−n≤0)=P(Yn≤n)=e−nk=0∑nk!nk→2π1∫−∞0e−x2/2dx=21.
评论
支持 Markdown 和 LaTeX 数学公式。