文档介绍：该【有关大气污染问题的数学建模论文 】是由【guoxiachuanyue009】上传分享，文档一共【29】页，该文档可以免费在线阅读，需要了解更多关于【有关大气污染问题的数学建模论文 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。1/26
问题重述
大气是指包围在地球外围的空气层,是地球自然环境的重要组成部分之一。人类生活在大气里,洁净大气是人类赖于生存的必要条件。一个人在五个星期内不吃饭或5天内不喝水,尚能维持生命,但超过5分钟不呼吸空气,便会死亡。随着地球上人口的急剧增加,人类经济增长的急速增大,地球上的大气污染日趋严重,其影响也日趋深刻,如由于一些有害气体的大量排放,不仅造成局部地区大气的污染,而且影响到全球性的气候变化。因此,加强大气质量的监测和预报是非常必要。目前对大气质量的监测主要是监测大气中SO、NO、悬浮颗粒物22
(主要为PM10)等的浓度,研究表明,城市空气质量好坏与季节与气象条件的关系十分密切。
附件给出城市A、B、C、D、E、F从2003年3月1日至2010年9月14日测量的污染物含量与气象参数的数据。
请运用数学建模的方法对下列问题作出回答:
找出各个城市SO、NO、PM10之间的特点,并将几个城市的空气质量进行
22
排序。
对未来一周即2010年9月15日至9月21日各个城市的SO、NO、PM10
22以与各气象参数作出预测。
分析空气质量与气象参数之间的关系。
就空气质量的控制对相关部门提出你的建议。
问题分析本题为生活中的实际问题,层层递进式提出四个问题,分别需要对空气污染因素以与气象参数进行分析求解。第一问为评价性问题,先从城市内部个污染物特点出发,再到城市之间空气质量进行比较。第二问是预测性问题,通过对给出的数据进行分析,预测各项参数之后的趋势。第三问是寻找关联性问题,要求找出空气质量与气象参数之间的关系。第四问为开放型问题,可通过之前得出的结论或者相关文章与模型提出建议。
问题1
通过查阅资料,运用已有的API对各个城市的各项污染指标进行计算,得出各个污染指数API月平均的折线图,观察,得出各城市各项指标的特点。鉴于求解城市API时有一定的误差,故选择综合评价模型,对数据进行标准化处理
之后,确定动态加权函数,对模型进行求解,排名。检验模型后确定结论的合理
2/26
性。
问题2预测模型主要有灰色预测,时间序列等模型。由所给数据以与问题可知该预测模型为时间序列。随机选取气象参数之一气温(tem)为例进行分析,先通过SPSS软件得到其时序图,观察其走势,对其做平稳化处理。然后以最小BIC为标准,构造模型,进一步应用SPSS软件求解,得出各项参数,并预测出2010年9月15日至2010年9月21日的数据。其余各城市各污染物浓度以与气象参数应用类似方法进行求解。最后,由于F城市所提供数据与需要预测日期相隔较远,故只做出定性的分析预测。
问题3
空气污染物与气象要素关系密切,研究的方向多为相关性分析与回归分析或从理论上描述气象要素对污染物迁移扩散的影响。但是回归分析应用于处理不相关变量之间关系,而典型相关性分析能很好地解决由于变量之间相关而导致回归准确性降低的问题。并且观察原始数据发现,其中只有一组气象参数,故猜测气象参数是在其中某一个城市所采集。现应用典型相关性分析分别分析A、B、C、三城市空气污染物s°2、no2、PM10与气象要素这两组数据间的关系。求出不同季节的相关系数,判定气象参数最有可能是属于哪一城市的。再对该城市进行偏相关性分析,最终得出污染物与气象参数之间的关系。该过程由SPSS直接完成。
问题4依据第三问所求得的气象参数和与其对应城市之间的关系,分析影响各污染物浓度的主要因素,依此对有关部门提出合理的建议,以提高该城市的空气质量。
数据处理
对附件中数据整体浏览,将不合理的数据进行删除:2005年11月7日的
,,依据常识,该两组数
据均为记录错误,故删去不予考虑
问题1
。
问题2
基于数据的不完整性,只选择具有连续性的数据(2010年1月20日至2010年9月20日)对问题二进行分析预测。
,4
将一年分为春季季风季(3-5月)和冬季采暖季(11-2月)两部分,分别进行分析。
3/26
模型基本假设
1、各组数据真实可信,且是在同一地点同一时间采集,不考虑人为因素,具有统计、预测意义。
2、假设A、B、C、D、E、F六个城市的发展状况一样,即发展速度没有明显
差异。
3、API指标真实可靠,所给数据具有参考统计意义。
4、月API平均值能很好的代表该月空气质量,具有比较意义。
5、对F城市进行定性预测时,A、F城市发展状况基本一样,有比较价值。
6、第三问中,灌输数据的对应关系,假设气象参数是在A、B、C三城市中某一个城市所采集。
符号说明
API
空气污染指数
I
某污染物的污染指数
C
该污染物的浓度
CC
大小
在API分区表中最接近C值得两个值
II
大小
在API分区表中最接近I值得两个值
x.(i1,2,3)
i
三项污染指标标准化值
m
J
区间最小值
4/26
Mj
区间最大值
a(j),b(j)
区间边界
kk
w.(x)(j=1,2…,m)j
权重
X
被评价指标的综合评价值
S
i
城市
B(S)
ji
每个城市的Borda数
{x}
t
时间序列
P
模型阶数
d
差分阶数
B
延迟算子
G
均方差
R
复相关系数
模型的建立与求解
问题1
问题一第一部分
通过查阅资料,可以找到API,即空气质量污染指数标准,由此计算每个城市各项指标的月API平均值,对各项数值进行比较,得出各个城市三项指标特点。
各项指标月API平均值求解
根据表一:
AFIFollullonInd盟的英文缩写}是空吒污染抬数,我国城市空气质量日报占FI分级杨准如表1:
5/26
表1空气W亏染J旨数对应的污粢物浓度限值
污染指数
污染物浓度(毫豆7立方米》
API
〔日均值)
叫
<日均值)
P怖{日均值〉
co
〔小时均苗■)
°3(小时均道、
□D
5
100
10
?CC
60
300
50
4aa
120
5Q0
<
150
计算各项指标的API值:
设I为某污染物的污染指数,C为该污染物的浓度。则:
I-1
I=亠小(C-C)+1
C-C小小
大小
式中:CC:在API分区表中最接近C值得两个值
大小
II:在API分区表中最接近I值得两个值
大小
AFIFollullonInd盟的英文缩写}是空吒污染抬数,我国城市空气质量日报占FI分级杨准如表1:
6/26
U城市
D城弔
E城市
F城市
7/26
图一:各城市各项指标月API平均值折线图
整体分析图表可以看出A、B、C、D、E五个城市SO2、NO2、PM10等污染物浓度均呈现波动性并且有缓慢下降趋势。
分析A城市数据,发现A城市PM10浓度与B城市差别并不显著,但是观察发现A城市PM10的值在2010年8月后有所回升,这一点也可由数据得到验证。
B城市SO2波动性很强,但是下降的趋势并不是非常的明显,说明B城市可能有一些周期性的污染源需要治理。而B城市的PM10波动性强有明显下降趋势,这说明B城市很有可能在2010年采取过一些相应的积极措施,使得该城市PM10浓度在短期内大幅度下降。
而A、B两城市的SO2和PM10数值均明显高于NO2的数值,且两城市污染物的波动方式相似,可粗略认为AB两城市有部分工业或者结构上的相似。
分析C、D两城市可知SO2、NO2、PM10浓度较平稳波动,只有PM10在个别时段有较大的起伏,而在其他时间序列内均趋于平缓变化。C、D两城市的PM10曲线在同一时间明显偏高,可推论在那一段时间有某些外界因素使得两个城市的PM10数值共同上升。
分析E城市空气污染物浓度可知,E城市SO2、NO2、PM10浓度均在一定X围内平稳变化,说明该城市在所选时间段内空气质量比较平稳。
由于F城市数据严重不足,只有从2004年9月1日到2009年12月27日的采集数据,故在F城市数具有统计意义的前提下,由图可知观看出F城月平均污染物浓度大致呈现平稳趋势。
问题一第二部分
根据问题对API分析发现,对于城市API值计算中,原理为取三项指标的最大值,这会造成相应的误差,故在分析第一问的第二部分时,只参考API的划分标准,应用综合评价模型。
数据的标准化处理
对所给的空气污染标准(API)进行标准化处理,记三项指标:SO、NO、
22
PM10的数值分别为x,x,x。
123
三项指标的数据均为极小型指标(即指标值越小越好),对其指标Xj做标准化处理,即令:
-m
匚(1<j<m)
-m
jj
其中m=min{x},M=max{x}。则相应的指标值变为{x'}e[0,1],即为无量j1<i<n可j1<i<nijij
纲的标准化指标,对应的分类区间[a(j),b(j))也随之相应的变化,在这里为了方便
kk
仍记为[a(j),b(j))(k=1,2…K;1<j<m)。
kk
(1)SO2的标准化
8/26
取m=0,M=,x、=二,则其标准化数据为:x(k)'e[0,1]
对应的分类区间为:
(0,],(,],(,],(,],(,],(,1],(1妙)
(2)NO2的标准化
取m=0,M=,x'二务,则其标准化数据为:x(k)'e[0,1]
对应的分类区间为:
(0,],(,],(,],(,],(,],(,1],(1,^)
PM10的标准化
取m=0,M=,x'二厶,则其标准化数据为:x(k)'e[0,1]
对应的分类区间为:
(0,],(,],(,],(,],
(,],(,1],(1,a)
动态加权函数确定
根据这一实际问题,通过对SO空、NO2、PM1O三项指标的变化关于空气质
量的分析,可得其变化的规律为:先是缓慢增长,中间有一个快速增长的过程,最后平缓增加趋于最大值。此增长规律可取动态加权函数为偏大型正态分布函数,即:
x—Q
w(x)
j
1-e"Q「,当x>a时
0,当x<a时
j
二一(b(j)—a(j)),Q
211j
其中Q不妨取指标x的第一类空气质量标准的中间值,即«
jj
由w(q(j))=(1<j<m)确定。j4
综合评价模型的构建
根据标准化后的评价值,不妨仍用x表示,以与相应的动态加权函数
i
w.(x)(j=1,2…,m),建立综合评价模型来对被评价的6个城市的空气质量进行
j
9/26
评价,在此,取综合评价模型为个评价指标的动态加权和,即:
X=£w(x)-x
jjj
j=1
其函数值X为被评价对象的综合指标值。
求出权后,可将6个城市的三项指标求期望,定量地得出每个城市中三项
指标的权值。
利用附件中给出的31个月的较为完整的数据,计算可得ABCDE五个城市
的空气质量评价性指标,即可得到一个综合评价矩阵(X),其结果如下:
ij5x31
A
B
C
D
E
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
表二:ABCDE五个城市综合评价矩阵
同时,利用附件中给出的4个月(,12)的数据,经计算可得
ABCDF六个城市的空气质量评价指标,得到矩阵(X),结果如下。
ij5x4
10/26
AB
C
D
F
表三:ABCDF五个城市综合评价矩阵
根据上表和表中的数据,根据其大小(即反映空气质量的高低程度)进行排
序,数值越大,说明其空气质量越差。编写C语言程序,对其进行排序。排序结果见附录
利用决策分析中的Borda函数方法来确定综合排序方法,记在第j个排序方案中排在第i个城市S后面的站点个数为B(S),则城市S的Borda函数为
ijii
B(S)=£B(S)(i=1,2・・・5)
iji
j=1
经计算,各城市Borda数与总排名如下
10/26
10/26