你会么MatLAB.docVIP

（2） 因为R是实对称矩阵（即rij=rji），所以只需计算上三角元素或下三角元素即可。 = 2 \* GB3 ② 计算特征值与特征向量 首先解特征方程，通常用雅可比法（Jacobi）求出特征值，并使其按大小顺序排列，即；然后分别求出对应于特征值的特征向量。这里要求=1，即，其中表示向量的第j个分量。 2.1程序结构 Cwprint.m 主函数 Cwprint.m 子函数 Cwfac.mCwscore.mCwstd.m Cwfac.m Cwscore.m Cwstd.m Cwprint.m——读入数据文件；调用以上三个函数并输出结果 3.源程序 3.1 cwstd.m %cwstd.m,用总和标准化法标准化矩阵 %cwfac.m function result=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) for z=1:length(y) newy(z)=y(length(y)+1-z); end newi(length(y)+1-k)=i(k); if sumrate0.85 break; end end %记下累积贡献率大85%的特征值的序号放入newi中 fprintf(主成分数：%g\n\n,length(newi)); fprintf(主成分载荷：\n) for p=1:length(newi) for q=1:length(y) result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p)); end end %计算载荷 disp(result) 3.3 cwscore.m %cwscore.m,计算得分 function score=cwscore(vector1,vector2); sco=vector1*vector2; csum=sum(sco,2); [newcsum,i]=sort(-1*csum); [newi,j]=sort(i); fprintf(计算得分：\n) score=[sco,csum,j] %得分矩阵：sco为各主成分得分；csum为综合得分；j为排序结果 3.4 cwprint.m %cwprint.m function print=cwprint(filename,a,b); %filename为文本文件文件名，a为矩阵行数(样本数)，b为矩阵列数(变量指标数) fid=fopen(filename,r) vector=fscanf(fid,%g,[a b]); fprintf(标准化结果如下：\n) v1=cwstd(vector) result=cwfac(v1); cwscore(v1,result); 4.程序测试 4.1原始数据 中国大陆35个大城市某年的10项社会经济统计指标数据见下表。 城 市 名 称 年底 总人口 (万人) 非农业 人口比(%) 农 业 总产值 (万元) 工业 总产值 (万元) 客运总量 (万人) 货运总量 (万吨) 地方财政 预算内收入(万元) 城乡居民年底储蓄余额 (万元) 在岗职工人数(万人) 在岗职工工资总额 (万元) 北 京 1 249.90 0.597 8 1 843 427 19 999 706 20 323 45 562 2 790 863 26 806 646 410.80 5 773 301 天 津 910.17 0.580 9 1 501 136 22 645 502 3 259 26 317 1 128 073 11 301 931 202.68 2 254 343 石 家 庄 875.40 0.233 2 2 918 680 6 885 768 2 929 1 911 352 348 7 095 875 95.60 758 877 太 原 299.92 0.656 3 236 038 2 737 750