多维高斯分布---【2】

1.一维高斯分布

\qquad 在介绍二维高斯分布之前我们先介绍一下一维高斯分布的函数图像,如下所示:
f ( x ) = 1 2 π ⋅ δ ⋅ e − ( x − μ ) 2 2 δ 2 f(x)= \frac{1}{\sqrt{2 \pi}\cdot \delta}\cdot e^{-\frac{{(x-\mu)}^2}{2\delta^2}} f(x)=2π δ1e2δ2(xμ)2
其中 μ \mu μ代表均值, δ \delta δ代表标准差,其图案如下所示:
在这里插入图片描述

2.二维高斯分布

\qquad 上边我们介绍了一维的高斯分布,接下来我们介绍一下二维高斯分布,首先贴出其函数表达式,为了简单起见,以下的推导我假设所有变量都是相对独立的,且都服从高斯分布,也就是对于概率分布函数 f ( x 0 , x 1 , ⋅ ⋅ ⋅ ⋅ , x n ) f(x_0,x_1,\cdot\cdot\cdot\cdot,x_n) f(x0x1,,xn)而言,存在以下等式:
f ( x 0 , x 1 , ⋅ ⋅ ⋅ ⋅ , x n ) = f ( x 0 ) ⋅ f ( x 1 ) ⋅ ⋅ ⋅ ⋅ f ( x n ) f(x_0,x_1,\cdot\cdot\cdot\cdot,x_n)=f(x_0)\cdot f(x_1)\cdot\cdot\cdot\cdot f(x_n) f(x0x1,,xn)=f(x0)f(x1)f(xn)
式中的 f ( x i ) f(x_i) f(xi)服从一维的高斯分布
f ( x i ) = 1 2 π ⋅ δ i ⋅ e − ( x − μ i ) 2 2 δ i 2 f(x_i)= \frac{1}{\sqrt{2 \pi}\cdot \delta_i}\cdot e^{-\frac{{(x-\mu_i)}^2}{2\delta_i^2}} f(xi)=2π δi1e2δi2(xμi)2
\qquad 其中 δ i \delta_i δi μ i \mu_i μi是第i个变量的标准差和均值。那当我们来描述二维的高斯分布的话,此处我们的 n n n取值为2。因为 x 1 x_1 x1 x 2 x_2 x2是相互独立的,所以二维的高斯分布函数可以表示为:
f ( x 1 , x 2 ) = f ( x 1 ) ⋅ f ( x 2 ) = 1 2 π ⋅ δ 1 ⋅ e − ( x − μ 1 ) 2 2 δ 1 2 ⋅ 1 2 π ⋅ δ 2 ⋅ e − ( x − μ 2 ) 2 2 δ 2 2 = 1 2 π ⋅ δ 1 δ 2 ⋅ e − [ δ 2 2 ( x 1 − μ 1 ) 2 + δ 1 2 ( x 2 − μ 2 ) 2 ] 2 δ 1 2 δ 2 2 f(x_1,x_2)=f(x_1)\cdot f(x_2)\\ \qquad\qquad\qquad\quad\quad\qquad\qquad\qquad=\frac{1}{\sqrt{2 \pi}\cdot \delta_1}\cdot e^{-\frac{{(x-\mu_1)}^2}{2\delta_1^2}}\cdot \frac{1}{\sqrt{2 \pi}\cdot \delta_2}\cdot e^{-\frac{{(x-\mu_2)}^2}{2\delta_2^2}}\\ \quad\\\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad=\frac{1}{{2 \pi}\cdot \delta_1 \delta_2}\cdot e^{-\frac{[\delta_2^2{(x_1-\mu_1)}^2+\delta_1^2{(x_2-\mu_2)}^2]}{2\delta_1^2\delta_2^2}} f(x1,x2)=f(x1)f(x2)=2π δ11e2δ12(xμ1)22π δ21e2δ22(xμ2)2=2πδ1δ21e2δ12δ22[δ22(x1μ1)2+δ12(x2μ2)2]
\qquad 其二维的高斯分布的图像如下所示:
在这里插入图片描述

3.多维高斯分布

\qquad 前边我们介绍了一维的和二维的高斯分布,并且画出了其图像,想必大家也在其它的资料中看到过,通过一维和二维扩展出来的多维的高斯分布的表达公式,就像这样:
N ( X ⃗ ∣ μ ⃗ , Σ ) = 1 ( 2 π ) D 2 ⋅ ∣ Σ ∣ 1 2 ⋅ e − ( X ⃗ − μ ⃗ ) T ⋅ Σ − 1 ⋅ ( X ⃗ − μ ⃗ ) 2 N(\vec{X}\mid\vec{\mu},{\Sigma})=\frac{1}{{(2\pi)}^{\frac{D}{2}}\cdot {\mid \Sigma\mid}^{\frac{1}{2}}}\cdot e^{-\frac{({\vec{X}-\vec{\mu})}^T\cdot{\Sigma^{-1}}\cdot{({\vec{X}-\vec{\mu})}}}{2}} N(X μ ,Σ)=(2π)2DΣ211e2(X μ )TΣ1(X μ )
\qquad 式中各个参数代表的意思如下:

  1. X ⃗ \vec{X} X 表示维度为D的向量,
  2. μ ⃗ \vec{\mu} μ 是由多个变量 x 1 , x 2 , x 3 ⋅ ⋅ ⋅ x n x_1,x_2,x_3\cdot\cdot\cdot x_n x1,x2,x3xn各自的均值 u i u_i ui组成的向量,
  3. Σ \Sigma Σ代表所有向量的协方差矩阵,是一个n维n列的矩阵,
  4. Σ − 1 \Sigma^{-1} Σ1代表协方差矩阵的逆,也是一个n维n列的矩阵。

\qquad 讲到这里感觉对于这个多维的公式中的理解还是不是特别透彻,那么我们可以这么想,既然超过二维的原理上来讲已经算是多维的高斯分布了,那么我们能不能通过第三节的公式推导出第二节的二维高斯分布的函数呢?答案当然是可以:
\qquad 因为要证明的是二维的情况,所以我们根据上边列出的各个向量分别得到:
X ⃗ = [ x 1 x 2 ]    μ ⃗ = [ μ 1 μ 2 ] \vec{X}=\left[ \begin{matrix} x_1 \\ x_2 \end{matrix} \right]\ \ \qquad\vec{\mu}=\left[ \begin{matrix} \mu_1 \\ \mu_2 \end{matrix} \right] X =[x1x2]  μ =[μ1μ2]
\qquad 式中 x 1 x_1 x1 x 2 x_2 x2均是长度为n的序列,因为是二维的情况下,所以我们不妨先去求一下协方差矩阵 ∑ \sum ,根据已学知识,我们知道其表达式如下:
Σ = [ δ 11 δ 12 δ 21 δ 22 ] = [ D X 1 C o v ( X 1 , X 2 ) C o v ( X 2 , X 1 ) D X 2 ] = [ δ 1 2 δ 12 δ 21 δ 2 2 ] \Sigma=\left[ \begin{matrix} \delta_{11} & \delta_{12}\\ \delta_{21} & \delta_{22} \end{matrix} \right] =\left[ \begin{matrix} DX_1 & Cov(X_1,X_2)\\ Cov(X_2,X_1)& DX_2 \end{matrix} \right] =\left[ \begin{matrix} \delta_{1} ^2& \delta_{12}\\ \delta_{21} & \delta_{2}^2 \end{matrix} \right] Σ=[δ11δ21δ12δ22]=[DX1Cov(X2,X1)Cov(X1,X2)DX2]=[δ12δ21δ12δ22]
\qquad 又因为我们假设 x 1 x_1 x1 x 2 x_2 x2是相互独立的,也就是有 δ 12 = δ 21 = 0 \delta_{12}=\delta_{21}=0 δ12=δ21=0,于是协方差矩阵即可简化为;
Σ = [ δ 1 2 0 0 δ 2 2 ] \Sigma=\left[ \begin{matrix} \delta_{1} ^2& 0\\ 0& \delta_{2}^2 \end{matrix} \right] Σ=[δ1200δ22]
\qquad 据此我们即可算得协方差的行列式 ∣ Σ ∣ = δ 1 2 δ 2 2 \mid \Sigma\mid={\delta_1}^2{\delta_2}^2 Σ=δ12δ22,又因为我们选取得是二维高斯分布,所以D取值为2,于是我们将其带入多维高斯分布得左半部分可得;
1 ( 2 π ) D 2 ⋅ ∣ Σ ∣ 1 2 = 1 ( 2 π ) ⋅ ∣ δ 1 2 δ 2 2 ∣ 1 2 = 1 2 π ⋅ δ 1 δ 2 \frac{1}{{(2\pi)}^{\frac{D}{2}}\cdot {\mid \Sigma\mid}^{\frac{1}{2}}}=\frac{1}{{(2\pi)}\cdot {\mid{\delta_1}^2{\delta_2}^2\mid}^{\frac{1}{2}}}=\frac{1}{{2 \pi}\cdot \delta_1 \delta_2} (2π)2DΣ211=(2π)δ12δ22211=2πδ1δ21
\qquad 可见左半部分是和二维分布得函数吻合得的,接下来就进行验证多维分布的右半部分是否与二维的分布相同。
\qquad 通过上边的部分,我们很轻易即可得到协方差矩阵的逆( Σ − 1 \Sigma^{-1} Σ1),其值如下:
Σ − 1 = 1 δ 1 2 δ 2 2 ⋅ [ δ 2 2 0 0 δ 1 2 ] {\Sigma}^{-1}=\frac{1}{\delta_1^2\delta_2^2}\cdot\left[ \begin{matrix} \delta_{2} ^2& 0\\ 0& \delta_{1}^2 \end{matrix} \right] Σ1=δ12δ221[δ2200δ12]
\qquad 于是我们将我们已知的变量带入即可得到:
e − ( X ⃗ − μ ⃗ ) T ⋅ Σ − 1 ⋅ ( X ⃗ − μ ⃗ ) 2 = e − 1 2 ⋅ [ x 1 − μ 1 x 2 − μ 2 ] ⋅ 1 δ 1 2 δ 2 2 ⋅ [ δ 2 2 0 0 δ 1 2 ] ⋅ [ x 1 − μ 1 x 2 − μ 2 ] = e − 1 2 δ 1 2 δ 2 2 ⋅ [ x 1 − μ 1 x 2 − μ 2 ] ⋅ [ δ 2 2 0 0 δ 1 2 ] ⋅ [ x 1 − μ 1 x 2 − μ 2 ] = e − 1 2 δ 1 2 δ 2 2 ⋅ [ δ 2 2 ⋅ ( x 1 − μ 1 ) δ 1 2 ⋅ ( x 2 − μ 2 ) ] ⋅ [ x 1 − μ 1 x 2 − μ 2 ] = e − 1 2 δ 1 2 δ 2 2 ⋅ [ δ 2 2 ⋅ ( x 1 − μ 1 ) 2 + δ 1 2 ⋅ ( x 2 − μ 2 ) 2 ] e^{-\frac{({\vec{X}-\vec{\mu})}^T\cdot{\Sigma^{-1}}\cdot{({\vec{X}-\vec{\mu})}}}{2}}=e^{-\frac{1}{2}\cdot \left[ \begin{matrix} x_1-\mu_1& x_2-\mu_2\end{matrix} \right]\cdot\frac{1}{\delta_1^2\delta_2^2}\cdot\left[ \begin{matrix} \delta_{2} ^2& 0\\ 0& \delta_{1}^2 \end{matrix} \right]\cdot\left[ \begin{matrix} x_1-\mu_1\\ x_2-\mu_2\end{matrix} \right]}\\ \qquad \qquad\\\qquad\qquad\qquad\quad=e^{-\frac{1}{2\delta_1^2\delta_2^2}\cdot \left[ \begin{matrix} x_1-\mu_1& x_2-\mu_2\end{matrix} \right]\cdot\left[ \begin{matrix} \delta_{2} ^2& 0\\ 0& \delta_{1}^2 \end{matrix} \right]\cdot\left[ \begin{matrix} x_1-\mu_1\\ x_2-\mu_2\end{matrix} \right]}\\ \qquad \qquad\\\qquad\qquad\qquad\quad=e^{-\frac{1}{2\delta_1^2\delta_2^2}\cdot \left[ \begin{matrix} \delta_2^2\cdot(x_1-\mu_1)& \delta_1^2\cdot(x_2-\mu_2)\end{matrix} \right]\cdot\left[ \begin{matrix} x_1-\mu_1\\ x_2-\mu_2\end{matrix} \right]}\\ \quad \\= e^{-\frac{1}{2\delta_1^2\delta_2^2}\cdot[ \delta_2^2\cdot(x_1-\mu_1)^2+\delta_1^2\cdot(x_2-\mu_2)^2]} e2(X μ )TΣ1(X μ )=e21[x1μ1x2μ2]δ12δ221[δ2200δ12][x1μ1x2μ2]=e2δ12δ221[x1μ1x2μ2][δ2200δ12][x1μ1x2μ2]=e2δ12δ221[δ22(x1μ1)δ12(x2μ2)][x1μ1x2μ2]=e2δ12δ221[δ22(x1μ1)2+δ12(x2μ2)2]
\qquad 好了通过上边公式的推导,我们可以看到通过多维情况下我们可以完整的得出二维情况下的高斯分布表达式。

4.心声

\qquad 原理虽简单,但还是在纸上详细推导一遍比较好,不要眼高手低~~
加油呀,各位~。

  • 7
    点赞
  • 25
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

王延凯的博客

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值