协方差与协方差矩阵的定义及计算

谁借莪１个温暖的怀抱￠ 2022-12-20 08:44 209阅读 0赞

这篇博客主要介绍协方差与协方差矩阵的定义及其计算

首先给出常用的一些数理统计的计算公式：  
均值：  
![在这里插入图片描述][20201111160217961.png_pic_center]

标准差：  
![在这里插入图片描述][20201111160057310.png_pic_center]  
方差：  
![在这里插入图片描述][20201111160156298.png_pic_center]  
这几个公式大家肯定都是比较熟悉的，这里需要注意的地方是标准差以及方差计算时：**分母是n-1,这是因为这样能使我们以较小的样本集更好的逼近总体的标准差，即统计上所谓的“无偏估计”**。  
均值描述的是样本集合的中间点，标准差给我们描述的则是样本集合的各个样本点到均值的距离之平均，也就是样本的散布度。

标准差和方差一般是用来描述一维数据的，但现实生活我们常常遇到含有多维数据的数据集，面对这样的数据集，我们当然可以按照每一维独立的计算其方差，但是通常我们还想了解更多的内在联系，协方差就是这样一种用来度量两个随机变量关系的统计量，其结果表示两个随机变量 X 与 Y 之间相互关系，其定义如下：  
![在这里插入图片描述][20201111161055798.png_pic_center]

从协方差的定义上我们也可以看出一些显而易见的性质，如：  
![在这里插入图片描述][20201111162413112.png_pic_center]  
![在这里插入图片描述][20201111162423262.png_pic_center]  
那么协方差是如何表示两个随机变量 X 与 Y 之间相互关系的呢，通过下述三张图片来进行讲解：  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center]  
当 X, Y 的联合分布像上图那样时，我们可以看出，大致上有： X 越大 Y 也越大， X 越小 Y 也越小，这种情况，我们称为“正相关”。  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center 1]  
当X, Y 的联合分布像上图那样时，我们可以看出，大致上有：X 越大Y 反而越小，X 越小 Y 反而越大，这种情况，我们称为“负相关”。  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center 2]  
当X, Y 的联合分布像上图那样时，我们可以看出：既不是X 越大Y 也越大，也不是 X 越大 Y 反而越小，这种情况我们称为“ 不相关”。

看上述3张图，在每张图的  
区域（1）中有 X>EX ，Y-EY>0 ，所以(X-EX)(Y-EY)>0；  
区域（2）中，有 X<EX ，Y-EY>0 ，所以(X-EX)(Y-EY)<0；  
区域（3）中，有 X<EX ，Y-EY<0 ，所以(X-EX)(Y-EY)>0；  
区域（4）中，有 X>EX ，Y-EY<0 ，所以(X-EX)(Y-EY)<0。

对于第一张图，当X 与Y 正相关时，它们的分布大部分在区域（1）和（3）中，小部分在区域（2）和（4）中，所以平均来说，有E(X-EX)(Y-EY)>0 。

对于第二张图，当 X与 Y负相关时，它们的分布大部分在区域（2）和（4）中，小部分在区域（1）和 (3)中，所以平均来说，有(X-EX)(Y-EY)<0 。

对于第三张图，当 X与 Y不相关时，它们在区域（1）和（3）中的分布，与在区域（2）和（4）中的分布几乎一样多，所以平均来说，有(X-EX)(Y-EY)=0 。

用协方差来表征这种关系：  
当 cov(X, Y)>0时，表明 X与Y 正相关；  
当 cov(X, Y)<0时，表明X与Y负相关；  
当 cov(X, Y)=0时，表明X与Y不相关。

经过上述讲解，各位应该对协方差有所了解，而对于协方差矩阵而言，协方差矩阵的每个元素是各个向量元素之间的协方差，是从标量随机变量到高维度随机向量的自然推广。  
![在这里插入图片描述][watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center 3]

在具体实现时，可以通过先对样本进行中心化（ x i = x i − 1 / m ∑ i = 1 m x i x\_i = x\_i - 1/m\\sum\_\{i=1\}^mx\_i xi=xi−1/m∑i=1mxi），然后利用 X X T XX^T XXT计算样本的协方差矩阵,注意此时是没有添加协方差前面的系数的，如果需要添加系数，再处理一下就行也就是除以 m − 1 m-1 m−1.

参考文献：

[协方差矩阵][Link 1]  
[协方差的意义][Link 2]  
[协方差、协方差矩阵定义与计算][Link 3]

[20201111160217961.png_pic_center]: /images/20221120/8460dc9fa4b741a29e202c77e6b11a49.png
[20201111160057310.png_pic_center]: https://img-blog.csdnimg.cn/20201111160057310.png#pic_center
[20201111160156298.png_pic_center]: https://img-blog.csdnimg.cn/20201111160156298.png#pic_center
[20201111161055798.png_pic_center]: https://img-blog.csdnimg.cn/20201111161055798.png#pic_center
[20201111162413112.png_pic_center]: https://img-blog.csdnimg.cn/20201111162413112.png#pic_center
[20201111162423262.png_pic_center]: https://img-blog.csdnimg.cn/20201111162423262.png#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center]: https://img-blog.csdnimg.cn/20201111161219349.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center 1]: https://img-blog.csdnimg.cn/20201111161332140.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center 2]: https://img-blog.csdnimg.cn/20201111161347832.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2,size_16,color_FFFFFF,t_70#pic_center
[watermark_type_ZmFuZ3poZW5naGVpdGk_shadow_10_text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2_size_16_color_FFFFFF_t_70_pic_center 3]: https://img-blog.csdnimg.cn/20201111164219194.png?x-oss-process=image/watermark,type_ZmFuZ3poZW5naGVpdGk,shadow_10,text_aHR0cHM6Ly9ibG9nLmNzZG4ubmV0L3FxXzM4MDQ4NzU2,size_16,color_FFFFFF,t_70#pic_center
[Link 1]: https://baike.baidu.com/item/%E5%8D%8F%E6%96%B9%E5%B7%AE%E7%9F%A9%E9%98%B5/9822183?fr=aladdin
[Link 2]: https://blog.csdn.net/wuhzossibility/article/details/8087863
[Link 3]: https://blog.csdn.net/xw20084898/article/details/42077141