function [V,S,E]=princa(X)[m,n]=size(X); %计算矩阵的行m和列n%-------------第一步：标准化矩阵-----------------%mv=mean(X); %计算各变量的均值st=std(X); %计算各变量的标准差X=(X-repmat(mv,m,1))./repmat(st,m,1); %标准化矩阵X%-------------第二步：计算相关系数矩阵-----------------%% R1=X‘*X/(m-1); %方法一：协方差矩阵计算公式% R2=cov(X);     %方法二：协方差矩阵计算函数R=corrcoef(X); %方法三：相关系数矩阵函数%-------------第三步：计算特征向量和特征值-----------------%[V,D]=eig(R);       %计算矩阵R的特征向量矩阵V和特征值矩阵D,特征值由小到大V=(rot90(V))‘;      %将特征向量矩阵V从大到小排序D=rot90(rot90(D));  %将特征值矩阵由大到小排序E=diag(D);          %将特征值矩阵转换为特征值向量%-------------第四步：计算贡献率和累计贡献率-----------------%ratio=0; %累计贡献率for k=1:n    r=E(k)/sum(E);   %第k主成份贡献率    ratio=ratio+r;  %累计贡献率    if(ratio>=0.9)  %取累计贡献率大于等于90%的主成分        break;    endend%-------------第五步：计算得分-----------------%S=X*V;