数学建模—大气污染预报问题.docxVIP

学生数学建模竞赛第一次预选赛 一、（必做题） （1）油罐的体积（本题 10 分） 一平放的椭圆柱体形状的油罐，长度为 L，椭圆的长半轴为 a，短半轴为 b，油的密度为 ρ，问当油罐中油的高度为 h 时油量是多少？ 解：由题意可话画出画出几何图形如图 1 所示 2b h 2a L 图 1.1 x a cost 椭圆方程为 y bsin t 如图 2，设阴影部分面积为 S/2 ，则油桶的底面积为 S。 Y b b-h a X 图 2 下面将会利用 mathematics 5.0 软件进行求解 , 求解的程序如下： Integrate[2*a*b*Cos[t]^2,{t,ArcSin[1-h/b],Pi/2}] 解得结果为： S a( (2b h)h ( b h) b arccos(1 h)) b2 b 当 h b 时，由椭圆对称性， A 中的 h 用 2b h 代替得到： (2b h)h h S ab a( (b h) b arccos( 1)) b 2 b 所以油液质量 M为： 若 h ，则 M (2b h)h h) h M SL b2 b ( 2b h)h h 若 h ，则 M abL aL b2 b （2）光的反射定律（本题 10 分） 费马原理：光总是沿用时最短的光程传播。试根据这一原理利用极值的有关知识证明光的反射定律：入射角等于反射角。 解：由于光在同一介质中的速度为常数，所以在同一介质中光总是沿直线传播。 如图 3，现假设有两种介质 1、 2 相接，光线在介质 1 中的传播速度为 v, 取两介质的分界线上的一条直线为 X 轴，设有一束光线从介质 1 中的 A(0, a) 点经 X 轴上的 P(x,0) 点反射，并沿直线方向行进到 B( d, b) 点。设直线 AP与 X 轴法线的夹角为 1 ，PB直线与 X 轴法线的夹角为 2 ，下面，根据最短时间效应来推导出光学中的反射定理。 P Y A B（ d,b） 1 2 0 P X 图三 光线由 A 点传到 P 点所需的时间为 : AP a2 x2 t1 v v 1 光线由 P 点传到 B 点所需的时间为 : t 2 PB b2 (d x)2 v v 故光线由 A 传到 B 所需的总时间为： t t1 t 2 a 2 x2 b2 (d x) 2 v 根据费马定理，最短时间效应对应的优化问题为： min t a2 x2 b2 ( d x) 2 (x [ 0, d ]) v 令 dt 1 [ x x2 d x ] 0 dx v a 2 b2 ( d x)2 于是可以得到： x d x a 2 x2 b2 ( d x)2 又由于 x sin 1 ， d x sin 2 a2 (d x)2 x2 b2 所以有： sin 1 sin 2 这就是光学中的反射定理。 证毕 2 大气污染预报问题 摘要 本文通过对四个城市的空气质量的排名以及城市 A的空气质量，利用 C语言、Excel 、 Mathematics 和 MATLAB等工具，分别建立了 层次模型 、多元线性回归预测模型 进行了合 理地分析。最后，我得到了以下一些比较满意的结果。 问题（1）：通过对问题（ 1）问题的分析，得出了这是一个比较典型的层次模型， 目标层是空气质量的排名 ，因素是三种污染物的浓度情况 ，对象是题目给出的 4 个城市 。查找资料后，我找到了一个非常关键的东西—— 空气污染指数的计算方法 ，于是利用 C 语言的编程知识我很快求出了这些城市的污染情况，结合层次模型的相关知识，建立层 次模型后很快得到了我们所期望的答案： 总权重 ： A-0.1374,B-0.1301 ， C-0.1028 ， D-0.6298. 于是城市环境排名 也就解决了，由优到劣的排名情况： C、B、 A、D。而且层次的模型的 一次性检验 也顺利通过。 问题（2）：问题（2）要求我们找出空气质量与气象因素之间的关系，一开始查阅了很多资料， 本想借助灰色预测模型进行求解，可是灰色预测模型的使用条件和咱们这个题目的要求似乎没有什么关联， 后来在网上浏览一片文章的时候，我找到了问题的突破口，便是利用多元回归预测模型进行求解 。然后根据这个模型的所要的处理数据，利用MATLAB、Excel 等工具，求解到本题的回归系数。得到了三种污染物与气象因子之间的 关系： SO2的浓度与气象因素的关系： c1 1.0173 0.0015p 0.0019t 0.0004f 0.0044v NO2的浓度与气象因素的关系： c2 0.2022 0.0003p 0.0017t 0.0001f 0.0077v PM10的浓度与气象因素的关系： c3 0.9329 0.0014p 0.0033t