数学建模第一题城市表层土壤重金属污染分析答案matlab.docx

PAGE PAGE 1 / 3 A 题 xx 表层土壤重金属污染分析 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。 按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等，分别记为 1 类区、2 类区、……、5 类区，不同的区域环境受人类活动影响的程度不同。 现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距 1 公里左右的网格子区域，按照每平方公里 1 个采样点对表层土（0~10 厘米深度）进行取样、编号，并用 GPS 记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元素的浓度数据。另一方面，按照2 公里的间距在那些远离人群及工业活动的自然区取样，将其作为该城区表层土壤中元素的背景值。 附件 1 列出了采样点的位置、海拔高度及其所属功能区等信息，附件2 列出了 8 种主要重金属元素在采样点处的浓度，附件 3 列出了 8 种主要重金属元素的背景值。 现要求你们通过数学建模来完成以下任务： 给出 8 种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度。 通过数据分析，说明重金属污染的主要原因。 分析重金属污染物的传播特征，由此建立模型，确定污染源的位置。(4)分析你所建立模型的优缺点，为更好地研究城市地质环境的演变模式， 还应收集什么信息？有了这些信息，如何建立模型解决问题？ 代码： rand(seed,0) x = rand(100,1)*4-2; y = rand(100,1)*4-2; z = x.*exp(-x.^2-y.^2); ti = -2:.25:2; [XI,YI] = meshgrid(ti,ti); ZI = griddata(x,y,z,XI,YI); mesh(XI,YI,ZI), hold plot3(x,y,z,o), hold off d = [-71560 -70 5 1 74 5 1 75 57 9 0.5 18 15 0.5 53 9 -40 11 -16 -40 45 -9 -70 32 -4 35 11 10 39 45 5 -29 18 -5 -52 37 -3 52 25 6 -18 53 -5 15 15 10 -57 9 3 65 41 2 -56 52 3 15 5 9 19 33 10 -19 6 -6 -45 23 -12 -63 17 -2 57 17 12]; [X Y] = meshgrid(-80:80, -80:80); Z = griddata(d(: ,1), d(: ,2), d(: ,3), X, Y); contour(X, Y, Z)