数学建模论文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 符号说明p IAQI 污染物项目P 的空气质量分指数；P 污染物项目P 的质量浓度值；Hi BP 表1中与p C 相近的污染物浓度限值的高值位；Lo BP 表1中与p C 相近的污染物浓度限值的低值位；Hi IAQI 表1中与Hi BP 对应的空气质量分指数；Lo IAQI 表1中与Lo BP 对应的空气质量分指数；IAQI空气质量分数n污染物项目S半集均方差；hP大于中位数半集的分指数；im大于中位数半集的分指数个数；P平均数；P加权后得到的综合指数；Q对各个i，S修正后的半集均方差；h4 问题分析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）算术均值（亦即一阶原点矩）是典型的也是基本的集中趋势测度，其定义式为：n i n X X i i,...2,1;==∑式中， i X 为样本观测值；n 为观测值个数。

算术均值不仅是度量数据离散程度的有效基准,还可以反映样本的总体水平。

因此, 算术均值是空气质量评价中的一个重要指标。

（2）样本观测值的离散程度, 最基本的度量是用标准差（或称均方差）,它是数据能变性测度，这里取定义式为：()n X X i i ∑-=2σ；n i ,...2,1=由标准差的定义式可以推论：i X 中最大值的增大或最小值的减小,均可导致标准差值的增大。

联系到空气质量评价, 也就是把远大于算术均值的分指数和远小于算术均值的分指数置于等同地位, 这显然与它们的环境效应相矛盾。

因为分指数中最大值对环境质量的影响是明显的, 而最小值的影响较小甚至可以忽略。

因此, 在空气质量评价中若简单地按定义式引入标准差, 尽管可以说明分指数的离散程度, 却未能真实地反映出各自对水质的影响，因而不能满足空气质量评价的需要。

因此,有必要从能够充分反映的分指数对空气质量影响的特点出发, 寻求满足于空气质量评价要求的数字特征。

（3）在研究各种指标集对某物影响的过程中，不仅指标集中的最大值具有最重要的作用，次大值等的作用也不容忽视，甚至具有与最大值类似的影响。

因此，为了在评价体系中引入次大值等的影响。

提出了用中位数来划分整个评价指令集，对大于中位数的分指数半集给予较大的权重。

联系到在不同的分布形态该半集数列的疏密程度不同,可通过计算它们对算术均值的均方差,以反映大于中位数半集的离散程度。

据此引出一个新的数字特征,可将它定义为半集均方差, 其定义式为：()m P P S m i ih ∑=-=12（4）根据上述分析, 空气质量指数AQI 作为反映大气环境质量的数字特征,应在反映分指数集合集中趋势的基础上,用大于中位数半集的离散程度加以调整。

据此可建立新的空气质量指数（AQI ）的数学模型为： h S P AQI +=该模型可简称为半集均方差模型，该模型既通过算术均值考虑了每项分指数对空气质量的影响, 也通过半集均方差对分指数中的大值给予较大的权重，据此, 该模型是按各参数对大气环境影响程度的不同,综合出大气环境质量状况的数量描述。

5.2 问题二的模型建立与求解由问题一我们得到新的空气质量指数计算模型，但是在实际应用的过程中又发现了半集均方差模型的缺陷。

由图1中的折线图可知：中国城市里PM2.5或PM10的AQI 值远超于其他污染物的AQI 值。