数学建模论文
A题 空气质量评价
摘要
本文主要研究空气质量评价的相关问题，为突出改进之后的模型中的实时特性而对数据做了必要的省略处理，然后在现有的国家最新空气污染物监测标准（HJ633-2012环境空气质量指数(AQI)技术规定）的基础上利用半集均方差原理对现有空气质量计算模型进行改进。在论证修正后模型可行性的基础上再对模型加以优化，最后利用优化后的模型对附表二中的各项监测结果得出其空气质量指数。
针对问题一，由于目标模型十分强调实时性，于是把附表一中臭氧8小时平均值﹑细颗粒物24小时平均值﹑可吸入颗粒物24小时平均值做了必要的省略处理。联系实际分析论证了现有模型的局限性，并在此基础上采用半集均方差原理对现有模型进行改进，结果顺利得到优化后的计算模型。
针对问题二，考虑到优化后的计算模型并没有对不同的污染物的危害做出差异化的评价，而是直接取表中所有污染物的AQI平均值进行分析。所以引入层次分析法根据污染物的危害性对不同的污染物赋予相应的权重，对半集均方差公式进行合理修正，最后得到修正后的空气质量计算模型。再代入附表二中的数据即得到各个观测点的空气质量指数。详细的matlab实现程序见附录二。
【关键词】 一维插值 半集均方差 层次分析加权法 优化后的半集均方差
1 问题重述
空气质量指数（AQI）是定量描述空气质量状况的无量纲指数。其数值越大、级别和类别越高，说明空气污染状况越严重，对人体的健康危害也就越大。
空气质量指数实时报一般是发布每个每一整点时刻的空气质量指数。 实时报的指标包括二氧化硫(SO2)、氧化碳(CO)、二氧化氮 (NO2)、臭氧(O3)1小时平均值、臭氧(O3)8小时平均值、一颗粒物(粒径小于等于10μm)、细颗粒物(粒径小于等于2.5μm)的1小时平均值和24小时平均值共计9个指标。福建1中列出了某地区11个城市过去7个时刻的空质量指标取值和相应的空气质量指数。
(1) 建立一种新的空气质量指数计算模型，并比较与现有计算模型的区别。 (2) 利用新的计算模型计算附件2中各个观测点的空气质量指数。
2 基本假设
（1）附表一和附表二中的数据是利用统一的污染物监测仪器并按照统一的测量方法测量得到的。
（2）附表一中的原有的空气质量指数（AQI）是按照国家最新出台的统一标准（HJ633-2012环境空气质量指数(AQI)技术规定）进行计算的。
（3）由于国家最新出台的标准中并没有PM2.5和PM10一小时平均浓度限值，所以计算时采用PM2.5和PM10二十四小时平均浓度限值近似代替。
（4）观测点的测量仪器所测量的不同种污染物浓度之间相互独立，互不影响。
（5）所测量的各个观测点附近的空气污染程度在测量的时刻较为稳定，不发生剧烈变化。
（6）在研究各种指标集对某物影响的过程中，不仅指标集中的最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。
（7）大气中各种污染物对环境和人类的危害程度是不一样的。
3 符号说明
IAQIp 污染物项目P的空气质量分指数；
P 污染物项目P的质量浓度值；
BPHi 表1中与Cp相近的污染物浓度限值的高值位； BPLo 表1中与Cp相近的污染物浓度限值的低值位； IAQIHi 表1中与BPHi对应的空气质量分指数； IAQILo 表1中与BPLo对应的空气质量分指数；
IAQI 空气质量分数 n 污染物项目 Sh 半集均方差；
Pi 大于中位数半集的分指数；
m 大于中位数半集的分指数个数；
平均数；
Q 对各个Pi加权后得到的综合指数；
Sh 修正后的半集均方差；
，
4 问题分析
4.1 问题一的分析
要建立新的空气质量指数计算模型，就必须先明确现有计算模型的不足之处。我国现有空气质量计算模型主要依据《HJ633-2012环境空气质量指数(AQI)技术规定》里的一维插值分析法。
我们以附表一中第一组数据为例进行分析：
（1）先对照下表中各种污染物浓度限值得到相关数值
表1 空气质量分指数及对应的污染物项目浓度限值
（2）把得到的数据代入以下的一维插值公式，分别进行计算
（3）得到空气中各污染物的AQI，并把11个城市的各种污染物AQI绘制成折线统计图。
如下图所示：
图1 某时刻11个城市各污染物所对应的AQI的折线图
由图可知，一般而言，对于某个特定的城市，它的PM10或PM2.5的AQI值远超过其余的各种污染物的AQI值。
而根据我国的最新标准，观测点取各种污染物所对应AQI的最大值来评价整体的空气质量。公式如下：
这就出现了一个问题：在现阶段的测量中，实际上只有PM10或PM2.5的AQI对最终的空气质量AQI有影响。
而在实际生活中，在研究各种指标集对某物影响的过程中，不仅指标集中的
最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。
现有的评价模型在实际应用中正是没有考虑到污染物次大值等的对空气总体质量的影响，不免有些以偏概全。为了解决这一问题，我们引入了半集均方差原理来构建新的空气质量指数计算模型。
4.2 问题二的分析
我们先查找出各种大气污染物对环境和人类有哪些方面的的危害：
（1）二氧化硫的主要危害：形成工业烟雾，高浓度时使人呼吸困难，是著名的伦敦烟雾事件的元凶；进入大气层后，氧化为硫酸（SO4）在云中形成酸雨，对建筑、森林、湖泊、土壤危害大；形成悬浮颗粒物，又称气溶胶，随着人的呼吸进入肺部，对肺有直接损伤作用。
（2）一氧化碳的主要危害：极易与血液中运载氧的血红蛋白结合， 结合速度比氧气快250倍，因此，在极低浓度时就能使人或动物遭到缺氧性伤害。轻者眩晕，头疼，重者脑细胞受到永久性损伤，甚至窒息死亡；对心脏病、贫血和呼吸道疾病的患者伤害性大；引起胎儿生长受损和智力低下。
（3）臭氧的主要危害：低空臭氧是一种最强的氧化剂， 能够与几乎所有的生物物质产生反应，浓度很低时就能损坏橡胶、油漆、织物等材料；臭氧对植物的影响很大。浓度很低时就能减缓植物生长，高浓度时杀死叶片组织，致使整个叶片枯死，最终引起植物死亡，比如高速公路沿线的树木死亡就被分析与臭氧有关；臭氧对于动物和人类有多种伤害作用，特别是伤害眼睛和呼吸系统，加重哮喘类过敏症。
（4）二氧化氮的主要危害：刺激人的眼，鼻，喉和肺， 增加病毒感染的发病率，例如引起导致支气管炎和肺炎的流行性感冒，诱发肺细胞癌变；形成城市的烟雾，影响可见度；破坏树叶的组织，抑制植物生长；在空中形成硝酸小滴，产生酸雨。
（5）PM10和PM2.5的主要危害：随呼吸进入肺，可沉积于肺，引起呼吸系统的疾病。颗粒物上容易附着多种有害物质，有些有致癌性，有些会诱发花粉过敏症；沉积在绿色植物叶面，干扰植物吸收阳光和二氧化碳和放出氧气和水分的过程， 从而影响植物的健康和生长；厚重的颗粒物浓度会影响动物的呼吸系统；杀伤微生物，引起食物链改变，进而影响整个生态系统；遮挡阳光而可能改变气候，这也会影响生态系统。
根据以上资料，可知各种污染物对生态系统危害的原理和方式是不尽相同的，那么，各种大气污染物对环境和人类的危害程度也是不一样的。因此各种污染物AQI的值在整个空气质量评价体系中所占的权重也是不一样的，并不能简单的利用平均值来反映污染的总体水平。
于是我们对问题一中得到的新的空气质量指数计算模型进行必要的加权修正，在比较了专家咨询加权法、相关系数加权法、主成分分析加权法和层次分析加权法的优缺点之后，我们选择了层次分析加权法得到各个污染物AQI所占的权重。并计算得到了优于平均值的加权指数来反映污染的总体水平。从而也得到了“优化后的半集均方差”（将公式中的平均数用加权指数替代）公式。
5 模型的建立与求解
5.1 问题一的模型建立与求解
首先进行数据处理，由于要建立的新空气质量指数模型十分强调实时性，于是把附表一和附表二中臭氧8小时平均值﹑细颗粒物24小时平均值﹑可吸入颗粒物24小时平均值做了必要的省略处理。且国家最新出台的标准中并没有PM2.5和PM10一小时平均浓度限值，所以计算时采用PM2.5和PM10二十四小时平均浓度限值近似代替。
按照数理统计的基本概念，用来评价空气质量的分指数集合，构成了容量为n的样本。而空气质量指数则可认为是反映空气质量特性的样本数字特征，它主要应包括样本观测值的集中趋势及其离散程度两个方面。
（1）算术均值（亦即一阶原点矩）是典型的也是基本的集中趋势测度，其定义式为：
=
∑X
i
i
n
;i=1,2,...n
式中， Xi为样本观测值；n为观测值个数。
算术均值不仅是度量数据离散程度的有效基准,还可以反映样本的总体水平。因此, 算术均值是空气质量评价中的一个重要指标。
（2）样本观测值的离散程度, 最基本的度量是用标准差（或称均方差）,它是数据能变性测度，这里取定义式为：
σ=
∑X
i
i
-2
n
；
i=1,2,...n
由标准差的定义式可以推论：Xi中最大值的增大或最小值的减小,均可导致标准差值的增大。联系到空气质量评价, 也就是把远大于算术均值的分指数和远小于算术均值的分指数置于等同地位, 这显然与它们的环境效应相矛盾。因为分指数中最大值对环境质量的影响是明显的, 而最小值的影响较小甚至可以忽略。因此, 在空气质量评价中若简单地按定义式引入标准差, 尽管可以说明分指数的离散程度, 却未能真实地反映出各自对水质的影响，因而不能满足空气质量评价的需要。因此,有必要从能够充分反映的分指数对空气质量影响的特点出发, 寻求满足于空气质量评价要求的数字特征。
（3）在研究各种指标集对某物影响的过程中，不仅指标集中的最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。因此，为了在评价体系中引入次大值等的影响。提出了用中位数来划分整个评价指令集，对大于中位数的分指数半集给予较大的权重。
联系到在不同的分布形态该半集数列的疏密程度不同,可通过计算它们对算术均值的均方差,以反映大于中位数半集的离散程度。据此引出一个新的数字特征,可将它定义为半集均方差, 其定义式为：
Sh=
∑(P-)
i
i=1
m
2
m
（4）根据上述分析, 空气质量指数AQI作为反映大气环境质量的数字特征,应在反映分指数集合集中趋势的基础上,用大于中位数半集的离散程度加以调整。据此可建立新的空气质量指数（AQI）的数学模型为：
AQI=P+Sh
该模型可简称为半集均方差模型，该模型既通过算术均值考虑了每项分指数对空气质量的影响, 也通过半集均方差对分指数中的大值给予较大的权重，据此, 该模型是按各参数对大气环境影响程度的不同,综合出大气环境质量状况的数量描述。
5.2 问题二的模型建立与求解
由问题一我们得到新的空气质量指数计算模型，但是在实际应用的过程中又发现了半集均方差模型的缺陷。
由图1中的折线图可知：中国城市里PM2.5或PM10的AQI值远超于其他污染物的AQI值。所以导致了水平偏低，不能充分反映空气污染的程度。直观的体现在最后计算得到的总指数AQI明显偏低。如此再按照原有的国家标准进行对比，会发现空气污染都大大减轻了，与实际严重不符。
所以我们采用层次分析加权法得到各个污染物AQI所占的权重。并计算得到了优于平均值的加权指数来反映污染的总体水平。从而也得到了“优化后的半集均方差”（将公式中的平均数用加权指数替代）公式。 （1）构造判断矩阵
通过对指标之间两两重要程度进行比较和分析判断，构造判断矩阵。层次分析法在对指标的相对重要程度进行测量时，我们引入了九分位的相对重要的比例标度。令A为判断矩阵，用以表示同一层次各个指标的相对重要性的判断值。根据心理学家提出的“人区分信息等级能力为7±2”的研究结论，有如下评分规
则：
表2 权重的评分规则
注：取8,6,4,2,1/2,1/4,1/6,1/8为上述评价值的中间值。
由此通过两两比较便可以构造出判断矩阵A，如下表所示：
表3 打分矩阵
（2）对判断矩阵进行一致性检验 ①计算判断矩阵的最大特征根：
1m(Aω)i
λ=∑ max
mi=1ωi
⎡a11⎢a⎢21
Aω=⎢.
⎢⎢.⎢⎣am1
a12a22
..am2
..........
式中Aω为判断矩阵A与特征向量ω的乘积，即为：
a1m⎤⎡ω1⎤
⎢ω⎥a2m⎥⎥⎢2⎥.⎥⎢.⎥⎥⎢⎥ .⎥⎢.⎥amm⎥⎦⎢⎣ωn⎥⎦
②计算判断矩阵的一致性指标：
λmax-m
CI=
m-1
③计算判断矩阵随机一致性比率：
CICR=≤0.10
RI
最后得到结果如下：
表4 各指标权重和一致性指标值
然后把各污染物的AQI与其相应的权重相乘得到综合指数Q。
由此对应的修正后的半集均方差公式：
∑m
(Pi-Q)2
S,
i=1h=
m
修正后的空气质量指数（AQI）的模型为：
AQI=Q+S,h
代入附表二中的数据，得到：
表5空气质量指数 AQI
SO2
CO
NO2
O3
O3
PM10
臭氧
臭氧细颗粒二氧一氧二氧1小8小物（1小化硫
化碳
化氮
时平时平时平均均
均
值）
观测点1 0.13 1.53 0.03 0.01 0.01 0.31 观测点2 0.07 1.38 0.08 0.07 0.04 0.19 观测点3 0.01 3.97 0.02 0.01 0.01 0.27 观测点4 0.03 0.76 0.07 0.03 0.04 0.15 观测点5 0.09 0.81 0.06 0.02 0.02 0.17 观测点6 0.01 0.01 0.05 0.01 0.01 0.16 观测点7 0.04 0.69 0.04 0.04 0.05 0.15 观测点8 0.03 1.05 0.07 0.01 0.04 0.15 观测点9 0.04 1.52 0.02 0.02 0.03 0.11 观测点10 0.02 1.18 0.03 0.03 0.03 0.2 观测点11 0.03
0.5
0.02
0.06
0.06
0.14
观测点1 0.14 1.97 0.08 0 0.01 0.36 观测点2 0.13 2.36 0.08 0.03 0.02 0.18 观测点3 0.1 4.24 0.08 0.02 0.01 0.31 观测点4
0.09 0.57 0.02 0.03 0.05 0.15 PM10
细颗粒物（24小时平均值）
0.31 0.18 0.21 0.49 0.18 0.15 0.14 0.13 0.11 0.17 0.16
0.28 0.15 0.22 0.19 PM2.5 PM2.5 可吸
可吸入颗入颗粒物粒物(1小（24时平小时均值) 平均值）
0.19 0.16 0.14 0.11 0.18 0.1 0.12 0.11 0.13 0.13 0.14 0.09 0.08 0.07 0.07 0.07
0.12 0.11 0.15 0.12 0.12
0.11 0.23 0.16 0.15 0.11 0.21 0.11 0.12
0.11 AQI 空气质量指数 215 159 203 137 152 153 107 95
136 171 135
263 166 234 136
观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 0.02 0.04 0.05 0.05 0.04 0.03
2.9 0.75 1.12 1.58 1.19 0.5
0.04 0.04 0.08 0.02 0.02 0.02
0.04 0.05 0.01 0.02 0.03 0.06
0.04 0.05 0.04 0.03 0.03 0.06
0.2 0.15 0.16 0.14 0.21 0.14
0.19 0.15 0.13 0.11 0.17 0.17
0.12 0.07 0.1 0.14 0.15 0.1
0.1 149 0.07
99
0.07 119 0.11 154 0.12 181 0.11 120
观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2
0.17 2.35 0.13 1.58 0 4.89 0.1 0.57 0.11 2.24 0.02 7.5 0.05 0.8 0.05 1.06 0.07 1.66 0.06 1.11 0.04
0.5
0.2 1.78 0.09 1.27 0.01 4.29 0.1 0.55 0.17 3.19 0.02 7.61 0.09 2.8 0.05 0.86 0.06 1.62 1.62 0.07 0.02
0.07
0.01 0.73 0.05
0.69
0.08 0 0 0.08 0.05 0.04 0.02 0.01 0.01 0.02 0.05 0.05 0.06 0.02 0.01 0.04 0.03 0.03 0.04 0.06 0.05 0.07 0.04 0.04 0.03 0.02 0.02 0.03 0.03 0.03 0.02
0.07
0.06
0.04 0 0 0.06 0.07 0.04 0.02 0.01 0.01 0.02 0.08 0.05 0.07 0.02 0.01 0.04 0.03 0.03 0.1 0.01 0.04 0.02 0.1 0.04 0.02 0.06 0.03 1.21 0.03 0.03 0.5
0.02
0.09
0.06 0.03 0.01 0.04
0.06
0.08
0.54 0.29 0.16 0.13 0.3 0.21 0.16 0.17 0.18 0.19 0.18 0.17 0.22 0.16 0.19 0.13 0.13 0.11 0.2 0.18 0.15
0.16
0.53 0.22 0.2 0.18 0.34 0.21 0.2 0.19 0.18 0.19 0.18 0.17 0.26 0.2 0.22 0.14 0.17
0.12
0.03 .206.177 0.06
0.19
0.05 0.09 0.13
0.1
0.25 0.16 0.15 0.11 0.21 0.11 0.12 0.11 0.13 0.13 0.13 0.12 0.09 0.07 0.09 0.07 0.15 0.11 0.15 0.12 0.11
0.11
0.25 0.16 0.11 0.1 0.2 0.12 0.14 0.11 0.17 0.13 0.11 0.12 0.06 0.05 0.1 0.07 0.14 0.12 0.14 0.12 0.16 0.1
0.04 0.04 0.06
0.05
396 169 225 141 153 137 132 116 164 171 128
371 142 222 165 185 126 127 135 166 196 171
55 86
观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11
0.03 0.89 0.01 0.51 0.01 2.15 0.06 1.12 0.03 0.8 0.01 1.94 0.01 0.75 0.01
1.24
0.02 0.89 0.06 0.86 0.02 0.66 0.03 0.68 0.01 0.5 0.02 3.13 0.06 1.27 0.04 0.79 0.02 1.57 0.02 0.72 0.03
1.2
0.02 0.87 0.07 0.8 0.03 0.84 0.03 0.75 0.01 0.5 0.02 2.81 0.06 1.18 0.03 0.73 0.02 1.64 0.02 0.72 0.02
1.17
0.06 0.02 0.02 0.03 0.03 0.03 0.03 0.05 0.06 0.04 0.05 0.06 0.03 0.07 0.08 0.02 0.02 0.04 0.02 0.03 0.03 0.04
0.02
0.02
0.07 0.01 0.04 0.03 0.09 0.06 0.03 0.06 0.06 0.04 0.03 0.03 0.04 0.03 0.04 0.03 0.06 0.08 0.04 0.07 0.05 0.02 0.08 0.06 0.01 0.06 0.05 0.02 0.04 0.04 0.04
0.03
0.06
0.07 0.01 0 0.03 0.1 0.08 0.03 0.06 0.06 0.05 0.02 0.03 0.03 0.03 0.04 0.03 0.06 0.06 0.05 0.06 0.06 0.02 0.08 0.06 0.01 0.06 0.06 0.02 0.03 0.04 0.04
0.03
0.05
0.04 0.09 0.04 0.07 0.05 0.09 0.1 0.09 0.08 0.06 0.02 0.04 0.06 0.11 0.07
0.07
0.05 0.1 0.06 0.1 0.05 0.05 0.09 0.1 0.05 0.09 0.06 0.09 0.1 0.08 0.1 0.04 0.01 0.05 0.06 0.1 0.05
0.08
0.05 0.1 0.05 0.1 0.07 0.05 0.06 0.1 0.04 0.09 0.06 0.09 0.09 0.08 0.09 0.05 0.03 0.05 0.07 0.1 0.06
0.08
0.03 0.04 0.04 0.04 0.04 0.04 0.03 0.02 0.04 0.02 0.05 0.05 0.04 0.05 0.02
0.02
0.05 0.03 0.02 0.04 0.02 0.02 0.06 0.04 0.04 0.04 0.03 0.04 0.03 0.02 0.03 0.02 0.03 0.04 0.04 0.04 0.02
0.02
0.08 0.04 0 0.04 0.02 0.02 0.05 0.04 0.04 0.04 0.04 0.04 0.02 0.02 0.03 0.02 0.03 0.04 0.05 0.04 0.03
0.02
43 47 50 65 58 50 59 49
57 46 40 75 53 50 62 61 36 55 40
82 39 47 61 51 50 60 58 39 63 46
观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11
0.08 0.02 0.02 0.01 0.01 0.04 0.03 0.01 0.01 0.01
1.35 1.05 1.02 0.46 2.3 1.08 0.69 1.78 0.73 1.3
0.05 0.04 0.07 0.03 0.03 0.05 0.02 0.02 0.03 0.05
0.05 0.06 0.01 0.04 0.06 0.05 0.08 0.04 0.03 0.01
0.08 0.06 0.03 0.04 0.06 0.06 0.08 0.06 0.04 0.04
0.07 0.09 0.04 0.04 0.05 0.09 0.08 0.03 0.07 0.05
0.1 0.05 0.09 0.08 0.09 0.09 0.05 0.05 0.11 0.08
0.03 0.03 0.04 0.04 0.04 0.02 0.04 0.04 0.05 0.03
0.04 0.02 0.04 0.04 0.04 0.02 0.02 0.04 0.04 0.02
56 60 51 49 50 57 58 45 65 42
6 模型的优缺点及推广
在涉及到城市空气质量状况评价时，用层次分析法能综合考虑各种不同因素，通过计算出各个因素的权重值来决定其在最后排序中所占的重要性。运用matlab软件计算出特征根与特征向量，在计算出权重值的大小。该方法能有效避免排序时顾此失彼所产生的差错，能准确反映出各个城市的空气质量状况排名。
从数理统计分析、实际环境影响以及计算结果可知，修正后的半集均方差模型基本上满足了实际要求，与所讨论的原有模型比较更为合理.但由于我国地域辽阔，空气状况差异甚大，该模型还需通过更广泛的应用加以多方面的验证,所以还有很大的优化空间。
7 参考文献
[1] 许鸿楷. 半集均方差水质评价模式 [J] .环境科学，1985（4）.
[2] 王斌会. 多元统计分析及R语言建模 [M] .广州：暨南大学出版社，2010.1.240-243.
[3] 沈恒范，概率论与数理统计教程（第五版) [M] .北京：高等教育出版社，2011.6.213-216. [4] MATLAB中文论坛，MATLAB神经网络30个案例分析 [M] .北京：北京航空航天大学出版社，2010.4.133-135.
[5] 环境空气质量指数(AQI)技术规定（试行）HJ633-2012 [Z].北京：中国环境科学出版社
[6] 环境空气质量标准GB3095-2012 [Z].北京：中国环境科学出版社
附录一
附表一 某地区11个城市的7个时刻的空气质量指标和空气质量指数
SO2
CO
NO2
O3
O3
PM10
PM10
PM2.5 可吸入颗粒物(1小时平均PM2.5 可吸入颗粒物（24小时平均空气质量指数 AQI
二氧化硫
一氧化碳
二氧化氮
臭氧1小时平均
臭氧8小时平均
细颗粒物（1小时平
细颗粒物（24小时平
均值） 均值）
城市1 0.027 城市2 0.221 城市3 0.025 城市4 0.09 城市5 0.089 城市6 0.031 城市7 0.087 城市8 0.009 城市9 0.045 城市10 0.035 城市11 0.048
城市1 0.017 城市2 0.146 城市3 0.021 城市4 0.096 城市5 0.08 城市6 0.016 城市7 0.077 城市8 0.006 城市9 0.047 城市10 0.019 城市11 0.031
城市1 0.02 城市2 0.06 城市3 0.018 城市4 0.031 城市5 0.015 城市6 0.014 城市7 0.048 城市8 0.04 城市9
0.017 0.516 0.025 3.861 0.101 0.639 0.036 0.637 0.027 0.372 0.036 2.817 0.026 0.826 0.04 0.305 0.034 1.486 0.026 0.811 0.031 0.897
0.028
0.608 0.035 2.255 0.082 596 0.029 0.687 0.036 0.411 0.038 2.713 0.023 1.084 0.046 0.255 0.025 1.508 0.023 0.599 0.021 0.885
0.029
1.138 0.038 0.878 0.031 0.555 0.028 0.641 0.033 0.474 0.033 2.794 0.025 1.34 0.034 0.844 0.021 1.569
0.012
0.065 0.055 0.012 0.03 0.043 0.035 0.045 0.071 0.039 0.052 0.082 0.067 0.043 0.063 0.035 0.035 0.036 0.047 0.028 0.019 0.066
0.06
0.048 0.054 0.029 0.026 0.053 0.039 0.033 0.066 0.04 0.05 0.084 0.069 0.037 0.059 0.049 0.036 0.039 0.048 0.033 0.021 0.066
0.059
0.039 0.061 0.091 0.053 0.061 0.053 0.038 0.132 0.034 0.039 0.063 0.075 0.072 0.051 0.077 0.54 0.06
0.051
0.039 0.031 0.22 0.088 0.045 0.051 0.066 0.147 0.174 0.129 0.042 0.031 0.085 0.046 0.044 0.034 0.045 0.057 0.126 0.112 0.069
0.12
0.079 0.118 0.219 0.094 0.049 0.051 0.104 0.136 0.111 0.128 0.075 0.104 0.134 0.049 0.063 0.035 0.056 0.056 0.103 0.111 0.068 0.115
0.034 0.085 0.057 0.101 0.048 0.049 0.036 0.104 0.062 0.094 0.059 0.095 0.123 0.078 0.089 0.042 0.025
0.051 值) 值）
0.108 0.122 0.09 0.034 0.019 0.016 0.031 0.039 0.045 0.034 0.092 0.106 0.018 0.016 0.027 0.013 0.051 0.034 0.054 0.034 0.02
0.021
0.039 0.03 0.083 0.037 0.018 0.016 0.04 0.038 0.04 0.033 0.036 0.031 0.024 0.016 0.024 0.014 0.055 0.035 0.037 0.034 0.015
0.019
0.019 0.018 0.022 0.041 0.013 0.016 0.054 0.037 0.037 0.035 0.035 0.036 0.029 0.018 0.021 0.016 0.033 0.037
79 135 46 59 113 53 68 44 71 88 60
65 135 50 78 81 63 92 57 75 77 70
58 54 48 91 56 55 87 70 48
城市11
0.02
1.138
0.038
0.039
0.061
0.034
0.085
0.019
0.018
34 城市1 0.07 城市2 0.123 城市3 0.03 城市4 0.043 城市5 0.067 城市6 0.025 城市7 0.063 城市8 0.026 城市9 0.057 城市10 0.031 城市11 0.039
城市1 0.138 城市2 0.092 城市3 0.117 城市4 0.041 城市5 0.077 城市6 0.026 城市7 0.049 城市8 0.019 城市9 0.053 城市10 0.034 城市11 0.034
城市1 0.168 城市2 0.113 城市3 0.105 城市4 0.041 城市5 0.076 城市6 0.024 城市7 0.06 城市8 0.031 城市9 0.073 城市10 0.041 城市11 0.033
城市1 0.135 城市2 0.132 城市3 0.097 城市4 0.092 城市5
0.13
0.874 0.04 1.918 0.083 1.116 0.054 0.967 0.041 0.823 0.058 2.312 0.044 1.208 0.058 0.741 0.042 1.72 0.023 0.806 0.04 0.438
0.034
1.542 0.067 1.375 0.062 4.919 0.062 0.895 0.03 1.218 0.058 2.769 0.042 1.329 0.049 1.061 0.059 1.885 0.018 1.009 0.033 0.5
0.023
1.987 0.068 1.578 0.066 4.409 0.064 0.947 0.031 1.27 0.055 2.904 0.039 1.642 0.065 1.091 0.067 1.956 0.017 1.078 0.033 0.5
0.022
1.97 0.077 2.357 0.075 4.243 0.078 0.567 0.019 1.211
0.046
0.076 0.087 0.041 0.095 0.065 0.117 0.082 0.104 0.049 0.053 0.078 0.08 0.093 0.109 0.099 0.114 0.052 0.068 0.042 0.047 0.088
0.125
0.003 0.01 0.041 0.025 0.01 0.012 0.034 0.046 0.009 0.013 0.034 0.041 0.043 0.08 0.014 0.034 0.019 0.031 0.021 0.022 0.057
0.057
0.004 0.006 0.044 0.03 0.019 0.012 0.036 0.044 0.012 0.012 0.035 0.038 0.028 0.068 0.016 0.032 0.026 0.029 0.022 0.022 0.059
0.057
0.003 0.005 0.032 0.02 0.024 0.008 0.034 0.052 0.008
0.003
0.17 0.272 0.144 0.141 0.259 0.26 0.383 0.263 0.147 0.221 0.117 0.175 0.148 0.093 0.151 0.092 0.096 0.069 0.141 0.153 0.163
0.136
0.292 0.278 0.183 0.155 0.246 0.193 0.138 0.253 0.186 0.205 0.204 0.186 0.203 0.164 0.154 0.132 0.158 0.117 0.19 0.177 0.144
0.164
0.352 0.281 0.179 0.156 0.279 0.195 0.152 0.25 0.201 0.199 0.202 0.188 0.203 0.17 0.157 0.135 0.174 0.12 0.197 0.178 0.143
0.165
0.364 0.278 0.176 0.152 0.307 0.221 0.15 0.19 0.246
0.214
0.092 0.142 0.082 0.092 0.058 0.104 0.131 0.114 0.102 0.136 0.062 0.093 0.033 0.049 0.068 0.034 0.1 0.089 0.08 0.101 0.1
0.093
0.178 0.153 0.121 0.097 0.149 0.093 0.112 0.113 0.113 0.128 0.125 0.098 0.057 0.042 0.096 0.065 0.152 0.108 0.138 0.116 0.118
0.106
0.196 0.155 0.122 0.098 0.17 0.096 0.111 0.112 0.129 0.125 0.124 0.1 0.061 0.044 0.112 0.067 0.146 0.11 0.144 0.118 0.1
0.105
0.23 0.162 0.145 0.113 0.211 0.107 0.119 0.11 0.122
0.111
122 110 155 248 135 84 100 101 132 106 132 229 159 199 147 148 165 127 127 202 183 154 247 161 221 146 171 164 127 147 194 192 132 280 193 262 156 180
城市7 城市8 城市9 城市10 城市11
0.041 0.053 0.049 0.036 0.033
0.748 1.117 1.584 1.193 0.5
0.04 0.082 0.024 0.024 0.022
0.046 0.013 0.017 0.031 0.059
0.051 0.04 0.025 0.027 0.057
0.15 0.158 0.142 0.21 0.143
0.149 0.132 0.107 0.174 0.165
0.068 0.101 0.136 0.152 0.1
0.067 0.069 0.108 0.123 0.105
100 134 180 220 132
附表二某地区11个观测点的8个整点时刻的空气质量指标
SO2
二氧化硫
观测点1 0.129 观测点2 0.067 观测点3 0.005 观测点4 0.034 观测点5 0.09 观测点6 0.005 观测点7 0.039 观测点8 0.031 观测点9 0.043
观测点10 0.024
观测点11 0.034
观测点1 0.135 观测点2 0.132 观测点3 0.097 观测点4 0.092 观测点5 0.13 观测点6 0.024 观测点7 0.041 观测点8 0.053 观测点9 0.049 观测点10 0.036 观测点
0.033
CO
NO2
一氧化二氧化碳
氮
1.533 0.031 1.383 0.078 3.968 0.02 0.756 0.067 0.811 0.059 0.006 0.054 0.692 0.037 1.054 0.065 1.517
0.024
1.184
0.03
0.5
0.023
1.97 0.077 2.357 0.075 4.243 0.078 0.567 0.019 1.211 0.046 2.904 0.039 0.748 0.04 1.117 0.082 1.584 0.024 1.193 0.024 0.5
0.022
O3
O3
臭氧1小臭氧8时平均
小时平均
0.005 0.007 0.069 0.038 0.007 0.007 0.033 0.043 0.017 0.016 0.0111 0.011 0.042 0.054 0.014 0.043 0.015
0.027
0.026
0.026
0.057
0.057 0.003 0.005 0.032 0.02 0.024 0.008 0.034 0.052 0.008 0.003 0.035 0.038 0.046 0.051 0.013 0.04 0.017 0.025
0.031 0.027 0.059
0.057
PM10
PM10
细颗粒物（1小细颗粒物时平均（24小时值）
平均值）
0.305 0.305 0.194 0.183 0.272 0.206 0.152 0.485 0.17 0.181 0.162 0.146 0.149 0.144 0.151 0.131 0.105
0.105
0.195
0.173
0.144
0.164
0.364 0.278 0.176 0.152 0.307 0.221 0.15 0.19 0.246 0.214 0.202 0.188 0.15 0.149 0.158 0.132 0.142 0.107 0.21 0.174 0.143
0.165
PM2.5
PM2.5 可吸入可吸入颗粒物颗粒物(1小时（24小平均值)
时平均值）
0.186 0.158 0.136 0.107 0.184 0.103 0.122 0.114 0.13 0.128 0.137 0.088 0.08 0.066 0.068 0.068 0.123 0.106 0.145 0.121 0.118 0.106
0.23 0.162 0.145 0.113 0.211 0.107 0.119 0.11 0.122 0.111 0.124 0.1 0.068 0.067 0.101 0.069 0.136 0.108 0.152 0.123 0.1
0.105
11
观测点1 0.174 观测点2 0.134 观测点3 0.004 观测点4 0.103 观测点5 0.109 观测点6 0.016 观测点7 0.048 观测点8 0.046 观测点9 0.065
观测点10 0.064
观测点11 0.044
观测点1 0.198 观测点2 0.093 观测点3 0.005 观测点4 0.101 观测点5 0.167 观测点6 0.02 观测点7 0.09 观测点8 0.045 观测点9 0.058 观测点10 1.62
观测点11 0.023
观测点1 0.013 观测点2 0.049 观测点3 0.023 观测点4 0.028 观测点5 0.009 观测点6 0.01 观测点7 0.055 观测点8 0.025 观测点9 0.013 观测点10 0.009
观测点11 0.012
2.354 0.075 1.577 0.082 4.889 0.023 0.574 0.019 2.235 0.056 7.502 0.037 0.798 0.035 1.062 0.074 1.66
0.026
1.109
0.031
0.5
0.021
1.777 0.0387 1.272 0.055 4.291 0.024 0.553 0.02 3.192 0.072 7.614 0.039 2.803 0.095 0.855 0.021 1.62 0.023 0.068
1.212
0.068
0.5
0.729 0.06 0.69 0.037 0.959 0.037 0.893 0.059 0.514 0.034 2.149 0.029 1.115 0.043 0.804 0.032 1.936 0.02 0.75
0.023
1.238
0.036
0.003 0.002 0.052 0.038 0.008 0.007 0.054 0.049 0.019 0.01 0.027 0.03 0.057 0.052 0.036 0.038 0.017
0.023
0.027
0.026
0.071
0.058 0.004 0.002 0.069 0.044 0.008 0.007 0.079 0.051 0.023 0.011 0.027 0.029 0.012 0.041 0.103 0.044 0.058 0.03
0.034
0.029
0.024
0.087 0.028 0.013 0.064 0.075 0.042 0.056 0.018 0.021 0.031 0.034 0.051 0.058 0.054 0.059 0.069 0.076 0.023 0.043
0.03
0.033
0.016
0.023
0.539 0.285 0.162 0.134 0.304 0.21 0.157 0.167 0.181 0.194 0.183 0.172 0.22 0.156 0.188 0.134 0.129
0.107
0.196
0.176
0.153
0.164
0.527 0.2179 0.195 0.182 0.337 0.214 0.199 0.186 0.18 0.19 0.179 0.172 0.258 0.195 0.221 0.137 0.172 0.116 0.026 .206.177
0.062
0.19
0.054 0.086 0.125 0.101 0.053 0.054 0.042 0.09 0.039 0.07 0.049 0.085 0.103 0.094 0.078 0.058 0.022 0.042 0.057
0.106
0.069
0.071
0.254 0.164 0.149 0.113 0.213 0.112 0.124 0.109 0.13 0.128 0.126 0.119 0.088 0.07 0.086 0.07 0.148
0.109
0.145
0.115
0.109
0.105
0.245 0.161 0.114 0.104 0.195 0.115 0.141 0.11 0.17 0.128 0.114 0.119 0.057 0.052 0.101 0.072 0.144 0.124 0.14
0.116
0.164
0.097
0.044 0.042 0.057 0.047 0.027 0.02 0.033 0.04 0.038 0.04 0.038 0.043 0.028 0.022 0.037 0.023 0.046 0.045 0.044
0.048
0.024
0.02
观测点1 0.015 观测点2 0.061 观测点3 0.021 观测点4 0.027 观测点5 0.011 观测点6 0.015 观测点7 0.056 观测点8 0.035 观测点9 0.018 观测点10 0.017 观测点11 0.028
观测点1 0.015 观测点2 0.068 观测点3 0.025 观测点4 0.027 观测点5 0.01 观测点6 0.016 观测点7 0.059 观测点8 0.03 观测点9 0.017 观测点10 0.015 观测点11 0.018
观测点1 0.015 观测点2 0.084 观测点3 0.023 观测点4 0.024 观测点5 0.008 观测点6 0.013 观测点7 0.044 观测点8 0.031 观测点9 0.014 观测点10 0.013 观测点11
0.009
0.891 0.065 0.863 0.032 0.658 0.03 0.68 0.038 0.501 0.036 3.134 0.025 1.27 0.039 0.794 0.021 1.572 0.014 0.722 0.017 1.202
0.042
0.867 0.071 0.799 0.031 0.837 0.031 0.753 0.045 0.497 0.033 2.81 0.03 1.184 0.05 0.729 0.021 1.636 0.013 0.719 0.02 1.17
0.04
1.023 0.073 1.3504 0.054 1.054 0.035 1.017 0.067 0.459 0.029 2.295 0.031 1.076 0.052 0.687 0.023 1.781 0.018 0.731 0.026 1.298
0.049
0.014 0.039 0.093 0.063 0.06 0.056 0.032 0.033 0.034 0.037 0.062 0.075 0.069 0.053 0.08 0.062 0.06 0.052 0.037 0.035 0.034
0.057
0.009 0.0032 0.096 0.075 0.062 0.059 0.021 0.03 0.033 0.036 0.06 0.063 0.059 0.055 0.08 0.064 0.059 0.055 0.034 0.036 0.031 0.052 0.008 0.023 0.054 0.082 0.056 0.059 0.008 0.027 0.038 0.036 0.06 0.061 0.052 0.06 0.08 0.075 0.04 0.056 0.032 0.035 0.013
0.039
0.05 0.097 0.061 0.102 0.049 0.048 0.091 0.1 0.046 0.09 0.063 0.094 0.095 0.081 0.1 0.044 0.014 0.049 0.061 0.102 0.046
0.083
0.053 0.096 0.049 0.101 0.065 0.049 0.061 0.096 0.044 0.089 0.056 0.093 0.093 0.084 0.093 0.047 0.026 0.048 0.067 0.103 0.057 0.081
0.055 0.093 0.069 0.103 0.086 0.052 0.04 0.089 0.035 0.078 0.048 0.091 0.087 0.089 0.083 0.052 0.027 0.046 0.072 0.106 0.048
0.078
0.048 0.034 0.021 0.042 0.015 0.017 0.056 0.038 0.043 0.036 0.034 0.037 0.028 0.019 0.027 0.017 0.032 0.038 0.036 0.04 0.024
0.018
0.077 0.037 0.002 0.042 0.02 0.017 0.048 0.038 0.041 0.037 0.036 0.038 0.024 0.02 0.027 0.018 0.034 0.039 0.045 0.041 0.026
0.019
0.044 0.04 0.034 0.044 0.032 0.019 0.044 0.04 0.042 0.038 0.039 0.04 0.024 0.02 0.035 0.019 0.04 0.043 0.047 0.044 0.027
0.02
附录二
程序一：利用一维插值计算各污染物的AQI
%计算so2的AQI的函数 function so2=f1(a) f=[0,150,500,650,800]; t=[0,50,100,150,200]; so2=interp1(f,t,a);
%计算co的AQI的函数 function co=f2(a)
f=[0,5,10,35,60,90,120,150]; t=[0,50,100,150,200,300,400,500]; co=interp1(f,t,a);
%计算no2的AQI的函数 function no2=f3(a)
f=[0,100,200,700,1200,2340,3090,3840]; t=[0,50,100,150,200,300,400,500]; no2=interp1(f,t,a);
%计算o3的AQI的函数 function o31=f4(a)
f=[0,160,200,300,400,800,1000,1200]; t=[0,50,100,150,200,300,400,500]; o31=interp1(f,t,a);
%计算pm10的AQI的函数 function pm10=f6(a)
f=[0,50,150,250,350,420,500,600]; t=[0,50,100,150,200,300,400,500]; pm10=interp1(f,t,a);
%计算pm2.5的AQI的函数 function pm25=f7(a)
f=[0,35,75,115,150,250,350,500];
t=[0,50,100,150,200,300,400,500]; pm25=interp1(f,t,a);
程序二：绘制某一时刻11个城市的各污染物的AQI折线图
%绘制某一时刻11个城市的各污染物的AQI x=[1,2,3,4,5,6]; y1=sum(1,1:6); y2=sum(2,1:6); y3=sum(3,1:6); y4=sum(4,1:6); y5=sum(5,1:6); y6=sum(6,1:6); y7=sum(7,1:6); y8=sum(8,1:6); y9=sum(9,1:6); y10=sum(10,1:6); y11=sum(11,1:6);
plot(x,y1,x,y2,x,y3,x,y4,x,y5,x,y6,x,y7,x,y8,x,y9,x,y10,x,y11) ylabel('AQI')
xlabel('1,2,3,4,5,6·Ö±ð±íʾso2,co,no2,o3,pm10,pm2.5') title('ijһʱ¿Ì¸÷³ÇÊеÄIAQI')
程序三：半集均方差
%某一时刻某一城市各污染物AQI的中位数，平均数 zw=median(sum,2) pj=mean(sum,2)
%筛选某一时刻某一城市各污染物AQI比中位数大的AQI b=zeros(11,3); k=0; for j=1:11 k=0; for i=1:6 if(sum(j,i)>zw(j,1)) k=k+1;
b(j,k)=sum(j,i);
end
end
if(k
b(j,3)=zw(j,1);
end
end
%半集均方差公式
for c=1:11
e=0;
for d=1:3
e=e+(b(c,d)-pj(c,1))^2;
end
s(c,1)=(e/3)^0.5;
end
程序四：“优化后的半集均方差”
%利用层次分析加权法计算各污染物的权重
A=[1 5 5/3 5 5/7 5/7
1/5 1 1/3 1 1/7 1/7
3/5 3 1 3 3/7 3/7
1/5 1 1/3 1 1/7 1/7
7/5 7 7/3 7 1 1
7/5 7 7/3 7 1 1];
[D,X]=eig(A);
ZD=max(max(X));
[m,m]=size(D);
p=max(X);
for i=1:m
if p(i)==ZD
q=i;
end
end
z=sum(D);
for i=1:m
for j=1:m
D(i,j)=D(i,j)/z(j);
end
end
B=(D(:,q))'
ZD
RI=[0 0 0.58 0.9 1.12 1.24 1.32 1.41 1.45 1.49 1.51]; CI=(ZD-m)/(m-1)
CR=CI/RI(m)
%加权指数计算
Q=zeros(11,1);
for i=1:11
Q(i,1)=0;
for j=1:6
Q(i,1)=Q(i,1)+sum(i,j)*B(1,j);
end
end
%“优化后的半集均方差”公式
for c=1:11
e=0;
for d=1:3
e=e+(b(c,d)-p(c,1))^2;
end
s1(c,1)=(e/3)^0.5;
end
%AQI最终的计算式
AQI=s1+Q
数学建模论文
A题 空气质量评价
摘要
本文主要研究空气质量评价的相关问题，为突出改进之后的模型中的实时特性而对数据做了必要的省略处理，然后在现有的国家最新空气污染物监测标准（HJ633-2012环境空气质量指数(AQI)技术规定）的基础上利用半集均方差原理对现有空气质量计算模型进行改进。在论证修正后模型可行性的基础上再对模型加以优化，最后利用优化后的模型对附表二中的各项监测结果得出其空气质量指数。
针对问题一，由于目标模型十分强调实时性，于是把附表一中臭氧8小时平均值﹑细颗粒物24小时平均值﹑可吸入颗粒物24小时平均值做了必要的省略处理。联系实际分析论证了现有模型的局限性，并在此基础上采用半集均方差原理对现有模型进行改进，结果顺利得到优化后的计算模型。
针对问题二，考虑到优化后的计算模型并没有对不同的污染物的危害做出差异化的评价，而是直接取表中所有污染物的AQI平均值进行分析。所以引入层次分析法根据污染物的危害性对不同的污染物赋予相应的权重，对半集均方差公式进行合理修正，最后得到修正后的空气质量计算模型。再代入附表二中的数据即得到各个观测点的空气质量指数。详细的matlab实现程序见附录二。
【关键词】 一维插值 半集均方差 层次分析加权法 优化后的半集均方差
1 问题重述
空气质量指数（AQI）是定量描述空气质量状况的无量纲指数。其数值越大、级别和类别越高，说明空气污染状况越严重，对人体的健康危害也就越大。
空气质量指数实时报一般是发布每个每一整点时刻的空气质量指数。 实时报的指标包括二氧化硫(SO2)、氧化碳(CO)、二氧化氮 (NO2)、臭氧(O3)1小时平均值、臭氧(O3)8小时平均值、一颗粒物(粒径小于等于10μm)、细颗粒物(粒径小于等于2.5μm)的1小时平均值和24小时平均值共计9个指标。福建1中列出了某地区11个城市过去7个时刻的空质量指标取值和相应的空气质量指数。
(1) 建立一种新的空气质量指数计算模型，并比较与现有计算模型的区别。 (2) 利用新的计算模型计算附件2中各个观测点的空气质量指数。
2 基本假设
（1）附表一和附表二中的数据是利用统一的污染物监测仪器并按照统一的测量方法测量得到的。
（2）附表一中的原有的空气质量指数（AQI）是按照国家最新出台的统一标准（HJ633-2012环境空气质量指数(AQI)技术规定）进行计算的。
（3）由于国家最新出台的标准中并没有PM2.5和PM10一小时平均浓度限值，所以计算时采用PM2.5和PM10二十四小时平均浓度限值近似代替。
（4）观测点的测量仪器所测量的不同种污染物浓度之间相互独立，互不影响。
（5）所测量的各个观测点附近的空气污染程度在测量的时刻较为稳定，不发生剧烈变化。
（6）在研究各种指标集对某物影响的过程中，不仅指标集中的最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。
（7）大气中各种污染物对环境和人类的危害程度是不一样的。
3 符号说明
IAQIp 污染物项目P的空气质量分指数；
P 污染物项目P的质量浓度值；
BPHi 表1中与Cp相近的污染物浓度限值的高值位； BPLo 表1中与Cp相近的污染物浓度限值的低值位； IAQIHi 表1中与BPHi对应的空气质量分指数； IAQILo 表1中与BPLo对应的空气质量分指数；
IAQI 空气质量分数 n 污染物项目 Sh 半集均方差；
Pi 大于中位数半集的分指数；
m 大于中位数半集的分指数个数；
平均数；
Q 对各个Pi加权后得到的综合指数；
Sh 修正后的半集均方差；
，
4 问题分析
4.1 问题一的分析
要建立新的空气质量指数计算模型，就必须先明确现有计算模型的不足之处。我国现有空气质量计算模型主要依据《HJ633-2012环境空气质量指数(AQI)技术规定》里的一维插值分析法。
我们以附表一中第一组数据为例进行分析：
（1）先对照下表中各种污染物浓度限值得到相关数值
表1 空气质量分指数及对应的污染物项目浓度限值
（2）把得到的数据代入以下的一维插值公式，分别进行计算
（3）得到空气中各污染物的AQI，并把11个城市的各种污染物AQI绘制成折线统计图。
如下图所示：
图1 某时刻11个城市各污染物所对应的AQI的折线图
由图可知，一般而言，对于某个特定的城市，它的PM10或PM2.5的AQI值远超过其余的各种污染物的AQI值。
而根据我国的最新标准，观测点取各种污染物所对应AQI的最大值来评价整体的空气质量。公式如下：
这就出现了一个问题：在现阶段的测量中，实际上只有PM10或PM2.5的AQI对最终的空气质量AQI有影响。
而在实际生活中，在研究各种指标集对某物影响的过程中，不仅指标集中的
最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。
现有的评价模型在实际应用中正是没有考虑到污染物次大值等的对空气总体质量的影响，不免有些以偏概全。为了解决这一问题，我们引入了半集均方差原理来构建新的空气质量指数计算模型。
4.2 问题二的分析
我们先查找出各种大气污染物对环境和人类有哪些方面的的危害：
（1）二氧化硫的主要危害：形成工业烟雾，高浓度时使人呼吸困难，是著名的伦敦烟雾事件的元凶；进入大气层后，氧化为硫酸（SO4）在云中形成酸雨，对建筑、森林、湖泊、土壤危害大；形成悬浮颗粒物，又称气溶胶，随着人的呼吸进入肺部，对肺有直接损伤作用。
（2）一氧化碳的主要危害：极易与血液中运载氧的血红蛋白结合， 结合速度比氧气快250倍，因此，在极低浓度时就能使人或动物遭到缺氧性伤害。轻者眩晕，头疼，重者脑细胞受到永久性损伤，甚至窒息死亡；对心脏病、贫血和呼吸道疾病的患者伤害性大；引起胎儿生长受损和智力低下。
（3）臭氧的主要危害：低空臭氧是一种最强的氧化剂， 能够与几乎所有的生物物质产生反应，浓度很低时就能损坏橡胶、油漆、织物等材料；臭氧对植物的影响很大。浓度很低时就能减缓植物生长，高浓度时杀死叶片组织，致使整个叶片枯死，最终引起植物死亡，比如高速公路沿线的树木死亡就被分析与臭氧有关；臭氧对于动物和人类有多种伤害作用，特别是伤害眼睛和呼吸系统，加重哮喘类过敏症。
（4）二氧化氮的主要危害：刺激人的眼，鼻，喉和肺， 增加病毒感染的发病率，例如引起导致支气管炎和肺炎的流行性感冒，诱发肺细胞癌变；形成城市的烟雾，影响可见度；破坏树叶的组织，抑制植物生长；在空中形成硝酸小滴，产生酸雨。
（5）PM10和PM2.5的主要危害：随呼吸进入肺，可沉积于肺，引起呼吸系统的疾病。颗粒物上容易附着多种有害物质，有些有致癌性，有些会诱发花粉过敏症；沉积在绿色植物叶面，干扰植物吸收阳光和二氧化碳和放出氧气和水分的过程， 从而影响植物的健康和生长；厚重的颗粒物浓度会影响动物的呼吸系统；杀伤微生物，引起食物链改变，进而影响整个生态系统；遮挡阳光而可能改变气候，这也会影响生态系统。
根据以上资料，可知各种污染物对生态系统危害的原理和方式是不尽相同的，那么，各种大气污染物对环境和人类的危害程度也是不一样的。因此各种污染物AQI的值在整个空气质量评价体系中所占的权重也是不一样的，并不能简单的利用平均值来反映污染的总体水平。
于是我们对问题一中得到的新的空气质量指数计算模型进行必要的加权修正，在比较了专家咨询加权法、相关系数加权法、主成分分析加权法和层次分析加权法的优缺点之后，我们选择了层次分析加权法得到各个污染物AQI所占的权重。并计算得到了优于平均值的加权指数来反映污染的总体水平。从而也得到了“优化后的半集均方差”（将公式中的平均数用加权指数替代）公式。
5 模型的建立与求解
5.1 问题一的模型建立与求解
首先进行数据处理，由于要建立的新空气质量指数模型十分强调实时性，于是把附表一和附表二中臭氧8小时平均值﹑细颗粒物24小时平均值﹑可吸入颗粒物24小时平均值做了必要的省略处理。且国家最新出台的标准中并没有PM2.5和PM10一小时平均浓度限值，所以计算时采用PM2.5和PM10二十四小时平均浓度限值近似代替。
按照数理统计的基本概念，用来评价空气质量的分指数集合，构成了容量为n的样本。而空气质量指数则可认为是反映空气质量特性的样本数字特征，它主要应包括样本观测值的集中趋势及其离散程度两个方面。
（1）算术均值（亦即一阶原点矩）是典型的也是基本的集中趋势测度，其定义式为：
=
∑X
i
i
n
;i=1,2,...n
式中， Xi为样本观测值；n为观测值个数。
算术均值不仅是度量数据离散程度的有效基准,还可以反映样本的总体水平。因此, 算术均值是空气质量评价中的一个重要指标。
（2）样本观测值的离散程度, 最基本的度量是用标准差（或称均方差）,它是数据能变性测度，这里取定义式为：
σ=
∑X
i
i
-2
n
；
i=1,2,...n
由标准差的定义式可以推论：Xi中最大值的增大或最小值的减小,均可导致标准差值的增大。联系到空气质量评价, 也就是把远大于算术均值的分指数和远小于算术均值的分指数置于等同地位, 这显然与它们的环境效应相矛盾。因为分指数中最大值对环境质量的影响是明显的, 而最小值的影响较小甚至可以忽略。因此, 在空气质量评价中若简单地按定义式引入标准差, 尽管可以说明分指数的离散程度, 却未能真实地反映出各自对水质的影响，因而不能满足空气质量评价的需要。因此,有必要从能够充分反映的分指数对空气质量影响的特点出发, 寻求满足于空气质量评价要求的数字特征。
（3）在研究各种指标集对某物影响的过程中，不仅指标集中的最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。因此，为了在评价体系中引入次大值等的影响。提出了用中位数来划分整个评价指令集，对大于中位数的分指数半集给予较大的权重。
联系到在不同的分布形态该半集数列的疏密程度不同,可通过计算它们对算术均值的均方差,以反映大于中位数半集的离散程度。据此引出一个新的数字特征,可将它定义为半集均方差, 其定义式为：
Sh=
∑(P-)
i
i=1
m
2
m
（4）根据上述分析, 空气质量指数AQI作为反映大气环境质量的数字特征,应在反映分指数集合集中趋势的基础上,用大于中位数半集的离散程度加以调整。据此可建立新的空气质量指数（AQI）的数学模型为：
AQI=P+Sh
该模型可简称为半集均方差模型，该模型既通过算术均值考虑了每项分指数对空气质量的影响, 也通过半集均方差对分指数中的大值给予较大的权重，据此, 该模型是按各参数对大气环境影响程度的不同,综合出大气环境质量状况的数量描述。
5.2 问题二的模型建立与求解
由问题一我们得到新的空气质量指数计算模型，但是在实际应用的过程中又发现了半集均方差模型的缺陷。
由图1中的折线图可知：中国城市里PM2.5或PM10的AQI值远超于其他污染物的AQI值。所以导致了水平偏低，不能充分反映空气污染的程度。直观的体现在最后计算得到的总指数AQI明显偏低。如此再按照原有的国家标准进行对比，会发现空气污染都大大减轻了，与实际严重不符。
所以我们采用层次分析加权法得到各个污染物AQI所占的权重。并计算得到了优于平均值的加权指数来反映污染的总体水平。从而也得到了“优化后的半集均方差”（将公式中的平均数用加权指数替代）公式。 （1）构造判断矩阵
通过对指标之间两两重要程度进行比较和分析判断，构造判断矩阵。层次分析法在对指标的相对重要程度进行测量时，我们引入了九分位的相对重要的比例标度。令A为判断矩阵，用以表示同一层次各个指标的相对重要性的判断值。根据心理学家提出的“人区分信息等级能力为7±2”的研究结论，有如下评分规
则：
表2 权重的评分规则
注：取8,6,4,2,1/2,1/4,1/6,1/8为上述评价值的中间值。
由此通过两两比较便可以构造出判断矩阵A，如下表所示：
表3 打分矩阵
（2）对判断矩阵进行一致性检验 ①计算判断矩阵的最大特征根：
1m(Aω)i
λ=∑ max
mi=1ωi
⎡a11⎢a⎢21
Aω=⎢.
⎢⎢.⎢⎣am1
a12a22
..am2
..........
式中Aω为判断矩阵A与特征向量ω的乘积，即为：
a1m⎤⎡ω1⎤
⎢ω⎥a2m⎥⎥⎢2⎥.⎥⎢.⎥⎥⎢⎥ .⎥⎢.⎥amm⎥⎦⎢⎣ωn⎥⎦
②计算判断矩阵的一致性指标：
λmax-m
CI=
m-1
③计算判断矩阵随机一致性比率：
CICR=≤0.10
RI
最后得到结果如下：
表4 各指标权重和一致性指标值
然后把各污染物的AQI与其相应的权重相乘得到综合指数Q。
由此对应的修正后的半集均方差公式：
∑m
(Pi-Q)2
S,
i=1h=
m
修正后的空气质量指数（AQI）的模型为：
AQI=Q+S,h
代入附表二中的数据，得到：
表5空气质量指数 AQI
SO2
CO
NO2
O3
O3
PM10
臭氧
臭氧细颗粒二氧一氧二氧1小8小物（1小化硫
化碳
化氮
时平时平时平均均
均
值）
观测点1 0.13 1.53 0.03 0.01 0.01 0.31 观测点2 0.07 1.38 0.08 0.07 0.04 0.19 观测点3 0.01 3.97 0.02 0.01 0.01 0.27 观测点4 0.03 0.76 0.07 0.03 0.04 0.15 观测点5 0.09 0.81 0.06 0.02 0.02 0.17 观测点6 0.01 0.01 0.05 0.01 0.01 0.16 观测点7 0.04 0.69 0.04 0.04 0.05 0.15 观测点8 0.03 1.05 0.07 0.01 0.04 0.15 观测点9 0.04 1.52 0.02 0.02 0.03 0.11 观测点10 0.02 1.18 0.03 0.03 0.03 0.2 观测点11 0.03
0.5
0.02
0.06
0.06
0.14
观测点1 0.14 1.97 0.08 0 0.01 0.36 观测点2 0.13 2.36 0.08 0.03 0.02 0.18 观测点3 0.1 4.24 0.08 0.02 0.01 0.31 观测点4
0.09 0.57 0.02 0.03 0.05 0.15 PM10
细颗粒物（24小时平均值）
0.31 0.18 0.21 0.49 0.18 0.15 0.14 0.13 0.11 0.17 0.16
0.28 0.15 0.22 0.19 PM2.5 PM2.5 可吸
可吸入颗入颗粒物粒物(1小（24时平小时均值) 平均值）
0.19 0.16 0.14 0.11 0.18 0.1 0.12 0.11 0.13 0.13 0.14 0.09 0.08 0.07 0.07 0.07
0.12 0.11 0.15 0.12 0.12
0.11 0.23 0.16 0.15 0.11 0.21 0.11 0.12
0.11 AQI 空气质量指数 215 159 203 137 152 153 107 95
136 171 135
263 166 234 136
观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 0.02 0.04 0.05 0.05 0.04 0.03
2.9 0.75 1.12 1.58 1.19 0.5
0.04 0.04 0.08 0.02 0.02 0.02
0.04 0.05 0.01 0.02 0.03 0.06
0.04 0.05 0.04 0.03 0.03 0.06
0.2 0.15 0.16 0.14 0.21 0.14
0.19 0.15 0.13 0.11 0.17 0.17
0.12 0.07 0.1 0.14 0.15 0.1
0.1 149 0.07
99
0.07 119 0.11 154 0.12 181 0.11 120
观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2
0.17 2.35 0.13 1.58 0 4.89 0.1 0.57 0.11 2.24 0.02 7.5 0.05 0.8 0.05 1.06 0.07 1.66 0.06 1.11 0.04
0.5
0.2 1.78 0.09 1.27 0.01 4.29 0.1 0.55 0.17 3.19 0.02 7.61 0.09 2.8 0.05 0.86 0.06 1.62 1.62 0.07 0.02
0.07
0.01 0.73 0.05
0.69
0.08 0 0 0.08 0.05 0.04 0.02 0.01 0.01 0.02 0.05 0.05 0.06 0.02 0.01 0.04 0.03 0.03 0.04 0.06 0.05 0.07 0.04 0.04 0.03 0.02 0.02 0.03 0.03 0.03 0.02
0.07
0.06
0.04 0 0 0.06 0.07 0.04 0.02 0.01 0.01 0.02 0.08 0.05 0.07 0.02 0.01 0.04 0.03 0.03 0.1 0.01 0.04 0.02 0.1 0.04 0.02 0.06 0.03 1.21 0.03 0.03 0.5
0.02
0.09
0.06 0.03 0.01 0.04
0.06
0.08
0.54 0.29 0.16 0.13 0.3 0.21 0.16 0.17 0.18 0.19 0.18 0.17 0.22 0.16 0.19 0.13 0.13 0.11 0.2 0.18 0.15
0.16
0.53 0.22 0.2 0.18 0.34 0.21 0.2 0.19 0.18 0.19 0.18 0.17 0.26 0.2 0.22 0.14 0.17
0.12
0.03 .206.177 0.06
0.19
0.05 0.09 0.13
0.1
0.25 0.16 0.15 0.11 0.21 0.11 0.12 0.11 0.13 0.13 0.13 0.12 0.09 0.07 0.09 0.07 0.15 0.11 0.15 0.12 0.11
0.11
0.25 0.16 0.11 0.1 0.2 0.12 0.14 0.11 0.17 0.13 0.11 0.12 0.06 0.05 0.1 0.07 0.14 0.12 0.14 0.12 0.16 0.1
0.04 0.04 0.06
0.05
396 169 225 141 153 137 132 116 164 171 128
371 142 222 165 185 126 127 135 166 196 171
55 86
观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11 观测点1 观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11
0.03 0.89 0.01 0.51 0.01 2.15 0.06 1.12 0.03 0.8 0.01 1.94 0.01 0.75 0.01
1.24
0.02 0.89 0.06 0.86 0.02 0.66 0.03 0.68 0.01 0.5 0.02 3.13 0.06 1.27 0.04 0.79 0.02 1.57 0.02 0.72 0.03
1.2
0.02 0.87 0.07 0.8 0.03 0.84 0.03 0.75 0.01 0.5 0.02 2.81 0.06 1.18 0.03 0.73 0.02 1.64 0.02 0.72 0.02
1.17
0.06 0.02 0.02 0.03 0.03 0.03 0.03 0.05 0.06 0.04 0.05 0.06 0.03 0.07 0.08 0.02 0.02 0.04 0.02 0.03 0.03 0.04
0.02
0.02
0.07 0.01 0.04 0.03 0.09 0.06 0.03 0.06 0.06 0.04 0.03 0.03 0.04 0.03 0.04 0.03 0.06 0.08 0.04 0.07 0.05 0.02 0.08 0.06 0.01 0.06 0.05 0.02 0.04 0.04 0.04
0.03
0.06
0.07 0.01 0 0.03 0.1 0.08 0.03 0.06 0.06 0.05 0.02 0.03 0.03 0.03 0.04 0.03 0.06 0.06 0.05 0.06 0.06 0.02 0.08 0.06 0.01 0.06 0.06 0.02 0.03 0.04 0.04
0.03
0.05
0.04 0.09 0.04 0.07 0.05 0.09 0.1 0.09 0.08 0.06 0.02 0.04 0.06 0.11 0.07
0.07
0.05 0.1 0.06 0.1 0.05 0.05 0.09 0.1 0.05 0.09 0.06 0.09 0.1 0.08 0.1 0.04 0.01 0.05 0.06 0.1 0.05
0.08
0.05 0.1 0.05 0.1 0.07 0.05 0.06 0.1 0.04 0.09 0.06 0.09 0.09 0.08 0.09 0.05 0.03 0.05 0.07 0.1 0.06
0.08
0.03 0.04 0.04 0.04 0.04 0.04 0.03 0.02 0.04 0.02 0.05 0.05 0.04 0.05 0.02
0.02
0.05 0.03 0.02 0.04 0.02 0.02 0.06 0.04 0.04 0.04 0.03 0.04 0.03 0.02 0.03 0.02 0.03 0.04 0.04 0.04 0.02
0.02
0.08 0.04 0 0.04 0.02 0.02 0.05 0.04 0.04 0.04 0.04 0.04 0.02 0.02 0.03 0.02 0.03 0.04 0.05 0.04 0.03
0.02
43 47 50 65 58 50 59 49
57 46 40 75 53 50 62 61 36 55 40
82 39 47 61 51 50 60 58 39 63 46
观测点2 观测点3 观测点4 观测点5 观测点6 观测点7 观测点8 观测点9 观测点10 观测点11
0.08 0.02 0.02 0.01 0.01 0.04 0.03 0.01 0.01 0.01
1.35 1.05 1.02 0.46 2.3 1.08 0.69 1.78 0.73 1.3
0.05 0.04 0.07 0.03 0.03 0.05 0.02 0.02 0.03 0.05
0.05 0.06 0.01 0.04 0.06 0.05 0.08 0.04 0.03 0.01
0.08 0.06 0.03 0.04 0.06 0.06 0.08 0.06 0.04 0.04
0.07 0.09 0.04 0.04 0.05 0.09 0.08 0.03 0.07 0.05
0.1 0.05 0.09 0.08 0.09 0.09 0.05 0.05 0.11 0.08
0.03 0.03 0.04 0.04 0.04 0.02 0.04 0.04 0.05 0.03
0.04 0.02 0.04 0.04 0.04 0.02 0.02 0.04 0.04 0.02
56 60 51 49 50 57 58 45 65 42
6 模型的优缺点及推广
在涉及到城市空气质量状况评价时，用层次分析法能综合考虑各种不同因素，通过计算出各个因素的权重值来决定其在最后排序中所占的重要性。运用matlab软件计算出特征根与特征向量，在计算出权重值的大小。该方法能有效避免排序时顾此失彼所产生的差错，能准确反映出各个城市的空气质量状况排名。
从数理统计分析、实际环境影响以及计算结果可知，修正后的半集均方差模型基本上满足了实际要求，与所讨论的原有模型比较更为合理.但由于我国地域辽阔，空气状况差异甚大，该模型还需通过更广泛的应用加以多方面的验证,所以还有很大的优化空间。
7 参考文献
[1] 许鸿楷. 半集均方差水质评价模式 [J] .环境科学，1985（4）.
[2] 王斌会. 多元统计分析及R语言建模 [M] .广州：暨南大学出版社，2010.1.240-243.
[3] 沈恒范，概率论与数理统计教程（第五版) [M] .北京：高等教育出版社，2011.6.213-216. [4] MATLAB中文论坛，MATLAB神经网络30个案例分析 [M] .北京：北京航空航天大学出版社，2010.4.133-135.
[5] 环境空气质量指数(AQI)技术规定（试行）HJ633-2012 [Z].北京：中国环境科学出版社
[6] 环境空气质量标准GB3095-2012 [Z].北京：中国环境科学出版社
附录一
附表一 某地区11个城市的7个时刻的空气质量指标和空气质量指数
SO2
CO
NO2
O3
O3
PM10
PM10
PM2.5 可吸入颗粒物(1小时平均PM2.5 可吸入颗粒物（24小时平均空气质量指数 AQI
二氧化硫
一氧化碳
二氧化氮
臭氧1小时平均
臭氧8小时平均
细颗粒物（1小时平
细颗粒物（24小时平
均值） 均值）
城市1 0.027 城市2 0.221 城市3 0.025 城市4 0.09 城市5 0.089 城市6 0.031 城市7 0.087 城市8 0.009 城市9 0.045 城市10 0.035 城市11 0.048
城市1 0.017 城市2 0.146 城市3 0.021 城市4 0.096 城市5 0.08 城市6 0.016 城市7 0.077 城市8 0.006 城市9 0.047 城市10 0.019 城市11 0.031
城市1 0.02 城市2 0.06 城市3 0.018 城市4 0.031 城市5 0.015 城市6 0.014 城市7 0.048 城市8 0.04 城市9
0.017 0.516 0.025 3.861 0.101 0.639 0.036 0.637 0.027 0.372 0.036 2.817 0.026 0.826 0.04 0.305 0.034 1.486 0.026 0.811 0.031 0.897
0.028
0.608 0.035 2.255 0.082 596 0.029 0.687 0.036 0.411 0.038 2.713 0.023 1.084 0.046 0.255 0.025 1.508 0.023 0.599 0.021 0.885
0.029
1.138 0.038 0.878 0.031 0.555 0.028 0.641 0.033 0.474 0.033 2.794 0.025 1.34 0.034 0.844 0.021 1.569
0.012
0.065 0.055 0.012 0.03 0.043 0.035 0.045 0.071 0.039 0.052 0.082 0.067 0.043 0.063 0.035 0.035 0.036 0.047 0.028 0.019 0.066
0.06
0.048 0.054 0.029 0.026 0.053 0.039 0.033 0.066 0.04 0.05 0.084 0.069 0.037 0.059 0.049 0.036 0.039 0.048 0.033 0.021 0.066
0.059
0.039 0.061 0.091 0.053 0.061 0.053 0.038 0.132 0.034 0.039 0.063 0.075 0.072 0.051 0.077 0.54 0.06
0.051
0.039 0.031 0.22 0.088 0.045 0.051 0.066 0.147 0.174 0.129 0.042 0.031 0.085 0.046 0.044 0.034 0.045 0.057 0.126 0.112 0.069
0.12
0.079 0.118 0.219 0.094 0.049 0.051 0.104 0.136 0.111 0.128 0.075 0.104 0.134 0.049 0.063 0.035 0.056 0.056 0.103 0.111 0.068 0.115
0.034 0.085 0.057 0.101 0.048 0.049 0.036 0.104 0.062 0.094 0.059 0.095 0.123 0.078 0.089 0.042 0.025
0.051 值) 值）
0.108 0.122 0.09 0.034 0.019 0.016 0.031 0.039 0.045 0.034 0.092 0.106 0.018 0.016 0.027 0.013 0.051 0.034 0.054 0.034 0.02
0.021
0.039 0.03 0.083 0.037 0.018 0.016 0.04 0.038 0.04 0.033 0.036 0.031 0.024 0.016 0.024 0.014 0.055 0.035 0.037 0.034 0.015
0.019
0.019 0.018 0.022 0.041 0.013 0.016 0.054 0.037 0.037 0.035 0.035 0.036 0.029 0.018 0.021 0.016 0.033 0.037
79 135 46 59 113 53 68 44 71 88 60
65 135 50 78 81 63 92 57 75 77 70
58 54 48 91 56 55 87 70 48
城市11
0.02
1.138
0.038
0.039
0.061
0.034
0.085
0.019
0.018
34 城市1 0.07 城市2 0.123 城市3 0.03 城市4 0.043 城市5 0.067 城市6 0.025 城市7 0.063 城市8 0.026 城市9 0.057 城市10 0.031 城市11 0.039
城市1 0.138 城市2 0.092 城市3 0.117 城市4 0.041 城市5 0.077 城市6 0.026 城市7 0.049 城市8 0.019 城市9 0.053 城市10 0.034 城市11 0.034
城市1 0.168 城市2 0.113 城市3 0.105 城市4 0.041 城市5 0.076 城市6 0.024 城市7 0.06 城市8 0.031 城市9 0.073 城市10 0.041 城市11 0.033
城市1 0.135 城市2 0.132 城市3 0.097 城市4 0.092 城市5
0.13
0.874 0.04 1.918 0.083 1.116 0.054 0.967 0.041 0.823 0.058 2.312 0.044 1.208 0.058 0.741 0.042 1.72 0.023 0.806 0.04 0.438
0.034
1.542 0.067 1.375 0.062 4.919 0.062 0.895 0.03 1.218 0.058 2.769 0.042 1.329 0.049 1.061 0.059 1.885 0.018 1.009 0.033 0.5
0.023
1.987 0.068 1.578 0.066 4.409 0.064 0.947 0.031 1.27 0.055 2.904 0.039 1.642 0.065 1.091 0.067 1.956 0.017 1.078 0.033 0.5
0.022
1.97 0.077 2.357 0.075 4.243 0.078 0.567 0.019 1.211
0.046
0.076 0.087 0.041 0.095 0.065 0.117 0.082 0.104 0.049 0.053 0.078 0.08 0.093 0.109 0.099 0.114 0.052 0.068 0.042 0.047 0.088
0.125
0.003 0.01 0.041 0.025 0.01 0.012 0.034 0.046 0.009 0.013 0.034 0.041 0.043 0.08 0.014 0.034 0.019 0.031 0.021 0.022 0.057
0.057
0.004 0.006 0.044 0.03 0.019 0.012 0.036 0.044 0.012 0.012 0.035 0.038 0.028 0.068 0.016 0.032 0.026 0.029 0.022 0.022 0.059
0.057
0.003 0.005 0.032 0.02 0.024 0.008 0.034 0.052 0.008
0.003
0.17 0.272 0.144 0.141 0.259 0.26 0.383 0.263 0.147 0.221 0.117 0.175 0.148 0.093 0.151 0.092 0.096 0.069 0.141 0.153 0.163
0.136
0.292 0.278 0.183 0.155 0.246 0.193 0.138 0.253 0.186 0.205 0.204 0.186 0.203 0.164 0.154 0.132 0.158 0.117 0.19 0.177 0.144
0.164
0.352 0.281 0.179 0.156 0.279 0.195 0.152 0.25 0.201 0.199 0.202 0.188 0.203 0.17 0.157 0.135 0.174 0.12 0.197 0.178 0.143
0.165
0.364 0.278 0.176 0.152 0.307 0.221 0.15 0.19 0.246
0.214
0.092 0.142 0.082 0.092 0.058 0.104 0.131 0.114 0.102 0.136 0.062 0.093 0.033 0.049 0.068 0.034 0.1 0.089 0.08 0.101 0.1
0.093
0.178 0.153 0.121 0.097 0.149 0.093 0.112 0.113 0.113 0.128 0.125 0.098 0.057 0.042 0.096 0.065 0.152 0.108 0.138 0.116 0.118
0.106
0.196 0.155 0.122 0.098 0.17 0.096 0.111 0.112 0.129 0.125 0.124 0.1 0.061 0.044 0.112 0.067 0.146 0.11 0.144 0.118 0.1
0.105
0.23 0.162 0.145 0.113 0.211 0.107 0.119 0.11 0.122
0.111
122 110 155 248 135 84 100 101 132 106 132 229 159 199 147 148 165 127 127 202 183 154 247 161 221 146 171 164 127 147 194 192 132 280 193 262 156 180
城市7 城市8 城市9 城市10 城市11
0.041 0.053 0.049 0.036 0.033
0.748 1.117 1.584 1.193 0.5
0.04 0.082 0.024 0.024 0.022
0.046 0.013 0.017 0.031 0.059
0.051 0.04 0.025 0.027 0.057
0.15 0.158 0.142 0.21 0.143
0.149 0.132 0.107 0.174 0.165
0.068 0.101 0.136 0.152 0.1
0.067 0.069 0.108 0.123 0.105
100 134 180 220 132
附表二某地区11个观测点的8个整点时刻的空气质量指标
SO2
二氧化硫
观测点1 0.129 观测点2 0.067 观测点3 0.005 观测点4 0.034 观测点5 0.09 观测点6 0.005 观测点7 0.039 观测点8 0.031 观测点9 0.043
观测点10 0.024
观测点11 0.034
观测点1 0.135 观测点2 0.132 观测点3 0.097 观测点4 0.092 观测点5 0.13 观测点6 0.024 观测点7 0.041 观测点8 0.053 观测点9 0.049 观测点10 0.036 观测点
0.033
CO
NO2
一氧化二氧化碳
氮
1.533 0.031 1.383 0.078 3.968 0.02 0.756 0.067 0.811 0.059 0.006 0.054 0.692 0.037 1.054 0.065 1.517
0.024
1.184
0.03
0.5
0.023
1.97 0.077 2.357 0.075 4.243 0.078 0.567 0.019 1.211 0.046 2.904 0.039 0.748 0.04 1.117 0.082 1.584 0.024 1.193 0.024 0.5
0.022
O3
O3
臭氧1小臭氧8时平均
小时平均
0.005 0.007 0.069 0.038 0.007 0.007 0.033 0.043 0.017 0.016 0.0111 0.011 0.042 0.054 0.014 0.043 0.015
0.027
0.026
0.026
0.057
0.057 0.003 0.005 0.032 0.02 0.024 0.008 0.034 0.052 0.008 0.003 0.035 0.038 0.046 0.051 0.013 0.04 0.017 0.025
0.031 0.027 0.059
0.057
PM10
PM10
细颗粒物（1小细颗粒物时平均（24小时值）
平均值）
0.305 0.305 0.194 0.183 0.272 0.206 0.152 0.485 0.17 0.181 0.162 0.146 0.149 0.144 0.151 0.131 0.105
0.105
0.195
0.173
0.144
0.164
0.364 0.278 0.176 0.152 0.307 0.221 0.15 0.19 0.246 0.214 0.202 0.188 0.15 0.149 0.158 0.132 0.142 0.107 0.21 0.174 0.143
0.165
PM2.5
PM2.5 可吸入可吸入颗粒物颗粒物(1小时（24小平均值)
时平均值）
0.186 0.158 0.136 0.107 0.184 0.103 0.122 0.114 0.13 0.128 0.137 0.088 0.08 0.066 0.068 0.068 0.123 0.106 0.145 0.121 0.118 0.106
0.23 0.162 0.145 0.113 0.211 0.107 0.119 0.11 0.122 0.111 0.124 0.1 0.068 0.067 0.101 0.069 0.136 0.108 0.152 0.123 0.1
0.105
11
观测点1 0.174 观测点2 0.134 观测点3 0.004 观测点4 0.103 观测点5 0.109 观测点6 0.016 观测点7 0.048 观测点8 0.046 观测点9 0.065
观测点10 0.064
观测点11 0.044
观测点1 0.198 观测点2 0.093 观测点3 0.005 观测点4 0.101 观测点5 0.167 观测点6 0.02 观测点7 0.09 观测点8 0.045 观测点9 0.058 观测点10 1.62
观测点11 0.023
观测点1 0.013 观测点2 0.049 观测点3 0.023 观测点4 0.028 观测点5 0.009 观测点6 0.01 观测点7 0.055 观测点8 0.025 观测点9 0.013 观测点10 0.009
观测点11 0.012
2.354 0.075 1.577 0.082 4.889 0.023 0.574 0.019 2.235 0.056 7.502 0.037 0.798 0.035 1.062 0.074 1.66
0.026
1.109
0.031
0.5
0.021
1.777 0.0387 1.272 0.055 4.291 0.024 0.553 0.02 3.192 0.072 7.614 0.039 2.803 0.095 0.855 0.021 1.62 0.023 0.068
1.212
0.068
0.5
0.729 0.06 0.69 0.037 0.959 0.037 0.893 0.059 0.514 0.034 2.149 0.029 1.115 0.043 0.804 0.032 1.936 0.02 0.75
0.023
1.238
0.036
0.003 0.002 0.052 0.038 0.008 0.007 0.054 0.049 0.019 0.01 0.027 0.03 0.057 0.052 0.036 0.038 0.017
0.023
0.027
0.026
0.071
0.058 0.004 0.002 0.069 0.044 0.008 0.007 0.079 0.051 0.023 0.011 0.027 0.029 0.012 0.041 0.103 0.044 0.058 0.03
0.034
0.029
0.024
0.087 0.028 0.013 0.064 0.075 0.042 0.056 0.018 0.021 0.031 0.034 0.051 0.058 0.054 0.059 0.069 0.076 0.023 0.043
0.03
0.033
0.016
0.023
0.539 0.285 0.162 0.134 0.304 0.21 0.157 0.167 0.181 0.194 0.183 0.172 0.22 0.156 0.188 0.134 0.129
0.107
0.196
0.176
0.153
0.164
0.527 0.2179 0.195 0.182 0.337 0.214 0.199 0.186 0.18 0.19 0.179 0.172 0.258 0.195 0.221 0.137 0.172 0.116 0.026 .206.177
0.062
0.19
0.054 0.086 0.125 0.101 0.053 0.054 0.042 0.09 0.039 0.07 0.049 0.085 0.103 0.094 0.078 0.058 0.022 0.042 0.057
0.106
0.069
0.071
0.254 0.164 0.149 0.113 0.213 0.112 0.124 0.109 0.13 0.128 0.126 0.119 0.088 0.07 0.086 0.07 0.148
0.109
0.145
0.115
0.109
0.105
0.245 0.161 0.114 0.104 0.195 0.115 0.141 0.11 0.17 0.128 0.114 0.119 0.057 0.052 0.101 0.072 0.144 0.124 0.14
0.116
0.164
0.097
0.044 0.042 0.057 0.047 0.027 0.02 0.033 0.04 0.038 0.04 0.038 0.043 0.028 0.022 0.037 0.023 0.046 0.045 0.044
0.048
0.024
0.02
观测点1 0.015 观测点2 0.061 观测点3 0.021 观测点4 0.027 观测点5 0.011 观测点6 0.015 观测点7 0.056 观测点8 0.035 观测点9 0.018 观测点10 0.017 观测点11 0.028
观测点1 0.015 观测点2 0.068 观测点3 0.025 观测点4 0.027 观测点5 0.01 观测点6 0.016 观测点7 0.059 观测点8 0.03 观测点9 0.017 观测点10 0.015 观测点11 0.018
观测点1 0.015 观测点2 0.084 观测点3 0.023 观测点4 0.024 观测点5 0.008 观测点6 0.013 观测点7 0.044 观测点8 0.031 观测点9 0.014 观测点10 0.013 观测点11
0.009
0.891 0.065 0.863 0.032 0.658 0.03 0.68 0.038 0.501 0.036 3.134 0.025 1.27 0.039 0.794 0.021 1.572 0.014 0.722 0.017 1.202
0.042
0.867 0.071 0.799 0.031 0.837 0.031 0.753 0.045 0.497 0.033 2.81 0.03 1.184 0.05 0.729 0.021 1.636 0.013 0.719 0.02 1.17
0.04
1.023 0.073 1.3504 0.054 1.054 0.035 1.017 0.067 0.459 0.029 2.295 0.031 1.076 0.052 0.687 0.023 1.781 0.018 0.731 0.026 1.298
0.049
0.014 0.039 0.093 0.063 0.06 0.056 0.032 0.033 0.034 0.037 0.062 0.075 0.069 0.053 0.08 0.062 0.06 0.052 0.037 0.035 0.034
0.057
0.009 0.0032 0.096 0.075 0.062 0.059 0.021 0.03 0.033 0.036 0.06 0.063 0.059 0.055 0.08 0.064 0.059 0.055 0.034 0.036 0.031 0.052 0.008 0.023 0.054 0.082 0.056 0.059 0.008 0.027 0.038 0.036 0.06 0.061 0.052 0.06 0.08 0.075 0.04 0.056 0.032 0.035 0.013
0.039
0.05 0.097 0.061 0.102 0.049 0.048 0.091 0.1 0.046 0.09 0.063 0.094 0.095 0.081 0.1 0.044 0.014 0.049 0.061 0.102 0.046
0.083
0.053 0.096 0.049 0.101 0.065 0.049 0.061 0.096 0.044 0.089 0.056 0.093 0.093 0.084 0.093 0.047 0.026 0.048 0.067 0.103 0.057 0.081
0.055 0.093 0.069 0.103 0.086 0.052 0.04 0.089 0.035 0.078 0.048 0.091 0.087 0.089 0.083 0.052 0.027 0.046 0.072 0.106 0.048
0.078
0.048 0.034 0.021 0.042 0.015 0.017 0.056 0.038 0.043 0.036 0.034 0.037 0.028 0.019 0.027 0.017 0.032 0.038 0.036 0.04 0.024
0.018
0.077 0.037 0.002 0.042 0.02 0.017 0.048 0.038 0.041 0.037 0.036 0.038 0.024 0.02 0.027 0.018 0.034 0.039 0.045 0.041 0.026
0.019
0.044 0.04 0.034 0.044 0.032 0.019 0.044 0.04 0.042 0.038 0.039 0.04 0.024 0.02 0.035 0.019 0.04 0.043 0.047 0.044 0.027
0.02
附录二
程序一：利用一维插值计算各污染物的AQI
%计算so2的AQI的函数 function so2=f1(a) f=[0,150,500,650,800]; t=[0,50,100,150,200]; so2=interp1(f,t,a);
%计算co的AQI的函数 function co=f2(a)
f=[0,5,10,35,60,90,120,150]; t=[0,50,100,150,200,300,400,500]; co=interp1(f,t,a);
%计算no2的AQI的函数 function no2=f3(a)
f=[0,100,200,700,1200,2340,3090,3840]; t=[0,50,100,150,200,300,400,500]; no2=interp1(f,t,a);
%计算o3的AQI的函数 function o31=f4(a)
f=[0,160,200,300,400,800,1000,1200]; t=[0,50,100,150,200,300,400,500]; o31=interp1(f,t,a);
%计算pm10的AQI的函数 function pm10=f6(a)
f=[0,50,150,250,350,420,500,600]; t=[0,50,100,150,200,300,400,500]; pm10=interp1(f,t,a);
%计算pm2.5的AQI的函数 function pm25=f7(a)
f=[0,35,75,115,150,250,350,500];
t=[0,50,100,150,200,300,400,500]; pm25=interp1(f,t,a);
程序二：绘制某一时刻11个城市的各污染物的AQI折线图
%绘制某一时刻11个城市的各污染物的AQI x=[1,2,3,4,5,6]; y1=sum(1,1:6); y2=sum(2,1:6); y3=sum(3,1:6); y4=sum(4,1:6); y5=sum(5,1:6); y6=sum(6,1:6); y7=sum(7,1:6); y8=sum(8,1:6); y9=sum(9,1:6); y10=sum(10,1:6); y11=sum(11,1:6);
plot(x,y1,x,y2,x,y3,x,y4,x,y5,x,y6,x,y7,x,y8,x,y9,x,y10,x,y11) ylabel('AQI')
xlabel('1,2,3,4,5,6·Ö±ð±íʾso2,co,no2,o3,pm10,pm2.5') title('ijһʱ¿Ì¸÷³ÇÊеÄIAQI')
程序三：半集均方差
%某一时刻某一城市各污染物AQI的中位数，平均数 zw=median(sum,2) pj=mean(sum,2)
%筛选某一时刻某一城市各污染物AQI比中位数大的AQI b=zeros(11,3); k=0; for j=1:11 k=0; for i=1:6 if(sum(j,i)>zw(j,1)) k=k+1;
b(j,k)=sum(j,i);
end
end
if(k
b(j,3)=zw(j,1);
end
end
%半集均方差公式
for c=1:11
e=0;
for d=1:3
e=e+(b(c,d)-pj(c,1))^2;
end
s(c,1)=(e/3)^0.5;
end
程序四：“优化后的半集均方差”
%利用层次分析加权法计算各污染物的权重
A=[1 5 5/3 5 5/7 5/7
1/5 1 1/3 1 1/7 1/7
3/5 3 1 3 3/7 3/7
1/5 1 1/3 1 1/7 1/7
7/5 7 7/3 7 1 1
7/5 7 7/3 7 1 1];
[D,X]=eig(A);
ZD=max(max(X));
[m,m]=size(D);
p=max(X);
for i=1:m
if p(i)==ZD
q=i;
end
end
z=sum(D);
for i=1:m
for j=1:m
D(i,j)=D(i,j)/z(j);
end
end
B=(D(:,q))'
ZD
RI=[0 0 0.58 0.9 1.12 1.24 1.32 1.41 1.45 1.49 1.51]; CI=(ZD-m)/(m-1)
CR=CI/RI(m)
%加权指数计算
Q=zeros(11,1);
for i=1:11
Q(i,1)=0;
for j=1:6
Q(i,1)=Q(i,1)+sum(i,j)*B(1,j);
end
end
%“优化后的半集均方差”公式
for c=1:11
e=0;
for d=1:3
e=e+(b(c,d)-p(c,1))^2;
end
s1(c,1)=(e/3)^0.5;
end
%AQI最终的计算式
AQI=s1+Q