数学建模： 城市空气质量评估及预测(省级优秀奖).docx

编号：时间：2021年x月x日书山有路勤为径，学海无涯苦作舟页码：第1页 共1页2010年西南交通大学“新秀杯”数学建模竞赛参赛论文论文题目：城市空气质量评估及预测 参加2011年建模创新团队选拔： 是 （“是” 或 “否”）2010年11月城市空气质量评估及预测 摘要：本文对我国的成都，杭州，北京，上海，广州，拉萨，乌鲁木齐，郑州，武汉，西安等10座城市的每日空气质量详细列表进行科学分析，利用层次分析法和指数平滑法等数学建模方法对其空气质量进行研究，综合考虑各种因素建立如下数学模型1、十个城市空气污染严重程度的科学排名本文采用层次分析法来对10个城市的环境污染情况进行比较用excel统计出各个城市不同级别的污染天数在采用层次分析法，首先建立层次结构模型然后构造成对比较阵，用matlab软件计算权向量并做一致性检验，计算组合权向量并做组合一致性检验最后根据权重的大小就可以比较出10个城市空气污染严重程度的科学排名2、成都市11月的空气质量状况预测本文采用指数平滑的方法来建立数学模型根据前面的统计数据，考虑到污染级别不同和首要污染物的种类两个因素来预测成都11月的空气质量状况。

3、分析影响城市空气污染程度的主要因素先通过excel统计出每个城市的各种首要污染物所占总天数的比例，然后再综合考虑各方面因素找出造成该污染物超标的以原因，以找出影响城市空气污染程度的主要因素最后本文就10个城市的空气污染严重程度的排名给出了相关的分析以及应对策略就成都11月份的空气质量状况预测给出出行和生活方面的建议并结合当下倡导建设环境友好型和资源节约型社会出发，就如何兼顾经济发展与环境保护给出指导性建议关键词：空气质量污染等级 层次分析法 指数平滑法 首要空气污染物一 问题的提出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 有如下关系，如表n12345678910RI000.580.901.121.241.321.411.451.49(5)权重计算方法计算矩阵A 的特征根及特征向量，将所求的特征向量单位化后得到的就是权重值。

2. 成都11月份空气质量进行预测 指数平滑法（Exponential Smoothing，ES）是布朗(Robert G..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/31 5/45/35/25四级4/154/514/324三级1/53/53/413/23二级2/152/51/22/314一级1/151/51/41/31/21由表中数据， 计算可知： λ(max) = 5.00 ， CI = 0.00 ， RI = 1.12 ，CR = 0.00 < 0.1。

因为CR = 0.00 < 0.1，所以此排序有满意的一致性3.给出对象层对准则层的各个因素的判断矩阵并进行分析由于各个城市只存在污染程度的不同，所以它们两者之间各因素之间的关系2010空气污染指数统计优良轻度污染中度污染严重污染成都492243300杭州502203002北京352015921上海12715614701广州951931000拉萨21291301乌鲁木齐192166651郑州152493611武汉522054802西安142444301成都杭州 北京上海广州拉萨乌鲁木齐郑州武汉西安权重成都149/507/549/12749/9549/21249/1949/1549/527/20.073杭州50/49110/750/12725/9525/10650/1910/325/2625/70.074北京5/77/10135/12735/9535/21235/197/335/525/20.052上海127/49127/50127/351127/95127/212127/19127/15127/52127/140.190广州95/4919/1019/795/127195/212519/395/5295/140.142拉萨212/49106/25212/35212/127212/951212/19212/1553/13106/70.317乌鲁木齐19/4919/5019/3519/1271/519/212119/1519/5219/140.028郑州15/493/103/715/1273/1915/21215/19115/5215/140.022武汉52/4926/2552/3552/12752/9513/5352/1952/15126/70.077西安2/77/252/514/12714/957/10614/1914/157/2610.021根据表中数据，类比（2）中方法，计算出各种不同污染等级对不同城市的权重空气质量一级对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.1 4. 模型求解进行层次总排序，方法： 将上面5个空气质量等级对10个城市的不同权重表单位化后作为列向。