matlab與對應(yīng)分析

1、第九章對應(yīng)分析第九章對應(yīng)分析安慶師范學(xué)院安慶師范學(xué)院胡云峰胡云峰習(xí)題習(xí)題用對應(yīng)分析法，對我國16個地區(qū)農(nóng)民1982年支出情況的6項指標(biāo)進行分類。16個地區(qū)農(nóng)民生活水平的調(diào)查數(shù)據(jù)（單位:元）食品衣著燃料住房生活用品及其他文化生活服務(wù)支出地區(qū)x1x2x3x4x5x6北京190.3343.779.7360.5449.019.04天津135.236.410.4744.1613.493.94河北95.2122.839.322.4422.812.8

2、山西104.7825.116.49.8918.173.25內(nèi)蒙128.4127.638.9412.5823.993.27遼寧145.6832.8317.7927.2939.093.47吉林159.3733.3818.3711.8125.295.22黑龍江116.2229.5713.2413.7621.756.04上海221.1138.6412.53115.6550.825.89江蘇144.9829.1211.6742.627.35.74

3、浙江169.9232.7512.7247.1234.355安徽153.1123.0915.6223.5418.186.39福建144.9221.2616.9619.5221.756.73江西140.5421.517.6419.1915.974.94山東115.8430.2612.233.6133.773.85河南101.1823.268.4620.220.54.3解第一步求出協(xié)方差矩陣A=Z’Z程序：x=[190.3343.779.73

4、60.5449.019.04135.2036.4010.4744.1613.493.94...95.2122.839.3022.4422.812.80104.7825.116.409.8918.173.25...128.4127.638.9412.5823.993.27145.6832.8317.7927.2939.093.47...159.3733.3818.3711.8125.295.22116.2229.5713.2413.762

5、1.756.04...221.1138.6412.53115.6550.825.89144.9829.1211.6742.6027.305.74...169.9232.7512.7247.1234.355.00153.1123.0915.6223.5418.186.39...144.9221.2616.9619.5221.756.73140.5421.517.6419.1915.974.94...115.8430.2612.2033.6

6、133.773.85101.1823.268.4620.220.54.3][np]=size(x)T=sum(sum(x))xliehe=zeros(1p)xhanghe=zeros(1n)fk=1:pxliehe(k)=sum(x(:k))end程序F=zeros(p2)ft=1:1:pF(t1)=(0.0357^0.5)(X(t6))F(t2)=(0.0069^0.5)(X(t5))endFG=ZF輸出結(jié)果F=0.05620.026

7、00.02320.02890.05820.01980.16740.01910.01610.06740.01830.0089G=0.00790.00210.00540.00240.00010.00130.00750.00110.00800.00120.00350.00240.01280.00010.00830.00070.02460.00120.00440.00060.00420.00000.00620.00270.00770.00180