原始代码片段
清除环境变量
清除命令窗口显示
定义变量x，存储一系列数值数据
x=[ 60 55 59 51 44 39 37 36 29 44 37 0 12 14 19 19 17 1 23 23 27 14 15 18 24 29 36 39 37 28 21 18 15 65 61 58 55 54 52 52 52 46 51 52 ];%这些数值代表x轴的坐标信息
定义变量y，存储一系列数值数据
y=[ 0 17 20 18 19 19 21 12 12 8 1 0 16 25 20 13 5 0 4 11 8 12 38 27 24 24 20 27 33 37 18 15 20 59 47 40 43 46 42 48 52 48 39 36 ];%这些数值代表y轴的坐标信息
定义变量z，存储一系列数值数据
z=[ 18 17 18 18 18 18 19 17 18 17 17 17 17 16 17 17 18 17 17 18 17 17 16 17 17 18 19 18 17 17 17 17 17 18 17 18 18 18 18 18 17 17 19 19 ];%这些数值代表z轴的坐标信息，即高程数据
使用meshgrid函数创建二维网格矩阵X和Y，通过设定起始值、步长和终止值来加密网格
[X,Y]=meshgrid(min(x):0.2:max(x),min(y):0.2:max(y)); %加密网格
利用griddata函数，根据原始数据点的坐标(x,y)和高程值(z)，对加密后的网格点进行插值计算，得到对应的高程值Z
Z=griddata(x,y,z,X,Y,’v4′);%在加密网格处取Z值
创建图形窗口1，用于绘制三维曲面图
figure(1);
使用surf函数绘制三维曲面图，将加密后的网格矩阵X、Y和高程矩阵Z作为输入参数
surf(X,Y,Z);%画三维图
设置x轴标签，指定字体大小、粗细和颜色
xlabel(‘X’,’fontsize’,12,’fontweight’,’b’,’color’,’r’);
设置y轴标签，指定字体大小、粗细和颜色
ylabel(‘Y’,’fontsize’,12,’fontweight’,’b’,’color’,’r’);
设置z轴标签，指定字体大小、粗细和颜色，并标注为”高程”
zlabel(‘高程’,’fontsize’,12,’fontweight’,’b’,’color’,’r’);
保持当前图形窗口，以便后续添加新的图形元素
hold on;
使用plot3函数，根据原始数据点的坐标(x,y)和高程值(z)，在三维图中用星号(*)标记原始数据点
plot3(x,y,z,’*’);%表示出原始数据点
释放图形窗口的保持状态
hold off;
创建图形窗口2，用于绘制等高线图
figure(2);
使用contour函数，根据加密后的网格矩阵X、Y和高程矩阵Z，绘制等高线图，并将返回值赋给变量C和h
[C,h] = contour(X,Y,Z);%画等高线图
设置等高线标签的显示属性，包括开启文本显示、调整文本步长等
set(h,’ShowText’,’on’,’TextStep’,get(h,’LevelStep’)*2)
应用”cool”色彩映射表，为等高线图设置颜色方案
colormap cool