注意,本文内容来自于吴恩达老师cs229课堂笔记的中文翻译项目:https://github.com/Kivy-CN/Stanford-CS-229-CN 中部分的内容进行翻译学习。
更多关于多元高斯分布
介绍
到目前为止的课堂上,多元高斯分布已经出现在许多应用中,比如线性回归的概率解释、高斯判别分析、高斯混合聚类,以及最近学习的因子分析。在本节的笔记中,我们试图揭开多元高斯函数在最近学习的因子分析课程中引入的一些奇特的性质。本节笔记的目的是让大家对这些性质的来源有一些直观的了解,这样你就可以在作业(提醒你写作业的线索!)中更加明确地使用这些性质。
1. 定义
我们称一个概率密度函数是一个均值为 μ ∈ R n \mu\in R^n μ∈Rn,协方差矩阵为 Σ ∈ S + + n \Sigma\in S_{++}^n Σ∈S++n的 1 ^1 1一个多元正态分布(或高斯分布)(multivariate normal (or Gaussian) distribution), 其随机变量是向量值 x ∈ R n x\in R^n x∈Rn,该概率密度函数可以通过下式表达:
上一小段上标1的说明(详情请点击本行)
1 复习一下线性代数章节中介绍的 S + + n S_{++}^n S++n是一个对称正定的 n × n n\times n n×n矩阵空间,定义为:
S + + n = { A ∈ R n × n : A = A T a n d x T A x > 0 f o r a l l x ∈ R n s u c h t h a t x ≠ 0 } S_{++}^n=\{A\in R^{n\times n}:A=A^T\quad and\quad x^TAx>0\quad for\quad all\quad x\in R^n\quad such\quad that\quad x\neq 0\} S++n={A∈Rn×n:A=ATandxTAx>0forallx∈Rnsuchthatx̸=0}
p ( x ; μ , Σ ) = 1 ( 2 π ) n / 2 ∣ Σ ∣ 1 / 2 exp ( − 1 2 ( x − μ ) T Σ − 1 ( x − μ ) ) p(x;\mu,\Sigma)=\frac{1}{(2\pi)^{n/2}|\Sigma|^{1/2}} \exp\left(-\frac{1}{2}(x-\mu)^T\Sigma^{-1}(x-\mu)\right) p(x;μ,Σ)=(2π)n/2∣Σ∣1/21exp(−21(x−μ)TΣ−1(x−μ))
我们可以写作 x ∼ N ( μ , Σ ) x\sim\mathcal{N}(\mu,\Sigma) x∼N(μ,Σ)。
2. 高斯分布的特点
多元高斯在实践中非常方便,因为其如下的特点:
-
特点 #1: 如果你知道以 x x x为随机变量的高斯分布的均值 μ \mu μ和协方差矩阵 Σ \Sigma Σ。则你可以直接写出关于 x x x的概率密度函数。
-
特点 #2: 下列高斯积分具有闭式解(closed-form solutions):
∫ x ∈ R n p ( x ; μ , Σ ) d x = ∫ − ∞ ∞ ​ ⋯ ∫ − ∞ ∞ p ( x ; μ , Σ ) d x 1 … d x 2 = 1 ∫ x ∈ R n x i p ( x ; μ , σ ) d x = μ i ∫ x ∈ R n ( x i − μ i ) ( x j − μ j ) p ( x ; μ , σ ) d x = Σ i j \begin{aligned} \int_{x\in R^n}p(x;\mu,\Sigma)dx &= \int_{-\infin}^{\infin}\dots\int_{-\infin}^{\infin}p(x;\mu,\Sigma)dx_1\dots dx_2=1 \\ \int_{x\in R^n}x_ip(x;\mu,\sigma)dx &= \mu_i \\ \int_{x\in R^n}(x_i-\mu_i)(x_j-\mu_j)p(x;\mu,\sigma)dx &=\Sigma_{ij} \end{aligned} ∫x∈Rnp(x;μ,Σ)dx∫x∈Rnxip(x;μ,σ)dx∫x∈Rn(xi−μi)(xj−μj)p(x;μ,σ)dx=∫−∞∞⋯∫−∞∞p(x;μ,Σ)dx1…dx2=1=μi=Σij
- 特点 #3: 高斯函数遵循一些封闭性质(closure properties:):
- 独立高斯随机变量的和是高斯分布。
- 联合高斯分布的边缘分布是高斯分布。
- 联合高斯分布的条件是高斯分布。
乍一看,这些事实中的一些结论,尤其是第 1 1 1和第 2 2 2条,似乎要么是直观上显而易见的,要么至少是可信的。然而,我们可能不太清楚的是为什么这些特点如此有用。在本文档中,我们将提供一些直观解释说明如何在平常操作处理多元高斯随机变量时使用这些特点。
3. 封闭性质
在本节中,我们将详细讨论前面描述的每个封闭属性,我们将使用特点#1和#2来证明属性,或者至少给出一些关于属性正确性的直觉。
下面是我们本节将要介绍的内容的路线图:
独立高斯分布的和 | 联合高斯分布的边缘分布 | 联合高斯分布的条件分布 | |
---|---|---|---|
为什么是高斯函数的解释 | 不介绍 | 介绍 | 介绍 |
概率密度函数的结果 | 介绍 | 介绍 | 介绍 |
3.1 独立高斯分布的和是高斯分布
本规则的正式表述为:
设有 y ∼ N ( μ , Σ ) y\sim\mathcal{N}(\mu,\Sigma) y∼N(μ,Σ)和 z ∼ N ( μ ′ , Σ ′ ) z\sim\mathcal{N}(\mu',\Sigma') z∼N(μ′,Σ′)为独立高斯分布,其中随机变量 μ , μ ′ ∈ R n \mu,\mu'\in R^n μ,μ′∈Rn且 Σ , Σ ′ ∈ S + + n \Sigma,\Sigma'\in S_{++}^n Σ,Σ′∈S++n。则它们的和也同样是高斯分布:
y + z ∼ N ( μ + μ ′ , Σ + Σ ′ ) y+z\sim\mathcal{N}(\mu+\mu',\Sigma+\Sigma') y+z∼N(μ+μ′,Σ+Σ′)
在我们证明上面的结论前,先给出一些直观结果:
- 首先要指出的是上述规则中独立假设的重要性。为了了解为什么这很重要,假设 y ∼ N ( μ , σ ) y\sim\mathcal{N}(\mu,\sigma) y∼N(μ,σ)是服从于均值 μ \mu μ方差 σ \sigma σ的多元高斯分布,并且假设 z = − y z=-y z=−y。很明显, z z z也是服从于与多元高斯分布(事实上, z ∼ N ( − μ , σ ) z\sim\mathcal{N}(-\mu,\sigma) z∼N(−μ,σ)),但是 y + z y+z y+z等于零(不是高斯分布)!
- 第二件需要指出的事情是许多学生感到困惑的一点:如果我们把两个高斯概率密度函数(多维空间中的“肿块(bumps)”)加在一起,我们会得到一些峰(即“双峰(two-humped)”的概率密度函数)么?在这里,我们要注意到随机变量 y + z y + z y+z的概率密度函数并不是简单的将两个单独的概率密度函数的随机变量 y y y和 z z z相加,而是会变成 y y y和 z z z的卷积的概率密度函数。 2 ^2 2 然而证明“两个高斯概率密度函数的卷积得到一个高斯概率密度函数”超出了这门课的范围。
上一小段上标2的说明(详情请点击本行)
2 例如,如果 y y y和 z z z是单变量高斯函数(即: y ∼ N ( μ , σ 2 ) , z ∼ N ( μ , σ ′ 2 ) y\sim\mathcal{N}(\mu,\sigma^2),z\sim\mathcal{N}(\mu,\sigma'^2) y∼N(μ,σ2),z∼N(μ,σ′2)),则它们的概率密度的卷积由下式给出:
p ( y + z ; μ , μ ′ , σ , σ ′ 2 ) = ∫ − ∞ ∞ p ( w ; μ , σ 2 ) p ( y + z − w ; μ ′ , σ ′ 2 ) d w = ∫ − ∞ ∞ 1 2 π σ exp ( − 1 2 σ 2 ( w − μ ) 2 ) ⋅ 1 2 π σ ′ exp ( − 1 2 σ ′ 2 ( y + z − w − μ ′ ) 2 ) d w \begin{aligned} p(y+z;\mu,\mu',\sigma,\sigma'^2) &=\int_{-\infin}^{\infin}p(w;\mu,\sigma^2)p(y+z-w;\mu',\sigma'^2)dw \\ &= \int_{-\infin}^{\infin}\frac 1{\sqrt{2\pi}\sigma}\exp\left(-\frac 1{2\sigma^2}(w-\mu)^2\right)\cdot \frac 1{\sqrt{2\pi}\sigma'}\exp\left(-\frac 1{2\sigma'^2}(y+z-w-\mu')^2\right)dw \end{aligned} p(y+z;μ,μ′,σ,σ′2)=∫−∞∞p(w;μ,σ2)p(y+z−w;μ′,σ′2)dw=∫−∞∞2πσ1exp(−2σ21(w−μ)2)⋅2πσ′1exp(−2σ′21(y+z−w−μ′)2)dw
转换一下思路,让我们用卷积给出高斯概率密度函数的观察结果,加上特点#1,来算出概率密度函数 p ( y + z ∣ μ , Σ ) p(y+z|\mu,\Sigma) p(y+z∣μ,Σ)的解析解。如果我们要计算卷积。我们该怎么做呢?回顾特点#1,高斯分布完全由它的均值向量和协方差矩阵指定。如果我们能确定这些值是什么,那么我们就能计算出其解析解了。
这很简单!对应期望而言,我们有:
E [ y i + z i ] = E [ y i ] + E [ z i ] = μ i + μ i ′ E[y_i+z_i]=E[y_i]+E[z_i]=\mu_i+\mu_i' E[yi+zi]=E[yi]+E[zi]=μi+μi′
上式的结果根据期望的线性性质。因此, y + z y + z y+z的均值可以简单的写作 μ + μ ′ \mu+\mu' μ+μ′。 同时,协方差矩阵的第 ( i , j ) (i, j) (i,j)项由下式给出:
E [ ( y i + z i ) ( y j + z j ) ] − E [ y i + z i ] E [ y j + z j ] = E [ y i y j + z i y j + y i z j + z i z j ] − ( E [ y i ] + E [ z i ] ) ( E [ y j ] + E [ z j ] ) = E [ y i y j ] + E [ z i y j ] + E [ y i z j ] + E [ z i z j ] − E [ y i ] E [ y j ] − E [ z i ] E [ y j ] − E [ y i ] E [ z j ] − E [ z i ] E [ z j ] = ( E [ y i y j ] − E [ y i ] E [ y j ] ) + ( E [ z i z j ] − E [ z i ] E [ z j ] ) + ( E [ z i y j ] − E [ z i ] E [ y j ] ) + ( E [ y i z j ] − E [ y i ] E [ z j ] ) \begin{aligned} &E[(y_i+z_i)(y_j+z_j)]-E[y_i+z_i]E[y_j+z_j] \\ &\qquad=E[y_iy_j+z_iy_j+y_iz_j+z_iz_j]-(E[y_i]+E[z_i])(E[y_j]+E[z_j]) \\ &\qquad=E[y_iy_j]+E[z_iy_j]+E[y_iz_j]+E[z_iz_j]-E[y_i]E[y_j]-E[z_i]E[y_j]-E[y_i]E[z_j]-E[z_i]E[z_j] \\ &\qquad=(E[y_iy_j]-E[y_i]E[y_j])+(E[z_iz_j]-E[z_i]E[z_j]) \\ &\qquad\qquad+(E[z_iy_j]-E[z_i]E[y_j])+(E[y_iz_j]-E[y_i]E[z_j]) \\ \end{aligned} E[(yi+zi)(yj+zj)]−E[yi+zi]E[yj+zj]=E[yiyj+ziyj+yizj+zizj]−(E[yi]+E[zi])(E[yj]+E[zj])=E[yiyj]+E[ziyj]+E[yizj]+E[zizj]−E[yi]E[yj]−E[zi]E[yj]−E[yi]E[zj]−E[zi]E[zj]=(E[yiyj]−E[yi]E[yj])+(E[zizj]−E[zi]E[zj])+(E[ziyj]−E[zi]E[yj])+(E[yizj]−E[yi]E[zj])
利用 y y y和 z z z相互独立的事实,我们得到 E [ z i y j ] = E [ z i ] E [ y j ] E[z_iy_j]=E[z_i]E[y_j] E[ziyj]=E[zi]E[yj]和 E [ y i z j ] = E [ y i ] E [ z j ] E[y_iz_j]=E[y_i]E[z_j] E[yizj]=E[yi]E[zj]。因此,最后两项消去了,剩下:
E [ ( y i + z i ) ( y j + z j ) ] − E [ y i + z i ] E [ y j + z j ] = ( E [ y i y j ] − E [ y i ] E [ y j ] ) + ( E [ z i z j ] − E [ z i ] E [ z j ] ) = Σ i j + Σ i j ′ \begin{aligned} &E[(y_i+z_i)(y_j+z_j)]-E[y_i+z_i]E[y_j+z_j] \\ &\qquad=(E[y_iy_j]-E[y_i]E[y_j])+(E[z_iz_j]-E[z_i]E[z_j]) \\ &\qquad=\Sigma_{ij}+\Sigma_{ij}' \end{aligned} E[(yi+zi)(yj+zj)]−E[yi+zi]E[yj+zj]=(E[yiyj]−E[yi]E[yj])+(E[zizj]−E[zi]E[zj])=Σij+Σij′
由此,我们可以得出 y + z y + z y+z的协方差矩阵可以简单的写作 Σ + Σ ′ \Sigma+\Sigma' Σ+Σ′。
此刻,让我们回顾一下刚刚我们做了什么?利用一些简单的期望和独立性的性质,我们计算出了 y + z y + z y+z的均值和协方差矩阵。根据特点#1,我们可以立即写出 y + z y + z y+z的概率密度函数,而不需要做卷积! 3 ^3 3
3 当然,我们首先需要知道 y + z y + z y+z是高斯分布。
3.2 联合高斯分布的边缘分布是高斯分布
本规则的正式表述为:
假设
[ x A x B ] ∼ N ( [ μ A μ B ] , [ Σ A A Σ A B Σ B A Σ B B ] ) \begin{bmatrix}x_A\\x_B\end{bmatrix}\sim\mathcal{N}\begin{pmatrix}\begin{bmatrix}\mu_A\\\mu_B\end{bmatrix},\begin{bmatrix}\Sigma_{AA}&\Sigma_{AB}\\\Sigma_{BA}&\Sigma_{BB}\end{bmatrix}\end{pmatrix} [xAxB]∼N([μAμB],[ΣAAΣBAΣABΣBB])
其中 x A ∈ R m , x B ∈ R n x_A\in R^m,x_B\in R^n xA∈Rm,xB∈Rn并选择均值向量和协方差矩阵子块的维数与 x A x_A xA和 x B x_B xB进行匹配。则边缘概率密度函数如下所示:
p ( x A ) = ∫ x B ∈ R n p ( x A , x B ; μ , Σ ) d x B p ( x B ) = ∫ x A ∈ R m p ( x A , x B ; μ , Σ ) d x A p(x_A)=\int_{x_B\in R^n}p(x_A,x_B;\mu,\Sigma)dx_B \\ p(x_B)=\int_{x_A\in R^m}p(x_A,x_B;\mu,\Sigma)dx_A p(xA)=∫xB∈Rnp(xA,xB;μ,Σ)dxBp(xB)=∫xA∈Rmp(xA,xB;μ,Σ)dxA
上面式子都是高斯分布:
x A ∼ N ( μ A , Σ A A ) x B ∼ N ( μ B , Σ B B ) x_A\sim\mathcal{N}(\mu_A,\Sigma_{AA}) \\ x_B\sim\mathcal{N}(\mu_B,\Sigma_{BB}) xA∼N(μA,ΣAA)xB∼N(μB,ΣBB)
为了证明这个规则,我们只关注变量 x A x_A xA的边缘分布。
4 一般来说,对于一个高斯分布的随机向量 x x x,只要我们对均值向量的项和协方差矩阵的行/列按对应的方式进行置换,则总是可以对 x x x的项进行置换。因此,只看 x A x_A xA就足够了, x B x_B xB的结果也立即得到了。
首先,请注意计算边缘分布的均值和协方差矩阵很简单:只需从联合概率密度函数的均值和协方差矩阵中提取相应的子块。为了确保这是绝对清楚的,我们来看看 x A , i x_{A,i} xA,i和 x A , j x_{A,j} xA,j( x A x_A xA的第 i i i个部分和 x A x_A xA的第 j j j个部分)之间的协方差。注意 x A , i x_{A,i} xA,i和 x A , j x_{A,j} xA,j同样也是下面式子的第 i i i个和第 j j j个部分:
[ x A x B ] \begin{bmatrix}x_A\\x_B\end{bmatrix} [xAxB]
(因为 x A x_A xA出现在这个向量的上部分)。要找到它们的协方差,我们只需简单的使用下面式子的那个协方差矩阵的第 ( i , j ) (i, j) (i,j)个元素即可:
[ Σ A A Σ A B Σ B A Σ B B ] \begin{bmatrix}\Sigma_{AA}&\Sigma_{AB}\\\Sigma_{BA}&\Sigma_{BB}\end{bmatrix} [ΣAAΣBAΣABΣBB]
第 ( i , j ) (i, j) (i,j)个元素在可以在 Σ A A \Sigma_{AA} ΣAA子块矩阵中找到。事实上就是 Σ A A , i j \Sigma_{AA,ij} ΣAA,ij。对所有的 i , j ∈ { 1 , … , m } i,j\in \{1,\dots,m\} i,j∈{1,…,m}使用这个参数,我们可以发现 x A x_A xA的协方差矩阵可以简化为 Σ A A \Sigma_{AA} ΣAA。类似的方法可以用来求 x A x_A xA的均值简化为 μ A \mu_A μA。因此,上面的论证告诉我们,如果我们知道 x A x_A xA的边缘分布是高斯分布,那么我们就可以用合适的均值子矩阵以及联合概率密度函数的协方差矩阵立即写出 x A x_A xA的概率密度函数。
上面的论证虽然简单,但多少有些不令人满意:我们如何才能真正确定 x A x_A xA是一个多元高斯分布?关于这一点的论述有点冗长,因此,与其节外生枝,不如先列出我们的推导过程:
- 明确写出边缘概率密度函数的积分形式。
- 通过对逆协方差矩阵进行分块来重写积分。
- 使用“平方和”参数来计算 x B x_B xB上的积分。
- 论述得到的概率密度函数是高斯的。
下面让我们分别研究一下上面提到的每一个步骤。
3.2.1 边缘概率密度函数的积分形式
假设我们想直接计算 x A x_A xA的密度函数。然后,我们需要计算积分:
p ( x A ) = ∫ x B ∈ R n p ( x A , x B ; μ , Σ ) d x B = 1 ( 2 π ) m + n 2 ∣ Σ A A Σ A B Σ B A Σ B B ∣ 1 / 2 ∫ x B ∈ R n exp ( − 1 2 [ x A − μ A x B − μ B ] T [ Σ A A Σ A B Σ B A Σ B B ] − 1 [ x A − μ A x B − μ B ] ) d x B \begin{aligned} p(x_A) &= \int_{x_B\in R^n}p(x_A,x_B;\mu,\Sigma)dx_B \\ &= \frac{1}{(2\pi)^{\frac{m+n}{2}} \begin{vmatrix}\Sigma_{AA}&\Sigma_{AB}\\\Sigma_{BA}&\Sigma_{BB}\end{vmatrix}^{1/2}}\int_{x_B\in R^n}\exp\left(-\frac12\begin{bmatrix}x_A-\mu_A\\x_B-\mu_B\end{bmatrix}^T\begin{bmatrix}\Sigma_{AA}&\Sigma_{AB}\\\Sigma_{BA}&\Sigma_{BB}\end{bmatrix}^{-1}\begin{bmatrix}x_A-\mu_A\\x_B-\mu_B\end{bmatrix}\right)dx_B \end{aligned} p(xA)=∫xB∈Rnp(xA,xB;μ,Σ)dxB=(2π)2m+n∣∣∣∣ΣAAΣBAΣABΣBB∣∣∣∣1/21∫xB∈Rnexp(−21[xA−μAxB−μB]T[ΣAAΣBAΣABΣBB]−1[xA−μAxB−μB])dxB
3.2.2 逆协方差矩阵的分块
为了进一步推导,我们需要把指数中的矩阵乘积写成稍微不同的形式。特别地,让我们定义下面这个矩阵:
V = [ V A A V A B V B A V B B ] = Σ − 1 V=\begin{bmatrix}V_{AA}&V_{AB}\\V_{BA}&V_{BB}\end{bmatrix}=\Sigma^{-1} V=[VAAVBAVABVBB]=Σ−1
这里我们可能会有下面这种诱人的推导想法:
V = [ V A A V A B V B A V B B ] = [ Σ A A Σ A B Σ B A Σ B B ] − 1 “ = ” [ Σ A A − 1 Σ A B − 1 Σ B A − 1 Σ B B − 1 ] V=\begin{bmatrix}V_{AA}&V_{AB}\\V_{BA}&V_{BB}\end{bmatrix}=\begin{bmatrix}\Sigma_{AA}&\Sigma_{AB}\\\Sigma_{BA}&\Sigma_{BB}\end{bmatrix}^{-1}“=”\begin{bmatrix}\Sigma_{AA}^{-1}&\Sigma_{AB}^{-1}\\\Sigma_{BA}^{-1}&\Sigma_{BB}^{-1}\end{bmatrix} V=[VAAVBAVABVBB]=[ΣAAΣBAΣABΣBB]−1“=”[ΣAA−1ΣBA−1ΣAB−1ΣBB−1]
然而,最右边的等号并不成立!我们将在稍后的步骤中讨论这个问题;不过,现在只要将 V V V定义为上述形式就足够了,而不必担心每个子矩阵的实际内容是什么。
利用 V V V的这个定义,积分扩展到下面的式子:
p ( x A ) = 1 Z ∫ x B ∈ R n exp ( − [ 1 2 ( x A − μ A ) T V A A ( x A − μ A ) + 1 2 ( x A − μ A ) T V A B ( x B − μ B ) + 1 2 ( x B − μ B ) T V B A ( x A − μ A ) + 1 2 ( x B − μ B ) T V B B ( x B − μ B ) ] ) d x B \begin{aligned} p(x_A)=\frac 1Z\int_{x_B\in R^n}\exp(-&[\frac 12(x_A-\mu_A)^TV_{AA}(x_A-\mu_A)+\frac 12(x_A-\mu_A)^TV_{AB}(x_B-\mu_B) \\ & +\frac 12(x_B-\mu_B)^TV_{BA}(x_A-\mu_A)+\frac 12(x_B-\mu_B)^TV_{BB}(x_B-\mu_B)])dx_B \end{aligned} p(xA)=Z1∫xB∈Rnexp(−[21(xA−μA)TVAA(xA−μA)+21(xA−μA)TVAB(xB−μB)+21(xB−μB)TVBA(xA−μA)+21(xB−μB)TVBB(xB−μB)])dxB
其中 Z Z Z是一个常数,不依赖于 x A x_A xA或 x B x_B xB,我们暂时忽略它。如果你以前没有使用过分块矩阵,那么上面的展开对你来说可能有点神奇。这类似于当定义一个二次形式基于某个矩阵 A A A时,则可得:
x T A x = ∑ i ∑ j A i j x i x j = x 1 A 11 x 1 + x 1 A 12 x 2 + x 2 A 21 x 1 + x 2 A 22 x 2 x^TAx=\sum_i\sum_jA_{ij}x_ix_j=x_1A_{11}x_1+x_1A_{12}x_2+x_2A_{21}x_1+x_2A_{22}x_2 xTAx=i∑j∑Aijxixj=x1A11x1+x1A12x2+x2A21x1+x2A22x2
花点时间自己研究一下,上面的矩阵推广也适用。
3.2.3 x B x_B xB上的积分
为了求积分,我们要对 x B x_B xB积分。然而,一般来说,高斯积分是很难手工计算的。我们能做些什么来节省计算时间吗?事实上,有许多高斯积分的答案是已知的(见特点#2)。那么,本节的基本思想是将上一节中的积分转换为一种形式,在这种形式中,我们可以应用特点#2中的一个结果,以便轻松地计算所需的积分。
这其中的关键是一个数学技巧,称为“配方法(completion of squares)”。考虑二次函数 。其中
1 2 x T A x + b T z + c = 1 2 ( z + A − 1 b ) T A ( z + A − 1 b ) + c − 1 2 b T A − 1 b \frac 12x^TAx+b^Tz+c=\frac 12(z+A^{-1}b)^TA(z+A^{-1}b)+c-\frac 12b^TA^{-1}b 21xTAx+bTz+c=21(z+A−1b)TA(z+A−1b)+c−21bTA−1b
下面使用单变量代数中的“配方法”来泛华的多元变量的等式:
1 2 a z 2 + b z + c = 1 2 a ( z + b z ) 2 + c − b 2 2 a \frac 12az^2+bz+c=\frac 12a(z+\frac bz)^2+c-\frac {b^2}{2a} 21az2+bz+c=21a(z+zb)2+c−2ab2
若要将配方法应用于上述情形,令
z = x B − μ B A = V B B b = V B A ( x A − μ A ) c = 1 2 ( x A − μ A ) T V A A ( x A − μ A ) \begin{aligned} z &= x_B-\mu_B \\ A &= V_{BB} \\ b &=V_{BA}(x_A-\mu_A) \\ c &=\frac 12(x_A-\mu_A)^TV_{AA}(x_A-\mu_A) \end{aligned} zAbc=xB−μB=VBB=VBA(xA−μA)=21(xA−μA)TVAA(xA−μA)
然后,这个积分可以重写为
p ( x A ) = 1 Z ∫ x B ∈ R n e x p ( − [ 1 2 ( x B − μ B ) T V A A ( x A − μ A ) + 1 2 ( x A − μ A ) T V A B ( x B − μ B ) + 1 2 ( x B − μ B ) T V B A ( x A − μ A ) + 1 2 ( x B − μ B ) T V B B ( x B − μ B ) ] ) d x B \begin{aligned} p(x_A)=\frac 1Z\int_{x_B\in R^n}exp(-&[\frac 12(x_B-\mu_B)^TV_{AA}(x_A-\mu_A)+\frac 12(x_A-\mu_A)^TV_{AB}(x_B-\mu_B) \\ & +\frac 12(x_B-\mu_B)^TV_{BA}(x_A-\mu_A)+\frac 12(x_B-\mu_B)^TV_{BB}(x_B-\mu_B)])dx_B \end{aligned} p(xA)=Z1∫xB∈Rnexp(−[21(xB−μB)TVAA(xA−μA)+21(xA−μA)TVAB(xB−μB)+21(xB−μB)TVBA(xA−μA)+21(xB−μB)TVBB(xB−μB)])dxB
我们可以提出不包括 x B x_B xB的项,
p ( x A ) = exp ( − 1 2 ( x A − μ A ) T V A A ( x A − μ A ) + 1 2 ( x A − μ A ) T V A B V B B − 1 V B A ( x A − μ A ) ) ⋅ 1 Z ∫ x B ∈ R n exp ( − 1 2 [ ( x B − μ B + V B B − 1 V B A ( x A − μ A ) ) T V B B ( x B − μ B + V B B − 1 V B A ( x A − μ A ) ) ] ) d x B \begin{aligned} p(x_{A})&=\exp\left(-\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A A}\left(x_{A}-\mu_{A}\right)+\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A B} V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right) \\ &\quad \cdot \frac{1}{Z} \int_{x_{B} \in \mathbb{R}^{n}} \exp \left(-\frac{1}{2}\left[\left(x_{B}-\mu_{B}+V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right)^{T} V_{B B}\left(x_{B}-\mu_{B}+V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right)\right]\right) d x_{B} \end{aligned} p(xA)=exp(−21(xA−μA)TVAA(xA−μA)+21(xA−μA)TVABVBB−1VBA(xA−μA))⋅Z1∫xB∈Rnexp(−21[(xB−μB+VBB−1VBA(xA−μA))TVBB(xB−μB+VBB−1VBA(xA−μA))])dxB
现在,我们可以应用特点#2。特别的,我们知道通常情况下随机变量为 x x x多元高斯分布,如果设均值 μ \mu μ,协方差矩阵 Σ \Sigma Σ,则概率密度函数可以得到如下式子:
1 ( 2 π ) n / 2 ∣ Σ ∣ 1 / 2 ∫ R n exp ( − 1 2 ( x − μ ) T Σ − 1 ( x − μ ) ) = 1 \frac{1}{(2 \pi)^{n / 2}|\Sigma|^{1 / 2}} \int_{\mathbf{R}^{n}} \exp \left(-\frac{1}{2}(x-\mu)^{T} \Sigma^{-1}(x-\mu)\right)=1 (2π)n/2∣Σ∣1/21∫Rnexp(−21(x−μ)TΣ−1(x−μ))=1
或等价与下式:
∫ R n exp ( − 1 2 ( x − μ ) T Σ − 1 ( x − μ ) ) = ( 2 π ) n / 2 ∣ Σ ∣ 1 / 2 \int_{R^{n}} \exp \left(-\frac{1}{2}(x-\mu)^{T} \Sigma^{-1}(x-\mu)\right)=(2 \pi)^{n / 2}|\Sigma|^{1 / 2} ∫Rnexp(−21(x−μ)TΣ−1(x−μ))=(2π)n/2∣Σ∣1/2
我们用这个事实来消去表达式中剩下的积分以得到 p ( x A ) p(x_A) p(xA):
p ( x A ) = 1 Z ⋅ ( 2 π ) n / 2 ∣ V B B ∣ 1 / 2 ⋅ exp ( − 1 2 ( x A − μ A ) T ( V A A − V A B V B B − 1 V B A ) ( x A − μ A ) ) p\left(x_{A}\right)=\frac{1}{Z} \cdot(2 \pi)^{n / 2}\left|V_{B B}\right|^{1 / 2} \cdot \exp \left(-\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T}\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)\left(x_{A}-\mu_{A}\right)\right) p(xA)=Z1⋅(2π)n/2∣VBB∣1/2⋅exp(−21(xA−μA)T(VAA−VABVBB−1VBA)(xA−μA))
3.2.4 论述得到的概率密度函数是高斯函数
这时我们几乎已经完成了全部计算!忽略前面的归一化常数,我们看到 x A x_A xA的概率密度函数是 x A x_A xA的二次形的指数。我们可以很快意识到概率密度函数就是均值向量为 μ A \mu_A μA,协方差矩阵为 ( V A A − V A B V B B − 1 V B A ) − 1 (V_{A A}-V_{A B} V_{B B}^{-1} V_{B A})^{-1} (VAA−VABVBB−1VBA)−1的高数分布。虽然协方差矩阵的形式看起来有点复杂,但是我们已经完成了我们开始想要展示的概念——即 x A x_A xA有一个边缘高斯分布。利用前面的逻辑,我们可以得出这个协方差矩阵必须以某种方式消去 Σ A A \Sigma_{AA} ΣAA。
但是,如果你好奇,也可以证明我们的推导与之前的证明是一致的。为此,我们对分块矩阵使用以下结果:
[ A B C D ] − 1 = [ M − 1 − M − 1 B D − 1 − D − 1 C M − 1 D − 1 + D − 1 C M − 1 B D − 1 ] \left[ \begin{array}{cc}{A} & {B} \\ {C} & {D}\end{array}\right]^{-1}=\left[ \begin{array}{cc}{M^{-1}} & {-M^{-1} B D^{-1}} \\ {-D^{-1} C M^{-1}} & {D^{-1}+D^{-1} C M^{-1} B D^{-1}}\end{array}\right] [ACBD]−1=[M−1−D−1CM−1−M−1BD−1D−1+D−1CM−1BD−1]
其中 M = A − B D − 1 C M=A-B D^{-1} C M=A−BD−1C。这个公式可以看作是 2 × 2 2\times 2 2×2矩阵显式逆矩阵的多变量推广:
[ a b c d ] − 1 = 1 a d − b c [ d − b − c a ] \left[ \begin{array}{ll}{a} & {b} \\ {c} & {d}\end{array}\right]^{-1}=\frac{1}{a d-b c} \left[ \begin{array}{cc}{d} & {-b} \\ {-c} & {a}\end{array}\right] [acbd]−1=ad−bc1[d−c−ba]
用这个公式,可以得出:
[ Σ A A Σ A B Σ B A Σ B B ] = [ V A A V A B V B A V B B ] − 1 = [ ( V A A − V A B V B B − 1 V B A ) − 1 − ( V A A − V A B V B B − 1 V B A ) − 1 V A B V B B − 1 − V B B − 1 V B A ( V A A − V A B V B B − 1 V B A ) − 1 ( V B B − V B A V A A − 1 V A B ) − 1 ] \begin{aligned} \left[ \begin{array}{cc}{\Sigma_{A A}} & {\Sigma_{A B}} \\ {\Sigma_{B A}} & {\Sigma_{B B}}\end{array}\right] &=\left[ \begin{array}{ll}{V_{A A}} & {V_{A B}} \\ {V_{B A}} & {V_{B B}}\end{array}\right]^{-1} \\ &=\left[ \begin{array}{cc}{\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1}} & {-\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1} V_{A B} V_{B B}^{-1}} \\ {-V_{B B}^{-1} V_{B A}\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1}} & {\left(V_{B B}-V_{B A} V_{A A}^{-1} V_{A B}\right)^{-1}}\end{array}\right] \end{aligned} [ΣAAΣBAΣABΣBB]=[VAAVBAVABVBB]−1=[(VAA−VABVBB−1VBA)−1−VBB−1VBA(VAA−VABVBB−1VBA)−1−(VAA−VABVBB−1VBA)−1VABVBB−1(VBB−VBAVAA−1VAB)−1]
正如我们所期望的那样,我们马上就能得出 ( V A A − V A B V B B − 1 V B A ) − 1 = Σ A A \left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1}=\Sigma_{A A} (VAA−VABVBB−1VBA)−1=ΣAA。
3.3 联合高斯分布的条件分布是高斯分布
本规则的正式表述为:
假设:
[ x A x B ] ∼ N ( [ μ A μ B ] , [ Σ A A Σ A B Σ B A Σ B B ] ) \left[ \begin{array}{l}{x_{A}} \\ {x_{B}}\end{array}\right]\sim\mathcal{N}\left(\left[ \begin{array}{l}{\mu_{A}} \\ {\mu_{B}}\end{array}\right], \left[ \begin{array}{cc}{\Sigma_{A A}} & {\Sigma_{A B}} \\ {\Sigma_{B A}} & {\Sigma_{B B}}\end{array}\right]\right) [xAxB]∼N([μAμB],[ΣAAΣBAΣABΣBB])
其中 x A ∈ R m , x B ∈ R n x_{A} \in \mathbf{R}^{m}, x_{B} \in \mathbf{R}^{n} xA∈Rm,xB∈Rn,并选择均值向量和协方差矩阵子块的维数来匹配 x A x_A xA和 x B x_B xB。则条件概率密度函数为:
p ( x A ∣ x B ) = p ( x A , x B ; μ , Σ ) ∫ x A ∈ R m p ( x A , x B ; μ , Σ ) d x A p ( x B ∣ x A ) = p ( x A , x B ; μ , Σ ) ∫ x B ∈ R n p ( x A , x B ; μ , Σ ) d x B \begin{aligned} p\left(x_{A} | x_{B}\right) &=\frac{p\left(x_{A}, x_{B} ; \mu, \Sigma\right)}{\int_{x_{A} \in \mathbb{R}^{m}} p\left(x_{A}, x_{B} ; \mu, \Sigma\right) d x_{A}} \\ p\left(x_{B} | x_{A}\right) &=\frac{p\left(x_{A}, x_{B} ; \mu, \Sigma\right)}{\int_{x_{B} \in \mathbb{R}^{n}} p\left(x_{A}, x_{B} ; \mu, \Sigma\right) d x_{B}} \end{aligned} p(xA∣xB)p(xB∣xA)=∫xA∈Rmp(xA,xB;μ,Σ)dxAp(xA,xB;μ,Σ)=∫xB∈Rnp(xA,xB;μ,Σ)dxBp(xA,xB;μ,Σ)
同样是高斯分布:
x A ∣ x B ∼ N ( μ A + Σ A B Σ B B − 1 ( x B − μ B ) , Σ A A − Σ A B Σ B B − 1 Σ B A ) x B ∣ x A ∼ N ( μ B + Σ B A Σ A A − 1 ( x A − μ A ) , Σ B B − Σ B A Σ A A − 1 Σ A B ) \begin{array}{l}{x_{A}\left|x_{B} \sim \mathcal{N}\left(\mu_{A}+\Sigma_{A B} \Sigma_{B B}^{-1}\left(x_{B}-\mu_{B}\right), \Sigma_{A A}-\Sigma_{A B} \Sigma_{B B}^{-1} \Sigma_{B A}\right)\right.} \\ {x_{B} | x_{A} \sim \mathcal{N}\left(\mu_{B}+\Sigma_{B A} \Sigma_{A A}^{-1}\left(x_{A}-\mu_{A}\right), \Sigma_{B B}-\Sigma_{B A} \Sigma_{A A}^{-1} \Sigma_{A B}\right)}\end{array} xA∣∣xB∼N(μA+ΣABΣBB−1(xB−μB),ΣAA−ΣABΣBB−1ΣBA)xB∣xA∼N(μB+ΣBAΣAA−1(xA−μA),ΣBB−ΣBAΣAA−1ΣAB)
和之前一样,我们只研究条件分布 x B ∣ x A x_B|x_A xB∣xA,另一个结果是对称的。我们的推导过程如下:
- 明确写出条件概率密度函数的表达式。
- 通过划分逆协方差矩阵重写表达式。
- 使用“平方和”参数。
- 论述得到的概率密度函数是高斯函数。
下面让我们分别研究一下上面提到的每一个步骤。
3.3.1 明确写出条件概率密度函数的表达式
假设我们想直接计算给定 x A x_A xA下 x B x_B xB的概率密度函数。则我们需要计算下式:
p ( x B ∣ x A ) = p ( x A , x B ; μ , Σ ) ∫ x B ∈ R m p ( x A , x B ; μ , Σ ) d x A = 1 Z ′ exp ( − 1 2 [ x A − μ A x B − μ B ] T [ Σ A A Σ A B Σ B A Σ B B ] − 1 [ x A − μ A x B − μ B ] ) \begin{aligned} p\left(x_{B} | x_{A}\right) &=\frac{p\left(x_{A}, x_{B} ; \mu, \Sigma\right)}{\int_{x_{B} \in R^m} p\left(x_{A}, x_{B} ; \mu, \Sigma\right) d x_{A}} \\ &=\frac{1}{Z^{\prime}} \exp \left(-\frac{1}{2} \left[ \begin{array}{c}{x_{A}-\mu_{A}} \\ {x_{B}-\mu_{B}}\end{array}\right]^{T} \left[ \begin{array}{cc}{\Sigma_{A A}} & {\Sigma_{A B}} \\ {\Sigma_{B A}} & {\Sigma_{B B}}\end{array}\right]^{-1} \left[ \begin{array}{c}{x_{A}-\mu_{A}} \\ {x_{B}-\mu_{B}}\end{array}\right]\right) \end{aligned} p(xB∣xA)=∫xB∈Rmp(xA,xB;μ,Σ)dxAp(xA,xB;μ,Σ)=Z′1exp(−21[xA−μAxB−μB]T[ΣAAΣBAΣABΣBB]−1[xA−μAxB−μB])
其中 Z ′ Z' Z′是一个归一化常数,我们用该常数表达不依赖于 x B x_B xB的因子。注意,这一次,我们甚至不需要计算任何积分——积分的值不依赖于 x B x_B xB,因此积分可以化简成归一化常数 Z ′ Z' Z′。
3.3.2 通过划分逆协方差矩阵重写表达式
和之前一样,我们用矩阵 V V V重新参数化概率密度函数,由此得到下式:
p ( x B ∣ x A ) = 1 Z ′ exp ( − 1 2 [ x A − μ A x B − μ B ] T [ V A A V A B V B A V B B ] [ x A − μ A x B − μ B ] ) = 1 Z ′ exp ( − [ 1 2 ( x A − μ A ) T V A A ( x A − μ A ) + 1 2 ( x A − μ A ) T V A B ( x B − μ B ) + 1 2 ( x B − μ B ) T V B A ( x A − μ A ) + 1 2 ( x B − μ B ) T V B B ( x B − μ B ) ] ) \begin{aligned} p\left(x_{B} | x_{A}\right) &=\frac{1}{Z^{\prime}} \exp \left(-\frac{1}{2} \left[ \begin{array}{c}{x_{A}-\mu_{A}} \\ {x_{B}-\mu_{B}}\end{array}\right]^{T} \left[ \begin{array}{cc}{V_{A A}} & {V_{A B}} \\ {V_{B A}} & {V_{B B}}\end{array}\right] \left[ \begin{array}{c}{x_{A}-\mu_{A}} \\ {x_{B}-\mu_{B}}\end{array}\right]\right) \\ &=\frac{1}{Z^{\prime}} \exp (-[\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A A}\left(x_{A}-\mu_{A}\right)+\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A B}\left(x_{B}-\mu_{B}\right) \\ &\qquad\qquad\qquad+\frac{1}{2}\left(x_{B}-\mu_{B}\right)^{T} V_{B A}\left(x_{A}-\mu_{A}\right)+\frac{1}{2}\left(x_{B}-\mu_{B}\right)^{T} V_{B B}\left(x_{B}-\mu_{B}\right) ] ) \end{aligned} p(xB∣xA)=Z′1exp(−21[xA−μAxB−μB]T[VAAVBAVABVBB][xA−μAxB−μB])=Z′1exp(−[21(xA−μA)TVAA(xA−μA)+21(xA−μA)TVAB(xB−μB)+21(xB−μB)TVBA(xA−μA)+21(xB−μB)TVBB(xB−μB)])
3.3.3 使用“平方和”参数
回忆下面这个式子:
1 2 z T A z + b T z + c = 1 2 ( z + A − 1 b ) T A ( z + A − 1 b ) + c − 1 2 b T A − 1 b \frac{1}{2} z^{T} A z+b^{T} z+c=\frac{1}{2}\left(z+A^{-1} b\right)^{T} A\left(z+A^{-1} b\right)+c-\frac{1}{2} b^{T} A^{-1} b 21zTAz+bTz+c=21(z+A−1b)TA(z+A−1b)+c−21bTA−1b
假设 A A A是一个对称的非奇异矩阵。如前所述,要将平方的补全应用于上述情况,令:
z = x B − μ B A = V B B b = V B A ( x A − μ A ) c = 1 2 ( x A − μ A ) T V A A ( x A − μ A ) \begin{aligned} z &=x_{B}-\mu_{B} \\ A &=V_{B B} \\ b &=V_{B A}\left(x_{A}-\mu_{A}\right) \\ c &=\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A A}\left(x_{A}-\mu_{A}\right) \end{aligned} zAbc=xB−μB=VBB=VBA(xA−μA)=21(xA−μA)TVAA(xA−μA)
然后,可以将 p ( x B ∣ x A ) p(x_B | x_A) p(xB∣xA)的表达式重写为:
p ( x B ∣ x A ) = 1 Z ′ exp ( − [ 1 2 ( x B − μ B + V B B − 1 V B A ( x A − μ A ) ) T V B B ( x B − μ B + V B B − 1 V B A ( x A − μ A ) ) + 1 2 ( x A − μ A ) T V A A ( x A − μ A ) − 1 2 ( x A − μ A ) T V A B V B B − 1 V B A ( x A − μ A ) ] ) \begin{array}{c}{p\left(x_{B} | x_{A}\right)=\frac{1}{Z^{\prime}} \exp \left(-\left[\frac{1}{2}\left(x_{B}-\mu_{B}+V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right)^{T} V_{B B}\left(x_{B}-\mu_{B}+V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right)\right.\right.} \\ {+\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A A}\left(x_{A}-\mu_{A}\right)-\frac{1}{2}\left(x_{A}-\mu_{A}\right)^{T} V_{A B} V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right) ] )}\end{array} p(xB∣xA)=Z′1exp(−[21(xB−μB+VBB−1VBA(xA−μA))TVBB(xB−μB+VBB−1VBA(xA−μA))+21(xA−μA)TVAA(xA−μA)−21(xA−μA)TVABVBB−1VBA(xA−μA)])
将不依赖于 x B x_B xB的指数部分化简到归一化常数中,得到:
p ( x B ∣ x A ) = 1 Z ′ ′ exp ( − 1 2 ( x B − μ B + V B B − 1 V B A ( x A − μ A ) ) T V B B ( x B − μ B + V B B − 1 V B A ( x A − μ A ) ) ) p\left(x_{B} | x_{A}\right)=\frac{1}{Z^{\prime \prime}} \exp \left(-\frac{1}{2}\left(x_{B}-\mu_{B}+V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right)^{T} V_{B B}\left(x_{B}-\mu_{B}+V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)\right)\right) p(xB∣xA)=Z′′1exp(−21(xB−μB+VBB−1VBA(xA−μA))TVBB(xB−μB+VBB−1VBA(xA−μA)))
3.3.4 论述得到的概率密度函数是高斯函数
看最后一个表达式,表达式 p ( x B ∣ x A ) p(x_B|x_A) p(xB∣xA)是均值为 μ B − V B B − 1 V B A ( x A − μ A ) \mu_B-V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right) μB−VBB−1VBA(xA−μA),协方差矩阵为 V B B − 1 V_{B B}^{-1} VBB−1的高斯概率密度函数。像往常一样,回忆一下矩阵等式:
[ Σ A A Σ A B Σ B A Σ B B ] = [ ( V A A − V A B V B B − 1 V B A ) − 1 − ( V A A − V A B V B B − 1 V B A ) − 1 V A B V B B − 1 − V B B − 1 V B A ( V A A − V A B V B B − 1 V B A ) − 1 ( V B B − V B A V A A − 1 V A B ) − 1 ] \left[ \begin{array}{cc}{\Sigma_{A A}} & {\Sigma_{A B}} \\ {\Sigma_{B A}} & {\Sigma_{B B}}\end{array}\right]= \left[ \begin{array}{c}{\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1}}&-\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1} V_{A B} V_{B B}^{-1} \\ {-V_{B B}^{-1} V_{B A}\left(V_{A A}-V_{A B} V_{B B}^{-1} V_{B A}\right)^{-1}}&\left(V_{B B}-V_{B A} V_{A A}^{-1} V_{A B}\right)^{-1}\end{array}\right] [ΣAAΣBAΣABΣBB]=[(VAA−VABVBB−1VBA)−1−VBB−1VBA(VAA−VABVBB−1VBA)−1−(VAA−VABVBB−1VBA)−1VABVBB−1(VBB−VBAVAA−1VAB)−1]
从上式可以推出:
μ B ∣ A = μ B − V B B − 1 V B A ( x A − μ A ) = μ B + Σ B A Σ A A − 1 ( x A − μ A ) \mu_{B | A}=\mu_{B}-V_{B B}^{-1} V_{B A}\left(x_{A}-\mu_{A}\right)=\mu_{B}+\Sigma_{B A} \Sigma_{A A}^{-1}\left(x_{A}-\mu_{A}\right) μB∣A=μB−VBB−1VBA(xA−μA)=μB+ΣBAΣAA−1(xA−μA)
反过来,我们也可以利用矩阵恒等式得到:
[ V A A V A B V B A V B B ] = [ ( Σ A A − Σ A B Σ B B − 1 Σ B A ) − 1 − ( Σ A A − Σ A B Σ A A − 1 Σ B B ) − 1 Σ A B Σ B B − 1 − Σ B B − 1 Σ B A ( Σ A A − Σ A B Σ B B − 1 Σ B A ) − 1 ( Σ B B − Σ B A Σ A A − 1 ∑ A B ) − 1 ] \left[ \begin{array}{cc}{V_{A A}} & {V_{A B}} \\ {V_{B A}} & {V_{B B}}\end{array}\right]= \left[ \begin{array}{c}{\left(\Sigma_{A A}-\Sigma_{A B} \Sigma_{B B}^{-1} \Sigma_{B A}\right)^{-1}}&-\left(\Sigma_{A A}-\Sigma_{A B} \Sigma_{A A}^{-1} \Sigma_{B B}\right)^{-1} \Sigma_{A B} \Sigma_{B B}^{-1} \\ {-\Sigma_{B B}^{-1} \Sigma_{B A}\left(\Sigma_{A A}-\Sigma_{A B} \Sigma_{B B}^{-1} \Sigma_{B A}\right)^{-1}}&\left(\Sigma_{B B}-\Sigma_{B A} \Sigma_{A A}^{-1} \sum_{A B}\right)^{-1}\end{array} \right] [VAAVBAVABVBB]=[(ΣAA−ΣABΣBB−1ΣBA)−1−ΣBB−1ΣBA(ΣAA−ΣABΣBB−1ΣBA)−1−(ΣAA−ΣABΣAA−1ΣBB)−1ΣABΣBB−1(ΣBB−ΣBAΣAA−1∑AB)−1]
由此推出:
Σ B ∣ A = V B B − 1 = Σ B B − Σ B A Σ A A − 1 Σ A B \Sigma_{B | A}=V_{B B}^{-1}=\Sigma_{B B}-\Sigma_{B A} \Sigma_{A A}^{-1} \Sigma_{A B} ΣB∣A=VBB−1=ΣBB−ΣBAΣAA−1ΣAB
我们完成了!
4. 总结
在本节的笔记中,我们使用了多元高斯的一些简单性质(加上一些矩阵代数技巧)来证明多元高斯分布满足许多封闭性质。一般来说,多元高斯分布是概率分布非常有用的表示形式,因为封闭性保证了这一点:即我们所希望的那样使用多元高斯分布执行的大多数类型的操作都可以以封闭形式完成。从分析的角度来看,涉及多元高斯的积分在实际应用中是往往是很好计算的,因为我们可以依赖于已知的高斯积分来避免自己进行积分。
5. 练习
理解题:令 A ∈ R n × n A\in R^{n\times n} A∈Rn×n是对称非奇异方阵, b ∈ R n , c b\in R^n,c b∈Rn,c,证明:
∫ x ∈ R n exp ( − 1 2 x T A x − x T b − c ) d x = ( 2 π ) n / 2 ∣ A ∣ 1 / 2 exp ( c − b T A − 1 b ) \int_{x \in \mathbf{R}^{n}} \exp \left(-\frac{1}{2} x^{T} A x-x^{T} b-c\right) d x=\frac{(2 \pi)^{n / 2}}{|A|^{1 / 2} \exp \left(c-b^{T} A^{-1} b\right)} ∫x∈Rnexp(−21xTAx−xTb−c)dx=∣A∣1/2exp(c−bTA−1b)(2π)n/2
参考资料
有关多元高斯的更多信息,请参见:
Bishop, Christopher M. Pattern Recognition and Machine Learning. Springer,2006.