主成分分析降维代码（直接调用版）.pdfVIP

clc

clearall

A=xlsread(这里换成你自己的数据根目录，例如：D:\资料库区\大三上\HUAWEI\MATLAB\主

成分分析.xls,B3:I17);

%得到的数据矩阵的行数和列数

a=size(A,1);

b=size(A,2);

%数据的标准化处理:得到标准化后的矩SA

fori=1:b

SA(:,i)=(A(:,i)-mean(A(:,i)))/std(A(:,i));

end

%计算系数矩:CM

CM=corrcoef(SA);

%计算CM的特征值和特征向量

[V,D]=eig(CM);

%将特征值按降序排列到DS中

forj=1:b

DS(j,1)=D(b+1-j,b+1-j);

end

%计算贡献率

fori=1:b

DS(i,2)=DS(i,1)/sum(DS(:,1));%单个贡献率

DS(i,3)=sum(DS(1:i,1))/sum(DS(:,1));%累计贡献率

end

%假定主成分的信息保留率

T=0.9;

fork=1:b

ifDS(k,3)=T

com_num=k;

break;

end

end

%提取主成分的特征向量

forj=1:com_n