概率论:多元高斯分布

一、多元高斯分布:

一元高斯分布的概率密度函数如下所示:
p(x)=1σ2πe12(xμσ)2(1)p(x)={1\over\sigma\sqrt{2\pi}}\centerdot e^{-{1\over2}({x-\mu\over\sigma})^2}\quad\quad\quad(1)

而如果我们对随机变量X进行标准化,用Z=XμσZ={X-\mu\over\sigma}对上式进行换元,可得:
x(z)=zσ+μx(z)=z\centerdot\sigma+\mu
p(x(z))=1σ2πe12(z)2p(x(z))={1\over\sigma\sqrt{2\pi}}\centerdot e^{-{1\over2}(z)^2}
1=+p(x(z))dx\therefore1=\int_{-\infty}^{+\infty}p(x(z))dx

=+1σ2πe12(z)2dx\quad\quad=\int_{-\infty}^{+\infty}{1\over\sigma\sqrt{2\pi}}\centerdot e^{-{1\over2}(z)^2}dx

=+12πe12(z)2dz\quad\quad=\int_{-\infty}^{+\infty}{1\over\sqrt{2\pi}}\centerdot e^{-{1\over2}(z)^2}dz

此时我们可以说随机变量ZN(0,1)Z\backsim N(0,1)服从一元标准高斯分布,其均值μ=0\mu=0,方差σ2=1\sigma^2=1,概率密度函数为:

p(z)=12πe12(z)2(2)p(z)={1\over\sqrt{2\pi}}\centerdot e^{-{1\over2}(z)^2}\quad\quad\quad(2)

1.1 多元标准高斯分布

多元标准高斯分布的概率密度函数是由(2)导出的
Z=[Z1,...,Zn]TZiN(0,1)Zi,Zj(i,j=1,...,nij)Zi2Z=[Z1,...,Zn]T假设有随机向量\overrightarrow{Z}=[Z_1,...,Z_n]^T,其中Z_i\backsim N(0,1)且Z_i,Z_j(i,j=1,...,n\bigwedge i\ne j)彼此独立,即随机向量中的每一个随机变量Z_i都服从标准高斯分布并且两两彼此独立,则由(2)与独立随机变量概率密度函数之间的关系,我们可得随机向量\overrightarrow{Z}=[Z_1,...,Z_n]^T的联合概率密度函数为:
p(z1,...,zn)=i=1n12πe12(zi)2=1(2π)n2e12(ZTZ)p(z_1,...,z_n)=\prod_{i=1}^n{1\over\sqrt{2\pi}}\centerdot e^{-{1\over 2}\centerdot(z_i)^2}={1\over(2\pi)^{n\over2}}\centerdot e^{-{1\over2}\centerdot(Z^TZ)}
且:1=+ ⁣+p(z1,...,zn)dz1...dzn1=\int_{-\infty}^{+\infty}\dots\int_{-\infty}^{+\infty}p(z_1,...,z_n)dz_1...dz_n

我们称随机向量ZN(0,I)\overrightarrow{Z}\backsim N(\overrightarrow{0},I),即随机向量服从均值为零向量,协方差矩阵为单位矩阵的高斯分布

1.2 多元高斯分布

对于普通的随机向量XN(μ,Σ)\overrightarrow{X}\backsim N(\mu,\Sigma),和其每个随机变量XiN(μi,σi2)(i=1,...,n)X_i\backsim N(\mu_i,\sigma_i^2)(i=1,...,n)Xi,Xj(i,j=1,...,n)X_i,X_j(i,j=1,...,n)彼此不独立的情况下,我们该怎么求随机向量X\overrightarrow{X}的联合概率密度函数呢,一个很自然的想法是:如果我们能通过线性变换,使得随机向量中的每个随机变量彼此独立,则我们也可以通过独立随机变量概率密度函数之间的关系求出其联合概率密度函数,事实上,我们有如下的定理可以完成这个工作:

1.2.1 定理1

XN(μ,Σ)μRnΣSn×nXBRn×n使Z=B1(Xμ)ZN(0,I)若存在随机向量\overrightarrow{X}\backsim N(\overrightarrow{\mu},\Sigma),其中\overrightarrow{\mu}\in R^n为均值向量,\Sigma\in S_{n\times n}为半正定实对称矩阵,是\overrightarrow{X}的协方差矩阵,则存在满秩矩阵B\in R^{n\times n},使得\overrightarrow{Z}=B^{-1}(\overrightarrow{X}-\overrightarrow{\mu}),而\overrightarrow{Z}\backsim N(\overrightarrow{0},I)

有了定理1,我们就可以对随机向量X\overrightarrow{X}作相应的线性变换,使其随机变量在线性变换后彼此独立,从而求出其联合概率密度函数,具体地:
Z=B1(Xμ),ZN(0,I)\because\overrightarrow{Z}=B^{-1}(\overrightarrow{X}-\overrightarrow{\mu}),Z\backsim N(\overrightarrow{0},I)

p(z1,...,zn)=1(2π)n2e12(ZTZ)\therefore p(z_1,...,z_n)={1\over(2\pi)^{n\over2}}\centerdot e^{-{1\over2}\centerdot(Z^TZ)}

p(z1(x1,...,xn),...)=1(2π)n2e12[(B1(Xμ))T(B1(Xμ))]\quad p(z_1(x_1,...,x_n),...) ={1\over(2\pi)^{n\over2}}\centerdot e^{-{1\over 2}[(B^{-1}(\overrightarrow{X}-\overrightarrow{\mu}))^T(B^{-1}(\overrightarrow{X}-\overrightarrow{\mu}))]}
=1(2π)n2e12[(Xμ)T(BBT)1(Xμ)](3)\quad\quad\quad\quad\quad\quad\quad\quad\quad={1\over(2\pi)^{n\over2}}\centerdot e^{-{1\over 2}[(\overrightarrow{X}-\overrightarrow{\mu})^T(BB^T)^{-1}(\overrightarrow{X}-\overrightarrow{\mu})]}\quad\quad\quad\quad (3)

https://zhuanlan.zhihu.com/p/58987388

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章