协方差矩阵的几何解释

Dear 丶 2022-07-13 04:44 195阅读 0赞

# A geometric interpretation of the covariance matrix #

http://www.visiondummy.com/2014/04/geometric-interpretation-covariance-matrix/

译文：http://demo.netfoucs.com/u010182633/article/details/45937051

## 介绍 ##

在本文中，我们通过探索线性变换与所得数据协方差之间的关系提供协方差矩阵一个直观的几何解释。大部分教科书基于协方差矩阵的概念解释数据的形状。相反，我们采取一个反向的方法，根据数据的形状来解释协方差矩阵的概念。

在《为什么样本方差除以N-1？》的文章中，我们会讨论方差的概念，并提供了众所周知的估算样本方差公式的推导和证明。这篇文章中使用的图1表明标准差（方差的平方根）提供了数据在特征空间上传播多少的量度。   
![这里写图片描述][20150523013736894]

我们发现，样本方差的无偏估计可由下式获得：   
![这里写图片描述][20150523013917027]

然而，方差只能用于解释平行于特征空间轴方向的数据传播。考虑图2所示的二维特征空间：   
![这里写图片描述][20150523014142950]

对于这个数据，我们可以计算出在x方向上的方差![这里写图片描述][20150523015110077]和y方向上的方差![这里写图片描述][20150523015132026]。然而，数据的水平传播和垂直传播不能解释明显的对角线关系。图2清楚地显示，平均而言，如果一个数据点的x值增加，则y值也将增加，这产生了正相关。这种相关性可以通过扩展方差概念到所谓的数据“协方差”捕捉到：   
![这里写图片描述][20150523015153507]

对于2D数据，我们得到![这里写图片描述][20150523015242351]，这些值可以用矩阵来表示，该矩阵叫做协方差矩阵：   
![这里写图片描述][20150523015344002]

如果x与y是正相关的，那么y和x也是正相关的。换句话说，![这里写图片描述][20150523020837615]。因此，协方差矩阵始终是一个对称矩阵，其对角线上是方差，非对角线上是协方差。二维正态分布数据由它的均值和2x2协方差矩阵就可以完全解释。同样，一个3x3协方差矩阵用于捕捉三维数据的传播，一个NxN协方差矩阵捕获N维数据的传播。

图3展示了数据的整体形状如何定义协方差矩阵：   
![这里写图片描述][20150523015824596]

## 协方差矩阵的特征值分解 ##

在下一节，我们将讨论协方差矩阵如何被解释为白色数据转换成我们观察到数据的线性操作。然而，在深入技术细节之前，对特征向量和特征值如何唯一地确定协方差矩阵（数据形状）有一个直观的认识是非常重要的。

正如我们在图3看到的，协方差矩阵定义了我们数据的传播（方差）和方向（协方差）。因此，如果我们想用一个向量和它的大小来表示协方差矩阵，我们应该简单地尝试找到指向数据最大传播方向上的向量，其大小等于这个方向上的传播（方差）。

如果我们定义这个向量为![这里写图片描述][20150523022152998]，那么我们数据D到这个向量上的映射为![这里写图片描述][20150523022351979]，映射数据的方差是![这里写图片描述][20150523022501357]。由于我们正在寻找指向最大方差方向的向量![这里写图片描述][20150523022300156]，所以我们应该选择它的成分，使得映射数据的协方差矩阵![这里写图片描述][20150523022524332]尽可能的大。最大化![这里写图片描述][20150523022313701]的形式为![这里写图片描述][20150523022639652]的任何函数，其中![这里写图片描述][20150523022211266]是归一化单位向量，可以用一个所谓的瑞利商表示。通过设置![这里写图片描述][20150523022358535]等于矩阵的最大特征特征向量![这里写图片描述][20150523022702643]可以获得这样瑞利商的最大值。

换句话说，协方差矩阵的最大特征向量总是指向数据最大方差的方向，并且该向量的幅度等于相应的特征值。第二大特征向量总是正交于最大特征向量，并指向第二[大数据][Link 1]的传播方向。

现在，让我们来看看一些例子。在文章《特征值和特征向量》中http://blog.csdn[.NET][]/u010182633/article/details/45921929，我们看到一个线性变换矩阵T完全由它的特征向量和特征值定义。应用到协方差矩阵，这意味着：   
![这里写图片描述][20150523023546860]   
![这里写图片描述][20150523023629620]

如果我们数据的协方差矩阵是对角矩阵，使得协方差是零，那么这意味着方差必须等于特征值λ。如图4所示，特征向量用绿色和品红色表示，特征值显然等于协方差矩阵的方差分量。   
![这里写图片描述][20150523023728563]

然而，如果协方差矩阵不是对角的，使得协方差不为零，那么情况稍微更复杂一些。特征值仍代表数据最大传播方向的方差大小，协方差矩阵的方差分量仍然表示x轴和y轴方向上的方差大小。但是，因为数据不是轴对齐的，所以这些值不再与图5所示的相同。   
![这里写图片描述][20150523024543701]

通过比较图5与图4，可以清楚地看到特征值表示沿特征向量方向数据的方差，而协方差矩阵的方差分量表示沿轴的传播。如果没有协方差，则这两个值是相等的。

## 协方差矩阵作为线性变换 ##

现在，让我们忘了协方差矩阵。图3的实例可以简单地认为是图6的一个线性变换实例：   
![这里写图片描述][20150523024820972]

图6所示的数据是D，则图3所示的每个实例可以通过线性变换D得到：![这里写图片描述][20150523025049967]

其中T是变换矩阵，包括一个旋转矩阵R和缩放矩阵S：   
![这里写图片描述][20150523025322514]

这些矩阵定义如下：   
![这里写图片描述][20150523025225564]   
其中![这里写图片描述][20150523025401052]是旋转角度。

![这里写图片描述][20150523025341349]   
![这里写图片描述][20150523025632943]分别是x方向和y方向的比例因子。

在下面的段落中，我们将讨论协方差矩阵![这里写图片描述][20150523025912734]与线性变换矩阵T= RS之间的关系。

让我们先从未缩放（缩放相当于1）和未旋转的数据开始。在统计中，这往往为“白数据’，因为它的样本是从标准正态分布引出的，因此对应于白（不相关）噪声：   
![这里写图片描述][20150523025801105]

这个“白色”数据的协方差矩阵等于单位矩阵，使得方差和标准差等于1，协方差等于零：   
![这里写图片描述][20150523025916438]

现在让我们用因子4在x方向缩放数据：   
![这里写图片描述][20150523030154210]

数据D’现在如下：   
![这里写图片描述][20150523030252570]

D’的协方差![这里写图片描述][20150523030337202]现在是：   
![这里写图片描述][20150523030418854]

D’的协方差![这里写图片描述][20150523030337202]与线性变换矩阵T有关系，D=TD，其中:   
![这里写图片描述][20150523030643872]

然而，虽然数据在x和y方向上缩放时等式（12）成立，但是应用旋转是否依然成立呢？为了调查一般情况下线性变换矩阵T和协方差矩阵![这里写图片描述][20150523030337202]之间的关系，我们试图分解协方差矩阵为旋转和缩放矩阵的乘积。

正如我们前面所看到的，我们可以用特征向量和特征值表示协方差矩阵：   
![这里写图片描述][20150523031549483]   
![这里写图片描述][20150523201002279]

等式（13）保存矩阵Σ的每个特征向量和特征值。在2D情况下，我们得到两个特征值和两个特征值。由公式（13）定义的两个等式可以有效地用矩阵符号来表示：   
![这里写图片描述][20150523201447706]   
其中V是矩阵，它的列是Σ的特征向量，L是对角矩阵，其非零元素对应特征值。

这意味着我们可以将协方差矩阵表示为特征向量和特征值的函数：   
![这里写图片描述][20150523201650899]

方程（15）就是所谓协方差矩阵特征值分解，并可以使用奇异值分解[算法][Link 2]来获得。而特征向量表示数据最大方差的方向，特征值表示那些方向方差的幅度。换言之，V表示旋转矩阵，而![这里写图片描述][20150523202059579]表示一个缩放矩阵。协方差矩阵可以进一步分解为：   
![这里写图片描述][20150523201857391]   
![这里写图片描述][20150523202004502]

在等式（6）中，我们定义了一个线性变换T= RS。由于S是对角缩放矩阵，所以S=ST。此外，由于R为正交矩阵，R-1=RT。因此，![这里写图片描述][20150523202406437]协方差矩阵可以写为：   
![这里写图片描述][20150523202355819]

换言之，如果我们应用由T=RS定义的线性变换到图7所示的原始白数据，我们得到了旋转和缩放的数据D’及协方差矩阵![这里写图片描述][20150523203047841]。这示于图10：   
![这里写图片描述][20150523202852276]  
图10的彩色箭头表示特征向量。最大特征向量，即与最大特征值对应的特征向量，总是指向数据最大方差的方向，并由此确定其方位。次特征向量总是正交于最大特征向量，因为旋转矩阵的正交性。

总结   
在本文中，我们表明观察到数据的协方差矩阵与白色不相关数据的线性变换有直接的关系。此线性变换完全由数据的特征向量和特征值确定。而特征向量表示旋转矩阵，特征值对应于每个维度上缩放因子的平方。

[20150523013736894]: /images/20220713/70ddc20b18c54750a220825bf0548b20.png
[20150523013917027]: /images/20220713/763e6881cf704c87a647f3979ea88e65.png
[20150523014142950]: /images/20220713/02c9ef01815242728c285a52298ed759.png
[20150523015110077]: https://img-blog.csdn.net/20150523015110077
[20150523015132026]: /images/20220713/1a69b776587d40e8b1c89a55c8844cc1.png
[20150523015153507]: /images/20220713/d72e2f2ae55847bead55a467b091a441.png
[20150523015242351]: /images/20220713/b162b7ec11134e9eb81ae33c3d3934f9.png
[20150523015344002]: /images/20220713/00b32fe13961402898a42fbdd86e4d31.png
[20150523020837615]: /images/20220713/e1cd7bd334b2492fb3fe5921282ca926.png
[20150523015824596]: /images/20220713/31b3469ca1d24b9c95e75b86eac397c7.png
[20150523022152998]: https://img-blog.csdn.net/20150523022152998
[20150523022351979]: https://img-blog.csdn.net/20150523022351979
[20150523022501357]: https://img-blog.csdn.net/20150523022501357
[20150523022300156]: https://img-blog.csdn.net/20150523022300156
[20150523022524332]: https://img-blog.csdn.net/20150523022524332
[20150523022313701]: https://img-blog.csdn.net/20150523022313701
[20150523022639652]: https://img-blog.csdn.net/20150523022639652
[20150523022211266]: https://img-blog.csdn.net/20150523022211266
[20150523022358535]: https://img-blog.csdn.net/20150523022358535
[20150523022702643]: /images/20220713/6783619b30524074a9472856650934a9.png
[Link 1]: http://lib.csdn.net/base/hadoop
[.NET]: http://lib.csdn.net/base/dotnet
[20150523023546860]: /images/20220713/cd75f93c36734b47bdec46ca9896470f.png
[20150523023629620]: /images/20220713/298ff1898d4c4be190d18f469929dd62.png
[20150523023728563]: /images/20220713/660c6e092bf6459694df960bf73bf87e.png
[20150523024543701]: /images/20220713/042a95be558f4e8cadbfd77e4ebdda7a.png
[20150523024820972]: /images/20220713/944d30fc6c1f4e048ee91f6a7044d132.png
[20150523025049967]: /images/20220713/cdd5fbb7c7364096989efc15cbf52df1.png
[20150523025322514]: /images/20220713/338c0a2738c440e5acbb736e319c21c0.png
[20150523025225564]: /images/20220713/5727f1574c75444eb9eb8f1abe36f3c3.png
[20150523025401052]: /images/20220713/0300a2cdded14a5e8d98cd495ab622ac.png
[20150523025341349]: /images/20220713/79506f1315ae47b28e397caf9f824c93.png
[20150523025632943]: /images/20220713/458a043db11a408ca81ff3388593acac.png
[20150523025912734]: /images/20220713/4f815b17fda144eca775cd92dc0ae1e7.png
[20150523025801105]: /images/20220713/755d063bb62b4b2d826ae15b90822dad.png
[20150523025916438]: /images/20220713/71b8c2a21e5742ffadae092c8d446fb3.png
[20150523030154210]: /images/20220713/04d417f6997f455da090e961b8464874.png
[20150523030252570]: /images/20220713/c8221729ef68436bb6e8fc808a87aa47.png
[20150523030337202]: /images/20220713/184f1cbfed164368b6d962c4ebf4b3ca.png
[20150523030418854]: /images/20220713/c700b87c7c65444db18b74888bc6c108.png
[20150523030643872]: /images/20220713/ad684b3854b5480985746d8be382cec8.png
[20150523031549483]: /images/20220713/de39b22e49c34e839d068213229ba506.png
[20150523201002279]: /images/20220713/3b711b8290b049f79c3864a5b249b5f7.png
[20150523201447706]: /images/20220713/2a278caa408840f48d758a8d4e012855.png
[20150523201650899]: /images/20220713/a6bef112dea74d989b7fd2bbd7894bfe.png
[Link 2]: http://lib.csdn.net/base/datastructure
[20150523202059579]: /images/20220713/07c1b78adcc3446088578c17519c59e0.png
[20150523201857391]: /images/20220713/3293c4d814d246ed857d1ab16a796b35.png
[20150523202004502]: /images/20220713/db62b3ee264b403eadb9a2ff88830bcf.png
[20150523202406437]: /images/20220713/b27d7fec998b4f79baad844f7f7d6e4a.png
[20150523202355819]: /images/20220713/495259fa7f3d4f1684a402ac0add6111.png
[20150523203047841]: /images/20220713/978ff0162ce947ceafcd6750d9ff6de2.png
[20150523202852276]: /images/20220713/7ee45c3a572a41548ba1cadd06bbe39b.png