主成分分析法matlab的實(shí)現(xiàn)-文庫(kù)吧在線文庫(kù)
【正文】 1219 251銀 川 9171 603661 2262 1061 19374 758814 103178 621烏魯木齊 478 5131 847 2412 6689 041254 8702 365 508517 622 cwprint(39。cwscore(v1,result)。 %filename為文本文件文件名，a為矩陣行數(shù)(樣本數(shù))，b為矩陣列數(shù)(變量指標(biāo)數(shù))fid=fopen(filename,39。 endend %計(jì)算載荷disp(result) %,計(jì)算得分function score=cwscore(vector1,vector2)。 newi(length(y)+1k)=i(k)。endfprintf(39。相關(guān)系數(shù)矩陣:\n39。② 計(jì)算特征值與特征向量首先解特征方程，通常用雅可比法（Jacobi）求出特征值，并使其按大小順序排列，即；然后分別求出對(duì)應(yīng)于特征值的特征向量。它起源于矩陣運(yùn)算，并已經(jīng)發(fā)展成一種高度集成的計(jì)算機(jī)語(yǔ)言。④ 計(jì)算主成分載荷其計(jì)算公式為 （3）得到各主成分的載荷以后，還可以按照（）式進(jìn)一步計(jì)算，得到各主成分的得分 （4）——用總和標(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)。)[vec,val]=eig(std) %求特征值(val)及特征向量(vec)newval=diag(val) 。fprintf(39。主成分?jǐn)?shù)：%g\n\n39。[newcsum,i]=sort(1*csum)。%g39。39。)v1=cwstd(vector)result=cwfac(v1)。)score=[sco,csum,j] %得分矩陣：sco為各主成分得分；csum為綜合得分；j為排序結(jié)果 %function print=cwprint(filename,a,b)。)for p=1:length(newi) for q=1:length(y) result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p))。for k=length(y):1:1 sumrate=sumrate+rate(k)。)for z=1:length(y) newy(z)=y(length(y)+1z)。fprintf(39。 ① 計(jì)算相關(guān)系數(shù)矩陣 （1）在（）式中，rij（i，j=1，2，…，p）為原變量的xi與xj之間的相關(guān)系數(shù)，其計(jì)算公式為 （2）因?yàn)镽是實(shí)對(duì)稱(chēng)矩陣（即rij=rji），所以只需計(jì)算上三角元素或下三角元素即可。它提供了強(qiáng)大的科學(xué)運(yùn)算、靈活的程序設(shè)計(jì)流程、高質(zhì)量的圖形可視化與界面設(shè)計(jì)、與其他程序和語(yǔ)言的便捷接口的功能。 %對(duì)列求和[a,b]=size(vector)。[y,i]=sort(newval) 。\n貢獻(xiàn)率：\n39。,length(newi))。
教學(xué)課件相關(guān)推薦
鄂ICP備17016276號(hào)-1