设 X 是一个随机变量,其分布已知,g(x) 是一个连续函数,如何求 Y=g(X) 的分布?
这是概率论中的基本问题之一,主要有以下几种方法。
分布函数法是求随机变量函数分布的最基本方法,步骤如下:
- 写出 Y 的分布函数 FY(y)=P(Y≤y)=P(g(X)≤y)
- 将不等式 g(X)≤y 转化为关于 X 的不等式
- 利用 X 的分布求出 FY(y)
- 对 FY(y) 求导得到 fY(y)(连续型)
若 g(x) 严格单调且可导,则 Y=g(X) 的密度函数为:
fY(y)={fX(g−1(y))⋅∣(g−1(y))′∣,0,y∈g 的值域其他
推导:设 g 严格单调递增,则
FY(y)=P(g(X)≤y)=P(X≤g−1(y))=FX(g−1(y))
fY(y)=FY′(y)=fX(g−1(y))⋅(g−1(y))′
若 g 严格单调递减,则
FY(y)=P(g(X)≤y)=P(X≥g−1(y))=1−FX(g−1(y))
fY(y)=fX(g−1(y))⋅∣(g−1(y))′∣
例题:设 X∼N(0,1),求 Y=X2 的分布。
解:
当 y≤0 时,FY(y)=0。
当 y>0 时:
FY(y)=P(X2≤y)=P(−y≤X≤y)=Φ(y)−Φ(−y)=2Φ(y)−1
fY(y)=FY′(y)=2φ(y)⋅2y1=y1⋅2π1e−y/2=2πy1e−y/2
即 Y∼χ2(1)。
当 g(x) 严格单调且可导时,可直接使用公式:
fY(y)=fX(h(y))⋅∣h′(y)∣
其中 h(y)=g−1(y) 是 g 的反函数。
线性变换:Y=aX+b(a=0)
fY(y)=∣a∣1fX(ay−b)
特别地,若 X∼N(μ,σ2),则 Y=aX+b∼N(aμ+b,a2σ2)。
幂变换:Y=X2(X>0)
fY(y)=2y1fX(y),y>0
指数变换:Y=eX
fY(y)=y1fX(lny),y>0
对数变换:Y=lnX(X>0)
fY(y)=eyfX(ey)
当 g(x) 不单调时,公式法不能直接使用,需要用分布函数法分段处理。
若 g(x) 在若干区间上分别单调,则在每个单调区间上分别使用公式法,然后求和。
例题:设 X∼N(0,1),求 Y=∣X∣ 的分布。
解:y=∣x∣ 在 (−∞,0) 上单调递减,在 (0,+∞) 上单调递增。
当 y>0 时:
FY(y)=P(∣X∣≤y)=P(−y≤X≤y)=Φ(y)−Φ(−y)=2Φ(y)−1
fY(y)=2φ(y)=2π2e−y2/2,y>0
对于一般的 g(x):
- 对每个 y,求出 {x:g(x)≤y} 的范围
- 利用 X 的分布计算 P(g(X)≤y)
- 对 FY(y) 求导
设 X 为离散型随机变量,分布律为 P(X=xk)=pk,Y=g(X)。
- 列出 Y 的所有可能取值 y1,y2,⋯
- 对每个 yj,计算 P(Y=yj)=∑g(xk)=yjpk
例题:设 X 的分布律为:
求 Y=X2 的分布律。
解:Y 的可能取值为 0,1,4。
P(Y=0)=P(X=0)=0.3
P(Y=1)=P(X=−1)+P(X=1)=0.2+0.1=0.3
P(Y=4)=P(X=2)=0.4
| 变换 | X 的分布 | Y 的分布 |
|---|
| Y=aX+b | N(μ,σ2) | N(aμ+b,a2σ2) |
| Y=X2 | N(0,1) | χ2(1) |
| Y=σX−μ | N(μ,σ2) | N(0,1) |
| Y=eX | N(μ,σ2) | 对数正态分布 |
| 变换 | X 的分布 | Y 的分布 |
|---|
| Y=aX+b | U(0,1) | U(b,a+b) |
| Y=−lnX | U(0,1) | Exp(1) |
| Y=F−1(X) | U(0,1) | F |
| 变换 | X 的分布 | Y 的分布 |
|---|
| Y=λX | Exp(λ) | Exp(1) |
| Y=eX | Exp(λ) | 帕累托分布 |