直观解释PCA

作者: 掰不开桃子的男人 | 来源:发表于2019-07-18 11:19 被阅读0次

    是干什么的?

    首先有一组数据蓝色点,PCA所谓的降维操作就是找到一个新的坐标系(旋转的两条直线式垂直的,我们可以用一组标准正交基来指示),然后减掉其中一些维度,使误差足够小。

    image.png

    基本思路

    假设我们有一个数据X_{n*m},其中n代表了特征的个数,m代表了样本数。首先对X的特征进行零均值化。

    协方差矩阵C_{n*n} = XX^T(这里应该除以m,不妨碍推导),C的对角线代表了特征自身的方差,而其他位置比如C_{i,j}代表了特征i和特征j之间的协方差。

    我们想要对原来的数据X做一个行变换A_{n*n},把原来的特征组合成新的特征,而新的特征之间没有关联,也就是新的矩阵Y_{n*m}的协方差矩阵里只有对角线元素,其他位置都是0.

    Want: YY^T=\Lambda\Lambda是一个对角阵

    带入Y=AX得到(AX)(AX)^T=\Lambda,展开得到AXX^TA^T=\Lambda,同志们这不就是对XX^T进行对角化吗?因为XX^T是实对称矩阵,一定可以对角化XX^T=Q{\Lambda}Q^T,原式中的A就是对角化后的特征向量矩阵Q^T

    其中特征值\Lambda大小代表了Y里面特征的方差大小,特征值越大,特征的方差越大,信息量越多。我们挑选主成分的话,只要挑选前面比较大的特征值对应的特征就可以了。

    计算PCA的步骤

    总结一下PCA的步骤:

    1. 将原始数据按列组成n行m列矩阵X

    2. 将X的每一行(代表一个属性字段)进行零均值化,即减去这一行的均值

    3. 求出协方差矩阵C=1/m * XX^T

    4. 求出协方差矩阵的特征值\Lambda及对应的特征向量Q

    5. 将特征向量按对应特征值大小从左到右按列排列成矩阵,取前k行组成矩阵Q

    6. Y=Q^TX即为降维到k维后的数据

    与SVD的关系

    以上就是先求协方差矩阵然后特征值分解的方法,那么可以不求协方差矩阵吗?可以,直接用SVD。

    SVD可以分解任意矩阵:X = U \Sigma V ^ T

    上面通过XX^T特征值分解其实得到的就是这里的U矩阵,通过成熟的SVD方法很容易求得这些矩阵。

    参考:

    Machine Learning — Singular Value Decomposition (SVD) & Principal Component Analysis (PCA)

    从PCA和SVD的关系拾遗

    PCA的数学原理

    相关文章

      网友评论

        本文标题:直观解释PCA

        本文链接:https://www.haomeiwen.com/subject/tiiklctx.html