多元正态分布

最新推荐文章于 2025-03-14 08:17:01 发布

omadesala

最新推荐文章于 2025-03-14 08:17:01 发布

阅读量2.3w

点赞数 20

分类专栏：基础概念文章标签：多元高斯分布

本文链接：https://blog.csdn.net/omade/article/details/28232047

版权

基础概念专栏收录该内容

13 篇文章

订阅专栏

多元正态分布

先定义一个d元随机向量，这里用列向量来表示，每一个元素都是一个一元随机变量，如 $X_i$

$\mathbf{X}=\begin{bmatrix}X_1\\ X_2\\ X_3\\ \vdots \\ X_i\\\vdots \\ X_d\\\end{bmatrix}$ ,其转置为 $\mathbf{X}^T=[X_1,X_2,X_3,...,X_i,...,X_d]^T$

其中表示这个多元随机变量的第i个分量，它是一个一维的随机变量。

高斯分布主要是用均值和方差来作为参数的分布，我们来看看随机向量的均值和方差

$E[\mathbf{X}]=\boldsymbol{\mu}^T=[\mu_1,\mu_2,\mu_3,...,\mu_i,...,\mu_d]^T$

关于方差，在多元分布里面，就是协方差矩阵

$\boldsymbol{\sum}=E\{(\mathbf{x}-\boldsymbol{\mu})(\mathbf{x}-\boldsymbol{\mu})^T \}$

$\boldsymbol{\sum}= \begin{bmatrix}E\{(X_1-\mu_1)(X_1-\mu_1)\} &E\{(X_1-\mu_1)(X_2-\mu_2)\} &\cdots &E\{(X_1-\mu_1)(X_j-\mu_j)\} & \cdots &E\{(X_1-\mu_1)(X_d-\mu_d) \}\\ E\{(X_2-\mu_2)(X_1-\mu_1)\}&E\{(X_2-\mu_2)(X_2-\mu_2)\} &\cdots &E\{(X_2-\mu_2)(X_j-\mu_j)\} &\cdots &E\{(X_2-\mu_2)(X_d-\mu_d)\} \\ \vdots&\vdots & \vdots & \vdots &\vdots &\vdots \\ E\{(X_i-\mu_i)(X_1-\mu_1)\}&E\{(X_j-\mu_j)(X_2-\mu_2)\} & \cdots&E\{(X_i-\mu_i)(X_j-\mu_j)\} &\cdots &E\{(X_i-\mu_i)(X_d-\mu_d)\} \\ \vdots&\vdots & \vdots & \vdots &\vdots &\vdots \\ E\{(X_d-\mu_d)(X_1-\mu_1)\} &E\{(X_d-\mu_d)(X_2-\mu_2)\} & \cdots &E\{(X_d-\mu_d)(X_j-\mu_j)\} &\cdots &E\{(X_d-\mu_d)(X_d-\mu_d)\} \end{bmatrix}$

$\boldsymbol{\sum}= \begin{bmatrix}\delta_{11} &\delta_{12} &\cdots &\delta_{1j} & \cdots &\delta_{1d}\\ \delta_{21}&\delta_{22} &\cdots &\delta_{2j} &\cdots &\delta_{2d} \\ \vdots&\vdots & \vdots & \vdots &\vdots &\vdots \\ \delta_{i1}&\delta_{i2} & \cdots&\delta_{ij} &\cdots &\delta_{id} \\ \vdots&\vdots & \vdots & \vdots &\vdots &\vdots \\\delta_{d1} &\delta_{d2} & \cdots &\delta_{dj} &\cdots &\delta_{dd} \end{bmatrix}$

其中 $\delta_{ij}=E\{ (X_i-\mu_i)(X_j-\mu_j) \}$

$\boldsymbol{\sum}^{-1}$ 是协方差矩阵的逆， $\left | \boldsymbol{\sum} \right |$ 是协方差矩阵的行列式

$\left | \boldsymbol{\sum} \right |$ 是非负矩阵，但是我们只考虑正定矩阵的情况

现在来看多元正态分布的函数表达式

$f(\mathbf{X})=\frac{1}{(2\pi)^{\frac{d}{2}}\left | \boldsymbol{\sum}\right |^{\frac{1}{2}}}e^{-\frac{1}{2}(\mathbf{X}-\boldsymbol{\mu})^T\boldsymbol{\sum}^{-1}(\mathbf{X}-\boldsymbol{\mu})}$

现在来看二元正态分布的情况

均值 $\mathbf{X}=\begin{bmatrix}X_1\\ X_2\\ \end{bmatrix}$

协方差矩阵 $\boldsymbol{\sum}= \begin{bmatrix}E\{(X_1-\mu_1)(X_1-\mu_1)\} &E\{(X_1-\mu_1)(X_2-\mu_2)\} \\ E\{(X_2-\mu_2)(X_1-\mu_1)\}&E\{(X_2-\mu_2)(X_2-\mu_2)\} \end{bmatrix}=\begin{bmatrix}\delta_{11} &\delta_{12} \\ \delta_{21}&\delta_{22} \end{bmatrix}$

在来计算二维正态分布的协方差矩阵的行列式

因为 $\delta_{12}=\delta_{21}$

$\left |\boldsymbol{\sum} \right |=\delta_{11}*\delta_{22}-\delta_{12}*\delta_{21}=\delta_{11}*\delta_{22}-\delta_{12}^2$

为了求出协方差矩阵的逆，首先看一下一些线性代数的概念，假设有矩阵

$A=\begin{vmatrix}A11 &A12 &A13 \\ A21&A22 &A23 \\ A31 &A32 &A33 \end{vmatrix}$

那么第i行，第j列的代数余子式为去掉A第i行，第j列之后的矩阵的行列式，记为 $M_{ij}$

那上面的矩阵A为例子，那么 $M_{11}=\begin{vmatrix}&A22 &A23 \\ &A32 &A33 \end{vmatrix}$

那么代数余子式为 $C_{ij}=(-1)^{i+j}M_{ij}$

那么我们可以定义：矩阵A的伴随矩阵是A的余子矩阵的转置矩阵

同样用矩阵A来做例子，有如下公式

$A^*=\begin{pmatrix}(-1)^{(1+1)}\begin{vmatrix} &A22 &A23 \\ &A32 &A33 \end{vmatrix} &(-1)^{(1+2)}\begin{vmatrix} &A21 &A23 \\ &A31 &A33 \end{vmatrix} &(-1)^{(1+3)}\begin{vmatrix} &A21 &A22 \\ &A31 &A32 \end{vmatrix} \\ (-1)^{(2+1)}\begin{vmatrix} &A12 &A13 \\ &A32 &A33 \end{vmatrix} & (-1)^{(2+2)}\begin{vmatrix} &A11 &A13 \\ &A31 &A33 \end{vmatrix} & (-1)^{(2+3)}\begin{vmatrix} &A11 &A12 \\ &A31 &A32 \end{vmatrix} \\ (-1)^{(3+1)}\begin{vmatrix} &A12 &A13 \\ &A22 &A23 \end{vmatrix} & (-1)^{(3+2)}\begin{vmatrix} &A11 &A13 \\ &A21 &A23 \end{vmatrix} &(-1)^{(3+3)}\begin{vmatrix} &A11 &A12 \\ &A21 &A22 \end{vmatrix} \end{pmatrix}^T$

那么矩阵A的逆矩阵就是

$A^{-1}=\frac{A^*}{\left |A \right |}$

现在来计算二元高斯函数的协方差矩阵的逆矩阵

$\boldsymbol{\sum}^{-1}=\frac{\boldsymbol{\sum}^*}{\left|\boldsymbol{\sum}\right |}$

那么协方差矩阵的伴随矩阵为：

$\boldsymbol{\sum}^*=\begin{pmatrix}(-1)^{1+1}\delta_{22} &(-1)^{1+2}\delta_{21} \\ (-1)^{2+1}\delta_{12} &(-1)^{2+2}\delta_{11} \end{pmatrix}^T=\begin{pmatrix}\delta_{22} &-\delta_{21} \\ -\delta_{12} &\delta_{11} \end{pmatrix}^T=\begin{pmatrix}\delta_{22} &-\delta_{12} \\ -\delta_{21} &\delta_{11} \end{pmatrix}$

那么协方差矩阵的逆矩阵为

$\boldsymbol{\sum}^{-1}=\frac{1}{(\delta_{11}*\delta_{22}-\delta_{12}^2)}\begin{pmatrix}\delta_{22} &-\delta_{12} \\ -\delta_{21} &\delta_{11} \end{pmatrix}$

指数部分

$(\mathbf{x}-\boldsymbol{\mu})^T \boldsymbol{{\sum}^{-1}}(\mathbf{x}-\boldsymbol{\mu})=\begin{bmatrix} X_1-\mu_1 , X_2-\mu_2\\ \end{bmatrix}\left (\frac{1}{\delta_{11}*\delta_{22}-\delta_{12}^2}\begin{pmatrix}\delta_{22} &-\delta_{12} \\ -\delta_{21} &\delta_{11} \end{pmatrix} \right )\begin{bmatrix} X_1-\mu_1 \\ X_2-\mu_2 \end{bmatrix}$

$=\frac{1}{\delta_{11}*\delta_{22}-\delta_{12}^2}\begin{bmatrix}X_1-\mu_1 , X_2-\mu_2\\ \end{bmatrix}\begin{pmatrix}\delta_{22} &-\delta_{12} \\ -\delta_{21} &\delta_{11} \end{pmatrix} \begin{bmatrix} X_1-\mu_1 \\ X_2-\mu_2 \end{bmatrix}$

$=\frac{1}{ \delta_{11}* \delta_{22}-\delta_{12}^2}\begin{bmatrix} (X_1-\mu_1) \delta_{22} + (X_2-\mu_2)(-\delta_{21})+(X_1-\mu_1)(-\delta_{12})+(X_2-\mu_2)(\delta_{11}) \end{bmatrix} \begin{bmatrix} X_1-\mu_1 \\ X_2-\mu_2 \end {bmatrix}$

$=\frac{1}{\delta_{11}*\delta_{22}-\delta_{12}^2} \begin{bmatrix} (X_1-\mu_1)\delta_{22}+(X_2-\mu_2)(-\delta_{21})+(X_1-\mu_1)(-\delta_{12})+(X_2-\mu_2)(\delta_{11})\end{bmatrix}\begin{bmatrix} X_1-\mu_1 \\ X_2-\mu_2 \end{bmatrix}$

$=\frac{1}{\delta_{11}*\delta_{22}-\delta_{12}^2}\begin{bmatrix} ((X_1-\mu_1)\delta_{22}+(X_2-\mu_2)(-\delta_{21}))(X_1-\mu_1)+((X_1-\mu_1)\(-\delta_{12})+(X_2-\mu_2)(\delta_{11}))( X_2-\mu_2 )\end{bmatrix}$

$=\frac{1}{\delta_{11}*\delta_{22}-\delta_{12}^2}\begin{bmatrix} (X_1-\mu_1)^2\delta_{22}-(X_1-\mu_1)(X_2-\mu_2)\delta_{21}-(X_1-\mu_1)( X_2-\mu_2)\delta_{12}+(X_2-\mu_2)^2\delta_{11}\end{bmatrix}$

考虑到这里 $\delta_{ij}=\delta_{ji}$ 得到下面的结论

$=\frac{1}{\delta_{11}*\delta_{22}-\delta_{12}^2}\begin{bmatrix} (X_1-\mu_1)^2\delta_{22}-2(X_1-\mu_1)(X_2-\mu_2)\delta_{21}+(X_2-\mu_2)^2\delta_{11}\end{bmatrix}$

那么二元正态分布的函数表达式为

$f(X_1,X_2)=\frac{1}{2\pi \sqrt{\delta_{11}*\delta_{22}-\delta_{12}^2}}e^{-\frac{1}{2(\delta_{11}*\delta_{22}-\delta_{12}^2)}\begin{bmatrix}(X_1-\mu_1)^2\delta_{22}-2(X_1-\mu_1)(X_2-\mu_2)\delta_{21}+(X_2-\mu_2)^2\delta_{11}\end{bmatrix}}$

如果我们引入相关性，可以得到另外一种形式，先看相关性的定义

$\rho_{X1,X2}=\frac{Cov(X_1,X_2)}{\delta_{X1}\delta_{X2}}$

<b

其中 $Cov(X_1,X_2)=E[(X1-\mu_1)(X2-\mu2)]$

$\delta_{X1}=\sqrt{E[(X1-\mu_1)(X1-\mu1)]}$

        $=\sqrt{E[(X1-E[X1])(X1-E[X1])]}=\sqrt{E[E[X1]^2-2E[X]X1+E[X]^2]]}$
        $=\sqrt{E[X1]^2-2E[X]^2+E[X]^2]}==\sqrt{E[X1]^2-E[X1]^2}$

同样

$\delta_{X2}=\sqrt{E[X2]^2-E[X2]^2}$

那么

$\rho_{X1,X2}=\frac{\delta_{12}}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}$

变形得到 $\delta_{12}=\rho_{X1,X2}\sqrt{\delta_{11}}\sqrt{\delta_{22}}$ ，将其代入上面得到的二元正态分布函数中，可以得到另外一种表达方式

$f(X_1,X_2)=\frac{1}{2\pi \sqrt{\delta_{11}*\delta_{22}-(\rho_{X1X2}\sqrt{\delta_{11}}\sqrt{\delta_{22}})^2}}e^{-\frac{1}{2\left(\delta_{11}*\delta_{22}-\(\rho_{X1X2}\sqrt{\delta_{11}}\sqrt{\delta_{22}}\right )^2)}\begin{bmatrix}(X_1-\mu_1)^2\delta_{22}-2(X_1-\mu_1)(X_2-\mu_2)(\rho_{X1X2}\sqrt{\delta_{11}}\sqrt{\delta_{22}})+(X_2-\mu_2)^2\delta_{11}\end{bmatrix}}$

整理得到

$f(X_1,X_2)=\frac{1}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left(\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)\right)}\begin{bmatrix}(X_1-\mu_1)^2\delta_{22}-2(X_1-\mu_1)(X_2-\mu_2)(\rho_{X1X2}\sqrt{\delta_{11}}\sqrt{\delta_{22}})+(X_2-\mu_2)^2\delta_{11}\end{bmatrix}}$ 继续整理

$f(X_1,X_2)=\frac{1}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_1-\mu_1)^2\delta_{22}}{\delta_{11}*\delta_{22}}-\frac{2(X_1-\mu_1)(X_2-\mu_2)(\rho_{X1X2}\sqrt{\delta_{11}}\sqrt{\delta_{22}})}{\delta_{11}*\delta_{22}}+\frac{(X_2-\mu_2)^2\delta_{11}}{\delta_{11}*\delta_{22}}\end{bmatrix}}$ 继续整理

得到

$f(X_1,X_2)=\frac{1}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_1-\mu_1)^2}{\delta_{11}}-\frac{2\rho_{X1X2}(X_1-\mu_1)(X_2-\mu_2)}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}+\frac{(X_2-\mu_2)^2}{\delta_{22}}\end{bmatrix}}$

这就是用相关系数来表示的二维正态分布的表现形式

下面看看二维正态分布的条件分布公式

$f(X_1)=\frac{1}{\sqrt{2\pi\delta_{11}}}e^{-\frac{(X_1-\mu_1)^2}{\delta_{11}}}$

$f(X_2)=\frac{1}{\sqrt{2\pi\delta_{22}}}e^{-\frac{(X_2-\mu_2)^2}{\delta_{22}}}$

$f_{X1|X2}(X_1|X_2)=\frac{f(X_1,X_2)}{f_{X_2}(X_2)}$

$=\frac{\sqrt{2\pi\delta_{22}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_1-\mu_1)^2}{\delta_{11}}-\frac{2\rho_{X1X2}(X_1-\mu_1)(X_2-\mu_2)}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}+\frac{(X_2-\mu_2)^2}{\delta_{22}}-\frac{(1-\rho_{X_1X_2}^2)(X_2-\mu_2)^2}{\delta_{22}}\end{bmatrix}}$

$=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_1-\mu_1)^2}{\delta_{11}}-\frac{2\rho_{X1X2}(X_1-\mu_1)(X_2-\mu_2)}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}+\frac{\rho_{X_1X_2}^2(X_2-\mu_2)^2}{\delta_{22}}\end{bmatrix}}\\=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}[\frac{X_1-\mu1}{\sqrt{\delta_{11}}}-\frac{\rho_{X_1X_2}(X_2-\mu_2)}{\sqrt{\delta_{22}}}]^2}\\=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}[\frac{(X_1-\mu1)\sqrt{\delta_{22}}}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}-\frac{\rho_{X_1X_2}(X_2-\mu_2)\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}\sqrt{\delta_{11}}}]^2}\\=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}[\frac{(X_1-\mu1)\sqrt{\delta_{22}}-\rho_{X_1X_2}(X_2-\mu_2)\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}\sqrt{\delta_{11}}}]^2}$

$=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}[\frac{(X_1-\mu1)\sqrt{\delta_{22}}-\rho_{X_1X_2}(X_2-\mu_2)\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}\sqrt{\delta_{11}}}]^2}\\=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}[\frac{(X_1-\mu1)-\rho_{X_1X_2}(X_2-\mu_2)\frac{\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}}}{\sqrt{\delta_{11}}}]^2}\\=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}[\frac{X_1-(\mu1+\rho_{X_1X_2}(X_2-\mu_2))\frac{\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}}}{\sqrt{\delta_{11}}}]^2}\\=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-[\frac{\left(X_1-(\mu1+\rho_{X_1X_2}(X_2-\mu_2))\frac{\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}}\right )^2}{2\left((1-\rho_{X1X2}^2)\right)\delta_{11}}]}$

$f_{X2|X1}(X_2|X_1)=\frac{f(X_1,X_2)}{f_{X_1}(X_1)}$

$=\frac{\sqrt{2\pi\delta_{11}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{\rho_{X1X2}^2(X_1-\mu_1)^2}{\delta_{11}}-\frac{2\rho_{X1X2}(X_1-\mu_1)(X_2-\mu_2)}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}+\frac{(X_2-\mu_2)^2}{\delta_{22}}\end{bmatrix}}\\=\frac{\sqrt{2\pi\delta_{11}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{X_2-\mu_2}{\sqrt{\delta_{22}}}-\frac{\rho_{X1X2}(X_1-\mu_1)}{\sqrt{\delta_{11}}}\end{bmatrix}^2}\\$

$=\frac{\sqrt{2\pi\delta_{11}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{X_2-\mu_2}{\sqrt{\delta_{22}}}-\frac{\rho_{X1X2}(X_1-\mu_1)}{\sqrt{\delta_{11}}}\end{bmatrix}^2}\\=\frac{\sqrt{2\pi\delta_{11}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_2-\mu_2)\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}\sqrt{\delta_{11}}}-\frac{\rho_{X1X2}(X_1-\mu_1)\sqrt{\delta_{22}}}{\sqrt{\delta_{11}}\sqrt{\delta_{22}}}\end{bmatrix}^2}\\=\frac{\sqrt{2\pi\delta_{11}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_2-\mu_2)\sqrt{\delta_{11}}-(\rho_{X1X2}(X_1-\mu_1)\sqrt{\delta_{22}})}{\sqrt{\delta_{22}}\sqrt{\delta_{11}}}\end{bmatrix}^2}$

$=\frac{\sqrt{2\pi\delta_{11}}}{2\pi \sqrt{\delta_{11}*\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_2-\mu_2)\sqrt{\delta_{11}}-(\rho_{X1X2}(X_1-\mu_1)\sqrt{\delta_{22}})}{\sqrt{\delta_{22}}\sqrt{\delta_{11}}}\end{bmatrix}^2}\\=\frac{1}{ \sqrt{2\pi\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{1}{2\left((1-\rho_{X1X2}^2)\right)}\begin{bmatrix}\frac{(X_2-\mu_2)-(\rho_{X1X2}(X_1-\mu_1)\frac{\sqrt{\delta_{22}}}{\sqrt{\delta_{11}}})}{\sqrt{\delta_{22}}}\end{bmatrix}^2}\\=\frac{1}{ \sqrt{2\pi\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{\left((X_2-\mu_2)-(\rho_{X1X2}(X_1-\mu_1)\frac{\sqrt{\delta_{22}}}{\sqrt{\delta_{11}}})\right )}{2\left((1-\rho_{X1X2}^2)\right)\delta_{22}}^2}\\=\frac{1}{ \sqrt{2\pi\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{\left(X_2-\left(\mu_2+(\rho_{X1X2}(X_1-\mu_1)\right)\frac{\sqrt{\delta_{22}}}{\sqrt{\delta_{11}}})\right )}{2\left((1-\rho_{X1X2}^2)\right)\delta_{22}}^2}$

整理一下得到,条件分布公式为

$f_{X1|X2}(X_1|X_2)=\frac{1}{ \sqrt{2\pi\delta_{11}(1-\rho_{X1X2}^2)}}e^{-\frac{\left(X_1-(\mu1+\rho_{X_1X_2}(X_2-\mu_2))\frac{\sqrt{\delta_{11}}}{\sqrt{\delta_{22}}}\right )^2}{2\left((1-\rho_{X1X2}^2)\right)\delta_{11}}}$

$f_{X2|X1}(X_2|X_1)=\frac{1}{ \sqrt{2\pi\delta_{22}(1-\rho_{X1X2}^2)}}e^{-\frac{\left(X_2-\left(\mu_2+(\rho_{X1X2}(X_1-\mu_1)\right)\frac{\sqrt{\delta_{22}}}{\sqrt{\delta_{11}}})\right )}{2\left((1-\rho_{X1X2}^2)\right)\delta_{22}}^2}$

可以看出，只要求出其中一个，另外一个的坐标做相应的调换就可以了