模式识别pca的设计

function [out1 out2 out3 out4]=pca（Xn）%out1相关系数矩阵特征值，
                      %out2相关系数矩阵，out3各主成分，
                      %out4累计贡献率，X数据矩阵，n前n个主成分；
clc;
X=xlsread（‘zb0506.xls‘）; 
m=size（X）;
a=sum（X）/m（1）;
for i=1:m（1）;
for j=i:m（2）;
R（ij）=sum（（X（:i）-a（i））.*（X（:j）-a（j）））./（sum（（X（:i）-a（i））.^2）*sum（（X（:j）-a（j））.^2））.^0.5;
end
end
for i=1:m（2）;
for j=1:i-1;
R