每日一贴,今天的内容关键字为特征方向
主成份分析（PCA） 在很多教程中做了分析，但是为何通过协方差矩阵的特征值分解可以得到数据的主成份？协方差矩阵和特征值为何如此奇异，我却一直没弄清。今天终究把整个进程整理出来，方便自己学习，也和大家交流。
- 提出背景
以二维特征为例，两个特征之间可能存在线性关系的（例如这两个特征分别是运动的时速和秒速度），这样就造成了第二维信息是冗余的。PCA的目标是为了发现这类特征之间的线性关系，检测出这些线性关系，并且去除这线性关系。
还是以二维特征为例，如下图。特征之间可能不存在完全的线性关系，可能只是强的正相关。如果把x-y坐标分解成u1-u2坐标，而u1轴线上反响了特征的主要变更（intrinsic），而u2的特征变更较小，其实可以完全理解为一些噪声的扰动而不去斟酌它。PCA的任务就是找到u1和u2。
- 预处理：
将每一维特征的均值中心化，方差归一化。
- PCA的数学目标：
特征的主方向，就是特征幅度变更最大的方向（“major axis of variation”）。这一点理解很主要。从背面理解，幅度变更最小的方向就是没有变更，或者非常非常小的变更（可以忽略的变更），相对来说可利用价值最小，最可以忽略。而为了找到特征变更最大的方向，假设单位方向矢量为u，则特征点x在u方向的投影点x’距离原点的距离为d=x。所有的样本点都在一个方向投影后，他们就都在统一条直线上了。而要比拟它们之间变更的水平，只要比拟d的方差就行。方差最大的u对应的方向就是我们要寻觅的主方向。因此，我们的目标函数就成为了：
（1）
其中x的上标i表示数据集中的第i个样本，m表示数据集中的样本总数。（因为x已中心化了，所以xu的均值也为0，因此xu的平方只和就是方差。）
括号中的一项非常熟悉，就是协方差矩阵Σ！终究晓得协方差矩阵是怎么来的了。再看一看上面的式子，协方差矩阵与投影的方向无关，之于数据集中的样本有关，因此协方差矩阵完全决定了数据的分布及变更情况（请和自相关矩阵区分）。
目标函数如下：
（2）
用拉格朗日乘数法求解上面的最大化问题，很容易得到：
（3）
瞥见没？！u就是Σ的特征向量，λ就是特征值。我们再把（3）代入（2），目标函数就变成了
（4）
可见，可以通过协方差矩阵的迹衡量方差的巨细。最大的特征值λ（以及对应的特征向量u）决定了数据变更最大的方向。u就是这个单位方向。因此PCA的求解进程就是对协方差矩阵进行特征值分解，并找到最大的几个特征值的进程。
再前面的进程就是由最大的k个特征值对应的特征向量组成一组新的基（basis），让原特征对这个新的基投影得到降维后的新特征。这个进程很多教程都分析的很清晰了，我就不描述了。
不过我还是想补充一下矩阵及其特征值的意思。矩阵应理解为一种空间变换（从一个空间到另一个空间的变换）。矩阵M是m×n维的，如果m=n则变换后空间维数不变，如果n<m则是降维了。我们斟酌m=n的方阵。对M进行特征值分解，得到一组新的基，这组基可以理解为m维空间下的一组新的基坐标（引入它是因为这组基坐标是一种更有效的表达Rm维空间的方法），而对应的特征值，就是投影在基坐标下相应维的响应水平。如果特征值为0，代表在这一维下没有相应，不管是多少，乘过来都是0。特征值如果很接近0，那么任何数在这一维下投影后都市变小很多。所以从降维的角度斟酌，可以将这一维忽略（也就是PCA中保留前k个特征值的目标）。奇异值分解也是一样的情理。
PCA实际上是最简单的降维方法之一了，很明显的劣势是它仅去除数据之间的线性相关性。对线性的改善常常通过kernel技术拓展到非线性的应用上。另外，PCA的这类降维不一定有助于分类，用于分类的降维方法之一就是LDA。从另一方面说，PCA是一种线性投影，保留了数据与数据之间的欧式距离，即本来欧式距离大的两点在降维后的空间中距离也应大（这样才好保障方差大）。而事实上数据有可能呈现某种流型结构，用PCA降维后数据将不能保持原有的流型结构。在这一方面常用的非线性降维方法是Locally linear embedding和Laplacian Eigenmaps，如下图所示：
PCA的另一种推导方法是最小化投影后的损失（把降维理解为压缩，压缩后还原所得到的误差最小）,在这篇文章中也有具体分析，我也未几说了。
写到这里，才发现我啥也没说，都是提供了各种文献的链接
。
--------------------
jiang1st2010
原文地址：http://blog.csdn.net/jiang1st2010/article/details/8935219
文章结束给大家分享下程序员的一些笑话语录：
《诺基亚投资手机浏览器UCWEB，资金不详或控股》杯具了，好不容易养大的闺女嫁外国。（心疼是你养的吗？中国创业型公司创业初期哪个从国有银行贷到过钱？）