一、联合密度函数的推导
1.1 协方差矩阵的基本性质
给定:
X=[X1X2]∼N2(μ,Σ)
其中:
μΣ=[μ1μ2],=[σ11σ21σ12σ22]=[σ12ρσ1σ2ρσ1σ2σ22]>0注记:Σ>0 表示 Σ 是正定矩阵。
1.2 计算协方差矩阵的行列式
∣Σ∣=det(Σ)=σ12ρσ1σ2ρσ1σ2σ22=σ12⋅σ22−(ρσ1σ2)2=σ12σ22(1−ρ2)
由于 Σ>0(正定),故 ∣Σ∣>0,即:
σ12σ22(1−ρ2)>0⟹1−ρ2>0⟹∣ρ∣<1
1.3 计算协方差矩阵的逆矩阵
对于 2×2 矩阵 A=[acbd],其逆矩阵为:
A−1=ad−bc1[d−c−ba]
伴随矩阵求逆
应用到 Σ 上:
Σ−1=∣Σ∣1[σ22−ρσ1σ2−ρσ1σ2σ12]=σ12σ22(1−ρ2)1[σ22−ρσ1σ2−ρσ1σ2σ12]
化简得:
Σ−1=1−ρ21[σ121−σ1σ2ρ−σ1σ2ρσ221]
1.4 二次型的展开
令 x=[x1x2],则:
x−μ=[x1−μ1x2−μ2]
计算二次型 (x−μ)′Σ−1(x−μ):
(x−μ)′Σ−1(x−μ)=[x1−μ1x2−μ2]⋅1−ρ21[σ121−σ1σ2ρ−σ1σ2ρσ221]⋅[x1−μ1x2−μ2]=1−ρ21[x1−μ1x2−μ2]⋅[σ12x1−μ1−σ1σ2ρ(x2−μ2)−σ1σ2ρ(x1−μ1)+σ22x2−μ2]=1−ρ21[(x1−μ1)(σ12x1−μ1−σ1σ2ρ(x2−μ2))+(x2−μ2)(−σ1σ2ρ(x1−μ1)+σ22x2−μ2)]=1−ρ21[σ12(x1−μ1)2−σ1σ2ρ(x1−μ1)(x2−μ2)−σ1σ2ρ(x1−μ1)(x2−μ2)+σ22(x2−μ2)2]=1−ρ21[σ12(x1−μ1)2−2ρσ1σ2(x1−μ1)(x2−μ2)+σ22(x2−μ2)2]1.5 联合密度函数
根据多元正态分布的密度函数公式:
f(x)=(2π)p/2∣Σ∣1/21exp[−21(x−μ)′Σ−1(x−μ)]
对于 p=2,有:
(2π)p/2=(2π)2/2=2π
∣Σ∣1/2=σ12σ22(1−ρ2)=σ1σ21−ρ2
因此,联合密度函数为:
f(x1,x2)=2πσ1σ21−ρ21exp[−2(1−ρ2)1(σ12(x1−μ1)2−2ρσ1σ2(x1−μ1)(x2−μ2)+σ22(x2−μ2)2)]
二、边缘密度函数的推导
2.1 X1 的边缘密度函数
边缘密度函数 fX1(x1) 通过对 x2 积分得到:
fX1(x1)=∫−∞+∞f(x1,x2)dx2
将联合密度函数代入:
fX1(x1)=∫−∞+∞2πσ1σ21−ρ21exp[−2(1−ρ2)1(σ12(x1−μ1)2−2ρσ1σ2(x1−μ1)(x2−μ2)+σ22(x2−μ2)2)]dx2
2.2 完成平方
将指数中的表达式关于 x2 完成平方:
令 A=2(1−ρ2)1,则指数部分为:
−σ22A(x2−μ2)2+σ1σ22Aρ(x1−μ1)(x2−μ2)−σ12A(x1−μ1)2=−σ22A[(x2−μ2)2−2ρσ1σ2(x1−μ1)(x2−μ2)]−σ12A(x1−μ1)2=−σ22A[(x2−μ2)2−2ρσ1σ2(x1−μ1)(x2−μ2)+ρ2σ12σ22(x1−μ1)2−ρ2σ12σ22(x1−μ1)2]−σ12A(x1−μ1)2=−σ22A[(x2−μ2)−ρσ1σ2(x1−μ1)]2+σ12Aρ2(x1−μ1)2−σ12A(x1−μ1)2=−σ22A[(x2−μ2)−ρσ1σ2(x1−μ1)]2−σ12A(1−ρ2)(x1−μ1)22.3 积分计算
代入 A=2(1−ρ2)1:
fX1(x1)=2πσ1σ21−ρ21exp[−21⋅σ12(x1−μ1)2]×∫−∞+∞exp[−2(1−ρ2)σ221((x2−μ2)−ρσ1σ2(x1−μ1))2]dx2令 y=(x2−μ2)−ρσ1σ2(x1−μ1),则 dy=dx2,积分变为:
∫−∞+∞exp[−2(1−ρ2)σ22y2]dy
这是一个高斯积分,其值为:
∫−∞+∞e−y2/(2σ2)dy=σ2π
这里 σ2=(1−ρ2)σ22,所以:
∫−∞+∞exp[−2(1−ρ2)σ22y2]dy=σ22π(1−ρ2)
2.4 最终结果
代入得:
fX1(x1)=2πσ1σ21−ρ21exp[−21⋅σ12(x1−μ1)2]⋅σ22π(1−ρ2)=2πσ11exp[−21⋅σ12(x1−μ1)2]因此,X1 的边缘密度函数为:
fX1(x1)=2πσ11exp[−2σ12(x1−μ1)2]
即 X1∼N(μ1,σ12)。
2.5 X2 的边缘密度函数
同理,X2 的边缘密度函数为:
fX2(x2)=2πσ21exp[−2σ22(x2−μ2)2]
即 X2∼N(μ2,σ22)。
四、总结
4.1 联合密度函数
f(x1,x2)=2πσ1σ21−ρ21exp[−2(1−ρ2)1(σ12(x1−μ1)2−2ρσ1σ2(x1−μ1)(x2−μ2)+σ22(x2−μ2)2)]
4.2 边缘密度函数
fX1(x1)=2πσ11exp[−2σ12(x1−μ1)2],X1∼N(μ1,σ12)
fX2(x2)=2πσ21exp[−2σ22(x2−μ2)2],X2∼N(μ2,σ22)
4.3 ρ 的统计意义
- 定义:ρ=σ1σ2Cov(X1,X2)
- 取值范围:∣ρ∣<1(非退化情况)
- 线性相关程度:衡量 X1 和 X2 之间的线性相关程度
- 条件分布:
- E[X2∣X1=x1]=μ2+ρσ1σ2(x1−μ1)
- Var(X2∣X1=x1)=σ22(1−ρ2)
- 决定系数:ρ2 表示解释比例
- 几何意义:决定联合密度函数等高线的形状和方向