《数学模型选读实报告第二次实验》由会员分享,可在线阅读,更多相关《数学模型选读实报告第二次实验(7页珍藏版)》请在金锄头文库上搜索。
1、分数:数学模型选读实验报告实验一Matlab入门实验 姓 名: 戴安宇 学科专业: 电子信息与电气工程 学号: 201181359 完成日期: 2012.3.23 大连理工大学Dalian University of Technology - I -目 录引言3模型4算法设计7结果与讨论8结论13References14附录-source code15数学建模实践(一)实验报告问题重述用自己的语言叙述一下要解决的问题城市按照功能划分,可分为生活区,山区,工业区,主干道路区及公园区,现在对某城市城区地质土壤环境进行调查,现在给出采样点的位置、海拔高度、及其所属功能区的信息,8种主要重金属元素在采
2、样点处的浓度,以及8种主要重金属元素的背景值。现在要求我们给出8种主要重金属元素在该城区的空间分布,并分析改城区不同区域重金属的污染程度。根据问题(1)的结果,结合该城区的海拔高度信息及功能区信息等,说明重金属污染的主要原因。相关函数命令列举在解决问题的过程中用到的Matlab命令,并举例说明每个命令的用法。clera all 直接输入即可,是清除屏幕中的其他数据的意思。xlsread 是读入一个excel文件的意思。Meshgrid是生成网格的意思TriScatteredInterp某种插值方法round四舍五入subplot在同一画面里画多张表格shading interp;消除网格线su
3、rf画三维曲面实现过程按照实验讲义中例子的形式给出实验的代码实现过程。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);
4、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,
5、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,ygri
6、d);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
7、interp;subplot(3,3,3);surf(xgrid,ygrid,high,cd);shading interp;subplot(3,3,4);surf(xgrid,ygrid,high,cr);shading interp;subplot(3,3,5);surf(xgrid,ygrid,high,cu);shading interp;subplot(3,3,6);surf(xgrid,ygrid,high,hg);shading interp;subplot(3,3,7);surf(xgrid,ygrid,high,ni);shading interp;subplot(3,3,8
8、);surf(xgrid,ygrid,high,pb);shading interp;subplot(3,3,9);surf(xgrid,ygrid,high,zn);shading interp;结果与讨论 用图像或表格的形式体现出实验结果,并对结果进行分析讨论。通过图表可以看出,第一个图为整体地形图,以下八个图分别为As,Cd,CrCu,Hg,Ni,Pb,Zn,通过他们在各个图中的亮点,可以知道污染程度,黄色越深,污染程度越重。As (g/g)Cd (ng/g)Cr (g/g)Cu (g/g)Hg (ng/g)Ni (g/g)Pb (g/g)Zn (g/g)6.27289.9669.024
9、9.493.0418.3469.11237.01生活区7.25393.1153.41127.54642.3619.8193.04277.93工业区4.04152.3238.9617.3240.9615.4536.5673.29山区5.16292.1351.8147.54314.1716.9154.72187.44交通区0280.5443.6430.19114.9915.2960.71154.24公园绿地区 工业区会产生很多废水废气,而这些废水废气中含有Hg,Cd,Zn,生活区因为有很多人生活,人口密度较大,所以Cd,Zn等含量较高,公园绿地区因为树木较多,环境较好,所以仅仅Cd含量较高,交通区
10、因为车辆较多,并且滞留时间较长,所以Cd,Hg,Zn等含量都较高,山区树木较多,并且人烟稀少,所以重金属污染现象不是很严重。总结与心得体会对实验过程及结论进行总结,并写出在实验过程中遇到的一些问题和解决办法等心得。每一次实验一开始给人的感觉都是很难,给人一种无从下手的感觉,但是其实每一道题都是和老师给的讲义种题的类型很相像,只要自己把讲义中的问题看懂,把程序看明白,那么你就发现他们并不是多么难。我们还要注意同学之间的合作,我们既不能一味的问别人,因为那样不利于我们的思维,但是我们又不能忽视同学之间的讨论与交流,因为那样会使我们的思路更加开阔。References 若有参考文献,请依次列出,并统一格式,并在文中引用该文献处注明。