欢迎您访问 最编程 本站为您分享编程语言代码,编程技术文章!
您现在的位置是: 首页

多变量正态分布 (MND)

最编程 2024-04-13 14:01:17
...

多元正态分布(Multivariate normal distribution)

前言

我们通常讨论正态分布都是在一元(univariate)的情况下,相信下面的定义大家都很熟悉了:假设随机变量 X X X服从正态分布,则 X X X具有概率密度函数:
f ( x ) = ( 2 π σ ) − 1 exp ( − ( x − μ ) 2 2 σ 2 ) f(x)=(\sqrt{2\pi}\sigma)^{-1}\text{exp}(-\frac{(x-\mu)^2}{2\sigma^2}) f(x)=(2π σ)1exp(2σ2(xμ)2)
其中 μ \mu μ表示 X X X的均值, σ 2 \sigma^2 σ2表示其方差。

有不少读者应该也看到过下面这个公式:
f ( x 1 , x 2 ) = ( 2 π σ 1 σ 2 1 − ρ 2 ) − 1 exp [ − 1 2 ( 1 − ρ 2 ) ( ( x 1 − μ 1 ) 2 σ 1 2 − 2 ρ ( x 1 − μ 1 ) ( x 2 − μ 2 ) σ 1 σ 2 + ( x 2 − μ 2 ) 2 σ 2 2 ) ] \begin{aligned} f(x_1,x_2)=&(2\pi\sigma_1\sigma_2\sqrt{1-\rho^2} )^{-1}\text{exp}[-\frac{1}{2(1-\rho^2)}(\frac{(x_1-\mu_1)^2}{\sigma_1^2}\\ &-\frac{2\rho(x_1-\mu_1)(x_2-\mu_2)}{\sigma_1\sigma_2}+\frac{(x_2-\mu_2)^2}{\sigma_2^2})] \end{aligned} f(x1,x2)=(2πσ1σ21ρ2 )1exp[2(1ρ2)1(σ12(x1μ1)2σ1σ22ρ(x1μ1)(x2μ2)+σ22(x2μ2)2)]
没错,这正是将正态分布拓展到二维的情况,即:
X = [ X 1 , X 2 ] T X=[X_1,X_2]^T X=[X1,X2]T
其中 X 1 X_1 X1, X 2 X_2 X2分别服从正态分布。

有不少读者应该和我一样,看到这个二维的公式就头痛了,这他娘的一堆是啥玩意儿啊?老实说把上面的公式准确的打出来还花费了我不少功夫,可见公式之复杂,如果再往三元以上,简直不敢想象了。

由于许多本文许多内容我是从wikipedia看的,现学现卖,自己也是似懂非懂,不敢误人子弟,只能把自己确定的一些心得写一写,以作备忘,如果可以,也能给一些同有此问的后来者一些帮助。

多元正态分布

假设 X = ( X 1 , X 2 , ⋯   , X k ) T X=(X_1,X_2,\cdots,X_k)^T X=(X1,X2,,Xk)T是一个 k k k维的列向量,服从多元正态分布,我们可以把它记做:
X ∼ N ( μ , Σ ) X\sim N(\mu,\Sigma) XN(μ,Σ)
其中,
μ = E ( X ) = ( μ 1 , μ 2 , ⋯   , μ k ) Σ i , j = C o v ( X i , X j ) \begin{aligned} &\mu=E(X)=(\mu_1,\mu_2,\cdots,\mu_k)\\ &\Sigma_{i,j}=Cov(X_i,X_j) \end{aligned} μ=E(X)=(μ1,μ2,,μk)Σi,j=Cov(Xi,Xj)
对于多元随机变量,我们最关心的是它的概率函数,当上述协方差矩阵是正定的(positive definite),分布才有概率密度函数,这种情况被称为“非退化的”(non-degenerate)。这里笔者亦不甚解,猜测大概和协方差矩阵 Σ \Sigma Σ是否可逆有关。

如果多元正态分布的概率密度函数存在,它被定义如下:
f ( x 1 , x 2 , ⋯   , x k ) = exp ( − 1 2 ( x − μ ) T Σ − 1 ( x − μ ) ) ( 2 π ) k ∣ Σ ∣ f(x_1,x_2,\cdots,x_k)=\frac{\text{exp}(-\frac{1}{2}(x-\mu)^T\Sigma^{-1}(x-\mu))}{\sqrt{(2\pi)^k|\Sigma|}} f(x1,x2,,xk)=(2π)kΣ exp(21(xμ)TΣ1(x

上一篇: 均匀分布和指数分布的密度和分布函数

下一篇: ELF 文件解析 (II):ELF 文件头详细信息

推荐阅读