重金属的空间分布的matlab程序.docx

重金属的空间分布的 matlab程 序 韩山师院数信系数学与应用数学专业 《数学实验》实验报告 班级 学号 姓名 成绩 试验 内容 Matlab三维作图程序设计 试验 类别 自选 试验 试验 时间 2012-6-4 试验问题：重金属的空间分布的matlab程序(以金属锌为例) 试验目的：① 掌握matlab画图程序设计方法； 用matlab作出重金属锌的空间分布图； 通过试验，认识matlab画图的妙处，掌握画图程序设计基本方法。 问题分析(可含问题的背景、相关知识、数学建模与求解的方法等 )： 问题背景： 随着城市经济的快速发展和城市人口的不断增加，人类活动对城市环境质量的影响日显突 出。对城市土壤地质环境异常的查证，以及如何应用查证获得的海量数据资料开展城市环 境质量评价，研究人类活动影响下城市地质环境的演变模式，日益成为人们关注的焦点。 按照功能划分，城区一般可分为生活区、工业区、山区、主干道路区及公园绿地区等， 分别记为1类区、2类区、……、5类区，不同的区域环境受人类活动影响的程度不同。 现对某城市城区土壤地质环境进行调查。为此，将所考察的城区划分为间距 1公里左 右的网格子区域，按照每平方公里1个采样点对表层土( 0~10厘米深度)进行取样、编号, 并用GPS记录采样点的位置。应用专门仪器测试分析，获得了每个样本所含的多种化学元 素的浓度数据。另一方面，按照 2公里的间距在那些远离人群及工业活动的自然区取样， 将其作为该城区表层土壤中元素的背景值。 数学建模与求解的方法： 1导入数据，由于数据过多，在此省去，详细见试验过程中 A的矩阵。并给xyz分别赋值： x=A(:,1);y=A(:,2);z=A(:,3); 2运用函数作图 scatter(x,y,5,z)% 运用scatter函数作散点图 figure 运用并进行分[X,Y,Z]=griddata(x,y,z,li nspace(1.486,4.271),li nspace(1.604,4.276),v4);% griddata函数进行插值 运用 并进行分 pcolor(X,Y,Z);shadi ng in terp% 运用函数 pcolor 作伪彩色图 figure,contourf(X,Y,Z) % 运用函数 contourf 作等高线图 figure,surf(X,Y,Z)% 运用函数surf作三维曲面 3与已用excel做出取样点位置及其所属功能区图进行对比，找出适用的图形, 析。 试验步骤（根据问题分析及试验目的所计划的试验步骤 ）： 1导入数据（见实验过程） 2运用函数作图 scatter（x,y,5,z）% 运用scatter函数作散点图 20 16 16 14 12 10 o O o O o V 0 o % ◎ cr° ° n o % 0 ° 0 0 n 6 0 Qe4 0 Q e oO°° ° oO°° ° 0 0 5 10 15 20 25 30 figure [X,Y,Z]=griddata(x,y,z,li nspace(1.486,4.271),li nspace(1.604,4.276),v4);% griddata函数进行插值 pcolor(X,Y,Z);shadi ng in terp% 运用函数 pcolor 作伪彩色图 15 2 2.5 3 3.5 4 figure,contourf(X,Y,Z) % 运用函数 contourf 作等高线图 15 2 2.5 3 3.5 figure,surf(X,Y,Z)% 运用函数surf作三维曲面 2000 1500 1000 500 0 3 3 2 2 500 5 1 1 3与已用excel做出取样点位置及其所属功能区图进行对比，找出适用的图形，并进行分 析。 取样点位置及其所属功能区*生活区■工业区 厶 取样点位置及其所属功能区 *生活区 ■工业区 厶山区 主干道路区 公园绿地区 x/km 对比可得，运用函数contourf作等高线图最为理想。具体分析见结果分析。 9.09 5.365 403.27 试验过程（含详细试验步骤、 8.049 5.439 369.80 8.077 6.401 294.69 8.017 7.21 377.14 程序清单及异常情况记录等） 6.869 7.286 215.10 8种重金属的空间分布的 7.056 8.348 210.00 matlab程序：（锌兀素为例） 7.747 8.26 572.96 A=[0.074 0.781 72.35 8.457 8.991 427.04 1.373 0.731 94.59 9.46 8.311 538.98 1.321 1.791 218.37 9.062 7.639 186.33 0 1.787 117