资源描述:
《主成分分析报告法matlab地实现》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、实用标准文案MATLAB结课作业指导老师:张肃班级:信管121姓名:桂亚东学号:201200654118精彩文档实用标准文案利用Matlab编程实现主成分分析概述Matlab语言是当今国际上科学界(尤其是自动控制领域)最具影响力、也是最有活力的软件。它起源于矩阵运算,并已经发展成一种高度集成的计算机语言。它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、与其他程序和语言的便捷接口的功能。Matlab语言在各国高校与研究单位起着重大的作用。主成分分析是把原来多个变量划为少数几个综合指标的一种统计分析方法,从数学角度来看,这是一种降维处理技术。1.1主成分分析
2、计算步骤①计算相关系数矩阵(1)在(3.5.3)式中,rij(i,j=1,2,…,p)为原变量的xi与xj之间的相关系数,其计算公式为(2)因为R是实对称矩阵(即rij=rji),所以只需计算上三角元素或下三角元素即可。②计算特征值与特征向量精彩文档实用标准文案首先解特征方程,通常用雅可比法(Jacobi)求出特征值,并使其按大小顺序排列,即;然后分别求出对应于特征值的特征向量。这里要求=1,即,其中表示向量的第j个分量。③计算主成分贡献率及累计贡献率主成分的贡献率为累计贡献率为一般取累计贡献率达85—95%的特征值所对应的第一、第二,…,第m(m≤p)个主成分。④计算主成分载荷
3、其计算公式为(3)得到各主成分的载荷以后,还可以按照(3.5.2)式进一步计算,得到各主成分的得分精彩文档实用标准文案(4)2.函数作用Cwstd.m——用总和标准化法标准化矩阵Cwfac.m——计算相关系数矩阵;计算特征值和特征向量;对主成分进行排序;计算各特征值贡献率;挑选主成分(累计贡献率大于85%),输出主成分个数;计算主成分载荷Cwscore.m——计算各主成分得分、综合得分并排序Cwprint.m——读入数据文件;调用以上三个函数并输出结果3.源程序3.1cwstd.m总和标准化法标准化矩阵%cwstd.m,用总和标准化法标准化矩阵functionstd=cwstd(
4、vector)cwsum=sum(vector,1);%对列求和[a,b]=size(vector);%矩阵大小,a为行数,b为列数fori=1:aforj=1:bstd(i,j)=vector(i,j)/cwsum(j);endend3.2cwfac.m计算相关系数矩阵精彩文档实用标准文案%cwfac.mfunctionresult=cwfac(vector);fprintf('相关系数矩阵:')std=CORRCOEF(vector)%计算相关系数矩阵fprintf('特征向量(vec)及特征值(val):')[vec,val]=eig(std)%求特征值(val)及
5、特征向量(vec)newval=diag(val);[y,i]=sort(newval);%对特征根进行排序,y为排序结果,i为索引fprintf('特征根排序:')forz=1:length(y)newy(z)=y(length(y)+1-z);endfprintf('%g',newy)rate=y/sum(y);fprintf('贡献率:')newrate=newy/sum(newy)sumrate=0;newi=[];fork=length(y):-1:1sumrate=sumrate+rate(k);newi(length(y)+1-k)=i(k);ifs
6、umrate>0.85break;endend%记下累积贡献率大85%的特征值的序号放入newi中fprintf('主成分数:%g',length(newi));fprintf('主成分载荷:')forp=1:length(newi)forq=1:length(y)result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p));endend%计算载荷disp(result)3.3cwscore.m%cwscore.m,计算得分functionscore=cwscore(vector1,vector2);sco=vector1*vect
7、or2;csum=sum(sco,2);[newcsum,i]=sort(-1*csum);[newi,j]=sort(i);精彩文档实用标准文案fprintf('计算得分:')score=[sco,csum,j]%得分矩阵:sco为各主成分得分;csum为综合得分;j为排序结果3.4cwprint.m%cwprint.mfunctionprint=cwprint(filename,a,b);%filename为文本文件文件名,a为矩阵行数(样本数),b为矩阵列数(变量指标数