grid
axis square
6
grid off
4. 再看：下面表格是某城市的采样点的具体位置及海拔，取前50个样本点（共319个取样点）。
附件1. 取样点位置及其所属功能区 编号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18
x(m) 74 1373 1321 0 1049 1647 2883 2383 2708 2933 4233 4043 2427 3526 5062 4777 5868 6534 y(m) 781 731 1791 1787 2127 2728 3617 3692 2295 1767 895 1895 3971 4357 4339 4897 4904 5641 海拔(m) 5 11 28 4 12 6 15 7 22 7 6 14 2 7 5 8 16 6 功能区 4 4 4 2 4 2 4 2 4 4 5 1 1 4 4 1 4 1 7
编号 x(m) 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 5635 5394 5291 4742 4948 5567 7004 7304 7048 8180 9328 9090 8049 8077 8017 6869 7056 7747 y(m) 海拔(m) 功能区 7965 8631 7349 7293 7293 6782 6226 5230 4600 4496 4311 5365 5439 6401 7210 7286 8348 8260 29 12 10 9 6 7 11 10 24 15 24 20 18 29 39 18 37 49 4 4 4 2 2 2 4 1 4 4 1 4 4 1 4 4 1 4 19 20 21 22 23 24 25 1 2 3 4 5 5481 4592 2486 3299 3573 4741 5375 生活区 工业区 山区 交通区 6004 4603 5999 6018 6213 6434 8643 0 6 2 4 5 5 15 4 1 1 4 1 5 1 44 45 46 47 48 49 50 8457 9460 9062 9319 10631 10685 10643 8991 8311 7639 6799 6472 5528 4472 21 45 45 49 57 34 45 4 4 4 4 4 4 4 功能区 公园绿地区 下表是测量所得50个样本点的Cd浓度 表二为前50个样本点的Cd的浓度：
附件2. 8种主要重金属元素的浓度
编号 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25
Cd (ng/g) 153.80 146.20 439.20 223.90 525.20 1092.90 269.80 1066.20 1123.90 267.10 201.40 287.00 193.70 359.50 516.40 1044.50 445.40 347.90 345.70 614.00 257.20 1213.50 325.80 212.10 90.50 编号 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 Cd (ng/g) 583.40 366.40 323.90 424.50 630.00 635.30 463.40 532.00 778.70 754.80 396.30 687.80 526.00 449.10 852.70 459.00 337.30 568.10 599.00 635.50 600.70 567.60 228.50 568.60 214.70 8
根据数据可知，取样点之间的地理距离是无关的，在考虑重金属污染元素的分布式，只需把海拔换成浓度。
在编写程序中发现这不是一般的函数曲面，而是不规则采样点的集合。调用linespace函数在采样点范围内生成非均匀分布的数据。联合调用griddata函数将数据内插到图中。 程序清单
我们给出地形图和Cd元素的分析程序代码： 程序1：
A=xlsread('附件_数据.xls',1,'B4:E53');
%先将数据放在电子表格中，再读取把第1页中数据读取,必须会!!! x1=A(:,1);y1=A(:,2); z1=A(:,3);
[X1,Y1,Z1]=griddata(x1,y1,z1,linspace(min(x1),max(x1),200)',linspace(min(y1),max(y1),200),'v4'); %插值，v4是一种插值算法 surf(X1,Y1,Z1);
2shading interp; %在flat的基础上进行色彩的插值处理，使色彩平滑过渡
title('含量分布三维曲面');
程序2：
A=xlsread('附件_数据.xls',1,'B4:E53'); B=xlsread('附件_数据.xls',2,'H4:I53'); D=[A B];
x1=D(:,1);y1=D(:,2); z1=D(:,6);
[X1,Y1,Z1]=griddata(x1,y1,z1,linspace(min(x1),max(x1),200)',linspace(min(y1),max(y1),200),'v4');
surf(X1,Y1,Z1);shading interp; title('Cd含量分布三维曲面'); 实现图形分析
第一个图是城市地形图即地形地貌，海拔越高越趋向于红。在分析重金属污染源时，研究人员可以根据流体相关知识结合本图及当地降水量等信息，进行数据分析，以促进研究的准确性。
第二个图是Cd元素在这个城市分布的三位立体图，视觉越高颜色越深则此重金属元素含量越高。第三个图是和第二个图相对性的二维，颜色越趋近于红色则所含浓度越高。
9
6. 要在某山区方圆大约27平方公里范围内修建一条公路，从山脚出发经过一个居民区，再到达一个矿区。横向纵向分别每隔400米测量一次，得到一些点的高程，下表为数据表，请利用matlab绘制该山区地貌图、等高线图和立体等高线图。 提示：meshz(x,y,z)绘制地貌图
contour(x,y,z)绘制等高线图
10