数学建模大气污染预报问题.docx

资源描述

数学建模大气污染预报问题.docx

《数学建模大气污染预报问题.docx》由会员分享，可在线阅读，更多相关《数学建模大气污染预报问题.docx（25页珍藏版）》请在冰点文库上搜索。

数学建模大气污染预报问题.docx

数学建模大气污染预报问题

学生数学建模竞赛第一次预选赛

一、（必做题）

（1）油罐的体积（本题10分）

一平放的椭圆柱体形状的油罐，长度为L,椭圆的长半轴为a,短半轴为b,油的密度为p，问当油罐中油的高度为h时油量是多少？

解：

由题意可话画出画出几何图形如图1所示

图1.1

x=acost

椭圆方程为』=bsint如图2,设阴影部分面积为S/2，则油桶的底面积为S

图2

下面将会利用mathematics5.0软件进行求解,求解的程序如下:

Integrate[2*a*b*Cos[tF2,{t,ArcSin[1-h/b],Pi/2}]

解得结果为：

S=a（.（2b；h）h（_b+h）+barccos（1_h））

Xbb

当hb时，由椭圆对称性，A中的h用2b-h代替得到：

所以油液质量M为：

（2）光的反射定律（本题10分）

费马原理：

光总是沿用时最短的光程传播。

试根据这一原理利用极值的有关知识证明光的反射定律：

入射角等于反射角。

解：

由于光在同一介质中的速度为常数，所以在同一介质中光总是沿直线传播。

如图3,现假设有两种介质1、2相接，光线在介质1中的传播速度为v,取两介质的分界线上的一条直线为X轴，设有一束光线从介质1中的A（O,a）点经X轴上的P（x,O）点反射，并沿直线方向行进到B（d,b）点。

设直线AP与X轴法线的夹角为比，PB直线与X轴法线的夹角为屯，下面，根据最短时间效应来推导出光学中的反射定理。

图三

光线由A点传到P点所需的时间为：

光线由P点传到B点所需的时间为：

故光线由A传到B所需的总时间为：

根据费马定理，最短时间效应对应的优化问题为:

于是可以得到:

J—X厂b2（d-x）2'0

a2x2

d-x

.b2（d-x）2

又由于

a2x2

d-x

.b2（d-x）2

=sinJ2

证毕

所以有：

这就是光学中的反射定理

大气污染预报问题

摘要

本文通过对四个城市的空气质量的排名以及城市A的空气质量，利用C语言、Excel、

Mathematics和MATLA等工具，分别建立了层次模型、多元线性回归预测模型进行了合理地分析。

最后，我得到了以下一些比较满意的结果。

问题

（1）：

通过对问题

（1）问题的分析，得出了这是一个比较典型的层次模型，目标层是空

气质量的排名，因素是三种污染物的浓度情况，对象是题目给出的4个城市。

查找资料后，我找到了一个非常关键的东西一一空气污染指数的计算方法，于是利用C语言的编程知识我很快求出了这些城市的污染情况，结合层次模型的相关知识，建立层次模型后很快得到了我们所期望的答案：

总

权重：

A-0.1374,B-0.1301，C-0.1028，D-0.6298.于是城市环境排名也就解决了，由优到劣的排名情况：

C、B、A、D。

而且层次的模型的一次性检验也顺利通过。

问题

（2）：

问题

（2）要求我们找出空气质量与气象因素之间的关系，一开始查阅了很多资料，本想借助灰色预测模型进行求解，可是灰色预测模型的使用条件和咱们这个题目的要求似乎没有什么关联，后来在网上浏览一片文章的时候，我找到了问题的突破口，便是利用多元回归预测模型进行求解。

然后根据这个模型的所要的处理数据，利用MATLABExcel等工具，求解到本题的回归系数。

得到了三种污染物与气象因子之间的关系：

SO2的浓度与气象因素的关系：

q=1.0173-0.0015p0.0019-0.0004f-0.0044V

NO2的浓度与气象因素的关系:

c2--0.20220.0003p0.0017t0.0001f-0.0077V

PM10的浓度与气象因素的关系：

C3=0.9329-0.0014p0.0033:

-0.0022f0.0019/

接着，我又利用了F检验和复相关系数R用来判别回归方程在统计上是否合理。

结果还是很让人满

意的，回归模型的拟合度还是很高的。

关键词：

层次模型多元线性回归预测模型空气污染指数

问题提出

大气污染预报问题

大气是指包围在地球外围的空气层，是地球自然环境的重要组成部分之一。

人类生活在大气里，洁净大气是人类赖于生存的必要条件。

一个人在五个星期内不吃饭或5天内不喝水，尚能维持生命，但超过5分钟不呼吸空气，便会死亡。

随着地球上人口的急剧增加，人类经济增长的急速增大，地

球上的大气污染日趋严重，其影响也日趋深刻，如由于一些有害气体的大量排放，不仅造成局部地区大气的污染，而且影响到全球性的气候变化。

因此，加强大气质量的监测和预报是非常必要。

目前对大气质量的监测主要是监测大气中SO2、NO2、悬浮颗粒物（主要为PM10等的浓度。

附件一给出了城市AB、C、D从2009年6月1日至2009年7月25日测量的污染物含量及城市A的气象参数的数据；附件二给出了城市A从2009年7月26日至2009年7月30日测量的污染物含量及气象参数的数据。

请解决下面两个问题：

（1）建立由污染物浓度评价空气质量的数学模型，然后利用附件一中的数据对四个城市的空气质量进行排序。

（2）分析城市A的空气质量（指SO2、NO2、PM10勺浓度）与气象参数之间的关系，并利用附件二中的数据进行检验。

二、基本假设

1、题目所给的四个城市的污染物含量及城市A的气象参数等数据都准确可靠。

2、根据内地空气污染指数（API）来划分为5个等级：

API值小于等于50,空气质量为优，相当于国家空气质量一级标准；API值大于50且小于等于100,表明空气质量良好，相当于达到国家质量二级标准；API值大于100且小于等于200,表明空气质量为轻度污染，相当于国家空气质量三级标准；API值大于200表明空气质量差，称之为中度污染，为国家空气质量四级标准；API大于300表明空气质量极差，已严重污染。

三、符号说明

符号

意义

备注

入（max）

对角矩阵的最大特征值

矩阵的阶数

层次模型的一致性指标

CI=（入-n）/（n-1）

随即一致性指标

一次性比率

CR=Ki

大气压强的数值

单位：

mmhg

温度的数值

空气湿度的数值

风速的数值

m/s

F检验统计量

预测模型的复相关系数

四、问题分析

1、问题

（1）的分析:

要对A、BCD四个城市的空气质量进行排序，可从题目的要求中获知利用污染物浓度来进

行四个城市的排名。

经过分析和查阅相关资料，这个问题应该属于典型的层次模型的运用。

下面将

层次模型的相关内容说明如下：

（1）最大特征值入（max）的MATLA计算方法：

[V,D]=eig（A）,其中A为待计算特征值的矩阵，D为对角矩阵，其对角元素为A的特征值，最大的即为入（max）0

（2）一致性指标CI计算方法：

CI=（入-n）/（n-1）;（其中入为矩阵A的最大特征值，n为矩阵的阶数）⑶随即一致性指标RI的计算方法：

RI与n有如下关系，如表

0.58

0.90

1.12

1.24

1.32

1.41

1.45

1.49

（4）权重计算方法

计算矩阵A的特征根及特征向量，将所求的特征向量单位化后得到的就是权重值。

2、问题

（2）的分析：

通过仔细分析题目的要求，得知题目要求我们找出空气质量与气象因素的之间的的关系。

于是

我首先想到了用预测模型去处理，但是由于变量太多而且，处理起来的拟合度太低了，达不到我们满意的要求。

后来，通过仔细阅读相关资料找到了以个比较好的模型一一多元线性回归预测模型，去处理。

五、模型的建立与求解

5.1问题一模型建立与求解

5.1.1问题一的分析

要得到城市的空气质量的排名，首先我们得找一个量去权衡它们的关系才能得出比较理想的结果，在前面的假设中我们便得到了，一个空气污染指数，我们可以以这个关键因素作为突破口求解。

5.1.2问题一模型的建立

1、将研究目标（Z）、因素（P）、对象（C）按相关关系分成目标层Z、准则层P、对象层Co层次结构图如图所示：

空气质量的排名

城市

2、给出空气质量一级，二级，三级两两成对比较的判断矩阵P

污染级别

一级

二级

三级

根据上图得出如下例两两成对比较的判断矩阵P

一级

二级

三级

权重

一级

1/2

1/3

0.500

二级

2/3

0.333

三级

3/2

0.167

由表中数据，计算可知：

入（max）=3.00，CI=0.00，RI=0.58，CR=0.00<0.1。

因为CR

=0.00<0.1，所以此排序有满意的一致性。

3、给出对象层对准则层的各个因素的判断矩阵并进行分析。

由于各个城市只存在污染程度的不同，所以它们两者之间各因素之间的关系。

在这里我们利用了C语言的相关知识求解出了这55天中各个城市的空气污染指数:

四个城市的空气污染指数统计（单位：

天）

优

良

轻微污染

轻度污染

中度污染

重污染

根据表中数据，类比

（2）中方法，计算出各种不同污染等级对不同城市的权重

权重

11/21

11/48

11/16

0.115

21/11

21/48

21/16

0.219

48/11

48/21

48/16

0.500

16/11

16/39

16/48

0.167

由表中数据，计算可知：

入（max）=4.00,CI=0.00,RI=0.90,CR=0.00<0.1

所以此排序有满意的一致性。

同理，可以计算出其余空气质量等级4个城市的不同权重。

计算方法类似，用MATLAB^件的计算过程详见本文附表。

空气质量“良”级对4个城市的不同权重表

城市

权重

0.355

0.281

0.058

0.306

由表中数据，计算可知：

入（max）=4.00，CI=0.00，RI=0.90，CR=0.00<0.1

空气质量“轻微污染”级对3个城市的不同权重表

城市

权重

0.000

1.000

由表中数据，计算可知：

入（max）=1.00，CI=0.00，RI=0.00，CR=0.00<0.1

由于其他的污染指数均为零，在这里不再考虑了。

5.1.3问题一模型的求解

进行层次总排序，方法：

将上面3个空气质量等级对4个城市的不同权重表单位化后作为列向量构成4X3矩阵，和空

气质量一级，二级，三级两两成对比较的判断矩阵P相乘，结果便得到4个城市的权重值。

根据上述问题的分析中的假设可知，权重值越大，表明空气污染情况越严重。

因此，将4个城市的权重值,

按照从小到大依次排序，得出的结果便是4个城市的空气污染严重程度的排名

最终结果如下表所示：

优（0.167）

良（0.333）

轻微污染（0.500）

总权重

0.115

0.350

0.000

0.1374

0.219

0.276

0.000

0.1301

0.500

0.057

0.000

0.1028

0.167

0.317

1.000

0.6298

根据4个城市的总权重值进行从小到大依次排序，空气污染严重程度的排名如下：

C、B、A、D

5.1.4问题二结果的分析及验证

总的一致性检验：

CR=0.167X0+0.333X0+0.500X0=0<<0.1。

此结果说明排序结有非常满意的一致性。

结论显示城市C的空气质量状况最好，而D的空气质量状况最差。

而且我们把A和D故个比较可以看出：

虽然D的“优级”天数比A还多，但是由于D出现了两个轻度污染而造成D的污染权重显著增加。

5.2问题二模型建立与求解

5.2.1问题二的分析

我们首先利用Excel对SO2与各个气象因子之间的关系图如下面四个图所示：

由于SO2的浓度高低并不是由单一因素决定的，而是由于大气压强、地面风速、温度以及湿度等气象因素共同影响的结果。

因此，可以建立多元回归预测模型，对其浓度变化规律进行分析和预测，从而实现对可吸入颗粒物（PM。

）浓度的最优控制

5.2.2问题二模型的建立

当前，对于大气污染物浓度预测所采取的方法主要是从污染物排放量高低为基础进行预测的，典型的预测模型有：

箱式模型、高斯扩散模式、多源扩散模式、线源扩散模式、面源扩散模式和总悬浮微粒扩散模式。

随着灰色系统、模糊数学和人工神经网络的发展，预测方法又出现了以污染物排放相关因素为基础的模型，如：

灰色预测模型（GM卜多元统计分析理论、模糊识别方法和人工神经网络预测方法。

本题研究主要采用多元统计的方法进行分析。

在许多实际问题中，影响结果y的因素往往不止一个，而是多个变量X1，X2，•…，xp与y之间存在着如下线性关系：

+匕冷+陽12+…+PpX1p+E

（1）

E；=0,D（；）=；

BBBB,、

其中：

01，2，...，p是回归系数；x1，x2，•…，xp是p个可以精确测量或控制的变量，

及回归因子；；是不可观测的随机误差，满足

E；=0,D（；）=；「＜：

（二未知）

一般地，我们称由公

（1）和

（2）确定的模型为多元线性回归模型，记为：

y-0：

1捲1•：

2捲2…n’ipXip-；

2E；=0,D（；）=

（3）

具体方法为：

（1）计算各变量的平均值：

-1n,.

XiXti,yyt

ntd

（4）

⑵根据公式⑸计算出矩阵Lj和矩阵Li:

Lj='（冷—xJ（Xtj-Xj）

v（i,j=12HI,P）

Liy二'（xti-x）（yt-y）

（5）

（3）

根据公式（6）求出回归系数的估计值：

=L「Liy,o=y-为1—X2：

2--Xp:

p（6）

即可求出回归模型：

根据本题的特点，可以得到这样一个模型：

5.2.3问题二模型的求解

根据多元线性回归法的基本理论，分别考虑大气压强、温度、湿度和地面平均风速4个自变量,自变量分别以p、t、f、v表示，变量用C1表示，即武汉城区吸入颗粒物（PM。

）浓度，mg/m。

贝可设数学模型为：

以环境空气质量自动监测子站监测的城区可吸入颗粒物（PM。

）浓度数据，和相应的地面平均风速、气温、相对湿度3个气象因子为原始数据，先根据公式（4）利用Excel计算出各变量的平均值：

再按公式⑸利用MATLA计算出Lj和Liy:

‘89004

86.94

-1950

2315

二0.4879

8694

46677

-11249

-6.21

1.2595

Liy=<

|-1950

-11249

19001

-1740

，

-6.3447

[23.15

-6.21

-1740

10.78

-0.0219

Lij二

最后根据公式（6）计算出回归系数的估计值：

故根据多元线性回归方法，建立的城市A的S02的浓度拟合模型为:

其中：

ci为S02的预测浓度，mg/m；p为大气压强，mmhgt为地面温度，C；f为近地面空气中的湿度，%v为地面平均风速，m/s。

利用上面类似的方法可以求到：

城市A的NO2的浓度的拟合模型为：

其中：

C2为N02浓度，mg/m;p为大气压强，mmhgt为地面温度，C；f为近地面空气中的湿度，%v为地面平均风速，m/s。

②城市A的PM1（的浓度的拟合模型为

其中：

c3为PM10测浓度，mg/m；p为大气压强，mmhgt为地面温度，°C；f为近地面空气中的湿度，%v为地面平均风速，m/s。

5.2.4问题二结果的分析及验证1、首先利用Excel做出预测值与实际值之间的折线图:

2、利用附表二中的数据结合Excel表格进行检验:

检验结果如下表所示：

SO2

NO2

PM10

预测值

实际值

预测值

实际值

预测值

实际值

0.03142

0.031

0.031941

0.037

0.062343

0.047

0.0311

0.021

0.02869

0.022

0.065379

0.030

0.01712

0.025

0.030037

0.034

0.046869

0.034

0.022236

0.024

0.032847

0.035

0.053687

0.035

0.012537

0.026

0.033138

0.033

0.05011

0.081

上面的见表格中：

我们可以看出预测值和实际值之间还是比较吻合的。

特别是SO2与NO2的预测值

和实际值之间还是很吻合的，只有PM10勺值稍差了一点。

3、下面分别用F检验和复相关系数R用来判别回归方程在统计上是否合理。

F检验统计量F的计算公式见式（7）:

n-m-1

其中，m为回归变量的自由度，

见公式（8）:

n为观察值的组数，回归平方和U和残差平方和Q的计算公式

U八（O2

（8）

Q八（％-？

）2

（9）

复相关系数R的计算公式见式（9）:

其中，回归平方和U和残差平方和Q的计算公式见公式（8）。

1）S02函数关系的检验：

选择所建预测模型的显著性水平为0.05，而F检验的统计了F=17.18>Fo.°5，预测模型在统计意义上是显著成立的。

预测模型的复相关系数R2为0.9357，表明SO2浓度与气象因子（p、t、f、v）之间的关系为高度正相关。

预测模型的标准误差由相关表达式计算得0.0139，因此，表明预测模型的拟合程度很高。

2）各个污染物与气象参数之间关系式的检测情况表：

三个个污染物的检验情况表