PCA matlab代码

function main（）
%*************主成份分析************
%
%see also http://www.matlabsky.com
%
%读入文件数据
X=load（‘data.txt‘）; 
%==========方法1:求标准化后的协差矩阵再求特征根和特征向量=================
%标准化处理
[pn]=size（X）;
for j=1:n
    mju（j）=mean（X（:j））;
    sigma（j）=sqrt（cov（X（:j）））;
end
for i=1:p
    for j=1:n
        Y（ij）=（X（ij）-mju（j））/sigma（j）; 
    end
end
sigmaY=cov（Y）;
%求X标准化的协差矩阵的特征根和特征向量
[Tlambda]=eig（sigmaY）;
disp（‘特征根（由小到大）:‘）;
disp（lambda）;
disp（‘特征向量:‘）;
disp（T）;
%方差贡献率;累计方差贡献率
Xsum=sum（sum（lambda2）1）;
for i=1:n
    fai（i）=lambda（ii）/Xsum;
end
for i=1:n
    psai（i）= sum（sum（lambda（1:i1:i）2）1）/Xsum;
end
disp（‘方差贡献率:‘）;
disp（fai）;
disp（‘累计方差贡献率:‘）;
disp（psai）;
%综合评价....略

%%
%
%
%============方法2:求X的相关系数矩阵再求特征根和特征向量================
%X的标准化的协方差矩阵就是X的