第 6 讲 · 一维随机变量的函数分布

已知 XX 的分布,要求 Y=g(X)Y=g(X) 的分布,关键不是套哪条公式,而是回答两件事:YY 能取什么值;一个 yy 可能由哪些 xx 映过来。

1. 离散型:合并所有同像点

若 XX 取值 xix_i,则 Y=g(X)Y=g(X) 的可能值是 g(xi)g(x_i)。不同 xix_i 可能映到同一个 yy,这些概率必须相加:

P(Y=y)=∑i:g(xi)=yP(X=xi).P(Y=y)=\sum_{i:g(x_i)=y}P(X=x_i).

例:平方会合并正负取值

设 XX 的分布为

xx−2-2−1-10011
P(X=x)P(X=x)0.10.10.20.20.30.30.40.4

令 Y=X2Y=X^2,则

P(Y=0)=0.3,P(Y=0)=0.3, P(Y=1)=0.2+0.4=0.6,P(Y=1)=0.2+0.4=0.6, P(Y=4)=0.1.P(Y=4)=0.1.

漏加负分支是这类题最常见的错误。

2. CDF 法:最通用也最稳

直接从定义出发:

FY(y)=P(Y≤y)=P(g(X)≤y).F_Y(y)=P(Y\le y)=P(g(X)\le y).

把不等式解成关于 XX 的区间,再用 XX 的分布计算。若 FYF_Y 可导,最后求

fY(y)=FY′(y).f_Y(y)=F_Y'(y).

CDF 法尤其适合绝对值、平方、最大值这类非单调变换。

3. 单调变换公式

若 gg 在 XX 的支持集上严格单调,反函数为 x=h(y)x=h(y),则

fY(y)=fX(h(y))∣h′(y)∣,f_Y(y)=f_X(h(y))|h'(y)|,

只在 YY 的值域内成立。

绝对值来自“概率不因坐标伸缩凭空改变”。小区间 dydy 对应原坐标长度约为 ∣h′(y)∣dy|h'(y)|dy。

例:对数正态分布

若 X∼N(μ,σ2)X\sim N(\mu,\sigma^2),Y=eXY=e^X,则 y>0y>0,反函数 x=ln⁡yx=\ln y,导数为 1/y1/y:

fY(y)=1yσ2πexp⁡[−(ln⁡y−μ)22σ2],y>0.f_Y(y)= \frac1{y\sigma\sqrt{2\pi}} \exp\left[-\frac{(\ln y-\mu)^2}{2\sigma^2}\right], \qquad y>0.

这称为对数正态分布。由 LOTUS 或正态矩生成公式还能得到

E(Y)=eμ+σ2/2,E(Y)=e^{\mu+\sigma^2/2}, Var⁡(Y)=e2μ+σ2(eσ2−1).\operatorname{Var}(Y) =e^{2\mu+\sigma^2}(e^{\sigma^2}-1).

4. 非单调变换:所有分支都要加

若方程 g(x)=yg(x)=y 有多个根 x1(y),…,xm(y)x_1(y),\ldots,x_m(y),且各分支可导,则

fY(y)=∑j=1mfX(xj(y))∣dxjdy∣.f_Y(y)= \sum_{j=1}^m f_X(x_j(y)) \left|\frac{dx_j}{dy}\right|.

例:Y=X2Y=X^2

当 y>0y>0 时有 x=yx=\sqrt y 与 x=−yx=-\sqrt y 两条分支:

fY(y)=fX(y)+fX(−y)2y,y>0.f_Y(y)= \frac{f_X(\sqrt y)+f_X(-\sqrt y)}{2\sqrt y}, \qquad y>0.

若 X∼N(0,1)X\sim N(0,1),则

fY(y)=12πye−y/2,y>0,f_Y(y)=\frac1{\sqrt{2\pi y}}e^{-y/2}, \qquad y>0,

即 Y∼χ2(1)Y\sim\chi^2(1)。

例:正弦变换

若 Y=sin⁡XY=\sin X,必须先看 XX 的支持集。若支持集跨越多个单调区间,同一 y∈(−1,1)y\in(-1,1) 可能对应多个 xx;应按每个落在支持集内的反解分支相加。只写 x=arcsin⁡yx=\arcsin y 通常不完整。

5. 线性变换与常见分布结论

若 Y=aX+bY=aX+b 且 a≠0a\ne0,

fY(y)=1∣a∣fX(y−ba).f_Y(y)=\frac1{|a|}f_X\left(\frac{y-b}{a}\right).

特别地,

X∼N(μ,σ2)⟹aX+b∼N(aμ+b,a2σ2).X\sim N(\mu,\sigma^2) \Longrightarrow aX+b\sim N(a\mu+b,a^2\sigma^2).

若 X∼Ga(α,λ)X\sim Ga(\alpha,\lambda) 且 c>0c>0,则在率参数约定下

cX∼Ga(α,λc).cX\sim Ga\left(\alpha,\frac\lambda c\right).

尺度变大 cc 倍,率参数相应除以 cc。

6. 概率积分变换

若 FXF_X 连续且严格单调,则

U=FX(X)∼U(0,1).U=F_X(X)\sim U(0,1).

证明很直接。对 0<u<10<u<1,

P(U≤u)=P(FX(X)≤u)=P(X≤FX−1(u))=u.P(U\le u) =P(F_X(X)\le u) =P(X\le F_X^{-1}(u)) =u.

反过来,若 U∼U(0,1)U\sim U(0,1),令

X=F−1(U),X=F^{-1}(U),

就能生成 CDF 为 FF 的随机变量。这是逆变换采样的理论基础。

一个特别漂亮的例子

若 X∼Exp(2)X\sim Exp(2),则 FX(X)=1−e−2X∼U(0,1)F_X(X)=1-e^{-2X}\sim U(0,1)。而

e−2X=1−FX(X)e^{-2X}=1-F_X(X)

也是 U(0,1)U(0,1),因为均匀变量的补仍均匀。

7. 固定解题流程

  1. 写出 XX 的支持集。
  2. 求 Y=g(X)Y=g(X) 的值域。
  3. 离散型列映射表;连续型判断是否单调。
  4. 单调可用反函数公式,非单调优先 CDF 法或多分支公式。
  5. 写清新密度的支持集。
  6. 检查非负性与积分是否为 11。

若最终密度积分不为 11,通常不是少了常数,而是漏了一条反函数分支或值域写错。

评论