如何计算参数的协方差矩阵(透彻理解协方差矩阵)

协方差及协方差矩阵有着特别广泛的应用,在多元高斯分布、高斯过程、卡尔曼滤波等算法中多有用到,本文从协方差、协方差矩阵讲起,并重点讲解协方差矩阵在高斯分布中的用法及意义,也是讲解高斯过程、贝叶斯优化的铺垫。

协方差(Covariance)

X、Y两个随机变量的协方差在和中用于衡量两个变量的总体。用来刻画两个随机变量之间的相关性:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(1)

假定我们不知道潜在的概率分布,我们取n个样本来计算:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(2)

分别计算这n样本的两个变量的均值,这两个变量的协方差可以用下式来计算:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(3)

由于变量都有量纲,如果消除各自量纲影响,将协方差除以两个变量的标准差,则可得相关系数:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(4)

协方差矩阵

随机向量:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(5)

我们计算所有元素的两两协方差,形成协方差矩阵:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(6)

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(7)

这是一个对称矩阵,对角线是每个变量的方差。如果是对角阵,

协方差矩阵形式如下:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(8)

协方差矩阵与多元高斯分布

多元高斯分布概率密度的推导

设多元高斯分布如下:均值向量为μ,协方差矩阵为Σ

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(9)

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(10)

与一元高斯分布对比,概率密度函数形式有所变化,这个变化是怎么来的,我们通过二元高斯分布来推导一下这个密度函数的由来。

对于二元高斯分布,我们设定:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(11)

现在我们推导两个变量的高斯分布的密度函数公式:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(12)

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(13)

这个联合概率密度函数是各自概率密度函数的乘积,这表明两个变量是独立的。这个独立性反映在我们的协方差矩阵中,就是只有对角线元素不为零,两个变量是独立的,所以联合概率密度可以表示为两个变量概率密度的乘积

二维高斯分布函数图像

我们看相互独立的两个变量的二维高斯分布图像在XoY平面投影的函数表达式

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(14)

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(15)

令:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(16)

得:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(17)

显然这是一个椭圆曲线的表达式。

我们看两种情况,一种协方差矩阵是对角阵(变量相互独立),另一种是协方差矩阵是非对角阵(变量有关联)

  • 如果高斯随机向量具有对角线的协方差矩阵(所有变量都是不相关的,那么概率密度函数曲面在X0Y投影的椭圆曲线的两个轴平行于坐标轴。

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(18)

  • 如果高斯随机向量的协方差矩阵是非对角阵(一些变量是相关的),那么概率密度函数曲面在X0Y投影的椭圆曲线的两个轴仍然是相互重直,但与坐标轴并不平行。

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(19)

我们用matlab来形象展示一下:

下图是两个变量的均值都是零,协方差矩阵为:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(20)

其三维曲面如下:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(21)

在XOY平面的投影如下:

如何计算参数的协方差矩阵(透彻理解协方差矩阵)(22)

本文主要讲解了协方差矩阵及其在高斯分布中意义和用法。协方差矩阵在高斯过程中有着非常重要的意义,如果不能很好的理解协方差矩阵,就不能很好的理解高斯过程。

,

免责声明:本文仅代表文章作者的个人观点,与本站无关。其原创性、真实性以及文中陈述文字和内容未经本站证实,对本文以及其中全部或者部分内容文字的真实性、完整性和原创性本站不作任何保证或承诺,请读者仅作参考,并自行核实相关内容。文章投诉邮箱:anhduc.ph@yahoo.com

    分享
    投诉
    首页