文档介绍：莇§(尤其是自动控制领域)最具影响力、也是最有活力的软件。它起源于矩阵运算,并已经发展成一种高度集成的计算机语言。它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、与其他程序和语言的便捷接口的功能。Matlab语言在各国高校与研究单位起着重大的作用。主成分分析是把原来多个变量划为少数几个综合指标的一种统计分析方法,从数学角度来看,这是一种降维处理技术。①计算相关系数矩阵蚅(1)莁在()式中,rij(i,j=1,2,…,p)为原变量的xi与xj之间的相关系数,其计算公式为薀(2)芅因为R是实对称矩阵(即rij=rji),所以只需计算上三角元素或下三角元素即可。蒂②计算特征值与特征向量蒀首先解特征方程,通常用雅可比法(Jacobi)求出特征值,并使其按大小顺序排列,即;然后分别求出对应于特征值的特征向量。这里要求=1,即,其中表示向量的第j个分量。罿③计算主成分贡献率及累计贡献率羅主成分的贡献率为蒄袂累计贡献率为荿螆一般取累计贡献率达85—95%的特征值所对应的第一、第二,…,第m(m≤p)个主成分。薅④计算主成分载荷羀其计算公式为螈(3)蒆得到各主成分的载荷以后,还可以按照()式进一步计算,得到各主成分的得分莂(4):一是通过编程来实现;二是直接调用Matlab种自带程序实现。下面主要主要介绍利用Matlab的矩阵计算功能编程实现主成分分析。————计算相关系数矩阵;计算特征值和特征向量;对主成分进行排序;计算各特征值贡献率;挑选主成分(累计贡献率大于85%),输出主成分个数;——计算各主成分得分、——读入数据文件;%,用总和标准化法标准化矩阵蒈functionstd=cwstd(vector)莅cwsum=sum(vector,1);%对列求和肃[a,b]=size(vector);%矩阵大小,a为行数,b为列数膄fori=1:a袀forj=1:b聿std(i,j)=vector(i,j)/cwsum(j);%=cwfac(vector);莁fprintf('相关系数矩阵:\n')袈std=CORRCOEF(vector)%计算相关系数矩阵芅fprintf('特征向量(vec)及特征值(val):\n')肄[vec,val]=eig(std)%求特征值(val)及特征向量(vec)葿newval=diag(val);莇[y,i]=sort(newval);%对特征根进行排序,y为排序结果,i为索引羅fprintf('特征根排序:\n')袁forz=1:length(y)袂newy(z)=y(length(y)+1-z);螆end螅fprintf('%g\n',newy)羃rate=y/sum(y);羀fprintf('\n贡献率:\n')膆newrate=newy/sum(newy)蒆sumrate=0;羄newi=[];肈fork=length(y):-1:1衿sumrate=sumrate+rate(k);芆newi(length(y)+1-k)=i(k);螁ifsumrate>;蒁end艿end%记下累积贡献率大85%的特征值的序号放入newi中羇fprintf('主成分数:%g\n\n',length(newi));袃fprintf('主成分载荷:\n')蕿forp=1:length(newi)螈forq=1:length(y)蒃result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p));袄end羂end%计算载荷***disp(result)%,计算得分薇functionscore=cwscore(vector1,vector2);羄sco=vector1*vector2;螃csum=sum(sco,2);膈[newcsum,i]=sort(-1*csum);羆[newi,j]=sort(i);蚄fprintf('计算得分:\n')袄score=[sco,csum