关于北京市气候变化与大气污染的研究中央民族大学马越龙、姚琨、朱爱丽摘要本文建立了一个关于分析北京市气候变化趋势和大气污染情况的统计模型。
本文主要研究以下三个问题：北京市气候变化的趋势，验证北京市是否存在气候变暖现象；大气污染情况及变化趋势；大气污染和气候变化之间的关系。
首先就问题一，根据主成分分析思想，我们选用降水量、气温和风速来表征北京市气候变化，并从《北京市统计年鉴2010》获得了1978——2009年的相关数据，利用Matlab、Excel对这些数据进行时间序列分析并根据最小二乘法原理将散点图拟合成趋势曲线，又根据问题实际需要和判定系数等对其进行筛选，得到我们需要的最优回归方程及曲线。
据此验证北京市存在气候变暖现象，并选取气温的气候倾向率（气温线性回归方程的一次项系数）来度量气候变暖的程度。
其次就问题二，根据主成分分析思想，我们选用SO2浓度、NO2浓度、可吸入颗粒物含量来表征北京市大气污染状况，并从《北京市统计年鉴2010》获得了2000——2009年的相关数据。
利用内地空气污染指数的计算方法，计算出大气污染的评价度量指标API。
并利用问题一中的方法分析各数据，得到SO2浓度、NO2浓度、可吸入颗粒物含量及API的回归模型，从而分析度量大气污染情况的变化趋势。
然后就问题三，利用SPSS软件计算分析API与温度、降水量、风速的相关性。
通过分析比较相关系数得出大气污染与气候变化之间的关系。
最后我们利用F检验来判断以上所建立回归模型的回归效果。
综上，我们利用SPSS、Excel、Matlab等数学软件很好地实现了数据统计分析，并结合主成分分析及时间序列分析思想给出了北京市气候变化和大气污染变化趋势的分析。
将数据与图形相结合，既充分发挥了图形的直观功能，又根据数据进行了相关分析，排除了主观经验的干扰，使所建立的数学模型能较好的解决上述问题。
关键词：气候变化大气污染变化趋势API 回归模型F检验一、问题的提出与分析1、问题的提出目前，气候变暖和大气污染成为当今人类社会亟待解决的两大问题，是人类必须面对的严峻挑战。
北京作为中国的首都，因其特殊的政治、经济地位，在国家乃至世界的发展中都占有特殊地位，因此研究北京的气候变化和大气污染及其之间的关系有着特别重要的意义。
气候是地球上某一地区多年时段大气的一般状态，是该时段各种天气过程的综合表现。
气象要素(温度、降水、风速等)的各种统计量(均值、极值、概率等)是表征气候的基本依据。
大气污染是指自然或人为原因使大气中某些成分超过正常含量或排入有毒有害的物质，对人类、生物和物体造成危害的现象[1]。
按中国《大气环境质量标准》规定的常规分析指标有总悬浮微粒、二氧化硫、氮氧化物、一氧化碳和光化学氧化剂。
现在我们所关心的问题就是：（1）对获得的大量数据怎样进行合理有效的筛选及分析，进而建立适当的数学模型，来分析北京市气候变化的趋势，验证北京市是否存在气候变暖现象并选择一个指标来度量气候变暖的程度。
（2）导致北京市大气污染的因素众多，如何选取适当的指标评价和度量大气污染情况及变化趋势。
（3）如何依据现有数据分析大气污染和气候变化之间的关系。
（4）如何选取适当的方法在模型建立并求解之后进行模型的检验。
2、问题的分析解决上述问题需要大量的数据，结合主成分分析思想筛选数据、掌握数据所代表的实际意义及变化趋势和数据间的联系至关重要。
为此我们需要结合时间序列分析思想对表征北京市气候和大气污染的要素做个别分析及联合分析。
根据问题的要求，我们分三个问题进行分析。
问题一：对气候变化趋势的分析，可以先搜集历年的气候数据，选取适当的指标作散点图，采用一元回归的的方法建立模型，观察其波动性及变化趋势。
问题二：对大气污染的度量和变化趋势的分析，可以先搜集历年相关污染物数据，并计算出适当的指标来度量大气污染情况，然后同问题一的方法进行变化情况分析。
问题三：可以研究度量大气污染情况的指标与表征气候变化的要素之间关系，并用其来代表大气污染与气候变化之间的关系[2]。
二、基本假设1.温度、降水量、风力可以完全表征北京市气候2.假设北京市每个区的气候差异不大，北京市气象站的所测得的数据可以代表北京市整体的气候状况3.在回归模型中，拟合的二次多项式可以较好模拟出数据的变化趋势4.在变暖程度分析中，拟合的一次多项式系数可以作为气温的气候倾向率来进行变暖程度的度量5.可吸入颗粒物、SO2和NO2可完全代表北京市大气污染的污染物6.API指标能正确评价和度量大气污染情况7.假设气温和API都是正态分布，或接近正态的单峰分布，从而可以用皮尔逊相关系数表示二变量间的相关关系三、符号说明1.API:空气污染指数，具体分级标准见附录五2.I：某污染物的实际污染指数3.C：该污染物在大气中的实际的浓度4.C大与C小：在API分级限值表（附录五）中最贴近C值的两个值，C大为大于C的限值，C小为小于C的限值5.I大与I小：在API分级限值表（附录三）中最贴近I值的两个值，I大为大于I的值，I小为小于I的值6.API＝max(I1,I2，…，I n)：取所有污染物污染指数的最大值作为API的实际值，这里只取n=37.2R：回归模型的判定系数8.x:气候特征值9.a:气候倾向率四、模型的建立与求解4.1利用实际数据分析北京市气候变化趋势，验证北京市存在气候变暖现象，并给出变暖程度的度量4.1.1数据总体分析，并画出散点图经查阅《北京市统计年鉴2010》，我们得到了1978——2009年的气候数据（见附录一），由模型假定，我们选取温度、降水量、风速表征气候情况，并给出了对应的描述统计量（见表1）及散点图（见附录二）。
1978年-2009年气候数据描述统计量表14.1.2气候变化回归模型的建立及气候变化趋势的分析为分析北京市气候变化趋势，我们根据最小二乘法原理，使用Excel作了气候各表征量的各种回归分析并作出了相应的趋势变化图，包括线性回归、对数回归、多项式回归、乘幂回归、指数回归。
根据判定系数2R的大小，我们发现多项式的拟合效果最佳且与阶次呈正相关。
又由于实际分析趋势的需求，我们不需要阶次>=3的曲线拟合，最终选取了二阶多项式非线性回归作为数学模型来近似表征气候变化趋势[3]。
相应变化趋势图及二阶多项式、2R见下图。
图1图2图3气候变化趋势分析：从趋势线看，北京市1978年以来年降水量减少了200毫米以上，并继续呈缓慢下降趋势；平均气温呈上升趋势，在31年中上升了1摄氏度左右，且上升速度趋缓；风速虽有波动，但总体平稳。
4.1.3验证北京市是否存在气候变暖现象，如果存在，给出变暖程度的度量根据图2可以看出北京确实存在气候变暖现象。
用气温的气候倾向率作为变暖程度的度量。
随着时间序列的增长, 对要素序列的各值平均而言, 或是增加或是减少, 形成序列在相当长时期里向上或向下的缓慢移动, 这种有一定规则的变化成为趋势。
序列的趋势变化[4]可以用一次线性回归方程表示:x=at+b式中: x 为年气候特征值; t 为年; a 则为气候倾向率, 用于定量描述序列的趋势变化特征。
经过对北京市1978— 2009 年气温的系列值回归分析得到一次线性回归方程：x=0.0593t-105.45并用Excel作出了趋势变化图，见下图。
图4由此得出, 北京市气温的气候倾向率为0.0593, 说明在这31年中北京市气温以0.0593摄氏度/年的速率缓慢上升。
4.2分析北京市大气污染情况，选取适当的指标评价和度量大气污染情况及其变化趋势4.2.1数据总体分析，并画出散点图2表4.2.2 API的计算API（Air Pollution Index的英文缩写）[5]是空气污染指数，就是将常规监测的几种空气污染物浓度简化成为单一的概念性指数值形式，并分级表征空气污染程度和空气质量状况，适合于表示城市的短期空气质量状况和变化趋势，也是目前世界上许多国家评估环境空气质量状况的一种方式。
根据我国空气污染的特点和污染防治重点，目前计入空气污染指数的项目暂定为二氧化硫、氮氧化物和可吸入颗粒物或总悬浮颗粒物。
我们选取其作为评价和度量大气污染情况及其变化趋势的指标。
空气污染指数的计算方法：API=max(I1,I2…Ii…In）基本计算式：设I为某污染物的污染指数，C为该污染物的浓度。
则：式中：C大与C小：在API分级限值表（见附表三）中最贴近C值的两个值，C大为大于C的限值，C小为小于C的限值。
I 大与I小：在API分级限值表中最贴近I值的两个值，I大为大于I的值，I小为小于I的值。
由附录三得到各种污染物的污染分指数见表3。
年份SO2NO2PM10（污染分指数）（污染分指数）（污染分指数）2000 60.5 44.375 1062001 57 44.375 107.52002 58.5 47.5 1082003 55.5 45 95.52004 52.5 44.375 99.52005 50 41.25 962006 51.5 41.25 105.52007 47 41.25 992008 36 30.625 862009 34 33.125 85.5表3根据API = max(I1,I2…Ii…In），得到API年均值见表4。
年份API2000 1062001 107.52002 1082003 95.52004 99.52005 96表44.2.3大气污染情况回归模型的建立及其变化趋势的分析及度量为分析北京市大气污染情况及其变化趋势，我们根据最小二乘法原理，使用Excel 做了大气污染各表征量的污染分指数及API的各种回归分析并作出了相应的趋势变化图，包括线性回归、对数回归、多项式回归、乘幂回归、指数回归。
根据判定系数2R的大小，我们发现多项式的拟合效果最佳且与阶次呈正相关。
又由于实际分析趋势的需求，我们不需要阶次>=3的曲线拟合，最终选取了二阶多项式非线性回归作为数学模型来近似表征大污染变化趋势势。
相应变化趋势图及二阶多项式、2R见下图。
图67图图图9、北京市大气污染情况及其变化趋势分析：从趋势线看，近10年北京市空气的SO2、可吸入颗粒物含量[6]均呈下降状，分别下降了约20、10、20个单位。
总体来看，NO2即从API分析与度量，API以2.1424单位/年的速度呈递减趋势，10年中下降20个单位以上，也就是说，北京市空气质量日益进步，大气污染程度降低，并根据“空气污染指数范围及相应的空气质量类别”（见附录五）可知北京市空气质量等级由“轻度污染”转为“良”，大气污染治理情况比较乐观。
4.3研究大气污染和气候变化之间的关系通过散点图可以判断两变量之间有无相关关系，并对变量间的关系形态作大致的描述，但散点图不能准确反映变量之间的关系强度。