function?main（）
%*************主成份分析************
%讀入文件數(shù)據(jù)
X=load（‘\\data.txt‘）;

%==========方法1:求標(biāo)準(zhǔn)化后的協(xié)差矩陣再求特征根和特征向量=================
%標(biāo)準(zhǔn)化處理
[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（