数学建模 低碳经济评估

#70.00000080.00000090.000000100.000000110.000000120.000000130.000000140.000000150.000000160.000000170.000000180.000000190.000000200.000000地区2008年2009年2010年2011年2012年英国21.2621.2621.2621.2631.72美国3.173.173.173.179.41西班牙22.2522.2522.2522.2534.85中国30.2730.2730.2730.2734.31model:sets:nian/1..5/;gj/1..4/;links(gj,nian):x,y;endsetsdata:k=10;x=21.2621.2621.2621.2631.723.173.173.173.179.4122.2522.2522.2522.2534.8530.2730.2730.2730.2734.31a=28;enddata@for(links(i,j):y(i,j)=(@abs(k-@abs(a-x(i,j)))));endFeasiblesolutionfound.Totalsolveriterations:0VariableValue

10.00000A28.000001,1)X(21.260001,2)X(21.260001,3)X(21.260001,4)X(21.260001,5)X(31.720002,1)X(3.1700002,2)X(3.1700002,3)X(3.1700002,4)X(3.1700002,5)X(9.4100003,1)X(22.250003,2)X(22.250003,3)X(22.250003,4)X(22.250003,5)X(34.850004,1)X(30.270004,2)X(30.270004,3)X(30.270004,4)X(30.270004,5)X(34.310001,1)Y(3.2600001,2)Y(3.2600001,3)Y(3.2600001,4)Y(3.2600001,5)Y(3.1400002,1)Y(3.8300002,2)Y(3.8300002,3)Y(3.8300002,4)Y(3.8300002,5)Y(2.9690003,1)Y(4.2500003,2)Y(4.2500003,3)Y(4.2500003,4)Y(4.2500003,5)Y(3.1500004,1)Y(4.7300004,2)Y(4.7300004,3)Y(4.7300004,4)Y(4.7300004,5)Y(3.690000

10.00000020.00000030.00000040.00000050.00000060.00000070.00000080.00000090.000000100.000000110.000000120.000000130.000000140.000000150.000000160.000000170.000000180.000000190.000000200.000000RowSlackor地区2008年2009年2010年2011年2012年英国34.5739.8832.4932.6641.35美国18.1219.0119.2920.4868.70西班牙12.6113.6017.7720.0722.00中国25.1319.4519.3319.0521.09Surplusmodel:Surplussets:nian/1..5/;gj/1..4/;links(gj,nian):x,y;endsetsdata:k=5;x=34.5739.8832.4932.6641.3518.1219.0119.2920.4868.7012.6113.6017.7720.0722.0025.1319.4519.3319.0521.09a=25;enddata@for(links(i,j):y(i,j)=(@abs(k-@abs(a-x(i,j)))));endFeasiblesolutionfound.Totalsolveriterations:VariableValue.0000005.000001,1)X(34.570001,2)X(39.880001,3)X(32.490001,4)X(32.660001,5)X(41.350002,1)X(18.120002,2)X(19.010002,3)X(19.290002,4)X(20.480002,5)X(68.700003,1)X(12.610003,2)X(13.600003,3)X(17.770003,4)X(20.070003,5)X(22.000004,1)X(25.130004,2)X(0.9650004,3)X(0.9650004,4)X(0.9650004,5)X(1.2090001,1)Y(0.5700001,2)Y(0.8800001,3)Y(0.4900001,4)Y(0.3600001,5)Y(1.0035002,1)Y(0.2800002,2)Y(0.9900002,3)Y(0.7100002,4)Y(0.480000Y(2,5)3,1)Y(0.790000.700003,2)Y(0.7400003,3)Y(0.230000Y(3,4)0.70000E-013,5)Y(0.5650004,1)Y(0.8700004,2)Y(0.5500004,3)Y(0.6700004,4)Y(0.9500004,5)Y(1.090000RowSlackorSurplus10.00000020.00000030.00000040.00000050.00000060.00000070.00000080.00000090.000000100.000000110.000000120.000000130.000000140.000000150.000000160.000000170.000000180.000000190.000000200.000000四、地区2008年2009年2010年2011年2012年英国4601647087487464840048667美国98422104206110282116335116964西班牙7485074757781607947489998中国7832098507103024103846108122model:sets:nian/1..5/;gj/1..4/;

links(gj,nian):x,y;endsetsdata:k=6;x=4.60164.70874.87464.84004.86679.842210.420611.028211.633511.69647.48507.47577.81607.94748.99987.83209.850710.302410.384610.8122a=10;enddata@for(links(i,j):y(i,j)=(@abs(k-@abs(a-x(i,j)))));endFeasiblesolutionfound.Totalsolveriterations:Feasiblesolutionfound.Totalsolveriterations:Value0VariableK6.000000A10.000001,1)X(4.6016001,2)X(4.7087001,3)X(4.8746001,4)X(4.8400001,5)X(4.8667002,1)X(9.8422002,2)X(10.420602,3)X(11.028202,4)X(11.633502,5)X(11.696403,1)X(7.4850003,2)X(7.4757003,3)X(7.8160003,4)X(7.9474003,5)X(8.9998004,1)X(7.8320004,2)X(9.8507004,3)X(10.302404,4)X(10.384604,5)X(10.812201,1)Y(1.6016001,2)Y(1.7087001,3)Y(1.8746001,4)Y(1.8400001,5)Y(1.8667002,1)Y(5.8422002,2)Y(5.5794002,3)Y(4.9718002,4)Y(4.3665002,5)Y(4.3036003,1)Y(3.4850003,2)Y(3.4757003,3)Y(3.8160003,4)Y(3.9474003,5)Y(4.9998004,1)Y(3.8320004,2)Y(5.8507004,3)Y(5.6976004,4)Y(5.6154004,5)Y(5.187800SurplusSurplusRowSlackor10.00000020.00000030.00000040.00000050.00000060.00000070.00000080.00000090.000000100.000000110.000000120.000000130.000000140.000000150.000000160.000000170.000000180.000000190.000000200.000000五、地区2008年2009年2010年2011年2012年英国1.53971.56861.35541.23631.1394美国1.81472.00331.84051.72881.6310西班牙2.25292.21212.09152.27222.0996中国2.33512.31852.22822.23332.1970model:sets:nian/1..5/;gj/1..4/;links(gj,nian):x,y;endsetsdata:k=1;x=1.53971.56861.35541.23631.13941.81472.00331.84051.72881.63102.25292.21212.09152.27222.09962.33512.31852.22822.23332.1970a=1.5;enddata@for(links(i,j):y(i,j)=(@abs(k-@abs(a-x(i,j)))));endFeasiblesolutionfound.Totalsolveriterations:0ValueVariableValueK1.000000A1.5000001,1)X(1.5397001,2)X(1.5686001,3)X(1.3554001,4)X(1.2363001,5)X(1.1394002,1)X(1.8147002,2)X(2.0033002,3)X(1.8405002,4)X(1.7288002,5)X(1.631000

3,1)X(2.2529003,2)X(2.2121003,3)X(2.0915003,4)X(2.2722003,5)X(2.0996004,1)X(2.3351004,2)X(2.3185004,3)X(2.2282004,4)X(2.2333004,5)X(2.1970001,1)Y(0.96030001,2)Y(0.93140001,3)Y(0.85540001,4)Y(0.73630001,5)Y(0.63940002,1)Y(0.68530002,2)Y(0.49670002,3)Y(0.65950002,4)Y(0.77120002,5)Y(0.86900003,1)Y(0.24710003,2)Y(0.28790003,3)Y(0.40850003,4)Y(0.22780003,5)Y(0.40040004,1)Y(0.16490004,2)Y(0.18150004,3)Y(0.27180004,4)Y(0.26670004,5)Y(0.3030000SurplusSurplusRowSlackor10.00000020.00000030.00000040.00000050.00000060.00000070.00000080.00000090.000000100.000000110.000000

120.000000130.000000140.000000150.000000160.000000170.000000180.000000190.000000200.000000附件二：附录：统一化数据4国家有关数据。表1人均二氧化碳排放单位：吨/（人•年）地区2008年2009年2010年2011年2012年英国1.97732.81362.48442.17592.0389美国3.21933.10123.26133.40213.5024西班牙1.79221.92951.64470.86030.6355中国1.10011.3791.63141.92351.7259表2森林覆盖率单位：%地区2008年2009年2010年2011年2012年英国3.263.263.263.263.14美国3.833.833.833.832.969西班牙4.254.254.254.253.15中国4.734.734.734.733.69表3人均城沛绿地面积单位:平方米/人地区2008年2009年2010年2011年2012年英国0.570.880.490.361.0035美国0.280.990.710.481.7西班牙0.790.740.230.70.565中国0.870.550.670.951.09表4单位建成区面积能耗单位：万吨煤/平方千米地区2008年2009年2010年2011年2012年英国1.60161.70871.87461.841.8667美国5.84225.57944.97184.36654.3036西班牙3.4853.47573.8163.94744.9998中国3.8325.85075.69765.61545.1878表5单位GDP二氧化碳排放单位：吨/万元地区2008年2009年2010年2011年2012年英国0.96030.93140.85540.73630.6394美国0.68530.49670.65950.77120.869中国0.24710.28790.40850.22780.40040.16490.18150.27180.26670.303西班牙附件三：(第一问的求解)%求A的最大特征根A=[11/32/51/53/5;314/33/58/5;5/21/211/27/5;51/43/418/3;5/31/55/321];[v,d]=eig(A)x=max(max(d)')%求A的归一化特征向量(权向量)A=[11/32/51/53/5;314/33/58/5;5/21/211/27/5;51/43/418/3;5/31/55/321];fori=1:5b(i)=0forj=1:5b(i)=A(i,j)+b(i)j+1endendm=0fork=1:5m=b(k)+mendmfori=1:5n(i)=b(i)/mendn附件四：(多项式拟合)y2=[9.438611.503511.893512.057312.0993];y3=[2.79223.07053.35534.13974.3645];y4=[2.10002.37902.63142.92353.2741];subplot(2,2,1);plot(x,y1,'r*');ylabel('人均二氧化碳排放量')；title('英国')；subplot(2,2,2)；plot(x,y2,'r*')；ylabel(人均二氧化碳排放量')；title('美国')；subplot(2,2,3)；plot(x,y3,'r*')；ylabel(人均二氧化碳排放量')；title('西班牙')；subplot(2,2,4);plot(x,y4,'r*');ylabel('人均二氧化碳排放量')；title('中国')；附件五：(灰色预测程序)clearsymsab;c=[ab]';Y1=[2.10002.37902.63142.92353.2741];B=cumsum(Y1);%原始数据累加n=length(Y1);fori=1:(n-1)C(i)=(B(i)+B(i+1))/2;%生成累加矩阵end%计算待定参数的值D=Y1;D(1)=[];D=D';E=[-C;ones(1,n-1)];c=inv(E*E')*E*D;c=c';a=c(1);b=c(2);%预测后续数据F=[];F(1)=Y1(1);fori=2:(n+8)F(i)=(Y1(1)-b/a)/exp(a*(i-1))+b/a;endG=[];G(1)=Y1(1);fori=2:(n+8)G(i)=F(i)-F(i-1);%得到预测出来的数据endt1=2008:2012;t2=2008