《2011年数学建模A题国一》由会员分享，可在线阅读，更多相关《2011年数学建模A题国一（18页珍藏版）》请在金锄头文库上搜索。
1、2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A我们的参赛报名号为（如果赛区设置报名号的话
2、）：所属学校（请填写完整的全名）：兰州交通大学参赛队员（打印并签名）：1.2.3.指导教师或指导教师组负责人（打印并签名）：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要本文主要针对城市表层土壤重金属污染进行分析，并根据其传播特征建立模型以确定污染源的位置，更好地研究城市地质环境的演变模式。文中采取的方法是首先通过单因子综合污染指数和内梅罗综
3、合污染指数相结合的方法确定污染源，然后再采取的多项式拟合方法验证。对于问题一：主要采用绘制切片图的方法对空间分布情况进行表述，首先根据已知数据，采用内插法，运用MATLAB绘制三维地形图。之后将海拔分为5个区间进行切片绘图和分析，最后，采用折线图的表示方法从宏观的角度加以说明。在对污染程度进行分析时，本文采用污染指数法评价模型，其中包括单因子污染指数法和综合污染指数法。对于问题二：要求找出重金属污染的主要原因，本文从分析重金属污染的主要来源入手，运用因子分析的方法，分别列出变量相关矩阵、特征值和累计贡献率、旋转前因子载荷矩阵、方差极大正交旋转后因子载荷矩阵。分析出交通区和工业区产生的废气和粉尘
4、是重金属污染的主要原因。对于问题三：要求确定污染源，我们可以在基本假设的基础上，定性分析污染源对周边采样点的传播特征。然后用单因子指数法和内梅罗综合污染指数法对各点的污染指数进行评定得出各点的综合污染程度，再用熵值法进行验证，得出（2708,2295,22）、（13694、2357、33）、（15248、9106、16）是受重金属污染最严重的几个点，再根据扩散原理可知重金属污染最严重的几个地区就是污染源。最后，首先用一个污染源和其周边的采样点通过数据拟合建立出污染程度与污染源距离的函数关系（此此处采取的是多项式拟合），然后再用其他污染源和其他采样点对该函数进行修正得出结论。同时我们还对污染程度
5、与海拔高度进行了拟合得出相应的关系以弥补我们未考虑高度影响引起的误差。对于问题四：为了更好的研究城市地质环境的演变模式还需要收集本地多年的重金属污染情况，只有如此才能将各年的情况进行比较得出演变过程。同时我们还应该考虑到本地本身的资源状况，采集点取样应采取不同土层综合取样才能更好的代表该地数据。我们仔细分析了模型的优缺点。为更好地研究城市地质环境的演变模式，还应收集该地区的每年生活、工业等重要污染源的垃圾排放量，以及每年的生物降解量，降雨量对重金属元素扩散的影响，空气污染也应该考虑进去。有了这些数据以后建立因子分析法，回归分析，曲线拟合等模型解决问题。本模型的优点主要是将数学知识与计算机知识相
6、结合且比较简洁，使用起来比较方便，使用面也比较广泛。缺点是有些条件没有考虑，对一些东西只能作出定性的分析，还需要进行优化。关键词单因子指数内梅罗综合污染指数熵值法扩散原理数据拟合_关键词：插值法、均值法、扩散模型、因子分析、回归分析。1问题的提出由于经济的快速发展和人口的不断增加，人类活动对城市环境质量的影响日显突出。环境问题与人类的生活息息相关，良好的环境可以促进人类的身体健康，而如果环境受到污染，则会使人类的健康受到危害。而在土壤的污染物中，重金属元素污染物具有毒性大、不易被土壤微生物等降解，易被生长在土壤中的植物吸收，并通过生物富集作用进入食物链中，对人类的健康产生较大的影响，因此土壤的
7、重金属污染问题受到越来越多的重视，这个问题的研究也正成为环境科学中的一个热点问题。题目中已给出某城市城区土壤地质环境的调查数据。本论文针对不同的功能区，不同的海拔高度和不同的重金属元素，利用已知的采样点数据进行分析，通过绘制切片图的方法，最终得到准确直观的空间分布图。考虑了等因素，通过建立的数学模型，借助于计算机的高速运算与逻辑判断能力，求出其污染源。为了保护人类自身的健康，应该充分研究重金属元素在土壤中分布、来源及传播特征，建立相应的重金属污染物分布、来源及传播特征的模型，通过掌握的污染现状来预测其发展趋势，可以为制定某城区相应的防范和应对措施提供科学依据，以防止这些污染物更大范围的扩散，避
8、免引起更大范围的污染。因此，研究重金属元素在土壤中的分布、来源及传播特征，并对其建立模型，进行仿真计算，具有攸关人类生存、健康与可持发展的现实意义。主要解决的问题如下：(1) 给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2) 通过数据分析，说明重金属污染的主要原因。(3) 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4) 分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？2问题的分析2.1 空间分布情况分析方法为了准确表达该城区内重金属的污染情况，在考虑地形和海拔的基础上
9、，运用切片绘图的方法，将海拔按照分成若干区间，绘制出每一区间不同重金属元素的分布图，并辅以采样点的分布情况。切片图的X轴坐标为附件1中的x(m)列的数据值，即为采样点的横坐标；切片图的y轴坐标为附件1中的y(m)列的数据值，即为采样点的纵坐标。图中线条的密集程度表示该元素在特定海拔区间的分布情况。2.2 污染程度评价方法与标准对土壤环境质量进行评价时所采取的评价方法多种多样，主要有污染指数法、污染程度法，T值分级法、基准分级法、密切值法和模糊数学综合评判法、灰色聚类法等。文中在对研究区原始测试数据分析评价的基础上，采用污染指数法对该地区土壤的环境质量进行现状评价。污染指数法具有一定的客观性和可
10、比性，且易于计算，己在环境质量评价中得到了广泛的应用。在采取污染指数法进行该区土壤环境质量评价时，参照中国土壤环境质量标准(GB156181995)来确定评价标准。所采用的污染指数法评价模型包括单因子污染指数法和综合污染指数法。本文主要以单因子污染指数法为依据进行分析。2.3 重金属污染的主要原因分析要想分析重金属污染的主要原因，就要从污染物的主要来源入手进行分析。本文中所研究城区按功能分为生活区，工业区，山区，交通区，公园绿地区。主要重金属污染来源有个方面，即：居民生活污染源、工业生产污染源、交通污染源。研究时采用多元统计数学方法之一的因子分析，它根据多个实测变量之间按的相互关系，运用数学变
11、换，将多个变量转变为少数几个现行不相关的综合指标，从而简化数据处理。其目的在于对大量观测数据，用较少的有代表性的因子来说明众多变量所提取的主要信息，提出多个变量间的因果关系。因子分析从变量的相关矩阵出发将一个m维的随机向量X分解成低于m个且有代表性的公因子和一个特殊的m维向量，使其公因子数取得最佳的个数，从而使对m维随机向量的研究转化成对较少个数的公因子的研究。设有n个样本，n个指标构成样本空间X，X=(Xj)mxn，i=I，2，n；j二1,2,，m(1)因子分析过程一般经过以下步骤：(1) 原始数据的标准化，标准化的公式为Xj=(Xj-Xj)/5j其中X.为第i个样本的第j个指标值，而X.和
12、5.分别为j指标的均值和标准ijjj差。标准化的目的在于消除不同变量的量纲的影响，而且标准化转化不会改变变量的相关系数。(2) 计算标准化数据的相关系数阵，求出相关系数矩阵的特征值和特征向量。(3) 进行正交变换，使用方差最大法。其目的是使因子载荷两极分化，而且旋转后的因子仍然正交。(4) 确定因子个数，计算因子得分，进行统计分析。因子分析只强调变量的离差(变化量)而不强调变量在样品中的比重(百分含量)。因子分析的数学模型中，通过正交的方差最大旋转法使每一个主因子只与最少个数的变量有相关关系，而使足够多的因子负荷均很小。变量或因子的重要程度都是以其方差大小来衡量的。因子旋转后每个变量因子负荷代
13、表着在系统中作用或重要性程度，以各个变量目标因子载荷平方与因子方差贡献率乘积作为变量的权重，构成一个判别污染来源的综合指标，而且因子分析是一个客观计算同主观思维相结合的过程。其它多元统计分析(如判别分析，回归分析)的计算结果基本上是一个最终结果，可以直接予以应用，但因子分析的计算结果(因子解)只能看作是一个中间结果，剩下的部分要求人们用自己的思维来完成，这就涉及环境地球化学知识、经验，甚至于思维方式和哲学思想。2.4 重金属污染物的传播特征及模型分析3.模型假设(1)(2)4.符号说明C为浓度(mg/kg)；Xj为j指标的均值；j为j指标的标准差；C为土壤中污染物i的实测数据；iS为污染物i的
14、评价标准iP为内梅罗综合污染指数；综合5.模型的建立及求解5.1重金属元素的空间分布图及结论5.1.1重金属元素空间分布图为了方便描述重金属元素空间分布情况，首先模拟出所研究城区的地形图。具体方法为，建立三维直角坐标系，按照附件一所给出的数据，运用内插法，借助MATLAB绘图功能，模拟出地形图如图1所示。400,300.200.出100-.图1在了解地形的基础上，我们运用MATLAB做出采样点的水平面散点图（图2）,作为分析空间分布情况的考虑要素。由图2可得，5类区的划分有交叉点的存在，功能区之间并不独立。生活区工业区山区克通区公园绿地旄图2在分析了地形和采样点分布情况之后，将用Excel处理
15、过的数据导入至软件Surfer8.0中绘制切片图。首先，将海拔分成5个区间，分别为040、4080、80120、120160,160230,单位为米，m。之后，x轴、y轴分别为x(m),y(m)的值，绘制切片坐标。图3表示的是As元素和Cd元素的空间分布情况。海拔在040m的范围内，该城区西南地区的Cd重金属污染浓度较大，主要集中在工业区和交通区。随着海拔的上升，在4080m的范围内，重金属污染浓度相对减轻，但仍是工业区和交通区较高。两图对比来看，Cd的浓度普遍高于As元素。图3图4表示的是Cr元素和Cu元素的空间分布情况。海拔在040m的范围内，该城区西南地区的重金属污染浓度相对较重，主要集中在工业区周围，属于交通区和生活区的采样点Cr、Cu浓度也相对较大。随着海拔的上升，在4080m的范围内，重金属污染浓度相对减轻。两图对比来看，Cu的浓度普遍高于Cr元素。图4图5表示的是Hg元素和Ni元素的空间分布情况。海拔在040m的范围内，Ni元素的浓度较低，Hg元素的浓度较高，主要集中在工业区、交通区和生活区。随着海拔的上升，在4080m的范围内，重金属污染浓度相对减轻，两图对比来看,Hg的浓度普遍高于Ni元素。图5图6示的是Pb元素和Zn元素的空间分布情况。海拔在040m的范围内，Pb的污染主要集中在西南
