function?[out1?out2?out3?out4]=pca（Xn）%out1相關系數矩陣特征值，out2相關系數矩陣，out3各主成分，out4累計貢獻率，X數據矩陣，n前n個主成分；
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（ij）=R（ji）?;?
end
end
[xd]=eig（