城市空气质量评估及预测数学建模论文文档

如不慎侵犯了你的权益，请联系我们告知！西南交通大学“新秀杯”数学建模竞赛参赛论文论文题目：城市空气质量评估及预测

如不慎侵犯了你的权益，请联系我们告知！城市空气质量评估及预测摘要：本文对我国的成都，杭州，北京，上海，广州，拉萨，乌鲁木齐，郑州，武汉，西安等10座城市的每日空气质量详细列表进行科学分析，利用层次分析法和指数平滑法等数学建模方法对其空气质量进行研究，综合考虑各种因素建立如下数学模型。1、十个城市空气污染严重程度的科学排名。本文采用层次分析法来对10个城市的环境污染情况进行比较。用excel统计出各个城市不同级别的污染天数。在采用层次分析法，首先建立层次结构模型然后构造成对比较阵，用matlab软件计算权向量并做一致性检验，计算组合权向量并做组合一致性检验。最后根据权重的大小就可以比较出10个城市空气污染严重程度的科学排名。2、成都市11月的空气质量状况预测。本文采用指数平滑的方法来建立数学模型。根据前面的统计数据，考虑到污染级别不同和首要污染物的种类两个因素来预测成都11月的空气质量状况。3、分析影响城市空气污染程度的首要污染物所占总天数的比例，然后再综合考虑各方面因素原因，以找出影响城市空气污染程度的主要因素。本文就10个城市的空气污染严重程度的排名给出了相关的分析以及应对策略。就成都11月份的空气质量状况预测给出出行和生活方面的建议。并结合当下倡导建设环境友好型和资源节约型社会出发，就如何兼顾经济发展与环境保护给出指导性建议。关键词：空气质量主要因素。先通过excel统计出每个城市的各种找出造成该污染物超标的以最后污染等级层次分析法指数平滑法首要空气污染物

如不慎侵犯了你的权益，请联系我们告知！一问题的提出1.1背景介绍随着科技的发展，工业的进步和全球人口急剧增多的因素的影响，人们赖以生存的环境遭到了很大的破坏，很多地区相继出现了酸雨、物种灭绝、土地沙化等环境问题，环境问题已经成为当今世界各国普遍关注的问题之一，也是21世纪人类面临的重大挑战。我国是一个人口大国，城市众多，人口密集。但由于工业的发展，我们的很多城市都受到了不同程度的污染，尤其是空气的污染，直接对我们造成伤。害空气中的污染物主要是可吸入颗粒、二氧化硫、二氧化氮等物质。1.2问题的解决我们通过统计不同城市的每天空气质量详细列表，运用层次分析法可以对不同城市的空气质量进行统计和比较。在预测成都11月的空气质量状况时采用指数平滑的方法来建立数学模型。二基本假设1.根据空气污染指数（ＡＰＩ）来划分为５个等级：API值小于等于50，空气质量为优，相当于国家空气质量一级标准；API值大于50且小于等于100，表明空气质量良好，相当于达到国家质量二级标准；API值大于100且小于等于200，表明空气质量为轻度污染，相当于国家空气质量三级标准；API值大于200表明空气质量差，称之为中度污染，为国家空气质量四级标准；API大于300表明空气质量极差，已严重污染。2．数据取自国家统计局，具有使用价值。三问题的分析1.不同城市的空气质量统计和比较(1)最大特征值λ(max)的Matlab计算方法：[V,D]=eig(A),其中A为待计算特征值的矩阵，D为对角矩阵，其对角元素为A的特征值，最大的即为λ(max)。(2)一致性指标CI计算方法：CI=(λ-n)/(n-1);（其中λ为矩阵A的最大特征值，n为矩阵的阶数）(3)随即一致性指标RI的计算方法：RI与n有如下关系，如表n12345678910

如不慎侵犯了你的权益，请联系我们告知！RI000.580.901.121.241.321.411.451.49(5)权重计算方法计算矩阵A的特征根及特征向量，将所求的特征向量单位化后得到的就是权重值。2.成都11月份空气质量进行预测指数平滑法（ExponentialSmoothing，ES）是布朗(RobertG..Brown)所提出，布朗、认为时间序列的态势具有稳定性或规则性，所以时间序列可被合理地顺势推延；他认为最近的过去态势，在某种程度上会持续的未来，所以将较大的权数放在最近的资料。指数平滑法是在移动平均法基础上发展起来的一种时间序列分析预测法，它是通过计算指数平滑值，配合一定的时间序列预测模型对现象的未来进行预测。其原理是任一期的指数平滑值都是本期实际观察值与前一期指数平滑值的加权平均。基本模型如下：Ft+1=a﹒Yt+(1-a)FtFt+1为第t+1时期的时间序列预测值Yt为第t时间的时间序列的Ft为第t时间的a为平滑系数,其取值范围为[0,1]；实际值时间序列的预测值四模型的建立与求解4.1建立数学模型给出十个城市空气污染严重程度的科学排名1.将研究目标（Z）、因素（P）、对象（C）按相关关系分成目标层Z、准则层P、对象层C。层次结构图如图所示：

如不慎侵犯了你的权益，请联系我们告知！空气质量空气质量一级空气质量二级空气质量三级空气质量四级空气质量五级成都杭州北京上海………2.给出空气质量一级，二级，三级，四级和五级两两成对比较的判断矩阵P污染级别一级二级三级四级五级12345根据上图得出如下例两两成对比较的判断矩阵P权重五级四级三级二级一级五级1/3四级4/15三级1/5二级2/15一级1/1515/415/34/315/22543414/53/52/51/53/41/21/43/212/31/31/2由表中数据，计算可知：λ(max)=5.00，CI=0.00，RI=1.12，CR=0.00<0.1。因为CR=0.00<0.1，所以此排序有满意的一致性。3.给出对象层对准则层的各个因素的判断矩阵并进行分析。由于各个城市只存在污染程度的不同，所以它们两者之间各因素之间的关系。2010空气污染指数统计优良轻度中度严重污染污染污染495035224220201333059002021成都杭州北京

如不慎侵犯了你的权益，请联系我们告知！127951561939114710300051001011121上海广州拉萨212乌鲁木齐1921624920524466364843155214郑州武汉西安根据表中数据，类比（2）中方法，计算出各种不同污染等级对不同城市的权重成都杭州北京上海广州拉萨乌鲁木郑州武汉西安权重齐成1都49/507/549/12749/9549/21249/1949/1549/527/20.073杭50/491州10/750/12725/9525/10650/1910/325/2625/70.07435/12735/9535/21235/197/335/525/20.052北5/77/101京上127/49127/50127/351海127/95127/212127/19127/15127/52127/140.190广95/4919/1019/795/1271州95/212519/395/5295/140.142212/19212/1553/13106/70.31719/212119/1519/5219/140.028拉212/49106/25212/35212/127212/951萨乌19/4919/5019/3519/1271/5鲁木齐郑15/493/103/715/1273/1915/21215/191州15/5215/140.02226/70.077武52/4926/2552/3552/12752/9513/5352/1952/151汉西2/77/252/514/12714/957/10614/1914/157/261安0.021空气质量一级对10个城市的不同权重表由表中数据，计算可知：λ(max)=10.00，CI=0.00，RI=1.49，CR=0.00<0.1同理，可以计算出其余空气质量等级对10个城市的不同权重。计算方法类似，用matlab软件的计算过程详见本文附表。空气质量二级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁郑州武汉西安木齐权重0.1120.1100.1000.0780.0960.0450.1080.1240.1030.122

如不慎侵犯了你的权益，请联系我们告知！由表中数据，计算可知：λ(max)=10.00，CI=0.00，RI=1.49，CR=0.00<0.1空气质量三级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁郑州武汉西安木齐权重0.0690.0630.1240.3090.0210.0060.1390.0760.1010.091由表中数据，计算可知：λ(max)=10.00，CI=0.00，RI=1.49，CR=0.00<0.1空气质量四级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁郑州武汉西安木齐权重000.250000.6250.12500由表中数据，计算可知：λ(max)=10.00，CI=0.00，RI=1.49，CR=0.00<0.1空气质量五级对10个城市的不同权重表城市成都杭州北京上海广州拉萨乌鲁郑州武汉西安木齐权重00.20.10.100.10.10.10.20.1由表中数据，计算可知：λ(max)=10.00，CI=0.00，RI=1.49，CR=0.00<0.14.模型求解进行层次总排序，方法：将上面5个空气质量等级对10个城市的不同权重表单位化后作为列向量构成10×5矩阵，和空气质量一级，二级，三级，四级和五级两两成对比较的判断矩阵P相乘，结果便得到10个城市的权重值。根据上述三问题的分析中的假设可知，权重值越大，表明空气污染情况越严重。因此，将10个城市的权重值按照从小到大依次排序，得出的结果便是10个城市的空气污染严重程度的排名。最终结果如下表所示：一级二级三级四级五级总权重(0.067)(0.133)(0.200)(0.267)(0.333)成都0.0730.1120.0690杭州0.0740.1100.0630北京0.0520.1000.1240.25上海0.1900.0780.3090广州0.1420.0960.021000.0335870.0987880.1416340.1182040.0264820.20.10.10

如不慎侵犯了你的权益，请联系我们告知！拉萨0.3170.0450.00600.10.0617240.2442150.0998410.1056580.069133乌鲁木齐0.0280.1080.1390.6250.1郑州0.0220.1240.0760.1250.1武汉0.0770.1030.1010西安0.0210.1220.0910总的一致性检验0.20.1CR=0.067×0+0.133×0+0.200×0+0.267×0+0.333×0=0<<0.1此结果说明排序结有非常满意的一致性。5.结果分析根据10个城市的总权重值进行从小到大依次排序，空气污染严重程度的排名如下：(1)广州(2)成都(3)拉萨(4)西安(5)杭州(6)郑州(7)武汉(8)上海(9)北京(10)乌鲁木齐结论显示广州的空气质量状况最好，而上海，北京，乌鲁木齐的空气质量状况较差。广州今年在筹备2010广州亚运会，因此在城市环境保护方面下大力气进行整治，空气质量状况得到明显改善。而北京上海作为我国的一线城市，城市人口密度大，机动保车有量大，交通拥堵，工业化程度高，因此城市空气污染情况严重。乌鲁木齐由于地理位置的因素，处在大西北，经常遭受沙尘暴等恶劣天气影响，因此空气污染情况严重，特别是可吸入颗粒明显超标。4.2建立模型对成都市11月的空气质量状况进行预测预测成都市11月的空气质量状况主要采用指数平滑的方法来建立数学模型。1.用excel统计出成都市2005年到2010年的空气污染指数，如下表所示：成都2005~2010空气污染指数统计优良轻微污染轻度污染中度污染中度重污重污染染2005200620072008200920102135312534492662512742912782246063564338281510204201002000040000071152.根据上表，做出统计表，然后做出分布图。如下所示：

如不慎侵犯了你的权益，请联系我们告知！0~100100~200大于200年份占全年百分比/%空气污染指数20052006200720082009201078.356278.356283.561686.575386.187889.215720.547920.000015.890413.698613.535910.78431.09591.64380.24790027620分布图如下所示：3.分析2005年到2010年的空气污染指数由这个分布图我们可以看出每年成都空气污染指数的大体情况，结合2010年每月11月份的空气污染指数进行预测。的空气污染指数，我们对成都运用指数平滑法进行计算。分析:为了预测成都11月份的空气污染指数，不仅要知道前面几月的实际空气污染指数，还要知道前面几个月的预测值。例如：我们已经知道了成都11月份前4天的空气污染指数。(这只是一个举例)F1=Y1，由于空气污染指数有波动，但长期趋势变化不大时，可选稍大的α值，常在0.1～0.4之间取值；我们取a=0.3进行计算，则F2=0.3Y1+0.7F1,F3=0.3Y2+0.7F2……日期空气污染指数移动平均法预预测偏差偏方平方值测值2010/11/12010/11/22010/11/37684607678.48-18.464338.6

如不慎侵犯了你的权益，请联系我们告知！72.88-4.882010/11/46823.8426.4合计解得2010/11/5的空气污染指数为F5=0.3Y4+0.7F4=0.3×68+0.7×72.8871第5日的预测偏差的平方值为=426.4/3=142.1不同的α取值对第5日的预测值和预测偏差的平方值的影响，公式写为Ft+1=a.Yt+(1-a)Ft=a.Yt+Ft-aFt=Ft+a(Yt-Ft)选取不同的值获得结果如下表所示:a第5日预测值第5日预测偏差平方估计值07674737170696867676668128121136142.1600.10.20.30.40.50.60.70.80.91.056168170188210235当α取0.5时，用指数平滑法求得的第5日的预测值为69，预测偏差的平方估计值最小为56。因此，0.5为α最适合的取值，69为第5日最精确的预测值。其它日子的空气污染指数计算方法如上。4.最终对成都市11月的空气质量状况的预测情况如下：成都市11月的空气质量状况预测详细列表空气质量级空气质量状况污染指数日期首要污染物别可吸入颗粒2010-11-302010-11-292010-11-282010-11-2774ⅡⅡⅡⅡ良物可吸入颗粒725168良良良物可吸入颗粒物可吸入颗粒物

如不慎侵犯了你的权益，请联系我们告知！Ⅰ2010-11-262010-11-2544––优可吸入颗粒8177Ⅱ良物可吸入颗粒物2010-11-24Ⅱ良2010-11-232010-11-2247––49––ⅠⅠ优优可吸入颗粒2010-11-212010-11-202010-11-192010-11-182010-11-178777695567ⅡⅡⅡⅡⅡ良良良良良物可吸入颗粒物可吸入颗粒物可吸入颗粒物可吸入颗粒物2010-11-162010-11-1550––47––ⅠⅠ优优可吸入颗粒物可吸入颗粒2010-11-142010-11-132010-11-12539198Ⅱ良Ⅱ良物可吸入颗粒Ⅱ良物可吸入颗粒物2010-11-11117Ⅲ1Ⅲ1Ⅲ1Ⅲ1Ⅱ轻微污染轻微污染轻微污染轻微污染良可吸入颗粒106物2010-11-10可吸入颗粒2010-11-91152010-11-8105物可吸入颗粒物可吸入颗粒物2010-11-72010-11-62010-11-52010-11-49485696860可吸入颗粒物可吸入颗粒Ⅱ良Ⅱ良物可吸入颗粒物Ⅱ良可吸入颗粒物2010-11-32010-11-2Ⅱ良良84可吸入颗粒Ⅱ

如不慎侵犯了你的权益，请联系我们告知！物

如不慎侵犯了你的权益，请联系我们告知！可吸入颗粒2010-11-176Ⅱ良物4.3建立模型分析影响城市空气污染程度的主要因素是什么1.通过对10个城市的空气污染情况分析，用excel统计出污染物种类与污染的天数，结果如下表所示：污染种类可吸入颗粒物二氧化硫2532472591671989300成都杭州3北京2上海3广州0拉萨24226225428844240乌鲁木齐郑州武汉0西安统计结果显示，污染物分为可吸入颗粒和二氧化硫两类。而且对比每个城市的数据，结果显示可吸入颗粒物为首要污染物的天数占绝大多数，远远大于二氧化硫。2.影响城市空气污染程度的主要因素首要污染物的来源可吸入颗粒物的来源有以下五个方面：a.工业污染（燃煤排放和生产中产生的尘及有机物）；b.机动车尾气排放；c.来自道路、施工、裸露地面的d.采暖小煤炉及采暖锅炉的民用燃煤排放；e.餐饮、干洗等其他排放。扬尘；二氧化硫的来源主要是煤炭的燃烧通过调查分析，我们即来源a和e可以看作是不量，即民用采暖、汽车尾气、扬尘这三项指标进行空气质量的预报监测。认为在一定时期内，每天的工业污染和餐饮干洗等其他排放，变的。我们只需要根据每天污染源数值的变化3．各个因素的具体分析(1)民用采暖，分为采暖季和非采暖季两类，得到如下表：

如不慎侵犯了你的权益，请联系我们告知！采暖季非采暖季每日平均可吸入SO2颗粒SO2可吸入颗粒污染指63.8104.5数28.399.2采暖季二氧化硫（SO2）每日平均污染指数是非采暖季的63.8/28.3=225.4%；多出35.5/28.3*100%=125.4%；多出63.8-28.3=35.5个指数；采暖季多出的35.5个指数所对应的就是来自采暖小煤炉及采暖锅炉等民用燃煤排放；这35.5个指数所对应的是占KLW来源的一部分：采暖小煤炉及采暖锅炉的民用燃煤排放，也是占采暖季二氧化硫（SO2）每日污染指数总数的质量时，应当分采暖季和非采暖季两个不同时间段。在同等或相似的天气条件下，如果设非采暖季二氧化硫（SO2）污染指数为X，可吸入颗粒物（KLW）污染指数35.5/63.8*100%=55.6%的部分。因此在预报某日空气为Y，则有：时间段SO2污染指数KLW污染指数非采暖季XY采暖季225.4%X110.8%Y(2)汽车尾气北京上海等国际大都市的机动车保有量大，机动车尾气排放污染已成为城市空气质量的主要来源之一。我们从北京市交管局指挥中心得到了一部分反映北京日平均车流量的几组典型数据：工作日日均流2005年线路周六日日均流量量二环路三环路四环路二环路三环路四环路141896125603157628128562120530141337134732126519150758129513120623138044161574513710012342113335212191711499111560912673412408411684912856112161112378814880174月6月9月二环路三环路四环路二环路三环路10月四环路总流量因此，北京地区周六日日均汽车总流量/工作日日均汽车总流量为：1488017/1615745*100%=92.1%即周六日每天的汽车尾气排放量是工作日每天的92.1%.

如不慎侵犯了你的权益，请联系我们告知！

如不慎侵犯了你的权益，请联系我们告知！按照月份划分的尾气排放从访问北京市公安交通管理局指挥中心和新闻中心获得的资料来看，北京市每年的3月份、4月份和9月份的车流量较其他几个月份要大，这三个月之间相比，车流量差异不明显；一年中其余的月份之间差异也不明显。下图的曲线说明2004、2005年北京6-11月份日均车流量的趋势是大体相同的。市公安交通管理局介绍说，从每年数据统计结果来看，尽管同一个月份日均车流量不同，但是月和月之间的比值，每年基本上是一样的。2004、2005年6-11月北京二环路日平均车流量及车速曲线2004、2005年6-11月北京二环路日平均车流量及车速曲线

如不慎侵犯了你的权益，请联系我们告知！(该曲线数据来源：北京市公安交通管理局新闻中心)北京市2005年全月份日均车流量数据表(部分月份)2005年二环路三环路四环路总计3882374160604160873915593801743853424036533837661月1275823月1354094月1404195月1306216月1267028月1266869月13327410月129212114179122501124951121620119309121026125867120943146476158150150717139318134163137630144512133611(上表数字来源：北京市公安交通管理局指挥中心)2005年1、5、6、8、10月五个月份的日平均车流量为：（388237+391559+380174+385342+383766）/5=385815.62005年3、4、9月三（416060+416087+403653）/3=411933.3由上述两个算式的结果进一步得出：385815.6/411933.3×100%=