重金属地空间分布地matlab程序.docx
- 文档编号:18031628
- 上传时间:2023-08-07
- 格式:DOCX
- 页数:20
- 大小:443.10KB
重金属地空间分布地matlab程序.docx
《重金属地空间分布地matlab程序.docx》由会员分享,可在线阅读,更多相关《重金属地空间分布地matlab程序.docx(20页珍藏版)》请在冰点文库上搜索。
重金属地空间分布地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,linspace(1.486,4.271)',linspace(1.604,4.276),'v4');%运用griddata函数进行插值
pcolor(X,Y,Z);shadinginterp%运用函数pcolor作伪彩色图
figure,contourf(X,Y,Z)%运用函数contourf作等高线图
figure,surf(X,Y,Z)%运用函数surf作三维曲面
3与已用excel做出取样点位置及其所属功能区图进行对比,找出适用的图形,并进行分析。
试验步骤(根据问题分析及试验目的所计划的试验步骤):
1导入数据(见实验过程)
2运用函数作图
scatter(x,y,5,z)%运用scatter函数作散点图
figure
[X,Y,Z]=griddata(x,y,z,linspace(1.486,4.271)',linspace(1.604,4.276),'v4');%运用griddata函数进行插值
pcolor(X,Y,Z);shadinginterp%运用函数pcolor作伪彩色图
figure,contourf(X,Y,Z)%运用函数contourf作等高线图
figure,surf(X,Y,Z)%运用函数surf作三维曲面
3与已用excel做出取样点位置及其所属功能区图进行对比,找出适用的图形,并进行分析。
对比可得,运用函数contourf作等高线图最为理想。
具体分析见结果分析。
试验过程(含详细试验步骤、
程序清单及异常情况记录等)
8种重金属的空间分布的
matlab程序:
(锌元素为例)
A=[0.0740.78172.35
1.3730.73194.59
1.3211.791218.37
01.787117.35
1.0492.127726.02
1.6472.728966.73
2.8833.617166.73
2.3833.6921417.86
2.7082.295926.84
2.9331.767100.41
4.2330.895102.65
4.0431.895223.16
2.4273.97189.08
3.5264.357853.98
5.0624.339494.80
4.7774.897602.04
5.8684.904389.80
6.5345.641307.24
5.4816.004380.92
4.5924.6031013.47
2.4865.999223.27
3.2996.0181818.47
3.5736.213429.29
4.7416.434186.22
5.3758.64346.84
5.6357.965311.02
5.3948.631182.65
5.2917.349253.16
4.7427.293175.71
4.9487.2931626.02
5.5676.782615.10
7.0046.226193.37
7.3045.23297.14
7.0484.6351.63
8.184.496323.37
9.3284.3112893.47
9.095.365403.27
8.0495.439369.80
8.0776.401294.69
8.0177.21377.14
6.8697.286215.10
7.0568.348210.00
7.7478.26572.96
8.4578.991427.04
9.468.311538.98
9.0627.639186.33
9.3196.799208.06
10.6316.472102.04
10.6855.528196.73
10.6434.472403.98
11.7024.4892.55
11.735.532275.82
11.4826.354278.37
10.78.184295.61
10.638.774196.33
11.6788.618242.04
11.9027.709157.35
13.2447.056141.02
12.7468.45146.22
12.8558.945102.86
13.7979.6213760.82
14.3258.66685.61
15.4678.658122.96
12.4424.329135.71
13.0934.339155.00
13.925.354218.27
14.8445.519125.92
16.5696.05582.96
16.3876.609221.22
16.0617.35278.98
15.6587.594138.06
14.2987.41862.24
14.1776.684118.16
15.0926.93682.86
12.7785.799148.88
17.04410.69132.86
17.08711.93384.49
17.07512.92447.76
17.96212.823178.98
18.41311.72189.08
19.00711.48861.94
18.73810.92163.27
17.81410.70747.76
18.13410.046182.04
17.1989.81105.00
17.1449.08184.08
18.3939.18394.69
19.7678.8140.92
21.0068.81994.80
21.0919.482163.27
22.8469.14991.73
23.6649.7963.88
22.30410.52795.41
21.41810.72182.76
21.43911.38357.55
20.55411.228104.90
20.10110.77463.98
21.07210.404157.24
20.2159.951170.71
18.99312.37167.24
19.96812.96148.78
21.76612.34889.49
22.67412.17389.08
22.53511.293100.10
25.2215.79571.43
26.4535.57743.37
26.4166.508114.29
27.8165.58155.71
25.3616.42361.53
24.0657.35358.47
25.9987.03243.88
27.1777.77164.08
26.4248.63982.96
26.0738.80752.45
24.6319.42261.94
24.7029.52278.37
25.4619.83468.16
24.81310.79957.86
26.08611.09469.39
26.01512.07892.24
27.711.60974.90
27.69611.62163.88
27.34613.33179.18
26.59113.71557.96
27.82314.73796.53
27.23214.48258.06
24.5813.31969.90
24.15312.4587.65
22.96513.53599.49
23.19813.523142.86
24.68514.27861.43
28.6548.75553.06
24.00315.28651.53
21.68413.10147.55
22.19312.185229.80
17.0795.89434.39
15.2555.1154.29
15.0075.53576.53
3.5182.571133.98
3.4692.30876.12
3.7622.17127.35
3.9272.11154.49
4.1532.2991389.39
3.2670.793118.98
4.6841.364110.92
5.4951.20549.90
5.6641.65384.59
5.5412.09365.61
5.4512.75763.57
4.022.99259.29
4.0263.913227.55
5.1014.08184.69
5.4383.994208.27
5.3823.012256.94
5.3142.06144.69
5.5031.127118.88
5.6360.133109.29
6.6050.37487.14
7.0931.381213.47
7.12.449138.88
6.8373.49399.90
7.9063.978334.39
8.0453.052540.00
8.3942.035172.65
8.4031.075171.84
8.0790171.12
9.6631.288165.10
9.4692.286103.88
9.1783.299115.41
9.0953.975230.41
10.2253.821104.29
10.212.78964.49
10.341.764185.41
11.5571.581315.92
11.4152.585110.61
11.6493.515147.35
12.7344.015237.45
12.6963.0241635.92
12.42.06186.53
12.5911.063177.76
13.7651.353221.94
13.6942.357349.80
13.8553.345142.55
14.8622.52484.80
14.8961.603106.53
15.3870.72969.80
15.812.30780.00
16.0323.06194.08
15.8013.966104.18
15.0873.51237.14
16.8722.798101.73
17.7343.62991.02
16.8234.207115.20
17.0084.77559.49
17.2036.21881.73
17.0057.21270.51
16.9477.487104.90
16.3018.29990.20
17.9048.28753.57
18.3037.38587.35
18.4386.53947.24
18.5565.58895.92
18.9544.87487.04
18.0124.41471.53
19.0728.51960.00
20.2828.5941.12
21.4758.5453.27
21.457.555106.43
20.2617.58677.45
19.5697.34841.63
19.4116.93461.12
19.5016.09159.90
20.5826.54851.12
19.9095.3101.84
21.0185.76488.16
22.1765.492108.47
23.3595.325130.61
23.2386.50282.14
22.6244.818105.31
21.7036.59169.80
5.0068.846293.16
5.7349.659358.27
6.39510.443457.96
7.40510.981303.06
8.44611.2241.63
7.61211.938192.65
7.91212.84170.10
8.86613.143170.00
9.29613.10251.63
9.47512145.51
9.21211.305177.86
8.62912.086178.57
7.77610.613100.41
8.62210.638204.08
9.2379.872247.04
8.3079.726254.80
7.1069.46756.33
6.4238.831208.37
7.4588.92230.11
8.9048.868447.83
10.5479.591153.27
10.39810.36204.81
10.39511.203211.47
11.52911.243184.61
11.56310.298118.70
11.6469.381209.50
12.6419.5690.07
148.9794.34
14.2079.98117.87
14.06510.98795.48
12.73410.34457.68
12.7277.691531.45
14.17311.941388.69
15.46712.0857.27
15.1411.10192.46
15.19810.181.01
15.2489.106269.89
16.4289.06995.17
16.28910.07279.24
16.26711.05881.42
16.4412.068193.67
16.4413.23284.86
15.41212.98297.67
14.26912.877162.85
13.27713.204118.91
13.17512.238188.88
12.15312.33698.81
11.95813.31370.18
10.813.282321.12
10.02212.204132.86
9.33314.631175.34
9.27716.148168.05
11.12116.432125.16
10.85614.727197.63
12.64414.943218.24
12.62516.259149.52
9.03617.53875.91
10.59917.9895.69
12.63217.949155.98
14.40518.032172.53
14.07416.516101.00
14.26215.12953.10
14.62414.00494.96
16.62914.48183.19
18.4714.411109.33
20.59113.549121.62
20.98315.86247.38
20.17717.64292.36
19.04115.76963.83
18.90616.34644.67
18.46717.00168.20
17.41415.47665.70
15.74815.728141.50
15.51717.03483.40
16.60717.36569.24
15.95218.39784.24
22.60514.30180.90
23.14615.38271.01
22.04617.63481.01
23.78517.64361.64
25.98118.05150.92
27.3818.20238.21
25.02116.2957.16
23.32516.70147.58
26.85216.11494.02
17.98118.449103.50
14.48212.692147.75
14.31813.56946.86
10.35217.133112.35
9.09516.414435.44
10.5115.314117.66
13.9545.61561.02
10.1421.662121.20
17.7653.56163.31
6.9245.696168.05
4.6783.76556.23
6.1822.00567.06
5.9852.56760.50
7.6531.95295.90
];
x=A(:
1);y=A(:
2);z=A(:
3);
scatter(x,y,5,z)
figure
[X,Y,Z]=griddata(x,y,z,linspace(1.486,4.271)',linspace(1.604,4.276),'v4');
pcolor(X,Y,Z);shadinginterp
figure,contourf(X,Y,Z)
figure,surf(X,Y,Z)
xlabel('x轴');ylabel('y轴');zlabel('浓度')
结果分析:
Zn元素Pb元素
Ni元素Hg元素
Cu元素Cr元素
Cd元素As元素
图一(a)
图一(b)
根据图一(b)取样点位置及其所属功能区,由图一(a)可以看出,各个功能区表层土壤中8种重金属含量有明显差异,工业区和主干道路区表层土壤重金属含量较大,公园绿地区和生活区含量较小,而山区的重金属含量是最少的。
对不同种类重金属而言,As在各功能区表层土壤中含量差异不大,而Cu、Pb、Cr、Zn、Cd、Hg和Ni含量差异较大。
总结体会:
Matlab是一个很好有的工具,特别是在作图上,只要有足够的三维数据,就可以做出图形,可以是二维的彩色图,也可以是三维的彩色图,日常生活中可以运用其做等高线图,等温线图,或其他分布图,可以形象直观的表现出数据所隐含的东西,有利于数学建模中模型的选择,更利于外行人一目了然的了解数据所隐含的信息。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 重金 属地 空间 分布 matlab 程序