数学实验论文-城市空气污染数据的真实性判别及分析研究精选文档

《数学实验论文-城市空气污染数据的真实性判别及分析研究精选文档》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、城市空气污染数据的真实性判别及分析研究1.摘要空气质量问题始终是政府、环境保护部门和全国人民关注的热点问题。2016年的两会上,全国政协常委、环境保护部副部长吴晓青表示，政府工作报告中提出的今后五年地级市及以上城市空气质量优良天数比率超过80％的目标必须完成。其中，优良天数比率指的是：区域内城镇空气质量优良以上的监测天数占全年监测总天数的比例。针对此，我们小组刚好找到了浙江舟山的相关数据，AQI指数与影响空气质量的污染物（PM2.5、PM10、CO、NO2、SO2）并用matlab对此进行分析。关键词：多元线性回归模型；散点图；层次分析；残差分析；相关性

2、分析2.模型建立2.1问题分析针对按照常识，空气污染指数随影响空气质量的污染物（PM2.5、PM10、CO、NO2、SO2）的上升而增大。空气质量指数AQI为y，PM2.5的含量为x1，PM10的含量为x2，CO的含量为x3，NO2的含量为x4，SO2的含量为x5。利用matlab分别作出y对x1、x2、x3、x4、x5的散点图，如后图所示，所以假设y与x1、x2、x3、x4、x5的关系为线性关系。图2-1-1Y对x1的散点图图2-1-2Y对x1的散点图图2-11Y对x2的散点图Y对x2的散点图2.2符号假设r:残差向量rint：r的置信区间b:参数估计

3、值bint：b的置信区间：复相关系数F:检验统计量P：F统计量对应的概率：剩余平方和State：回归模型的检验统计量β：回归系数向量β=（，，，，，）Y：空气质量指数向量X1:PM2.5的含量（直径小于2.5um可吸入颗粒）X2：PM10的含量(直径小于10um的可吸入颗粒)X3:CO的含量（CO为一氧化碳）X4：NO2的含量（NO2为二氧化氮）X5：SO2的含量（SO2为二氧化硫）ε：随机误差p：相关系数矩阵2.3模型建立模型一的建立空气质量指数AQI为y与PM2.5的含量为x1，PM10的含量为x2，CO的含量为x3，NO2的含量为x4，SO2的含量

4、为x5之间的多元线性回归模型为y=+x1+x2+x3+x4+x5+ε其中，，，，，是带估计的回归系数，ε是随机误差。3.模型计算3．1问题一的计算利用matlab命令：[b,bint,r,rint,s]=regress(y,X)，rcoplot(r,rint)得到残差分析表图2-1-3再由残差置信区间得出异常数据第8、25、26、31、33、36、37、38、39、40、50、51、56、57、79、223、229、276组得到回归模型的系数（置信区间95%）表3-1回归系数回归系数估计值回归系数置信区间16.0434[11.993020.0938]0.

5、0819[-0.03950.2033]0.6750[0.58780.7622]13.9480[7.877120.0189]-0.2308[-0.3666-0.0951]-0.3412[-0.5515-0.1309]R2=0.8221F=439.0060p<0.0001s2=171.8599有表中数据可知的置信区间包含零点，所以这个系数的解释是不可靠的。再用matlab中的corrcoef命令直接得到相关系数矩阵1.00000.93100.72220.72940.70940.85420.93101.00000.66620.77140.75430.89480

6、.72220.66621.00000.57310.51470.66700.72940.77140.57311.00000.71830.64120.70940.75430.51470.71831.00000.61940.85420.89480.66700.64120.61941.0000又表中数据可知X1与X2关系显著，所以X1多余，应当去掉。去掉X1后再用上述方法得到残差分析表图2-1-4由残差置信区间得出异常数据第8、25、26、31、33、36、37、38、39、40、50、51、56、57、79、223、229、276，321组得到回归模型的系数（

7、置信区间95%）表3-2回归系数回归系数估计值回归系数置信区间15.0105[11.257918.7631]0.7207[0.66580.7756]15.4721[9.833121.1111]-0.2314[-0.3672-0.0955]-0.3380[-0.5484-0.1276]R2=0.8214F=547.4489p<0.0001s2=172.1328得到线性回归方程y=15.0105+0.7207x2+15.4721x3-0.2314x4-0.3380x5数据真实性判定:表3-3组数AQI指数PM2.5PM10CONO2SO2异常点（为0异常）3

8、015531490.75174321555110.544303475527540