数学模型选读实报告第二次实验.doc

数学模型选读实验报告 实验一 Matlab入门实验 姓 名： 戴安宇 学科专业： 电子信息与电气工程 学号： 201181359 完成日期： 2012.3.23 大连理工大学 Dalian University of Technology 目 录 引言 3 模型 4 算法设计 7 结果与讨论 8 结论 13 References 14 附录-source code 15 问题重述 用自己的语言叙述一下要解决的问题 城市按照功能划分，可分为生活区，山区，工业区，主干道路区及公园区，现在对某城市城区地质土壤环境进行调查，现在给出采样点的位置、海拔高度、及其所属功能区的信息，8种主要重金属元素在采样点处的浓度，以及8种主要重金属元素的背景值。现在要求我们给出8种主要重金属元素在该城区的空间分布，并分析改城区不同区域重金属的污染程度。根据问题（1）的结果，结合该城区的海拔高度信息及功能区信息等，说明重金属污染的主要原因。 相关函数命令 列举在解决问题的过程中用到的Matlab命令，并举例说明每个命令的用法。 clera all 直接输入即可，是清除屏幕中的其他数据的意思。 xlsread 是读入一个excel文件的意思。 Meshgrid是生成网格的意思 TriScatteredInterp某种插值方法 round四舍五入 subplot在同一画面里画多张表格 shading interp;消除网格线 surf画三维曲面 实现过程 按照实验讲义中例子的形式给出实验的代码实现过程。 clear all; x0=xlsread(shiyan2.xls,附件1,B4:B322); y0=xlsread(shiyan2.xls,附件1,C4:C322); h0=xlsread(shiyan2.xls,附件1,D4:D322); z0=xlsread(shiyan2.xls,附件1,E4:E322); as0=xlsread(shiyan2.xls,附件2,b4:b322); cd0=xlsread(shiyan2.xls,附件2,c4:c322); cr0=xlsread(shiyan2.xls,附件2,d4:d322); cu0=xlsread(shiyan2.xls,附件2,e4:e322); hg0=xlsread(shiyan2.xls,附件2,f4:f322); ni0=xlsread(shiyan2.xls,附件2,g4:g322); pb0=xlsread(shiyan2.xls,附件2,h4:h322); zn0=xlsread(shiyan2.xls,附件2,i4:i322); x=0:50:30000; y=0:50:20000; [xgrid,ygrid]=meshgrid(x,y); thigh=TriScatteredInterp(x0,y0,h0); tzone=TriScatteredInterp(x0,y0,z0); tas=TriScatteredInterp(x0,y0,as0); tcd=TriScatteredInterp(x0,y0,cd0); tcr=TriScatteredInterp(x0,y0,cr0); tcu=TriScatteredInterp(x0,y0,cu0); thg=TriScatteredInterp(x0,y0,hg0); tni=TriScatteredInterp(x0,y0,ni0); tpb=TriScatteredInterp(x0,y0,pb0); tzn=TriScatteredInterp(x0,y0,zn0); high=thigh(xgrid,ygrid); zone=round(tzone(xgrid,ygrid)); as=tas(xgrid,ygrid); cd=tcd(xgrid,ygrid); cr=tcr(xgrid,ygrid); cu=tcu(xgrid,ygrid); hg=thg(xgrid,ygrid); ni=tni(xgrid,ygrid); pb=tpb(xgrid,ygrid); zn=tzn(xgrid,ygrid); subplot(3,3,1); surf(xgrid,ygrid,high,zone); shading interp; subplot(3,3,2); surf(xgrid,ygrid,high,as); shading interp