节能减排建模论文.doc

节能减排摘要本文针对节能减排对大气环境质量的改善作用，确定了合理的评价指标体系，运用灰色系统理论，建立了灰色gm（1,1）模型，并通过设计相关算法，最终较完整地解决了题中的所有问题。针对问题一：建立定量的综合评价模型。确立为污染大气环境质量的主要因素，通过层次分析法确定在综合指标体系中所占的权重，从而成功的建立评价模型:利用最小二乘法曲线拟合函数，预测出2012年全国各省会城市浓度值。运用程序求解出的各分指数，代入评价模型得出各省会的空气综合评价指数，并进行图表分析，得出兰州空气质量最差，福州、乌鲁木齐次之；南宁海口空气质量最好。针对问题二：建立灰色gm（1,1）模型：并检验模型的可行性，在不采取节能减排的情况下，对2007-2011年的二氧化硫、烟尘和工业粉尘的排放量进行预测，并对预测结果进行分析，得出2011年的大气中二氧化硫排放量比2006年增长47.9%，烟尘排放量比2006年增长15.8%，虽然工业粉尘排放量比2006年减少14.3%，但远远不及二氧化硫与烟尘的增加量，由此可知中国大气环境将加剧恶化。 针对问题三：根据问题二预测结果与节能减排后2007-2011年的实际值进行图表比较，得出2011年二氧化硫、烟尘排放量同比不采取节能减排分别减少42.9%、24.1%，工业粉尘排放量下降到500万吨以下，减排效果显著，大气环境得到了很好的改善。针对问题四：运用最小二乘法曲线拟合函数，预测出不采取节能减排的情况下各省会2007-2010年的浓度，并进行数据处理，且定义的降低率公式，并进行图表分析，得出大多数省会在实施节能减排政策后大气污染物的排放量明显减少，例如南京、上海、福州、成都等城市，但仍有少数几个省会在实施节能减排政策后大气污染物的排放量并未明显减少反而增加，例如石家庄、沈阳、西宁等工业较发达和人口密集的城市。针对问题五：本文提出了大力发展循环经济，落实节能减排任务等建议。关键字：最小二乘法 层次分析 灰色系统模型 降低率一、问题重述环境保护是重大民生问题，随着社会对环境保护的日益重视，人们越来越重视环境的改善，工业革命以来，世界各国尤其是西方国家经济的飞速发展是以大量消耗能源资源为代价的，并且造成了生态环境的日益恶化。节约能源资源，保护生态环境，已成为世界人民的广泛共识。我国从2007年8月起，中央财政开始实施节能减排工作，既是对人类社会发展规律认识的不断深化，也是积极应对全球气候变化的迫切需要。国民经济和社会发展第十一个五年规划纲要提出了“十一五”期间单位国内生产总值能耗降低20%左右，主要污染物排放总量减少10%的约束性指标。根据这两个指标，如中国gdp年均增长一成，五年内就需要节能六亿吨标准煤，减排二氧化硫六百二十多万吨、化学需氧量五百七十多万吨。试根据我国近年污染物总量减排和大气环境相关数据，并结合经济发展情况，根据附录中的数据，结合你们收集到的相关资料，建立数学模型，完成以下问题：1、建立模型对全国各省会城市的大气环境质量做出定量的综合评价，并对2012年各地区大气的污染状况进行分析比较。2、假如不采取节能减排，依照过去几年的主要统计数据，对我国大气环境的发展趋势做出预测分析，3、建立模型分析讨论节能减排对大气环境质量改善所起作用。4、建立模型对节能减排实施前后各省会城市大气环境质量改善情况进行科学分析。5、对下一步实施节能减排提出建议。二、问题分析针对问题一，由附录知影响大气环境质量的最主要因素为 ，收集全国各省会城市历年的浓度数据，根据附录中提供的空气质量分指数的计算公式得出的分指数，再对各分指数加权求和得出各省会城市的大气环境质量指数（aqi）的定量综合评价。由建立的综合评价模型预测出2012年各地区大气的污染状况数据。针对问题二，要求在不采取节能减排的情况下依照过去几年的主要统计数据，对我国大气环境的发展趋势做出预测分析。通过观察附录中的数据，根据国家规定的从2007年开始节能减排的要求，以2001年-2006年的大气中污染物的排放量为依据，利用灰色系统预测模型对我国大气环境的发展趋势做出预测分析。针对问题三，以问题二预测的2007年-2011年大气中污染物的排放量与节能减排后实际值进行比较。针对问题四要求建立模型分析讨论节能减排实施前后各省会城市大气环境质量改善情况。因此，本文通过分别比较分析各省会城市节能减排前后的数据，实现对本问题的求解。三、模型假设1、 影响空气质量的因素主要为等，其他影响因素忽略掉。2、 官方网站的各个省份的指标的统计数据信息真实可靠3、 评价空气质量的各指标间相互作用关系忽略不计四、符号说明 环境空气质量指数 第项影响污染物项目的质量浓度值 第项与相近的污染物浓度限值的低位值 第项与相近的污染物浓度限值的高位值 第项与对应的空气质量分指数 第项与对应的空气质量分指数 第项指标的权重 数据矩阵 数据向量 参数估计 关联度 相对误差 均方差 二氧化硫的降低率 二氧化氮的降低率 粉尘颗粒的降低率、 、平均浓度的预测值和实际值。五、模型建立与求解5.1.1问题一模型建立通过查阅中国统计年鉴，获得2003-2010年所需数据（见附表1），根据（附表1）数据利用最小二乘法曲线拟合函数 ，预测出2012年全国各省会城市浓度值。根据hj633-2012环境空气质量指数技术规定（试行），空气质量分指数的 成分的计算公式，建立如下模型：运用层次分析法对的分指数加权，得出大气环境质量指数本文以分别对环境大气能见度、水体富营养化、有毒物质、酸雨、光化学烟雾、破坏臭氧层、温室效应、呼吸道疾病八方面的危害轻重，建立权重 的层次结构模型（图1）图1 层次结构模型5.1.2问题一的模型求解准则层判断矩阵如表1a b 1b2b3b4b5b6b7b8b1111/81/71/51/71/21/9b2111/711/51/611/8b387156431/2b4711/512121/3b5551/61/2111/21/5b6761/411131/3b7211/31/221/311/6b898235361表1 准则层判断矩阵方案层的判断矩阵如表2b1c1c2c3b2c1c2c3b3c1c2c3b4c1c2c3c1111/8c111/58c1117c1158c2111/8c2511c2117c21/516c3881c31/811c31/71/71c31/81/61b5c1c2c3b6c1c2c3b7c1c2c3b8c1c2c3c111/71c1111/2c1121c111/21/5c2718c2111/3c21/211c2211/4c311/81表2 方案层的判断矩阵c3231c3111c3541用matlab软件求解权向量w（程序见附录1）得， 因此环境空气质量指数运用matlab软件预测出2012年全国各省会城市浓度值（程序见附录2，结果见附表1），将各浓度值代入上式中运用matlab软件得出各省会的空气综合质量指数aqi值（程序见附录3，结果见附表2）。并用excel绘制出以下图表：从图中可以看出空气质量最差的是兰州，质量指数超过100，其次福州、乌鲁木齐空气质量较差，南宁海口空气质量最好。空气质量指数级别达到一级的省会为41.9%，空气质量指数级别达到二级的省会为54.8%，只有兰州空气质量指数级别为三级轻度污染。5.1.3问题一结果检验利用最小二乘法曲线拟合函数预测出2011年全国各省会城市浓度值（见附表1），与中国环境保护部中给出的2011年各省会浓度范围相符。5.2.1问题二模型建立及求解利用节能减排前2001-2006年中全国的二氧化硫、烟尘和工业粉尘的排放量（相关数据见给出资料），对其进行量变分析，对其运用灰色理论模型gm（1，1），预测我国大气环境的发展趋势。序列共有6个观察值(1)累加生成序列(2) 的紧邻均值生成序列(3)构造数据矩阵b 及数据向量y (4)用最小二乘法估计参数u记 (5)建立灰微分方程: 求解得(6)累减生成得到利用matlab软计算(程序见附录4)结果如下表格5.2.2问题二的模型检验（运用matlab软件对一下各检验公式求解）(1)残差检验序列与序列的残差序列：相对残差序列： (2)关联度检验计算关联系数： 计算关联度： （n=6）(3)均方差比检验的均值、方差分别为 =的均值、方差分别为： =计算c：实际值预测值相对误差关联度均方差20011947.81947.800.57800.107420021926.61968.50.0217620032158.72119.60.0181120042254.92282.30.0121420052549.32457.40.0360320062588.82646.00.02211检验结果表4 二氧化硫排放量的预测值及检验精确合格精确该模型的二氧化硫排放量预测值与实际值相的几种检测均合格，该模型可进行预测。我们可以利用此模型对2007年-2011年二氧化硫排放量进行预测，结果如下表，烟尘排放量与工业粉尘排放量的预测方法与二氧化硫预测方法相同，预测结果见下表：二氧化硫排放量烟尘排放量工业粉尘排放量合计工业生活合计工业生活20072849.12645.238531163.4929.1234.4812.120083067.82934.9390.11191.6950.6241.2780.320093303.23256.2394.91220.5972.5248.2749.720103556.83612.7399.71250.0994.9255.4720.320113829.7表5 预测2007-2011年全国空气中污染物排放量4008.3404.81280.31017.9262.8692.1由上表可以看出，若不采取节能减排，2011年的大气中二氧化硫排放量为3829.7万吨，比2006年增长47.9%，烟尘排放量为1280.3万吨，比2006年增长15.8%，虽然工业粉尘排放量同比2006年减少14.3%，但远远不及二氧化硫与烟尘的增加量，大气污染会越来越严重。5.3.问题三求解根据第二问的预测数据与采取节能减排后的数据对比见下图由二氧化硫排放量对比图看出采取节能减排措施后，虽然每年二氧化硫的减排趋势缓慢，但是抑制了减排前逐年迅速增长的状况；2011年二氧化硫排放量同比不采取节能减排减少42.9%由烟尘排放量对比图看出烟尘减排效果明显，2011年尘年比2006年减排了248.7万吨，同比不采取节能减排减少24.1%由工业粉尘排放量对比图看出，未采取减排措施前，工业粉尘排放量有所下降，采取措施后，减排速率明显，2011年工业粉尘排放量下降到500万吨以下。综上可知，采取节能减排措施后，大气中各污染物得到了明显的控制，空气环境质量有了显著的改善。5.4.1问题四模型建立在问题一的结论基础上，采用各省会的2007-2010年的空气质量指数的实际数据与不实施节能减排的情况下的预测数据进行对比分析，并用excel软件对2008-2010年全国各省会城市的二氧化硫、二氧化氮、粉尘颗粒平均浓度的实际数据与预测数据（四年平均值）计算出各省会降低率，其计算式分别如下：5.4.2问题四模型求解用excel软件对每年全国各省会城市的空气质量指数的实际数据与预测数据计算出各省会二氧化硫、二氧化氮、粉尘颗粒平均浓度降低率，结果见（附表4），由降低率可以看出节能减排的实施对各省会城市大气环境质量起到了明显的改善作用。大多数省会在实施节能减排政策后污染物的排放量明显减少，例如南京、上海、福州、成都等城市，但仍有少数几个省会在实施节能减排政策后大气污染物的排放量并未明显减少反而增加，例如石家庄、沈阳、西宁等工业较发达和人口密集的城市。由降低率可以看出节能减排的实施对各省会城市大气环境质量起到了明显的改善作用。大多数省会在实施节能减排政策后污染物的排放量明显减少，例如北京、重庆、济南、乌鲁木齐等城市，但也有少数几个省会在实施节能减排政策后大气污染物的排放量并未明显减少，例如合肥、兰州、武汉等重工业城市仍需加大力度。由降低率图可明显看出节能减排的实施对各省会城市大气环境质量得到了很好的改善，且大多数城市都保持在较高水平，只有极少数城市不明显。综上可知节能减排的实施对大部分省会城市大气环境质量起到了非常明显的作用。尽管在大力发展经济的前提下，污染物的排放并没有随着工厂数的增多而增多，反而是在进一步改善这说明节能减排所带来的环境效益是极大的，我们应继续将节能减排实施下去。但有少数城市改善并不明显，还需要继续加强节能减排的方案实施。55施节能减排的措施第一，.政府需加大对节能减排的宣传力度，深化全民节能意识，改善其生活方式和思想意识，引导全民从小事做起、从自我做起实现节能减排。第二，鼓励企业主动去落实节能减排任务，全面提高能源的使用效率，彻底淘汰落后的设备工艺，从而实现节能减排 第三，依靠技术，鼓励发展高新科技。研究开发如风能、太阳能等低碳环保的新能源。加快节能减排技术产业化示范和推广，加快建立节能减排技术服务体系，推动环保产业健康发展，加强国际交流合作。第四，坚持贯彻落实科学发展观，大力发展“循环经济”。走可持续发展的新型工业道路。从“高投入高能耗高污染”的粗放型经济增长方式向集约型经济转变，建设资源节约型社会，环境友好型社会。六、模型的评价与推广模型评价问题一中，本文以分别对环境大气能见度、水体富营养化、有毒物质、酸雨、光化学烟雾、破坏臭氧层、温室效应、呼吸道疾病八方面的危害轻重，建立权重 的层次结构模型。将多个指标转换为一个指标，根据（附表1）数据利用最小二乘法曲线拟合函数，拟合函数预测出2011年全国各省会城市浓度值，通过excel来综合评价各省会城市的大气环境质量。从而能够比较客观地进行比较，达到综合评价的目的。问题二中，在对2012年各地区大气的污染状况进行分析时，由于数据较少，我们运用了灰色理论模型求解，使预测结果更加精确。模型推广问题一中的权重分析模型不仅可以科学地评价空气质量问题，还可以推广到水质评价问题以及经济、军事等领域的许多综合评价问题，具有广泛的应用价值。问题二采用的灰色理论预测法是目前最为广泛的预测方法之一，在数据较少时相比其他预测模型也具有一定优势。可以在工资预测、粮食作物的生产、经济管理、工程管理等方面进行应用。七、参考文献【1】 中华人民共和国国家统计局年鉴【2】中华人民共和国国家环境保护标准hj633-2012环境空气质量指数（aqi）技术规定（试行）【3】中国环境保护部，中国环境状况公报，2011【4】姜启源，数学模型（第四版），高等教育出版社，2011年1月【5】卓金武，matlab在数学建模中的应用，北京航空航天大学出版社，2011年4月【6】严哲，易纲，excel2007图表制作应用于技巧详解，北京邮电出版社，009年2月附录附录1层次分析求解程序fid=fopen(e:matlab 2010bworktxt3.txt,r);n1=8;n2=3;a=;for i=1:n1tmp=str2num(fgetl(fid);a=a;tmp; %读准则层判断矩阵endfor i=1:n1str1=char(b,int2str(i),=;);str2=char(b,int2str(i),=b,int2str(i),;tmp;);eval(str1);for j=1:n2tmp=str2num(fgetl(fid);eval(str2); %读方案层的判断矩阵endendri=0,0,0.58,0.90,1.12,1.24,1.32,1.41,1.45; %一致性指标x,y=eig(a);lamda=max(diag(y);num=find(diag(y)=lamda);w0=x(:,num)/sum(x(:,num);cr0=(lamda-n1)/(n1-1)/ri(n1)for i=1:n1x,y=eig(eval(char(b,int2str(i);lamda=max(diag(y);num=find(diag(y)=lamda);w1(:,i)=x(:,num)/sum(x(:,num);cr1(i)=(lamda-n2)/(n2-1)/ri(n2);endcr1, ts=w1*w0, cr=cr1*w0附录2最小二乘法曲线拟合预测a=0.0310.0500.116 0.050.0440.090 0.0630.0460.089 0.0660.0180.092 0.0470.0320.064 0.0620.0350.095 0.0320.0510.081 0.0430.0450.100 0.0180.0480.074 0.0230.0430.113 0.0190.0530.092 0.0190.0330.111 0.0050.2000.075 0.0550.0500.086 0.0420.0320.119 0.050.0460.111 0.0370.0570.098 0.0340.0490.077 0.0260.0490.064 0.0250.0210.063 0.0060.0180.037 0.0410.0350.099 0.020.0550.102 0.0580.0320.069 0.0360.0510.068 0.0080.0190.044 0.040.0470.114 0.0560.0490.182 0.0570.0240.142 0.0370.0310.092 0.090.0670.125 ;xda=1:1:4;for i=1:30yda=a(i,:);x0=0,0,0;x,resnorm=lsqcurvefit(nihehanshu,x0,xda,yda);x1=6,7;f=x(1)*exp(x(2)*x1)+x(3)end附表1：二氧化氮城市200320042005200620072008200920102011（检验）2012(预测)北 京0.0720.0710.0660.0660.0660.0490.0530.0570.0510.050天 津0.0520.0520.0470.0480.0430.0410.0400.0450.0430.044石 家 庄0.0440.0420.0410.0390.0350.0310.0350.0410.0420.046太 原0.0310.0220.0200.0250.0270.0210.0220.0200.0190.018呼和浩特0.0460.0380.0410.0480.0480.0450.0400.0340.0340.032沈 阳0.0360.0350.0360.0430.0360.0370.0370.0350.0350.035长 春0.0220.0320.0350.0390.0380.0380.0430.0440.0470.051哈 尔 滨0.0650.0600.0560.0490.0600.0550.0540.0480.0470.045上 海0.0570.0620.0610.0550.0540.0560.0530.0500.0500.048南 京0.0490.0550.0540.0520.0510.0530.0480.0460.0450.043杭 州0.0560.0550.0580.0570.0570.0530.0520.0560.0530.053合 肥0.0250.0170.0250.0320.0260.0250.0270.0300.0310.033福 州0.0340.0410.0420.0490.0550.0460.0400.0320.0260.200南 昌0.0340.0290.0310.0320.0340.0360.0370.0420.0450.050济 南0.0460.0380.0240.0210.0230.0220.0250.0270.0290.032郑 州0.0330.0370.0390.0440.0450.0470.0460.0460.0460.046武 汉0.0520.0540.0500.0490.0550.0540.0540.0570.0560.057长 沙0.0380.0330.0360.0390.0410.0430.0420.0460.0470.049广 州0.0720.0730.0680.0670.0650.0560.0560.0530.0510.049南 宁0.0320.0340.0380.0350.0480.0440.0280.0300.0240.021海 口0.0130.0130.0150.0120.0120.0170.0160.0150.0170.018重 庆0.0460.0670.0480.0470.0440.0430.0370.0390.0360.035成 都0.0460.0480.0520.0490.0490.0520.0550.0510.0540.055贵 阳0.0190.0240.0130.0170.0230.0230.0260.0270.0290.032昆 明0.0330.0400.0380.0440.0420.0390.0460.0460.0480.051拉 萨0.0290.0200.0250.0260.0250.0240.0210.0210.0200.019西 安0.0350.0330.0320.0420.0430.0440.0460.0450.0460.047兰 州0.0500.0450.0370.0520.0420.0540.0430.0480.0480.049西 宁0.0310.0270.0260.0290.0350.0300.0320.0260.0260.024银 川0.0370.0400.0250.0270.0250.0210.0310.0260.0300.031乌鲁木齐0.0550.0560.0640.0670.0650.0680.0670.0670.067二氧化硫城市200320042005200620072008200920102011（检验）2012(预测)北 京0.0610.0550.050.0520.0470.0360.0340.0320.0310.031天 津0.0740.0730.0760.0670.0620.0610.0560.0540.0520.05石 家 庄0.1520.0870.0540.0440.0430.0460.0450.0540.0570.063太 原0.0990.0870.0770.080.0760.0730.0750.0680.0680.066呼和浩特0.0390.0450.050.0540.0660.0490.0490.0460.0470.047沈 阳0.0520.0520.0540.0580.0540.0590.0590.0580.060.062长 春0.0120.0130.0260.0260.0300.030.0340.030.0320.032哈 尔 滨0.0430.0420.0420.0340.0480.0430.0460.0450.0440.043上 海0.0430.0550.0610.0510.0550.0510.0350.0290.0240.018南 京0.030.0450.0520.0630.0580.0540.0350.0360.0280.023杭 州0.0490.0490.060.0560.0600.0520.0410.0340.0270.019合 肥0.01230.0130.0180.0240.0230.0220.0230.020.020.019福 州0.00770.010.0160.020.0270.0230.0140.0090.0070.005南 昌0.0510.0570.050.0560.0540.050.0540.0550.0550.055济 南0.0640.0450.060.040.0560.0520.050.0450.0440.042郑 州0.050.0570.0590.060.0690.060.0530.0530.0510.05武 汉0.0490.0480.0540.0570.0610.0510.0440.0410.0380.037长 沙0.0810.0840.0810.0820.0650.0530.0390.040.0350.034广 州0.0590.0770.0530.0540.0510.0460.0390.0330.030.026南 宁0.0470.0610.0580.0590.0590.040.0320.0280.0260.025海 口0.0090.0070.0120.010.0090.0090.0070.0070.0070.006重 庆0.1150.1130.0730.0740.0650.0630.0530.0480.0450.041成 都0.0520.0670.0770.0650.0620.0490.0380.0310.0250.02贵 阳0.0890.0940.0630.0650.0550.0640.0580.0570.0580.058昆 明0.0450.0690.0550.0620.0680.0510.0410.040.0370.036拉 萨0.0020.0030.010.0090.0070.0050.0080.0070.0070.008西 安0.0570.0490.0440.0560.0530.050.0480.0430.0430.04兰 州0.0860.0710.0680.0570.0600.070.0590.0570.0570.056西 宁0.0310.0240.0290.0240.0280.0290.0420.0390.0490.057银 川0.0630.0540.0540.0480.0490.0490.0440.0390.0390.037乌鲁木齐0.09670.1020.1160.1130.0880.1050.0930.0890.0910.09可吸入颗粒物城市200320042005200620072008200920102011（检验）2012(预测)北 京0.1410.1490.1410.1620.1480.1230.1210.1210.1190.116天 津0.1330.1110.1060.1140.0940.0880.1010.0960.0920.090石 家 庄0.1750.1230.1320.1420.1280.1160.1040.0980.0930.089太 原0.1720.1750.1390.1420.1240.0940.1060.0890.0960.092呼和浩特0.1160.0800.0970.1020.0840.0700.0740.0680.0660.064沈 阳0.1350.1370.1180.1170.1190.1180.1100.1010.0990.095长 春0.0980.0850.0990.0990.0990.0960.0850.0890.0840.081哈 尔 滨0.1210.1130.1040.1040.1020.1020.1010.1010.1000.100上 海0.0970.0990.0880.0860.0880.0840.0810.0790.0760.074南 京0.1200.1210.1100.1090.1070.0980.1000.1140.1110.113杭 州0.1190.1100.1120.1110.1070.1100.0970.0980.0950.092合 肥0.1000.1100.0950.0990.1160.1340.1110.1150.1130.111福 州0.0800.0740.0720.0720.0650.0710.0640.0730.0730.075南 昌0.1000.0990.0890.0860.0830.0830.0790.0870.0850.086济 南0.1490.1490.1280.1140.1180.1260.1230.1170.1200.119郑 州0.1070.1110.1090.1110.1050.0940.0990.1110.1080.111武 汉0.1330.1300.1190.1210.1230.1130.1050.1080.1010.098长 沙0.1350.1400.1220.1110.1040.0970.0920.0830.0810.077广 州0.0990.0990.0880.0760.0770.0710.0700.0690.0660.064南 宁0.0720.0780.0670.0660.0640.0560.0500.0690.0620.063海 口0.0300.0330.0400.0410.0430.0430.0380.0400.0380.037重 庆0.1470.1420.1200.1110.1080.1060.1050.1020.1000.099成 都0.1180.1150.1250.1230.1110.1110.1110.1040.1040.102贵 阳0.1040.0830.0760.0830.0850.0820.0740.0750.0720.069昆 明0.0860.0850.0820.0910.0750.0670.0670.0720.0680.068拉 萨0.0650.0520.0700.0620.0570.0510.0500.0480.0460.044西 安0.1360.1420.1290.1330.1350.1130.1130.1260.1160.114兰 州0.1740.1720.1580.1920.1290.1320.1500.1550.1680.182西 宁0.1390.1270.1140.1350.1150.1180.1410.1240.1370.142银 川0.1320.1220.0900.0970.0920.0840.0900.0930.0920.092乌鲁木齐0.1270.1140.1140.1520.1360.1450.1400.1330.1280.125附录3综合评价指标求解程序function aqi=fun(s,m,n)load e:matlab 2010bwork03.txts=1000.*x03(:,1);m=1000.*x03(:,2);n=1000.*x03(:,3);z=length(s);for k=1:z if (s(k)=0)&(s(k)50)&(s(k)150)&(s(k)=0)&(m(k)40)&(m(k)80)&(m(k)180)&(m(k)280)&(m(k)=0)&(n(k)50)&(n(k)150)&(n(k)250)&(n(k)=350) c(k)=1/2*(n(k)-250)+250; end end endendendfor k=1:z aqi(k)=0.312*a(k)+0.336*b(k)+0.353*c(k)end附表2 各省会综合评价城市综合指标排名 兰 州101.20341 福 州 87.46252 乌鲁木齐 77.50753 西 宁60.664 郑 州 59.97655 武 汉58.2466 石 家 庄58.12157 南 昌 58.0248 西 安57.8069 沈 阳 57.764510 北 京 56.61111 济 南56.372512 哈 尔 滨 55.43113 天 津 55.4314 重 庆 53.790515 成 都52.80816 贵 阳 51.291517 长 春 51.165518 太 原 50.71919 南 京 50.645520 长 沙 50.243521 昆 明50.11922 杭 州 49.89123 银 川 49.62