xlRange = 'B1:X427';
[X,textdata] = xlsread('20180319.csv',xlRange);
XZ=zscore(X);
[COEFF,SCORE,latent,tsquare]=princomp(XZ)
explained = 100*latent/sum(latent);
[m,n]=size(X);
result1=cell(n+1,4);
result1(1,:)={'特征值','差值','贡献率','累积贡献率'};
result1(2:end,1)=num2cell(latent);
result1(2:end-1,2)=num2cell(-diff(latent));
result1(2:end,3:4)=num2cell([explained,cumsum(explained)])
varname=textdata(1,1:end)';
result2=cell(n+1,3);
result2(1,:)={'标准化变量','主成分Prin1','主成分Prin2'};
result2(2:end,1)=varname;
result2(2:end,2:end)=num2cell(COEFF(:,1:2))
