地形地貌图形的绘制

[x,y]=meshgrid(1:10);%构造测量网格
h=-3+8*rand(10,10);%随机生成10*10的-3到5范围的矩阵
[xi,yi]=meshgrid(1:0.1:10);%构造插值网格
hi=interp2(x,y,h,xi,yi,'spline');%spline插值方法为三次样条
s=surf(xi,yi,hi);%画出地貌图
s.EdgeColor='none';%隐藏边
xlabel('x'),ylabel('y'),zlabel('z');%添加坐标名称