主成分分析报告matlab程序

标准实用Matlab编程实现主成分分析.程序结构及函数作用在软件Matlab中实现主成分分析可以采取两种方式实现：一是通过编程来实现；二是直接调用Matlab种自带程序实现。下面主要主要介绍利用Matlab的矩阵计算功能编程实现主成分分析。1程序结构Cwprint.m主函数子函数Cwstd.mCwfac.mCwscore.m2函数作用Cwstd.m——用总和标准化法标准化矩阵Cwfac.m——计算相关系数矩阵；计算特征值和特征向量；对主成分进行排序；计算各特征值贡献率；挑选主成分（累计贡献率大于85%），输出主成分个数；计算主成分载荷Cwscore.m——计算各主成分得分、综合得分并排序Cwprint.m——读入数据文件；调用以上三个函数并输出结果文案大全标准实用3.源程序3.1cwstd.m总和标准化法标准化矩阵%cwstd.m,用总和标准化法标准化矩阵functionstd=cwstd(vector)cwsum=sum(vector,1);%对列求和[a,b]=size(vector);%矩阵大小,a为行数,b为列数fori=1:aforj=1:bstd(i,j)=vector(i,j)/cwsum(j);endend3.2cwfac.m计算相关系数矩阵%cwfac.mfunctionresult=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')forz=1:length(y)newy(z)=y(length(y)+1-z);endfprintf('%g\n',newy)rate=y/sum(y);fprintf('\n贡献率：\n')newrate=newy/sum(newy)sumrate=0;newi=[];fork=length(y):-1:1sumrate=sumrate+rate(k);newi(length(y)+1-k)=i(k);ifsumrate>0.85break;endend%记下累积贡献率大85%的特征值的序号放入newi中fprintf('主成分数：%g\n\n',length(newi));fprintf('主成分载荷：\n')forp=1:length(newi)文案大全标准实用forq=1:length(y)result(q,p)=sqrt(newval(newi(p)))*vec(q,newi(p));endend%计算载荷disp(result)3.3cwscore.m%cwscore.m,计算得分functionscore=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.4cwprint.m%cwprint.mfunctionprint=cwprint(filename,a,b);%filename为文本文件文件名，a为矩阵行数(样本数)，b为矩阵列数(变量指标数)fid=fopen(filename,'r')vector=fscanf(fid,'%g',[ab]);fprintf('标准化结果如下：\n')v1=cwstd(vector)result=cwfac(v1);cwscore(v1,result);4.程序测试例 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 4.1原始数据中国大陆35个大城市某年的10项社会经济统计指标数据见下表。文案大全标准实用非地方财城乡居客运货运在岗在岗职年底农业农业工业政民年底城市总量总量职工工工资总人口人口总产值总产值预算内储蓄余名称(万(万人数总额(万人)比(万元)(万元)收入(万额人)吨)(万人)(万元)(%)元)(万元)北10.59718431999920452790268065773410.80京249.908427706323562863646301天0.5801501226453261128113012254910.17202.68津9136502259317073931343石家0.23329186885217095875.4035234895.60758877庄2680768929911875太0.65627371113943299.9223603820327788.65654023原3750937895100呼和浩0.441221396207.7836534381645210578342.11309337特2351623588沈0.6291295582671590161152677.08567919135.45阳9418733782412998811大0.4941879842610197556545.3170922794.15965922连6739385780187796长0.40618535966494803691.23357096102.63884447春8210343810532744哈尔0.462266341866764501309927.09481443172.79滨7855123720520020151上10.7382069545296444318259715605336.84海313.124019098406485500200445南0.53413072141156801357537.44989199664299113.81京1737269193472861杭0.355141412000171174251180616.0544959396.90州6737796883684967947宁0.25414281062222105246538.4150172362.15824034波7235866215298350合0.3182514411622429.9562876423362847.27369577肥4125893517931福0.27321526555875030583.1346752469.59680607州3288351851190220厦0.4865751322108128.9933337441875846.93657484门5124728570331南0.3982305332640424.2068828916771462.08479,555昌8881674189460文案大全标准实用济0.408148662855114126557.6346069083.31756696南5302882915775970青0.36923821149213174978702.97658435103.52961704岛3320036408038045郑0.34252871065135615.3667742538725284.66696848州4601433768338武0.5861211750691557481314740.20604658149.20汉9291085793442055766长0.31011463098853461582.4732366069.57596986沙7367179706718244广0.62116002334822231761204013047685.00182.81州4738139007854499811594深0.7932036884184795191890119.8529966291.26圳1295754274908900338南0.4061149532190285.8772048614970045.09371809宁4691130293918海0.83552162654.384481571746111517419.01198138口4345356800重30.20641688585522590901606898,912223.73庆072.347780525441124969804成1193558944019747912000.335561189132.89都003.56590289140632684671贵0.45522471541787321.5036206119790855.28419681阳7934703143748昆0.38636055124127473.3979335652421688.11842321明5729604042900西0.40936651095863674.50739905408896114.01885169安4942311766980兰0.5442940142641287.5925944416954065.83550890州5884832749568西0.52211133.95658487113104913485505127.21219251宁7746469银0.5702195.381716036612267475881410323.72178621川9106193乌鲁木0.8241847292365158.927851325487055.27517622齐4241668041508文案大全标准实用4.2运行结果>>cwprint('cwbook.txt',35,10)fid=6数据标准化结果如下：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.00860.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.04880.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.05080.02230.02430.01750.02000.02220.01830.01640.00600.02900.00790.01950.01020.00630.01790.00930.0124文案大全标准实用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.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.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.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.04910.00190.00630.00730.02210.01090.01050.01460.0125相关系数矩阵:std=文案大全标准实用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.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.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-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.38120.4879-0.57070.12170.17610.09870.35500.3280-0.01390.00710.3832-0.7894-0.16280.19250.2510-0.04220.26940.03960.0456文案大全标准实用0.16680.3799特征值(val)val=0.003900000000000.024000000000000.030700000000000.099100000000000.123200000000000.256600000000000.320700000000000.530000000000002.351400000000006.2602特征根排序：6.260222.351380.5300470.3206990.2566390.1232410.09909150.03070880.02403550.00393387各主成分贡献率：newrate=0.62600.23510.05300.03210.02570.01230.00990.0031文案大全标准实用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.00000.38060.03860.41924.00000.1848-0.04330.141421.00000.11860.03110.149720.00000.05490.01150.066433.00000.22880.00700.23587.00000.2364-0.00810.228310.00000.1778-0.01670.161116.00000.2292-0.03370.195514.00000.83820.13390.97211.00000.22760.00640.23408.00000.2279-0.02220.205612.00000.1989-0.03820.160718.00000.0789-0.00610.072832.00000.1711-0.03170.139423.00000.09260.02660.119225.00000.0900-0.00000.089928.00000.1692-0.00820.161017.00000.2441-0.03180.212411.00000.1507-0.01080.139922.00000.23160.00120.23289.00000.1294-0.02110.108327.00000.47160.03280.50453.00000.27370.08340.35705.00000.0754-0.00130.074131.0000文案大全标准实用0.04480.03490.079730.00000.4759-0.20280.27316.00000.2907-0.08830.202413.00000.0944-0.01180.082629.00000.15460.00350.158119.00000.1718-0.00920.162615.00000.08650.02300.109526.00000.03490.02160.056635.00000.03430.02280.057234.00000.08890.04220.131024.0000文案大全