雾霾天气与城市工业污染的关系.docVIP

雾霾天气与城市工业污染的关系 课程名称 数学模型 学 院 机械工程学院 年级 13级 专业班 机自实验2班 学 生 姓 名 肖方杰 学 号 开 课 时 间 2014至 2015学年第1学期 题目 雾霾天气增多是否说明城市工业污染加重，请收集数据，建模来分析，并验证模型的合理性。 问题重述与分析 要研究雾霾天气与城市工业污染的关系，可以通过研究空气质量指数（AQI）与PM2.5（含量）与PM10（含量）、CO（含量）、SO2（含量）、NO2（含量）的关系来说明，因为雾霾天气主要由空气中PM2.5等污染物形成，AOI能根据PM2.5等污染物来评价大气质量，而城市工业污染又是造成PM2.5等污染物增加的原因，所以可以通过分析建立AQI 与PM2.5（含量）与PM10（含量）、CO（含量）、SO2（含量）、NO2（含量）的模型来说明雾霾天气与城市工业污染的关系。 一、数据整理与分析 因为北京市的雾霾天气较多，所以本文模型是以北京市10月份的天气数据来建立，从而具有较高准确性和指导意义。 (1)数据来源 本文数据来源于某天气统计网站（/aqi/beijing-201410.html） 写程序时要把各变量单位统一 变量说明 符号 符号说明 Y 空气质量指数（AQI） x1 PM2.5含量 x2 PM10含量 x3 CO含量 x4 NO2含量 x5 SO2含量 模型建立和求解 本题建立多元线性模型求解,多元线性模型表示为： Y=b0+b1*x1+b2*x2+b3*x3+b4*x4=b5*x5+ 回归系数b0,b1,b2,b3,b4,b5由数据估计，是随机误差。 MATLAB程序如下； y=[91,75,119,101,48,47,179,309,380,351,233,25,29,72,81,73,127,194,289,254,67,100,154,263,341,85,50,84,138,203,201]; x1=[66,54,89,75,29,25,139,258,331,300,190,6,11,43,43,33,93,15,238,203,33,71,116,212,291,58,24,61,104,155,152]; x2=[77,51,108,45,40,50,181,292,354,268,126,8,26,90,99,93,176,229,307,226,81,135,175,260,327,61,58,109,160,194,180]; x3=[970,950,1590,1130,700,460,1290,2090,2290,2360,1800,250,310,690,590,630,1210,2000,2390,1550,1210,1260,1560,2060,2450,670,630,1000,1290,1680,2120]; x4=[57,46,52,38,29,46,86,102,92,84,60,22,30,61,43,54,104,113,118,75,56,75,87,87,88,39,52,71,88,99,88]; x5=[4,4,4,2,2,2,12,9,5,12,6,2,4,7,4,7,11,17,12,4,19,10,12,10,6,2,7,11,12,16,10]; figure(1);plot(x1,y,o),title(AQI与PM2.5含量的关系’) hold on figure(2);plot(x2,y,+),title( AQI与PM10含量的关系) hold on figure(3);plot(x3,y,*),title( AQI与CO含量的关系) hold on figure(4);plot(x4,y,.),title( AQI与NO2含量的关系) hold on figure(5);plot(x5,y,^),title( AQI与SO2含量的关系) %作散点图 hold on n=31; m=5; X=[ones(n,1),x1,x2,x3,x4,x5]; [b,bint,r,rint,s]=regress(y,X); b,bint,s, figure(6);rcoplot(r,rint) %作残差分析图 运行结果如下所示： b = 9.5986 0.4558 0.4662 0.0416 -0.4182 -0.3846 bint = -12.8205 32.0177 0.2634 0.6483 0.2046 0.7278 0.0170 0.0662