当前位置:网站首页>浙大概率第四版,证明正态分布随机变量不相关等价于独立
浙大概率第四版,证明正态分布随机变量不相关等价于独立
2022-07-21 17:34:00 【蓝域小兵】
证明命题:
∀ 满 足 正 态 分 布 的 随 机 变 量 X , Y , 且 其 联 合 分 布 也 满 足 正 态 分 布 , 有 X , Y 不 相 关 ⇔ X , Y 独 立 \forall 满足正态分布的随机变量X,Y,且其联合分布也满足正态分布,有X,Y不相关\Leftrightarrow X,Y独立 ∀满足正态分布的随机变量X,Y,且其联合分布也满足正态分布,有X,Y不相关⇔X,Y独立
如果 X , Y X,Y X,Y独立,则必有 X , Y X,Y X,Y不相关
因为 X , Y X,Y X,Y独立,那么有 E ( X Y ) = E ( X ) E ( Y ) E(XY)=E(X)E(Y) E(XY)=E(X)E(Y),则其协方差 C o v ( X , Y ) = E ( X Y ) − E ( X ) E ( Y ) Cov(X,Y)=E(XY)-E(X)E(Y) Cov(X,Y)=E(XY)−E(X)E(Y)等于 0 0 0,相关系数 ρ = C o v ( X , Y ) D ( X ) D ( Y ) \rho = \frac{Cov(X,Y)}{\sqrt{D(X)}\sqrt{D(Y)}} ρ=D(X)D(Y)Cov(X,Y),所以其相关系数为 0 0 0,所以不相关。接下来就要证明正态分布的 X , Y X,Y X,Y不相关,则其独立
由多维正态分布随机变量的联合分布概率密度函数公式:
f ( N → ) = 1 ( 2 π ) n 2 ∣ Σ ∣ 1 2 e − 1 2 ( N → − μ ) T Σ − 1 ( N → − μ ) f\left( \overrightarrow {N}\right) =\dfrac {1}{\left( 2\pi \right) ^{\dfrac {n}{2}}\left| \Sigma \right| ^{\dfrac {1}{2}}}e^{-\dfrac {1}{2}}\left( \overrightarrow {N}-\mu \right) ^{T} \Sigma ^{-1}\left( \overrightarrow {N}-\mu \right) f(N)=(2π)2n∣Σ∣211e−21(N−μ)TΣ−1(N−μ)
其中
N → 是 多 个 随 机 变 量 组 成 的 随 机 矢 量 , n 行 1 列 Σ 是 随 机 矢 量 的 协 方 差 矩 阵 μ 是 N → 的 期 望 矢 量 , n 行 1 列 \begin{aligned} &\overrightarrow {N} 是多个随机变量组成的随机矢量,n行1列 \\ &\Sigma是随机矢量的协方差矩阵 \\ &\mu 是\overrightarrow {N}的期望矢量,n行1列 \end{aligned} N是多个随机变量组成的随机矢量,n行1列Σ是随机矢量的协方差矩阵μ是N的期望矢量,n行1列
在二维随机变量中,协方差矩阵是一个二维矩阵,不难计算它的逆,
协方差矩阵的公式:
Σ = E [ ( N − E [ N ] ) ( N − E [ N ] ) T ] N = [ X , Y ] T \Sigma=\mathrm{E} \left[ \left( \mathbf{N} - \mathrm{E}[\mathbf{N}] \right) \left( \mathbf{N} - \mathrm{E}[\mathbf{N}] \right)^{\rm T} \right] \\ N=[X,Y]^T Σ=E[(N−E[N])(N−E[N])T]N=[X,Y]T
可以按照公式
A − 1 = A ∗ ∣ A ∣ A^{-1} = \frac{A^*}{|A|} A−1=∣A∣A∗,其中 A ∗ A^* A∗是A的伴随矩阵,伴随矩阵是余因子矩阵的转置矩阵,余因子矩阵的每个元素是原矩阵 A A A对应元素的代数余子式,某元素的代数余子式即去掉该元素所在行列之后得到的矩阵的行列式。这一段话涉及了很多线性代数的概念。可以看着下面的公式去一步步推导。
A = [ a b c d ] A − 1 = 1 a d − b c [ d − b − c a ] \begin{aligned} A&=\begin{bmatrix} a & b \\ c & d \end{bmatrix} \\ A^{-1}&=\dfrac {1}{ad-bc}\begin{bmatrix} d & -b \\ -c & a \end{bmatrix} \end{aligned} AA−1=[acbd]=ad−bc1[d−c−ba]
所以有,二维随机变量(也即2行1列的随机矢量)的概率密度函数 f ( x , y ) f(x,y) f(x,y):
f ( x , y ) = 1 2 π σ 1 σ 2 1 − ρ 2 e q 其 中 : q = − 1 2 ( 1 − ρ 2 ) [ ( x − μ 1 ) 2 σ 1 2 − 2 ρ ( x − μ 1 ) ( y − μ 2 ) σ 1 σ 2 + ( y − μ 2 ) 2 σ 2 2 ] ρ = C o v ( X , Y ) σ 1 σ 2 σ 1 , σ 2 分 别 是 X , Y 的 标 准 差 , μ 1 , μ 2 分 别 是 X , Y 的 期 望 \begin{aligned} f\left( x,y\right) &=\dfrac {1}{2\pi \sigma _{1}\sigma _{2}\sqrt {1-\rho ^{2}}}e^{q} \\ 其中:\\ q&=\dfrac {-1}{2\left( 1-\rho ^{2}\right) }\left[ \dfrac {\left( x-\mu _{1}\right) ^{2}}{\sigma_1^{2}}-2\rho \dfrac {\left( x-\mu _{1}\right) \left( y-\mu _{2}\right) }{\sigma _{1}\sigma _{2}}+\dfrac {\left( y-\mu _2\right) ^{2}}{\sigma ^2_{2}}\right] \\ \rho &=\dfrac {Cov\left( X,Y\right) }{\sigma _{1}\sigma _{2}} \\ &\sigma_1,\sigma_2分别是X,Y的标准差,\mu_1,\mu_2分别是X,Y的期望 \end{aligned} f(x,y)其中:qρ=2πσ1σ21−ρ21eq=2(1−ρ2)−1[σ12(x−μ1)2−2ρσ1σ2(x−μ1)(y−μ2)+σ22(y−μ2)2]=σ1σ2Cov(X,Y)σ1,σ2分别是X,Y的标准差,μ1,μ2分别是X,Y的期望
而:
f ( x ) = 1 2 π σ 1 e − ( x − μ 1 ) 2 2 σ 1 2 f ( y ) = 1 2 π σ 2 e − ( y − μ 2 ) 2 2 σ 2 2 f ( x ) f ( y ) = 1 2 π σ 1 σ 2 e ( x − μ 1 ) 2 σ 1 2 + ( y − μ 2 ) 2 σ 2 2 \begin{aligned} f\left( x\right) &= \Large{\dfrac {1}{\sqrt {2\pi }\sigma_1} e}^{\small{-\dfrac {\left( x-\mu_1 \right) ^{2}}{2\sigma_1^{2}}}}\\ \\ f\left( y\right) &= \Large{\dfrac {1}{\sqrt {2\pi }\sigma_2 }e}^{\small{-\dfrac {\left( y-\mu_2 \right) ^{2}}{2\sigma_2^{2}}}}\\ \\ f(x)f(y) &= \Large{\dfrac {1}{2\pi \sigma _{1}\sigma _{2}}e}^{\small{ \dfrac {\left( x-\mu _{1}\right) ^{2}}{\sigma_1^{2}}+\dfrac {\left( y-\mu _2\right) ^{2}}{\sigma ^2_{2}}}} \\ \end{aligned}\\ f(x)f(y)f(x)f(y)=2πσ11e−2σ12(x−μ1)2=2πσ21e−2σ22(y−μ2)2=2πσ1σ21eσ12(x−μ1)2+σ22(y−μ2)2
可以看出,有 ρ = 0 ⇒ f ( x , y ) = f ( x ) f ( y ) \rho=0 \Rightarrow f(x,y)=f(x)f(y) ρ=0⇒f(x,y)=f(x)f(y),即不相关可以推导出独立。
综上1,2,命题左边可以推出右边,右边也可以推出左边,所以命题成立。
边栏推荐
- Error: a connection was successfully established with the server, but
- How to build a smart park with 5g and mobile edge computing servers
- js 对象深拷贝
- Token
- 三星6818LED驱动的编写
- Two methods of redis installation and configuration in Windows Environment
- JS时间和时间戳的转换
- sed,流编辑器以及vim编辑器。
- 自定义图片裁剪框
- 鼠标移入select options会触发mouseleave 事件处理方案
猜你喜欢
Cannot read property ‘type‘ of undefined Occurred while linting **\index.jsx:1
vscode 配置代码自动格式化加修复
Module build failed: Error: Plugin/Preset files are not allowed to export objects, only functions.
Zhang Chi Consulting: the fifth Six Sigma project of a packaging and printing group passed the expert review
OpenGL贴图保持原比例的一种方法
海思Hi3531||瑞芯微RK1109用rtsp客户端实现h264拉流
指针总结篇
Hisilicon [hi3531]onvif protocol is implemented based on gSOAP and OpenSSL_ Code article
La Fondation n'est pas assez solide pour secouer la montagne
OSPF 重发布
随机推荐
更新版本号和静态资源(web内嵌页面的存放地址)
La Fondation n'est pas assez solide pour secouer la montagne
RunTime+compiler和RunTime only版本的区别
scrollIntoView
Kernel timer interrupt based on [hisilic hi3516dv300]
ts 学习记录(一)sudo忘记密码(乌龙)Try changing the ‘lib’ compiler option to include ‘dom’.
How to build a smart park with 5g and mobile edge computing servers
一篇让理解你Mysql存储引擎
C#异常处理 | 连接成功但登陆异常/管道另一端上无任何进程Error: A connection was successfully established with the server, but
Two methods of redis installation and configuration in Windows Environment
OSPF 重发布
Must use destructuring props assignmenteslint
解决IDEA的插件中心连接不上网络
Formation and destruction of function stack frames (26 pictures help you understand function stack frames in depth)
File operation and related functions
去除tinymce中粘贴的样式
基于[三星6818]芯片超声波测距驱动编写
STM32-f40x FSMC信号配置并点亮LCD屏
RK3288关于LVDS信号配置和1080p视频信号的详解
鼠标移入select options会触发mouseleave 事件处理方案