主成成分分析——MATLAB实现

主代码：

%下面为主要成分分析的程序
clear;clc;
load A;     %要分析的数据
A_aver=mean(A);
A_bzc=std(A,0,1);
A=(A-A_aver)./A_bzc;   %将矩阵A标准化
A_xfc=cov(A);   %求出A的协方差矩阵
[COEFF,latent,explained] = pcacov(A_xfc);
zcf=find(latent>1)';    %找到特征值大于1的成分以作为我们的主成分
COEFF=COEFF(:,zcf);
W=COEFF./sqrt(latent(1:size(COEFF,2)))';    %转换为主成分分析的系数
F=A*W;    %降维后的数据
E=(sum((F.*explained(1:size(F,2))')'))'/100;   %综合主成分得分

%下面程序是碎石图的绘制
hold on
plot(latent,'-o')
plot(ones(1,length(latent)),'--')
legend('latent', '参考线')
title('碎石图')

A矩阵对应的数据：

运行结果：