乐湖华研题库
学生教师

例 3.4.9

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

题干

二维正态分布 N(μ1,μ2,σ12,σ22,ρ)N(\mu_1,\mu_2,\sigma_1^2,\sigma_2^2,\rho)N(μ1​,μ2​,σ12​,σ22​,ρ) 的相关系数就是 ρ\rhoρ。

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

解析

下面先求 Cov(X,Y)\mathrm{Cov}(X,Y)Cov(X,Y)。

Cov(X,Y)=E[(X−E(X))(Y−E(Y))]\mathrm{Cov}(X,Y) = E[(X-E(X))(Y-E(Y))]Cov(X,Y)=E[(X−E(X))(Y−E(Y))] =12πσ1σ21−ρ2∫−∞∞∫−∞∞(x−μ1)(y−μ2)⋅exp⁡{−12(1−ρ2)[(x−μ1)2σ12−2ρ(x−μ1)(y−μ2)σ1σ2+(y−μ2)2σ22]}dxdy.= \frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}(x-\mu_1)(y-\mu_2) \cdot \exp\left\{-\frac{1}{2(1-\rho^2)}\left[\frac{(x-\mu_1)^2}{\sigma_1^2} - 2\rho\frac{(x-\mu_1)(y-\mu_2)}{\sigma_1\sigma_2} + \frac{(y-\mu_2)^2}{\sigma_2^2}\right]\right\}\mathrm{d}x\mathrm{d}y.=2πσ1​σ2​1−ρ2​1​∫−∞∞​∫−∞∞​(x−μ1​)(y−μ2​)⋅exp{−2(1−ρ2)1​[σ12​(x−μ1​)2​−2ρσ1​σ2​(x−μ1​)(y−μ2​)​+σ22​(y−μ2​)2​]}dxdy.

先将上式中方括号内化成

(x−μ1σ1−ρy−μ2σ2)2+(1−ρ2 y−μ2σ2)2,\left(\frac{x-\mu_1}{\sigma_1} - \rho\frac{y-\mu_2}{\sigma_2}\right)^2 + \left(\sqrt{1-\rho^2}\,\frac{y-\mu_2}{\sigma_2}\right)^2,(σ1​x−μ1​​−ρσ2​y−μ2​​)2+(1−ρ2​σ2​y−μ2​​)2,

再作变量变换

{u=11−ρ2(x−μ1σ1−ρy−μ2σ2),v=y−μ2σ2,\begin{cases} u = \dfrac{1}{\sqrt{1-\rho^2}}\left(\dfrac{x-\mu_1}{\sigma_1} - \rho\dfrac{y-\mu_2}{\sigma_2}\right), \\[6pt] v = \dfrac{y-\mu_2}{\sigma_2}, \end{cases}⎩⎨⎧​u=1−ρ2​1​(σ1​x−μ1​​−ρσ2​y−μ2​​),v=σ2​y−μ2​​,​

则

{x−μ1=σ1(u1−ρ2+ρv),y−μ2=σ2v,\begin{cases} x-\mu_1 = \sigma_1(u\sqrt{1-\rho^2}+\rho v), \\ y-\mu_2 = \sigma_2 v, \end{cases}{x−μ1​=σ1​(u1−ρ2​+ρv),y−μ2​=σ2​v,​ dxdy=∣J∣ dudv=σ1σ21−ρ2 dudv.\mathrm{d}x\mathrm{d}y = |J|\,\mathrm{d}u\mathrm{d}v = \sigma_1\sigma_2\sqrt{1-\rho^2}\,\mathrm{d}u\mathrm{d}v.dxdy=∣J∣dudv=σ1​σ2​1−ρ2​dudv.

由此得

Cov(X,Y)=σ1σ22π∫−∞∞∫−∞∞(uv1−ρ2+ρv2)exp⁡{−12(u2+v2)}dudv.\mathrm{Cov}(X,Y) = \frac{\sigma_1\sigma_2}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}(uv\sqrt{1-\rho^2}+\rho v^2)\exp\left\{-\frac{1}{2}(u^2+v^2)\right\}\mathrm{d}u\mathrm{d}v.Cov(X,Y)=2πσ1​σ2​​∫−∞∞​∫−∞∞​(uv1−ρ2​+ρv2)exp{−21​(u2+v2)}dudv.

上式右端积分可以分为两个积分之和,其中

∫−∞∞∫−∞∞uvexp⁡{−12(u2+v2)}dudv=0,\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}uv\exp\left\{-\frac{1}{2}(u^2+v^2)\right\}\mathrm{d}u\mathrm{d}v = 0,∫−∞∞​∫−∞∞​uvexp{−21​(u2+v2)}dudv=0, ∫−∞∞∫−∞∞v2exp⁡{−12(u2+v2)}dudv=2π.\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}v^2\exp\left\{-\frac{1}{2}(u^2+v^2)\right\}\mathrm{d}u\mathrm{d}v = 2\pi.∫−∞∞​∫−∞∞​v2exp{−21​(u2+v2)}dudv=2π.

从而

Cov(X,Y)=σ1σ22π⋅ρ⋅2π=ρσ1σ2,\mathrm{Cov}(X,Y) = \frac{\sigma_1\sigma_2}{2\pi}\cdot\rho\cdot 2\pi = \rho\sigma_1\sigma_2,Cov(X,Y)=2πσ1​σ2​​⋅ρ⋅2π=ρσ1​σ2​, Corr(X,Y)=Cov(X,Y)σ1σ2=ρ.\mathrm{Corr}(X,Y) = \frac{\mathrm{Cov}(X,Y)}{\sigma_1\sigma_2} = \rho.Corr(X,Y)=σ1​σ2​Cov(X,Y)​=ρ.

评论 0

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