matlab主成分分析法-文庫吧資料
2024-08-18 01:20本頁面
【正文】 的功能。主成分分析是把原來多個(gè)變量劃為少數(shù)幾個(gè)綜合指標(biāo)的一種統(tǒng)計(jì)分析方法，從數(shù)學(xué)角度來看，這是一種降維處理技術(shù)。② 計(jì)算特征值與特征向量首先解特征方程，通常用雅可比法（Jacobi）求出特征值，并使其按大小順序排列，即；然后分別求出對應(yīng)于特征值的特征向量。③ 計(jì)算主成分貢獻(xiàn)率及累計(jì)貢獻(xiàn)率主成分的貢獻(xiàn)率為累計(jì)貢獻(xiàn)率為一般取累計(jì)貢獻(xiàn)率達(dá)85—95%的特征值所對應(yīng)的第一、第二，…，第m（m≤p）個(gè)主成分。下面主要主要介紹利用Matlab的矩陣計(jì)算功能編程實(shí)現(xiàn)主成分分析。 %對列求和[a,b]=size(vector)。 endend %function result=cwfac(vector)。相關(guān)系數(shù)矩陣:\n39。特征向量(vec)及特征值(val)：\n39。[y,i]=sort(newval) 。特征根排序：\n39。endfprintf(39。,newy)rate=y/sum(y)。\n貢獻(xiàn)率：\n39。newi=[]。 newi(length(y)+1k)=i(k)。 end end %記下累積貢獻(xiàn)率大85%的特征值的序號放入newi中fprintf(39。,length(newi))。主成分載荷：\n39。 endend %計(jì)算載荷disp(result) %,計(jì)算得分function score=cwscore(vector1,vector2)。csum=sum(sco,2)。[newi,j]=sort(i)。計(jì)算得分：\n39。 %filename為文本文件文件名，a為矩陣行數(shù)(樣本數(shù))，b為矩陣列數(shù)(變量指標(biāo)數(shù))fid=fopen(filename,39。)vector=fscanf(fid,39。,[a b])。標(biāo)準(zhǔn)化結(jié)果如下：\n39。cwscore(v1,result)。城 市名 稱年底總?cè)丝?萬人) 非農(nóng)業(yè)人口比(%) 農(nóng) 業(yè)總產(chǎn)值(萬元)工業(yè)總產(chǎn)值(萬元)客運(yùn)總量(萬人)貨運(yùn)總量(萬噸)地方財(cái)政預(yù)算內(nèi)收入(萬元)城鄉(xiāng)居民年底儲蓄余額(萬元)在崗職工人數(shù)(萬人)在崗職工工資總額(萬元)北 京1 81 843 42719 999 70620 32345 5622 790 86326 806 6465 773 301天 津 91 501 13622 645 5023 25926 3171 128 07311 301 9312 254 343石 家 莊 22 918 6806 885 7682 9291 911352 3487 095 875758 877太 原 3236 0382 737 7501 93711 895203 2773 943 100654 023呼和浩特 2365 343816 4522 3512 623105 7831 396 588309 337沈 陽 91 295 4185 826 7337 78215 412567 9199 016 9981 152 811大 連 61 879 7398 426 38510 78019 187709 2277 556 796965 922長 春 81 853 2105 966 3434 8109 532357 0964 803 744884 447
數(shù)學(xué)相關(guān)推薦
鄂ICP備17016276號-1