function?[out1?out2?out3?out4]=pca（Xn）%out1相關(guān)系數(shù)矩陣特征值，
??????????????????????%out2相關(guān)系數(shù)矩陣，out3各主成分，
??????????????????????%out4累計(jì)貢獻(xiàn)率，X數(shù)據(jù)矩陣，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