乐湖华研题库
学生教师

例 3.3.8

medium一级题目发布者: ai-batch

题干

设随机变量 X∼Ga(α1,λ)X \sim Ga(\alpha_1,\lambda)X∼Ga(α1​,λ),Y∼Ga(α2,λ)Y \sim Ga(\alpha_2,\lambda)Y∼Ga(α2​,λ),且 XXX 与 YYY 独立,证明 Z=X+Y∼Ga(α1+α2,λ)Z=X+Y \sim Ga(\alpha_1+\alpha_2,\lambda)Z=X+Y∼Ga(α1​+α2​,λ)。

答案点击展开后可查看解析

解析

首先指出 Z=X+YZ=X+YZ=X+Y 仍在 (0,∞)(0,\infty)(0,∞) 上取值,所以当 z≤0z \leq 0z≤0 时,pZ(z)=0p_Z(z)=0pZ​(z)=0。而当 z>0z>0z>0 时,可用卷积公式 (3.3.13),此时被积函数 pX(z−y)pY(y)p_X(z-y)p_Y(y)pX​(z−y)pY​(y) 的非零区域为 0<y<z0<y<z0<y<z,故

pZ(z)=λα1+α2Γ(α1)Γ(α2)∫0z(z−y)α1−1e−λ(z−y)⋅yα2−1e−λydyp_Z(z)=\frac{\lambda^{\alpha_1+\alpha_2}}{\Gamma(\alpha_1)\Gamma(\alpha_2)} \int_0^z (z-y)^{\alpha_1-1}e^{-\lambda(z-y)}\cdot y^{\alpha_2-1}e^{-\lambda y}\mathrm{d}ypZ​(z)=Γ(α1​)Γ(α2​)λα1​+α2​​∫0z​(z−y)α1​−1e−λ(z−y)⋅yα2​−1e−λydy =λα1+α2e−λzΓ(α1)Γ(α2)∫0z(z−y)α1−1yα2−1dy=\frac{\lambda^{\alpha_1+\alpha_2}e^{-\lambda z}}{\Gamma(\alpha_1)\Gamma(\alpha_2)} \int_0^z (z-y)^{\alpha_1-1}y^{\alpha_2-1}\mathrm{d}y=Γ(α1​)Γ(α2​)λα1​+α2​e−λz​∫0z​(z−y)α1​−1yα2​−1dy =λα1+α2e−λzΓ(α1)Γ(α2)zα1+α2−1∫01(1−t)α1−1tα2−1dt,=\frac{\lambda^{\alpha_1+\alpha_2}e^{-\lambda z}}{\Gamma(\alpha_1)\Gamma(\alpha_2)} z^{\alpha_1+\alpha_2-1} \int_0^1(1-t)^{\alpha_1-1}t^{\alpha_2-1}\mathrm{d}t,=Γ(α1​)Γ(α2​)λα1​+α2​e−λz​zα1​+α2​−1∫01​(1−t)α1​−1tα2​−1dt,

最后的积分是贝塔函数,它等于 Γ(α1)Γ(α2)/Γ(α1+α2)\Gamma(\alpha_1)\Gamma(\alpha_2)/\Gamma(\alpha_1+\alpha_2)Γ(α1​)Γ(α2​)/Γ(α1​+α2​)。代入上式得

pZ(z)=λα1+α2Γ(α1+α2)zα1+α2−1e−λz.p_Z(z)=\frac{\lambda^{\alpha_1+\alpha_2}}{\Gamma(\alpha_1+\alpha_2)} z^{\alpha_1+\alpha_2-1}e^{-\lambda z}.pZ​(z)=Γ(α1​+α2​)λα1​+α2​​zα1​+α2​−1e−λz.

这正是形状参数为 α1+α2\alpha_1+\alpha_2α1​+α2​,尺度参数仍为 λ\lambdaλ 的伽马分布。

这个结论表明:两个尺度参数相同的独立的伽马变量之和仍为伽马变量,其尺度参数不变,而形状参数相加,即

Ga(α1,λ)∗Ga(α2,λ)=Ga(α1+α2,λ).(3.3.16)Ga(\alpha_1,\lambda) * Ga(\alpha_2,\lambda) = Ga(\alpha_1+\alpha_2,\lambda). \tag{3.3.16}Ga(α1​,λ)∗Ga(α2​,λ)=Ga(α1​+α2​,λ).(3.3.16)

显然这个结论可推广到有限个尺度参数相同的独立伽马变量之和上。

另外由第二章中我们知道,伽马分布有两个常用的特例:指数分布和卡方分布,即

Exp(λ)=Ga(1,λ),χ2(n)=Ga(n2,12),Exp(\lambda)=Ga(1,\lambda), \quad \chi^2(n)=Ga\left(\frac{n}{2},\frac{1}{2}\right),Exp(λ)=Ga(1,λ),χ2(n)=Ga(2n​,21​),

由此又可以得到另外两个结论:

(1)mmm 个独立同分布的指数变量之和为伽马变量,即

Exp(λ)∗Exp(λ)∗⋯∗Exp(λ)⏟m个=Ga(m,λ).(3.3.17)\underbrace{Exp(\lambda) * Exp(\lambda) * \cdots * Exp(\lambda)}_{m个} =Ga(m,\lambda). \tag{3.3.17}m个Exp(λ)∗Exp(λ)∗⋯∗Exp(λ)​​=Ga(m,λ).(3.3.17)

(2)mmm 个独立的 χ2\chi^2χ2 变量之和为 χ2\chi^2χ2 变量(χ2\chi^2χ2 分布的可加性),即

χ2(n1)∗χ2(n2)∗⋯∗χ2(nm)=χ2(n1+n2+⋯+nm).(3.3.18)\chi^2(n_1) * \chi^2(n_2) * \cdots * \chi^2(n_m) =\chi^2(n_1+n_2+\cdots+n_m). \tag{3.3.18}χ2(n1​)∗χ2(n2​)∗⋯∗χ2(nm​)=χ2(n1​+n2​+⋯+nm​).(3.3.18)

评论 0

还没有评论,来说第一句话吧