城市表层土壤重金属污染分析大学生数学建模论文

PAGE第2页，共1页2011高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则.我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的,如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规则的行为，我们将受到严肃处理。我们参赛选择的题号是（从A/B/C/D中选择一项填写）：A 我们的参赛报名号为（如果赛区设置报名号的话）：所属学校（请填写完整的全名）：参赛队员(打印并签名)：1.2.3.指导教师或指导教师组负责人(打印并签名)：日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：2011高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：城市表层土壤重金属污染分析摘要：本文考虑浓度与坐标之间的关系，通过类比法建立模型。运用Matlab中的命令求解模型和寻找重金属源的位置，再运用图形数据分析建立对类比法的补充分析模型。在对重金属浓度污染程度进行分析时，本文采用单因子指数法，综合指数法（内梅罗指数法）进行分析。在对模型的改善上本文加入时间变量建立类似人口模型进行分析。对问题一，利用Matlab绘出了各种重金属的浓度与坐标的关系，基于直观的重金属的空间分布情况，分别采用单因子指数法（）和综合指数法（），对重金属浓度进行分析，最终得出污染程度的结果如下：区域生活区工业区山区主干道路区公园绿地区划分的等级中度污染重污染轻度污染重污染中度污染对问题二，综合问题一中所得到的结论，可以看出污染程度排在前三位的是工业区，主干交通区，居民区。通过建立了地质累计指数（）模型进行分析，得出重金属污染的主要原因是来自工矿企业污染源，交通污染，居民活动污染。问题三：在对于问题三的分析上，本文只考虑浓度与坐标之间的关系，运用类比法建立模型，再利用Matlab中命令的求解模型和重金属源的位置，最后本文运用图形数据分析建立对类比法的补充分析模型。对问题四，为了研究城市地质环境的演变模式，在问题三的模型上又引进了时间参数。本文认为重金属的浓度与时间的关系与人口增长模型相类似，以此来建立相应的数学模型进行分析。关键词：，类比模型，相似函数模拟，地质累计指数，单因子指数法，综合指数法，Matlab一、问题重述A题城市表层土壤重金属污染分析随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距1公里左右的网格子区域，按照每平方公里1个采样点对表层土（0~10厘米深度）进行取样、编号，并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。附件1列出了采样点的位置、海拔高度及其所属功能区等信息，附件2列出了8种主要重金属元素在采样点处的浓度，附件3列出了8种主要重金属元素的背景值。现要求你们通过数学建模来完成以下任务：(1)给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。(2)通过数据分析，说明重金属污染的主要原因。(3)分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式，还应收集什么信息？有了这些信息，如何建立模型解决问题？二、问题分析和建模思路在对本题进行分析时，为了建立适当的模型求解出重金属源的位置。在对数据进行分析时发现每种重金属的源都不止一个，因而本文对数据进行了相应的处理（本文将数据按重金属源进行了分类，见附表中）。本文分析到重金属源的污染主要是影响周围的区域，对于较远处就没有什么影响。得出重金属的传播特性应该是一个小范围的扩散问题。问题一：在对问题一的分析上，本文首先采用Matlab绘出了各种重金属的浓度与坐标的关系，然后分别采用单因子指数法（）和综合指数法（）对重金属浓度进行分析。从而判断每个区域相应重金属的污染程度。问题二：在对问题二的分析上，本文不但运用问题一中所得到的结论，发现工业区，主干交通区，居民区的污染程度在前三位。还建立了地质累计指数（）进行分析。最终得该城区内重金属污染的主要原因。问题三：在对于问题三的分析上，本文只考虑浓度与坐标之间的关系，运用类比法建立模型，再利用Matlab中的命令求解模型和重金属源的位置，最后分析图形和数据进行相应的改善。问题四：在对于问题四的分析上，为了研究城市地质环境的演变模式，在问题三的模型上又引进了时间参数。在现实生活中随着经济的发展，重金属的产生率开始一段时间会越来越大，但是当人们意识到问题时，就会采取相应的措施，这是重金属的产生率就会减少。本文认为随着时间增长，浓度是一个先增长率加大，然后增长率减小的过程。最总浓度达到一个最大值。在此本文采用人口增长模型进行分析。三、符号说明与基本假设3.1符号说明：是区域内土壤中重金属的单因子指数：是区域内土壤中测得的重金属的浓度：是土壤中重金属的评价标准（即是它的背景值）：是指区域的综合评价结果：是区域内土壤中单因子值最大的那个数：可以取得值是1,2,3,4,5在下标中分别代表生活区、工业区、山区、主干道路区及公园绿地区五个区域：可以取得值是1,2,3,4,5,6,7,8在下标中分别代表As，Cd，Cr，Cu，Hg，Ni，Pb，Zn：是地质累计指数：是指区域中元素的实测值（此处为平均值）：是元素背景值1.5：是指地质累计指数中的修正系数：表示重金属的横坐标：表示重金属的纵坐标：表示模型一中坐标是（，）点处的重金属浓度：表示重金属源的浓度：表示重金属源的横坐标：表示重金属源的纵坐标：表示重金属的沉淀系数：是一个常数：是一个常数：表示是在一定区域内的点到题目中所给数据（此处的数据也是这个区域内的数据，本文在后面对本题数据进行了相应的划分处理，分成了许多小的区域，见附表问题三数据）最大点的距离：表示时间：表示时刻坐标为（，）处重金属的浓度：表示坐标为（，）处重金属的浓度的最大值：表示坐标为（，）处重金属的浓度的初始值：表示重金属浓度随时间变化的增长率3.2基本假设1.本文假设在污染重金属的传播过程中，不考虑风等自然因数的影响。2.本文假设题目中所给的数据都是真实可靠的。3.本文假设坐标，山高，时间等因数相互之间是没有影响的。它们都各自单独影响重金属的浓度。4.假设浓度的分布与水平高度没有关系。5.假设重金属的传播只对污染源周围的浓度有较大的影响，对于较远处本文在此不加考虑。单因子指数法：是区域内土壤中重金属的单因子指数，是区域内土壤中测得的重金属的浓度，是土壤中重金属的评价标准（即是它的背景值）。小于1表示没受污染，大于1表示受到了污染。越大污染越严重。综合指数法（内梅罗指数法）：式中是指区域的综合评价结果。是区域内土壤中重金属的单因子指数，是区域内土壤中单因子值最大的那个数。当评价参考标准采用研究自然区自然背景值时，上述式中的计算结果分别定义为土壤中重金属的单一和综合富集指数。地质累计指数式中是地质累计指数，是指区域中元素的实测值（此处为平均值），是元素背景值，1.5是指修正系数。四、模型的建立与求解4.1问题一4.1.1各种重金属的空间分布图（Z坐标代表相应的浓度）（图一）（图二）（图三）（图四）（图五）（图六）（图七）（图八）4.1.2单因子指数法本文运用这个公式可以求得数据如下表：区域As(μg/g)Cd(ng/g)Cr(μg/g)Cu(μg/g)Hg(ng/g)Ni(μg/g)Pb(μg/g)Zn(μg/g)11.742.6621.492.33.4322.013.021.729.6718.351.613.004.0361.311.171.261.181.06441.592.771.884.7112.771.432.053.5251.742.161.411.733.291.2421.962.24（表一）4.1.3单因子指数法的等级划分单因子指数<=11<<=22<<=3>3划分的等级无污染轻污染中污染重污染（表二）可以得出每个区域的污染程度如下：生活区因数AsCdCrCuHgNiPbZn划分等级轻污染中污染中污染重污染中污染度污染中污染重污染（表三）工业区因数AsCdCrCuHgNiPbZn划分等级中污染重污染轻污染重污染重污染轻污染重污染重污染（表四）山区因数AsCdCrCuHgNiPbZn划分等级轻污染轻污染轻污染重污染轻污染轻污染轻污染轻污染（表五）主干道路区因数AsCdCrCuHgNiPbZn划分等级轻污染中污染轻污染重污染重污染轻污染中污染重污染（表六）公园绿地区因数AsCdCrCuHgNiPbZn划分等级轻污染中污染轻污染轻污染重污染轻污染轻污染中污染（表七）4.1.4综合指数法本文运用公式可以求得数据如下表：区域生活区工业区山区主干道路区公园绿地区2.4721.181.259.432.71（表八）4.1.5综合指数法的等级划分综合因子指数<=0.70.7<<=1.01.0<<=2.02.0<<=3.0>3.0划分等级清洁尚清洁轻度污染中度污染重污染（数据来自国家土壤质量评价技术规范）（表九）划分得到的等级：区域生活区工业区山区主干道路区公园绿地区划分的等级中度污染重污染轻度污染重污染中度污染（表十）4.2问题二4.2.1地质累计指数运用公式可求得相应的地质累计指数为As(μg/g)Cd(ng/g)Cr(μg/g)Cu(μg/g)Hg(ng/g)Ni(μg/g)Pb(μg/g)Zn(μg/g)10.210.5760.571.320.83-0.00840.571.2020.431.010.202.693.6110.101.001.433-0.42-0.36-0.26-0.19-0.36-0.26-0.35-0.5040.0800.880.321.653.09-0.0670.451.2350.210.52-0.0920.241.136-0.270.3850.58（表十一）4.2.2地质累计指数等级的划分地质累计指数的范围<00--11--22--33--44--5>5划分的等级无污染无污染到中度污染中度污染中度到强度污染强度污染强度污染到极强污染极强污染（表十二）As(μg/g)Cd(ng/g)Cr(μg/g)Cu(μg/g)Hg(ng/g)Ni(μg/g)Pb(μg/g)Zn(μg/g)1无污染到中度污染无污染到中度污染无污染到中度污染中度污染无污染到中度污染无污染无污染到中度污染中度污染2无污染到中度污染中度污染无污染到中度污染中度到强度污强度污染无污染到中度污染中度污染中度污染3无污染无污染无污染无污染无污染无污染无污染-无污染4无污染到中度污染无污染到中度污染无污染到中度污染中度污染强度污染无污染无污染到中度污染中度污染5无污染到中度污染无污染到中度污染无污染无污染到中度污染中度污染无污染无污染到中度污染无污染到中度污染（表十三）4.2.3观点的阐述从第一问得到的结论，可以总体上看出五个地区分别受到重金属的污染程度如下：工业区>主干道路区>公园绿地区>生活区>山区（21.18>9.43>2.71>2.47>1.25）。在五个区的分别讨论中也可以看出污染排名前三位的地方有工业区和主干道路区，居民区。工业区有Cu，Hg，Zn，主干道区有Cu，Hg，居民区有Cu，Zn，这些重金属元素的浓度非常大。在二问的的数据分析中本文也可以得出相识的结论。通过查找大量文献，我们发现土壤中的重金属主要有以下几个来源。第一，大气中重金属的沉淀。第二，农药，化肥，塑料袋的使用。第三，污泥施肥。第四，污水灌溉。第五，含有重金属的固体的堆积。第六，含有重金属的污水的污染。对于工业区，工业污染主要通过工业生产所产生的废渣、废水、废气排入环境，在人和动物、植物中富集。排放到大气中的重金属通过自然沉淀和雨淋沉降进入土壤圈。而废弃物中重金属含量一般较高，污染的范围一般以废物堆中心向四周扩散。金属矿山开采、冶金厂排弃的尾矿渣是矿山主要固体废弃物之一。金属矿山的开采、冶炼、重金属尾矿、冶炼废渣和矿渣堆放等，以及可以被酸溶的含重金属的矿山酸性废水，随着矿山排水和降雨使之带入水环境或直接进入土壤，都可以直接或者间接地造成土壤重金属污染。对于交通区，重金属污染主要来自汽车尾气排放、汽车轮胎磨损产生的大量金属的有害气体和粉尘。它们主要分布在工矿的周围及公路、铁路两侧。对于公园绿地区，重金属污染源主要来自污泥中所含有大量的有机质和氮、磷、钾等营养元素，但同时污泥中也含有大量重金属元素，随着市政府污水处理产生的大量污泥被施加与绿色公园区，使得重金属含量也不断增高。对于生活区，重金属污染程度比较清，重金属污染的主要来源于居民生活垃圾的排放。对于山区，由于在山区，居住人口比较少，其环境近似与测量背景值的标准环境。其重金属污染源较小程度的来自随大气沉淀如土壤中的重金属元素。综上所述，本文分析到重金属污染的主要原因是来自工矿企业污染源，交通污染和居民活动污染。4.3问题三4.3．1重金属污染物的传播特征的分析及模型的建立在对重金属污染物的传播特征进行描述时，本文认为某地重金属浓度的变大是来源于重金属在此地的沉淀累计作用，而且由于重金属的质量比较大，这种沉淀累计作用只在小范围内有作用。从而得出重金属的传播只是在一个小范围与扩散模型相似，因而本文在对污染源的位置进行分析时，考虑在整个区域内存在多个污染源的情况。本文认为在有浓度明显高于其他地方的一个区域内就会存在一个污染源（在污染源个数的确定上，本文在此还结合问题一中所描绘的图像与数据进行对比分析）。本文对于每个污染源都设计相应的污染半径，以此污染半径内的数据来分析污染源所在的位置。（关于相应数据的处理见附表）在此本文引入沉淀系数（沉淀系数与扩散方程中扩散系数类似）。结合扩散方程的理论，本文认为可建立如下模型（本模型考虑到某点重金属浓度与它和重金属源的距离成反比，与接触系数成正比。在本模型中由于每种重金属的重金属源都不止一个，因而对于同种重金属源而言，，，，的值也不唯一）。4.3.2模型一中重金属浓度与坐标的关系及重金属源的位置重金属As的关系及重金属源的位置重金属As污染源序号相应污染源的位置重金属As与坐标的关系1（909.7,357.48）2（468．01，689．31）3（6723．2，6826．9）4（12951，2871）5（18461，10014）重金属Cd的关系及重金属源的位置重金属Cd污染源序号相应污染源的位置重金属Cd与坐标的关系1（21375,11217）2（3996.9，5642.1）3（3242.8,1703.6）4（2216.7，2825.8）重金属Cr的关系及重金属源的位置重金属Cr污染源序号相应污染源的位置重金属Cr与坐标的关系1（3639,5604.3）2（10571，5353）重金属Cu的关系及重金属源的位置重金属Cu污染源序号相应污染源的位置重金属Cu与坐标的关系1（24599,16116）2（27404，16703）重金属Hg的关系及重金属源的位置1（2519,2729）2（13690，2363）3（15247,9099）重金属Ni的关系及重金属源的位置1（3163.6,5881.5）2（22463，12062）重金属Pb的关系及重金属源的位置1（4665.9,5145.5）2（1821，3101）重金属Zn的关系及重金属源的位置重金属Zn污染源序号相应污染源的位置重金属Zn与坐标的关系1（13775,9618）2（9274.6，4135.6）3（3754.7，5637.5）4（12710，3133）5（4812.7，7031.2）4．3．3模型的改进本文在此观察重金属Hg的平面分析图如下：（图中颜色的不同代表了重金属Hg浓度的不同，可以看出浓度不较大的区域有三个）在此将某个区域进行放大如下：从此图可以明显看出Hg元素在它的浓度最大区域内是成面状分布的，本文将它的分布简化为成圆状分布）为了更好的说明Hg元素是成圆状分布，在此处本文建立一个类似正态分布模型进行说明模型如下：（用此模型来拟合Hg元素的局部分布，取5）本文在此计算出Hg从它的局部区域题目浓度最大值点下降到最大值20%的点时，两点之间的直线距离（及的大小）。得出=119.5753。从这个数据可以看出Hg的下降速度是非常的快，这不满足扩散方程，因而此处Hg可以看做是堆积在这个以为半径，局部区域题目浓度最大值的坐标为圆心的圆中。从以上结论本文可以得出在这几种重金属中有些元素的某几个或一个重金属源是成圆状分布的。在此本文定义某元素从它的局部区域题目浓度最大值点下降到最大值20%的点时，两点之间的直线距离为半（此值必须小于200）本文在此对它们进行详细分析结果如下：元素圆心半径Hg（2708,2295）119.58Hg（13694,2357）121.08Hg（15248,9106）117.57因此，我们可以说当（x,y）满足，，时我们就称此地为污染源。综合问题三中的理论，本文可以得出Hg必须要用改善模型进行分析，其它重金属可以用原来的模型分析。4.4问题四4.4.1优点缺点评价优点：本模型能够准确的计算出在重金属源一定范围内各点的重金属浓度，准确求出各个重金属源的位置。模型的计算求解采用了专业的数学软件，计算结果可信度高。文章中的模型原创性强，这些模型是通过类比得出缺点：由于题目中没有给出重金属浓度随时间的变化关系，因而本文上述模型中并没有涉及时间变量，但是在实际情况中，浓度的变化必然与时间有关，而且在没有特定原因的情况下，浓度与时间会呈现先快速增长，然后增长减慢的趋势。4.4.2新模型的提出为更好地研究城市地质环境的演变模式，在此本文应该再收集浓度与时间之间的对应数据，将时间参数加入到模型之中。结合问题三中所建立的数学模型和人口增长模型，建立浓度与时间坐标都有关系的模型如下：模型一：（第一个方程表示的是某一点重金属的浓度随时间的改变，第二个方程表示的是某一时刻任意一点的重金属浓度）第一个方程表示的是某一点重金属的浓度随时间的改变，第二个方程表示的是某一时刻任意一点的重金属浓度。本模型专门用于Hg元素）。五、模型的检验本模型主要是运用了拟合的观点对题目中所给的数据进行了函数关系的拟合，对于模型准确性的检验，本文只需要检验拟合出的函数值与实际函数值之间的差距。本文以为例进行分析。比较结果如下：实际值14.089.628.945.935.72拟合值14.3310.898.155.594.80误差1.8%13.19%8.79%5.81%16.11%从得出的数据来看，本模型能够很好的与题目所给的数据相匹配，因而本模型的建立是合理的。六、模型的优缺点评价与改进6.1模型的优缺点优点：1.本模型能够准确的求出重金属源所在的位置。2.本模型能够计算出污染源周围所有点的重金属浓度。3.文章中的模型原创性强，这些模型是通过类比得出。4.模型的计算求解采用了专业的数学软件，计算结果可信度高缺点：1.模型考虑的变量过少，会导致一些误差；2.题中已知数据较少，模型在解决实际问题时可能需要改进；6.2模型的改进在此模型中本文只考虑了坐标对重金属分布的影响，这使得模型与实际情况之间存在一定的差距。本模型在改进上应该将高度，区域和不同重金属之间存在不同的传播方式等加入到模型的考虑之中去。建立较为完善的模型进行分析。七、模型的推广与相关建议不管黑猫白猫，能抓老鼠的就是好猫。本模型虽然简单，但是本模型能够用于大多数城市中重金属源位置的寻找，方便有关部门对于环境的检测和管理，还有利于公安部门对于违章行为的发现和处理。因而本模型能够很好的被环境和公安等部门所接受和利用。八、参考文献【1】张江华，赵阿宁，王仲复，柯海玲，陈华清，《内梅罗指数和地质累计指数在土壤重金属评价中的差异探讨》，第8期第31卷：43-46,2010年【2】刘申，刘凤枝，李晓华，《天津公园土壤重金属污染评价及其空间分布》，第19期第5卷：1097-1102,2010年【3】廖国礼，周音达，吴超，《尾矿区金属污染浓度预测模型及其应用》，第6期第65卷：1009-1013,2004年【4】廖启林，华明，金洋，《江苏省土壤重金属分布特征及污染源初步研究》，第5期第36卷：1163-1174,2009年九、附录9.1问题一中的程序第一个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('As的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('As浓度')第二个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Cd的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Cd浓度')第三个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Cr的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Cr浓度')第四个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Cu的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Cu浓度')第五个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Hg的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Hg浓度')第六个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Ni的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Ni浓度')第七个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Pb的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Pb浓度')第八个图xlin=linspace(0,30000,319);ylin=linspace(0,20000,319);[X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)>>title('Zn的污染分布')>>xlabel('x坐标')>>ylabel('y坐标')>>zlabel('Zn浓度')9.2问题三中的程序及结果模型的程序functionf=jianmo(x,xdata)f=(x(1)./((xdata(:,1)-x(2)).^2+(xdata(:,2)-x(3)).^2).^0.5)*x(4)%x(1)=k;x(2)=x0;x(3)=y0;x(4)=U(x0,y0)x0=[696.6627,1600,2750,14.50];f=@(x,xdata)(x(1)./((xdata(:,1)-x(2)).^2+(xdata(:,2)-x(3)).^2).^0.5)*x(4)[x,resnorm]=lsqcurvefit(f,x0,xdata,cdata)改进模型的程序functionf=fcurv(x,tdata)f=x(1).*exp(-(tdata).^2./(x(2).^2))+x(3);endf=@(x,tdata)x(1).*exp(-(tdata).^2./(x(2).^2))+x(3);x=lsqcurvefit(f,x0,tdata,cdata)模型As的结果：x=1.0e+003*0.69670.90973.57480.0231resnorm=17.8723x=1.0e+003*0.14164.68016.89310.0625resnorm=0.7329x=1.0e+003*0.09756.72326.82690.0836resnorm=13.0247x=1.0e+004*0.00901.29510.28710.0079resnorm=9.1717x=1.0e+004*0.01311.84671.00140.0077resnorm=2.1860模型Cd的结果x=1.0e+004*0.01762.13751.12170.1638resnorm=1.3109e+004x=1.0e+003*0.65763.99695.64211.4197resnorm=1.9675e+005x=1.0e+003*0.98573.24281.70360.9857resnorm=1.0725e+006x=1.0e+003*0.79622.21672.82580.7962resnorm4.4705e+005模型Cr的结果x=1.0e+003*0.47073.63905.60431.004resnorm=2.7728e+005x=1.0e+004*0.01801.05710.53530.0356resnorm=515.8698模型中Cu的结果x=1.0e+004*0.04482.45991.61160.2576resnorm=6.5586e+004x=1.0e+004*0.07252.74041.67030.1519resnorm=1.2779e+004模型Zn的结果x=1.0e+004*0.00221.37750.96180.3770resnorm=1.0274e+004x=1.0e+003*0.17969.27464.13562.9542resnorm=3.3278e+004x=1.0e+003*0.37233.75475.63752.8837resnorm=1.5491e+005x=1.0e+004*0.01091.27100.31330.1642resnorm=4.1042e+003x=1.0e+003*0.28924.81277.03121.6579resnorm=2.2450e+003模型Hg的结果x=1.0e+004*0.04320.25190.27291.6010resnorm=1.0177e+008x=1.0e+004*0.00071.36900.23631.4910resnorm=9.8607e+003x=1.0e+004*0.00071.52470.90991.3805resnorm=1.1949e+003模型中Ni的结果x=1.0e+003*0.01863.16365.88151.4501resnorm=1.0696e+003x=1.0e+004*0.01412.24631.20620.0152resnorm=392.5825模型中Pb的结果x=1.0e+003*0.24294.66595.14550.5274resnorm=4.7155e+003x=1.0e+003*0.32131.82103.10100.5269resnorm=5.7956e+004改进模型的结果x=1.0e+0041.57760.01000.0224r=100.*(-log((4000-224)./15776)).^0.5r=119.5753第二个x=1.0e+004*1.48220.01000.0079r=100.*(-log((3500-78.875)./14822.13)).^0.5r=121.0848第三个x=1.0e+004*1.37520.01000.0048r=100.*(-log((3500-47.92833)./13752.07)).^0.5r=117.56769.3问题三中处理过后的数据As：编号x(m)y(m)海拔(m)功能区As(μg/g)d距1距离6164727286214.080823833692729.621212.845401787426.561896.8637288336171548.941522.5041029331767748.721605.402213737311145.932015.711435264357744.492486.822139351825715957.411877.576142392721102755.722362.271编号x(m)y(m)海拔(m)功能区As(μg/g)d距离229474272939221.87030494872936218.3820626563579652946.141117.6011954816004047.411485.81427539486311245.311488.40525537586431514.091491.036编号x(m)y(m)海拔(m)功能区As(μg/g)d距离68697286416869728618416.58042705683483717.411078.33840801772103946.561150.51343774782604945.931311.3226563579652946.141408.47339807764012914.291497.49428529173491043.691579.257238642388314023.371608.086编号x(m)y(m)海拔(m)功能区As(μg/g)距离12696302417812696302427423.7201791240020601346.471008.4211761164935152716.051156.4131831385533457916.261202.6311751141525851245.621354.135511170244807145.311762.9441801259110631814.171963.8091901508735126953.372440.292185148961603421.612619.015编号x(m)y(m)海拔(m)功能区As(μg/g)距离181341004684181341004641430.13082187381092153310.991063.2228119007114888436.141685.6731991790482872515.831773.97325916289100724315.411845.1832581642890692014.381965.9519720101107744044.492097.39726016267110586044.382123.63710018993123717834.292478.61Cd编号x(m)y(m)海拔(m)功能区Cd(ng/g)距离9521439113834541619.800982107210404324386.601045.5291042253511293543233.201099.6891352219312185793205.501100.782932230410527404232.901216.947972010110774404180.601470.0771312468514278983109.301735.382101199681296142387.102157.296811900711488843100.302434.266251127341034432395.602638.003编号x(m)y(m)海拔(m)功能区Cd(ng/g)距离2232996018441213.5001647774897811044.501855.027304948729362630.002084.425195481600404345.702182.04515140263913134401.802227.006315567678272635.302393.2242852917349104323.902395.751728833617154269.802436.772编号x(m)y(m)海拔(m)功能区Cd(ng/g)距离9270822952241123.900616472728621092.901145.954823833692721066.201434.306143415322997351024.901445.006313211791284439.201475.73215040202990274332.601484.712132427397121193.701699.39321373731114146.202056.28831646783765405100.702458.007编号x(m)y(m)海拔(m)功能区Cd(ng/g)距离616472728621092.9009270822952241123.901145.954823833692721066.201212.845728833617154269.801522.504102933176774267.101605.40213935182571595222.401877.57640178742223.901896.86387183939183264194.002187.3717478154153.802503.026Cr：编号x(m)y(m)海拔(m)功能区Cr(μg/g)距离32996018223299601844920.840204592460361744.461916.787195481600404159.452182.04513242739712180.352224.9931514026391313484.942227.00631556767827264.032393.224285291734910435.142395.751编号x(m)y(m)海拔(m)功能区Cr(μg/g)距离10685552849106855528344306.0204810631647257449.27945.54325211730553254442.011045.0085111702448071447.241460.3417110225382119443.631767.8947512778579993439.502110.47217611649351527135.922231.925410700818450433.792656.042Cu：编号x(m)y(m)海拔(m)功能区Cu(μg/g)距离303250211629010442528.48031417765356185308.611212.845114260738807383123.271322.3742591628910072431151.641434.3062368307972614460.111657.797276126251625966146.011781.17220418012441420431.432053.80719316823420767537.112210.5792309475120009213.312296.161编号x(m)y(m)海拔(m)功能区Cu(μg/g)距离305268521611422531364.850114260738807383123.271676.44312139545615615136.971855.0272862059113549423130.551916.78731646783765405114.812084.425260162671105860471.032182.0452368307972614460.112227.0067415092693632534.662395.751Hg编号x(m)y(m)海拔(m)功能区Hg(ng/g)距离2708229592708229522416000.0008238336927213500.001434.306104510492127124800.001667.48463311423389565259.002070.17511315140263913134138.002086.875176145468413643751339.292184.3390314352643577477.002218.3254953164678376540542.002458.0073231555314206040428.002616.574287编号x(m)y(m)海拔(m)功能区Hg(ng/g)距离13694235718213694235733414901.000181137651353154162.001006.507327178126963024274160.001200.37202617912400206013496.001327.64641418012591106318158.001700.30732517712734401543257.001915.8716031871581023078141.002116.590655186153877298138.002348.75137189158013966115319.002651.099772编号x(m)y(m)海拔(m)功能区Hg(ng/g)距离15248910625715248910616413800.00025816428906920184.001180.5799426113797962118463.211539.6837347214298741836165.361936.967733260162671105860434.002201.968438250140651098725426.002222.08235776170441069193315.002395.379093Ni编号x(m)y(m)海拔(m)功能区Ni(μg/g)距离223299601844142.50023357362135131.50336.304920459246036132.801916.78729474272939227.801925.58419548160040418.102182.0451514026391313419.162227.006285291734910413.902395.7513164678376540511.502641.524编号x(m)y(m)海拔(m)功能区Ni(μg/g)距离135221931218579374.030102217661234867414.51457.053695214391138345415.411100.78294214181072135115.531555.148129229651353578315.531656.47896205541122843418.341897.938286205911354942318.602104.0298210721040432412.262104.424297226051430193310.202155.73797201011077440411.932523.368Pb编号x(m)y(m)海拔(m)功能区Pb(μg/g)距离164777489781472.480204592460361228.64347.3629155062433954121.72626.56921758684904164111.241091.0221514026391313496.761237.84419548160040489.341311.89423357362135190.901783.66818653456416170.821908.0321504020299027468.102051.75572883361715462.912285.9653470484600244110.202290.338315692456967540.182290.854编号x(m)y(m)海拔(m)功能区Pb(μg/g)距离61647272862434.80082383369272381.641212.84513242739712157.641467.46372883361715462.911522.50410293317677436.941605.4021403469230852537.641869.7824017874232.281896.8632137373111436.182015.711747815435.382503.026Zn编号x(m)y(m)海拔(m)功能区Zn(μg/g)距离137979621611379796211843760.820249142079980142117.87544.959624814000897014294.34681.91646214325866623385.611091.24224712641956011290.071157.60860128558945184102.861159.457250140651098725495.481392.0422992204617634171381.011480.6227214298741836162.242259.25编号x(m)y(m)海拔(m)功能区Zn(μg/g)距离36932843112412893.4703790905365204403.271080.53750106434472454403.981324.81916279063978224334.391460.473880495439184369.801705.3523373045230101297.142222.8673980776401291294.692435.79616483942035274172.652460.18948106316472574102.042523.436编号x(m)y(m)海拔(m)功能区Zn(μg/g)距离2232996018441818.470143526435774853.981676.442045924603611013.471916.787315567678272615.102393.224155062433954494.802434.5861758684904164389.802800.13515354383994102208.272944.80815040202990274259.293112.656编号x(m)y(m)海拔(m)功能区Zn(μg/g)距离1781269630242741635.920177127344015432237.45991.7283179124002060134186.531008.421176116493515271147.351156.413183138553345791142.551202.63164124424329654135.711329.489175114152585124110.611354.1355111702448071492.551762.94419015087351269537.142440.292编号x(m)y(m)海拔(m)功能区Zn(μg/g)距离3049487293621626.020315567678272615.10802.6718195481600404380.921394.8512215006884662293.161554.0834168697286184215.101921.01323864238831402208.372130.9783270046226114193.372316.3823156924569675168.052540.666基于C8051F单片机直流电动机反馈控制系统的设计与研究基于单片机的嵌入式Web服务器的研究MOTOROLA单片机MC68HC（8）05PV8/A内嵌EEPROM的工艺和制程方法及对良率的影响研究基于模糊控制的电阻钎焊单片机温度控制系统的研制基于MCS-51系列单片机的通用控制模块的研究基于单片机实现的供暖系统最佳启停自校正（STR）调节器单片机控制的二级倒立摆系统的研究基于增强型51系列单片机的TCP/IP协议栈的实现基于单片机的蓄电池自动监测系统基于32位嵌入式单片机系统的图像采集与处理技术的研究基于单片机的作物营养诊断专家系统的研究基于单片机的交流伺服电机运动控制系统研究与开发基于单片机的泵管内壁硬度测试仪的研制基于单片机的自动找平控制系统研究基于C8051F040单片机的嵌入式系统开发基于单片机的液压动力系统状态监测仪开发模糊Smith智能控制方法的研究及其单片机实现一种基于单片机的轴快流CO〈,2〉激光器的手持控制面板的研制基于双单片机冲床数控系统的研究基于CYGNAL单片机的在线间歇式浊度仪的研制基于单片机的喷油泵试验台控制器的研制基于单片机的软起动器的研究和设计基于单片机控制的高速快走丝电火花线切割机床短循环走丝方式研究基于单片机的机电产品控制系统开发基于PIC单片机的智能手机充电器基于单片机的实时内核设计及其应用研究基于单片机的远程抄表系统的设计与研究基于单片机的烟气二氧化硫浓度检测仪的研制基于微型光谱仪的单片机系统单片机系统软件构件开发的技术研究基于单片机的液体点滴速度自动检测仪的研制基于单片机系统的多功能温度测量仪的研制基于PIC单片机的电能采集终端的设计和应用基于单片机的光纤光栅解调仪的研制气压式线性摩擦焊机单片机控制系统的研制基于单片机的数字磁通门传感器基于单片机的旋转变压器-数字转换器的研究基于单片机的光纤Bragg光栅解调系统的研究单片机控制的便携式多功能乳腺治疗仪的研制基于C8051F020单片机的多生理信号检测仪基于单片机的电机运动控制系统设计Pico专用单片机核的可测性设计研究基于MCS-51单片机的热量计基于双单片机的智能遥测微型气象站MCS-51单片机构建机器人的实践研究基于单片机的轮轨力检测基于单片机的GPS定位仪的研究与实现基于单片机的电液伺服控制系统用于单片机系统的MMC卡文件系统研制基于单片机的时控和计数系统性能优化的研究基于单片机和CPLD的粗光栅位移测量系统研究单片机控制的后备式方波UPS提升高职学生单片机应用能力的探究基于单片机控制的自动低频减载装置研究基于单片机控制的水下焊接电源的研究基于单片机的多通道数据采集系统基于uPSD3234单片机的氚表面污染测量仪的研制基于单片机的红外测油仪的研究96系列单片机仿真器研究与设计基于单片机的单晶金刚石刀具刃磨设备的数控改造基于单片机的温度智能控制系统的设计与实现基于MSP430单片机的电梯门机控制器的研制基于单片机的气体测漏仪的研究基于三菱M16C/6N系列单片机的CAN/USB协议转换器基于单片机和DSP的变压器油色谱在线监测技术研究基于单片机的膛壁温度报警系统设计基于AVR单片机的低压无功补偿控制器的设计基于单片机船舶电力推进电机监测系统基于单片机网络的振动信号的采集系统基于单片机的大容量数据存储技术的应用研究基于单片机的叠图机研究与教学方法实践基于单片机嵌入式Web服务器技术的研究及实现基于AT89S52单片机的通用数据采集系统基于单片机的多道脉冲幅度分析仪研究机器人旋转电弧传感角焊缝跟踪单片机控制系