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

cic clear all A=xlsrea”这里换成你自己的数据根目录，例如：D:\资料库区大三上\HUAWEI\MATLAB\主 成分分析.xls?B3:ll7); %得到的数据矩阵的行数和列数 a=size(A,l); b=size(A,2); %数据的标准化处理:得到标准化后的矩阵SA for i=l:b SA(:J)=(A(:,i)-mean(A(:,i)))/std(A(:J)); end %计算系数矩阵:CM CM=corrcoef(SA); %计算CM的特征值和特征向量 [V,D]=eig(CM); %将特征值按降序排列到DS中 for j=l:b DS(jzl)=D(b+l-j,b+l-j)； end %计算贡献率 for i=l:b DS(i,2)=DS(i,l)/sum(DS(:,l));% 单个贡献率 DS(i,3)=sum(DS(l:i,l))/sum(DS(:,l));% 累计贡献率 end %假定主成分的信息保留率 T=0.9; for k=l:b if DS(kz3)=T com_num=k; break; end end %提取主成分的特征向最 for j=l:com_num PV(:,j)=V(:,b+l-j); end %计算主成分得分 new_score=SA* PV; for i=l:a total_score(i,l)=sum(new_score(i,:)); total_score(i/2)=i; end %强主成分得分与总分放到同一个矩阵中 result_report=[new_score/total_score]; %按总分降序排列 result_report=sortrows（result_report/-4）; %输出结果 dispC特征值、贡献率、累计贡献率：,） DS disp（信息保留率T对应的主成分数与特征向量：,） com_num PV disp（，主成分得分及排序（按第4列的总分进行降序排列，前3列为个各成分得分，第5列 为企业编号）,） result_report 可以用下方数据训练一下 企业序号 净利润率建 固定资产利涸率八 ,总产值利润率6 植售收入利润率小 产品成本利润率/? 物耗利润率d 人均利润（千元/人） 流动资金利?润率 1 40.4 24.7 7.2 6.1 8.3 8.7 2.442 20 2 25 12.7 11.2 11 12.9 20.2 3. 542 9.1 3 13.2 3.3 3.9 4.3 4.4 5.5 0. 578 3.6 4 22.3 6.7 5.6 3.7 6 7.4 0.176 7.3 5 34.3 11.8 7.1 7.1 8 8.9 1.726 27.5 6 35.6 12.5 16.4 16.7 22.8 29.3 3. 017 26.6 7 22 7.8 9.9 10.2 12.6 17.6 0. 847 10.6 8 48.4 13.4 10.9 9.9 10.9 13.9 1.772 17.8 9 40.6 19.1 19.8 19 29.7 39.6 2.449 35.8 10 24.8 8 9.8 8.9 11.9 16.2 0.789 13.7 11 12.5 9.7 4.2 4.2 4.6 6.5 0. 874 3.9 12 L8 0.6 0.7 0.7 0,8 1.1 0. 056 1 13 32.3 13.9 8.3 8.3 9.8 13.3 2.126 17.1 14 38.5 9.1 9.5 9.5 12.2 16.4 1.327 11.6 15 26.2 10.1 15.6 15.6 7.7 30.1 0.126 25.9 企业综合实力评价装 5 6 8 9 10 12 13 14 15 16