12.1 主成分分析原理(PCA)

主成分分析步骤如下:

设有nm维数据

1. 将原始数据按列组成nm列矩阵X

2. 将矩阵 X 的每一行进行零均值化;

3. 求出协方差矩阵C=\frac{1}{m}XX^{^{T}}

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

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

6. Y=PX即为降维到k维后的数据。

        如果对线性代数的实对成矩阵的相似对角化熟悉的人可以很好的了解上面的步骤,根据线性代数的求解步骤如下(也可以是代码步骤,python的np有求解方法):

1. 将原始数据按列组成nm列矩阵X

2. 将矩阵 X 的每一行进行零均值化;

    这里的每一行零均值化就是每一行的元素减去该行的均值,均值可以用np.mean()求解。

3.求出协方差矩阵C=\frac{1}{m}XX^{^{T}};这个矩阵实际上就是实对称矩阵。

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

       这里其实就是对实对称矩阵相似对角化,在相似对角化的过程中会得到特征值和对应的特征向量,得到的对角矩阵的对角线上的元素就是特征值,对角线以外的元素都为0。(当然,要对这组特征向量进行正交化,因为我们最终要的是正交向量,如果在这里不求,也可以在最后一步求k个向量的正交向量,计算量肯定比现在求小,因为现在要求n个向量,而最后只需要求k个)。

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

6.这里就是对特征值做从大到小的排序,可以用np.sorted()函数排序。

7.然后根据排序取前k行的向量组成一个新的矩阵。

你可能感兴趣的:(人工智能,算法)