matlab主成分分析法.docxVIP

（完整word fig）matlab主成分分析法 编辑整理: 尊敬的读者朋友们： 这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们 对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（（完整word版）matlab主 成分分析法）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈， 这将是我们进步的源泉，前进的动力。 本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步，以 下为（完整word版）matlab主成分分析法的全部内容。 §10?利用Matlab编程实现主成分分析 1o概述 Mat lab语言是当今国际上科学界（尤其是自动控制领域）最具影响力、也是 最有活力的软件。它起源于矩阵运算，并已经发展成一种高度集成的计算机语言。 它提供了强大的科学运算、灵活的程序设计流程、高质量的图形可视化与界面设计、 与其他程序和语言的便捷接口的功能。Mat lab语言在各国高校与研究单位起着重 大的作用。主成分分析是把原来多个变量划为少数几个综合指标的一种统计分析方 法，从数学角度来看，这是一种降维处理技术。 主成分分析计算步骤 ①计算相关系数矩阵 (1) 在（，5（2, 2,Q为原变量的灯与灯之间的相关系数，其计算公式 为 n 习不厂召）鴿-初 不厂*虫％ ―耳）‘ V =1 *=1 因为R是实对称矩阵（即eg ,所以只需计算上三角元素或下三角元素即可。 ②计算特征值与特征向量 功能编程实现主成分分析。 主函数Cwscore.m2. 1程序结构 主函数 Cwscore.m Cwpr i nt。m 子函数 Cwstd o Cwfaco m 2.2歯数作用 Cwstd. m—用总和标准化法标准化矩阵 Cwfac. m—计算相关系数矩阵；计算特征值和特征向量；对主成分进行排序；计 算各特征值贡献率；挑选主成分（累计贡献率大于85%）,输出主成分个数；计算 主成分载荷 Cwscore. m 计算各主成分得分、综合得分并排序 Cwpr i nt. l读入数据文件；调用以上三个函数并输出结果 3.源程序 3o 1 cwstdo m总和标准化法标准化矩阵 %cwstd. m,用总和标准化法标准化矩阵 function std二cwstd (vector) cwsum=sum(vector, 1) ； %对列求和 [a, b]二size (vector) ； %矩阵大小，a为行数,b为列数 for i=1: a for j=1: b std (i, j)= vector (i, j) /cwsum(j)； end end 3o 2 cwfaco m计算相关系数矩阵 %cwfac. m function resu11二cwfac (vector); fprintf C相关系数矩阵：\n) std二CORRCOEF (vector) %计算相关系数矩阵 fprintfC特征向量(vec)及特征值(val) : \nr ) [vec, val]=eig (std) %求特征值(va I)及特征向量(vec) newval=diag (val)； [y. i]=sort(newval) ； %对特征根进行排序，y为排序结果，i为索引 fprintf (特征根排序：\n*) for z=1:length (y) newy (z) =y (I ength (y) +1 —z); end fpr intf (%g\rT , newy) rate二y/sum(y); fpr intf ( \n 贡献率:\n) ne wra te 二 ne wy/sum (newy) sumrate二0; newi = []; for k=Iength(y)：-1:1 sumrate=sumrate+rate(k)； newi (length (y)+1—k) =i(k); if sumrate〉0.85 break; end end %记下累积贡献率大85%的特征值的序号放入newi中 fpr intf (主成分数:%g\n\n, Iength (newi))； fprintf (主成分载荷:\rT) for p=1: length(newi) for q=1:length (y) result (q, p)=sqrt (newvaI (newi (p) ))*vec (q, newi (p))； end end %计算载荷 d i sp(resuIt) 3.3 cwscoreo m %cwscoreo m,计算得分 function score二cwscore (vectorh vector2); sco二vector1*vector2; csum=sum (sco,2)； [ne