主成分解析總結(jié)報(bào)告matlab程序

1、標(biāo)準(zhǔn)合用Matlab編程實(shí)現(xiàn)主成分解析.程序結(jié)構(gòu)及函數(shù)作用在軟件Matlab中實(shí)現(xiàn)主成分解析可以采用兩種方式實(shí)現(xiàn)：一是經(jīng)過編程來實(shí)現(xiàn)；二是直接調(diào)用Matlab種自帶程序?qū)崿F(xiàn)。下面主要主要介紹利用Matlab的矩陣計(jì)算功能編程實(shí)現(xiàn)主成分解析。程序結(jié)構(gòu)Cwprint.m主函數(shù)子函數(shù)Cwstd.mCwfac.mCwscore.m函數(shù)作用Cwstd.m用總和標(biāo)準(zhǔn)化法標(biāo)準(zhǔn)化矩陣Cwfac.m計(jì)算相關(guān)系數(shù)矩陣；計(jì)算特色值和特色向量；對(duì)主成分進(jìn)行排序；計(jì)算各特色值貢獻(xiàn)率；優(yōu)選主成分（累計(jì)貢獻(xiàn)率大于85%），輸出主成分個(gè)數(shù)；計(jì)算主成分載荷Cwscore.m計(jì)算各主成分得分、綜合得分并排序Cwprint.m讀入

2、數(shù)據(jù)文件；調(diào)用以上三個(gè)函數(shù)并輸出結(jié)果文案大全標(biāo)準(zhǔn)合用源程序3.1cwstd.m總和標(biāo)準(zhǔn)化法標(biāo)準(zhǔn)化矩陣%cwstd.m,用總和標(biāo)準(zhǔn)化法標(biāo)準(zhǔn)化矩陣functionstd=cwstd(vector)cwsum=sum(vector,1);%對(duì)列求和a,b=size(vector);%矩陣大小,a為行數(shù),b為列數(shù)fori=1:aforj=1:bstd(i,j)=vector(i,j)/cwsum(j);endend3.2cwfac.m計(jì)算相關(guān)系數(shù)矩陣%cwfac.mfunctionresult=cwfac(vector);fprintf(相關(guān)系數(shù)矩陣:n)std=CORRCOEF(vector)%計(jì)算

3、相關(guān)系數(shù)矩陣fprintf(特色向量(vec)及特色值(val)：n)vec,val=eig(std)%求特色值(val)及特色向量(vec)newval=diag(val);y,i=sort(newval);%對(duì)特色根進(jìn)行排序，y為排序結(jié)果，i為索引fprintf(特色根排序：n)forz=1:length(y)newy(z)=y(length(y)+1-z);endfprintf(%gn,newy)rate=y/sum(y);fprintf(n貢獻(xiàn)率：n)newrate=newy/sum(newy)sumrate=0;newi=;fork=length(y):-1:1sumrate=sum

4、rate+rate(k);newi(length(y)+1-k)=i(k);ifsumrate0.85break;endend%記下累積貢獻(xiàn)率大85%的特色值的序號(hào)放入newi中fprintf(主成分?jǐn)?shù)：%gnn,length(newi);fprintf(主成分載荷：n)forp=1:length(newi)文案大全標(biāo)準(zhǔn)合用forq=1:length(y)result(q,p)=sqrt(newval(newi(p)*vec(q,newi(p);endend%計(jì)算載荷disp(result)3.3cwscore.m%cwscore.m,計(jì)算得分functionscore=cwscore(vec

5、tor1,vector2);sco=vector1*vector2;csum=sum(sco,2);newcsum,i=sort(-1*csum);newi,j=sort(i);fprintf(計(jì)算得分：n)score=sco,csum,j%得分矩陣：sco為各主成分得分；csum為綜合得分；j為排序結(jié)果3.4cwprint.m%cwprint.mfunctionprint=cwprint(filename,a,b);%filename為文本文件文件名，a為矩陣行數(shù)(樣本數(shù))，b為矩陣列數(shù)(變量指標(biāo)數(shù))fid=fopen(filename,r)vector=fscanf(fid,%g,ab);

6、fprintf(標(biāo)準(zhǔn)化結(jié)果以下：n)v1=cwstd(vector)result=cwfac(v1);cwscore(v1,result);程序測試?yán)}4.1原始數(shù)據(jù)中國大陸35個(gè)大城市某年的10項(xiàng)社會(huì)經(jīng)濟(jì)統(tǒng)計(jì)指標(biāo)數(shù)據(jù)見下表。文案大全標(biāo)準(zhǔn)合用非客運(yùn)貨運(yùn)地方財(cái)年終農(nóng)業(yè)農(nóng)業(yè)工業(yè)政城總量總量市人口總產(chǎn)值總產(chǎn)值估量內(nèi)名總?cè)丝?萬(萬稱比(萬元)(萬元)收入(萬(萬人)人)噸)(%)元)北10.59718431999920452790京249.908427706323562863天910.170.5801501226453261128津9136502259317073城鄉(xiāng)居在崗在崗職民年終職工工薪水存儲(chǔ)

7、余人數(shù)總數(shù)額(萬人)(萬元)(萬元)26806410.80577364630111301202.682254931343石家875.400.2332918688521352348709595.60758877莊2680768929911875太299.920.6562360382737111203277394388.65654023原3750937895100呼和浩207.780.44136534381645222105783139642.11309337特2351623588沈677.080.629129558267155679199016135.451152陽941873378241299

8、8811大545.310.494187984261019709227755694.15965922連6739385780187796長691.230.40618535966493570964803102.63884447春8210343810532744哈爾927.090.46226634186674814436450172.791309濱7855123720520020151上10.738206954529644431825971336.845605海313.124019098406485500200445南537.440.5349891991307214116642995680113.81

9、1357京1737269193472861杭616.050.3551414120001711449593742596.901180州6737796883684967947寧538.410.2541428106222210501723524662.15824034波7235866215298350合429.950.318628764251441233628162247.27369577肥4125893517931福583.130.2732152655587467524503069.59680607州3288351851190220廈128.990.486333374575132418758210

10、846.93657484門5124728570331南424.200.398688289230533167714264062.08479,555昌8881674189460文案大全標(biāo)準(zhǔn)合用濟(jì)557.630.40814866285511460690412683.31756696南5302882915775970青702.970.36923821149213176584354978103.52961704島3320036408038045鄭615.360.3426774255287106387252513584.66696848州4601433768338武740.200.586121175069

11、156046585748149.201314漢9291085793442055766長582.470.3101146309885323660346169.57596986沙7367179706718244廣685.000.6211600233482223176120401182.813047州4738139007854499811594深119.850.79329966220368841847951991.261890圳1295754274908900338南285.870.406720486114953149700219045.09371809寧4691130293918海54.380.83

12、54481571746152115174162619.01198138口4345356800重30.206416885855225898,9129090223.731606慶072.347780525441124969804成119355894401974791200都003.560.335590289140632561189684132.89671貴321.500.4553620612247154197908178755.28419681陽7934703143748昆473.390.3867933563605512524216412788.11842321明5729604042900西674

13、.500.40973990536651094088965863114.01885169安4942311766980蘭287.590.544259444294014169540264165.83550890州5884832749568西133.950.52265848711310114913485505127.21219251寧7746469銀95.380.570171603661226217475881410323.72178621川9106193烏魯木158.920.82478513184729254870236555.27517622齊4241668041508文案大全標(biāo)準(zhǔn)合用4.2運(yùn)行結(jié)

14、果cwprint(cwbook.txt,35,10)fid=6數(shù)據(jù)標(biāo)準(zhǔn)化結(jié)果以下：v1=0.05810.03560.04350.06800.05570.11120.11940.11840.10830.13920.04230.03460.03540.07700.00890.06420.04830.04990.05340.05440.04070.01390.06880.02340.00800.00470.01510.03140.02520.01830.01390.03910.00560.00930.00530.02900.00870.01740.02340.01580.00970.02630.00

15、860.00280.00640.00640.00450.00620.01110.00750.03150.03750.03050.01980.02130.03760.02430.03980.03570.02780.02530.02950.04430.02860.02950.04680.03040.03340.02480.02330.03210.02420.04370.02030.01320.02330.01530.02120.02700.02130.04310.02760.06280.01420.01840.01840.02060.02850.04550.03160.06100.04400.04

16、880.18530.01760.10860.18480.11480.08880.13520.02500.03180.02330.04440.03910.02730.02840.02510.03000.03270.02860.02120.03340.04080.04900.02850.01920.03280.02550.02850.02500.01520.03370.03610.06090.02510.02150.02320.01640.01990.02000.01900.01480.00850.01340.00370.01000.00720.01250.00890.02710.01630.05

17、080.02230.02430.01750.02000.02220.01830.01640.00600.02900.00790.01950.01020.00630.01790.00930.0124文案大全標(biāo)準(zhǔn)合用0.01590.01970.02370.01620.00780.01010.00780.00720.01170.01640.01160.02590.02430.03500.02140.01620.02870.01970.01820.02200.01820.03270.02200.05620.03910.03670.04160.02820.02200.02730.02320.02860.

18、02040.01600.01800.02860.01650.01660.02270.02230.01680.03440.03490.02860.02550.02680.03770.02590.02540.03930.03170.02710.01850.02700.01050.02390.01400.01390.01530.01830.01440.03180.03700.03770.07930.06030.05820.07540.09010.04820.07350.00560.04720.00710.06920.02400.01040.07910.04210.02400.04560.01330.

19、02420.01700.00390.01410.00800.00640.00970.01190.00900.00250.04970.00110.00240.01460.00570.00490.00720.00500.00480.14280.01230.09830.02920.14370.06130.03850.04020.05900.03870.04660.01990.04560.02000.11000.04790.02400.03310.03500.02900.01490.02710.00850.00760.04300.01010.00850.00790.01460.01010.02200.

20、02300.01870.01230.01540.02940.02240.01820.02320.02030.03130.02440.01740.01250.02830.02380.01750.02590.03000.02130.01340.03240.00610.01000.00500.01160.00730.01170.01730.01330.00620.03110.00160.00240.00480.00360.00210.00380.00720.00530.00440.03400.00400.00220.00580.00290.00320.00360.00630.00430.00740.

21、04910.00190.00630.00730.02210.01090.01050.01460.0125相關(guān)系數(shù)矩陣:std=文案大全標(biāo)準(zhǔn)合用1.0000-0.34440.84250.36030.73900.62150.40390.49670.67610.4689-0.34441.0000-0.47500.3096-0.35390.19710.35710.26000.15700.30900.8425-0.47501.00000.33580.58910.50560.32360.44560.55750.37420.36030.30960.33581.00000.15070.76640.94120.

22、84800.73200.86140.7390-0.35390.58910.15071.00000.42940.19710.31820.38930.25950.62150.19710.50560.76640.42941.00000.83160.89660.93020.90270.40390.35710.32360.94120.19710.83161.00000.92330.83760.95270.49670.26000.44560.84800.31820.89660.92331.00000.92010.97310.67610.15700.55750.73200.38930.93020.83760

23、.92011.00000.93960.46890.30900.37420.86140.25950.90270.95270.97310.93961.0000特色向量(vec)：vec=-0.13670.2282-0.26280.19390.6371-0.21630.3176-0.1312-0.41910.2758-0.0329-0.02170.00090.0446-0.1447-0.44370.4058-0.55620.54870.0593-0.0522-0.02800.2040-0.0492-0.5472-0.42250.34400.3188-0.44380.24010.0067-0.4176

24、-0.2856-0.23890.1926-0.4915-0.41890.27260.20650.34030.04040.14080.08960.0380-0.1969-0.0437-0.4888-0.6789-0.44050.1861-0.03430.23600.0640-0.82940.03770.26620.1356-0.12900.02780.37820.29810.47390.56850.23580.1465-0.1502-0.26310.12450.21520.36440.15670.3464-0.64850.2489-0.40430.2058-0.07040.04620.12140

25、.38120.4879-0.57070.12170.17610.09870.35500.3280-0.01390.00710.3832-0.7894-0.16280.19250.2510-0.04220.26940.03960.0456文案大全標(biāo)準(zhǔn)合用0.16680.3799特色值(val)val=0.003900000000000.024000000000000.030700000000000.099100000000000.123200000000000.256600000000000.320700000000000.530000000000002.351400000000006.2602

26、特色根排序：6.260222.351380.5300470.3206990.2566390.1232410.09909150.03070880.02403550.00393387各主成分貢獻(xiàn)率：newrate=0.62600.23510.05300.03210.02570.01230.00990.0031文案大全標(biāo)準(zhǔn)合用0.00240.0004第一、二主成分的載荷：0.6901-0.64270.14830.84140.6007-0.68050.85150.31670.4656-0.67540.94630.04260.91170.32990.95370.18620.95890.01090.95060.2558第一、二、三、四主成分的得分：score=0.71850.04990.76842.0000