matlab主成分分析法(參考版)
2024-08-16 01:20本頁面
【正文】 39。中國(guó)大陸35個(gè)大城市某年的10項(xiàng)社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)指標(biāo)數(shù)據(jù)見下表。)v1=cwstd(vector)result=cwfac(v1)。fprintf(39。%g39。r39。)score=[sco,csum,j] %得分矩陣：sco為各主成分得分；csum為綜合得分；j為排序結(jié)果 %function print=cwprint(filename,a,b)。fprintf(39。[newcsum,i]=sort(1*csum)。sco=vector1*vector2。)for p=1:length(newi) for q=1:length(y) result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p))。fprintf(39。主成分?jǐn)?shù)：%g\n\n39。 if sumrate break。for k=length(y):1:1 sumrate=sumrate+rate(k)。)newrate=newy/sum(newy)sumrate=0。fprintf(39。%g\n39。)for z=1:length(y) newy(z)=y(length(y)+1z)。 %對(duì)特征根進(jìn)行排序，y為排序結(jié)果，i為索引fprintf(39。)[vec,val]=eig(std) %求特征值(val)及特征向量(vec)newval=diag(val) 。)std=CORRCOEF(vector) %計(jì)算相關(guān)系數(shù)矩陣fprintf(39。fprintf(39。 %矩陣大小,a為行數(shù),b為列數(shù)for i=1:a for j=1:b std(i,j)= vector(i,j)/cwsum(j)。 主函數(shù) 子函數(shù)——用總和標(biāo)準(zhǔn)化法標(biāo)準(zhǔn)化矩陣——計(jì)算相關(guān)系數(shù)矩陣；計(jì)算特征值和特征向量；對(duì)主成分進(jìn)行排序；計(jì)算各特征值貢獻(xiàn)率；挑選主成分（累計(jì)貢獻(xiàn)率大于85%），輸出主成分個(gè)數(shù)；計(jì)算主成分載荷——計(jì)算各主成分得分、綜合得分并排序——讀入數(shù)據(jù)文件；調(diào)用以上三個(gè)函數(shù)并輸出結(jié)果 %,用總和標(biāo)準(zhǔn)化法標(biāo)準(zhǔn)化矩陣function std=cwstd(vector)cwsum=sum(vector,1)。④ 計(jì)算主成分載荷其計(jì)算公式為 （3）得到各主成分的載荷以后，還可以按照（，得到各主成分的得分 （4）在軟件Matlab中實(shí)現(xiàn)主成分分析可以采取兩種方式實(shí)現(xiàn)：一是通過編程來實(shí)現(xiàn)；二是直接調(diào)用Matlab種自帶程序?qū)崿F(xiàn)。這里要求=1，即，其中表示向量的第j個(gè)分量。 ① 計(jì)算相關(guān)系數(shù)矩陣 （1）在（，rij（i，j=1，2，…，p）為原變量的xi與xj之間的相關(guān)系數(shù)，其計(jì)算公式為 （2）因?yàn)镽是實(shí)對(duì)稱矩陣（即rij=rji），所以只需計(jì)算上三角元素或下三角元素即可。Matlab 語言在各國(guó)高校與研究單位起著重大的作用。它起源于矩陣運(yùn)算，并已經(jīng)發(fā)展成一種高度集成的計(jì)算機(jī)語言。167。 Matlab語言是當(dāng)今國(guó)際上科學(xué)界 (尤其是自動(dòng)控制領(lǐng)域) 最具影響力、也是最有活力的軟件。它提供了強(qiáng)大的科學(xué)運(yùn)算、靈活的程序設(shè)計(jì)流程、高質(zhì)量的圖形可視化與界面設(shè)計(jì)、與其他程序和語言的便捷接口
數(shù)學(xué)相關(guān)推薦
鄂ICP備17016276號(hào)-1