char4_5.m

matlab 文件

%在500米间距正方形网格系统上测得的海底深度（米）
%绘制海底深度图
x=0:0.5:4;  %0.5km为单位
y=0:0.5:6;  %0.5km为单位
z=[100  99 100  99 100  99  99  99 100
   100  99  99  99 100  99 100  99  99
    99  99  98  98 100  99 100 100 100
   100  98  97  97  99 100 100 100  99
   101 100  98  98 100 102 103 100 100
   102 103 101 100 102 106 104 101 100
    99 102 100 100 103 108 106 101  99
    97  99 100 100 102 105 103 101 100
   100 102 103 101 102 103 102 100  99
   100 102 103 102 101 101 100  99  99
   100 100 101 101 100 100 100  99  99
   100 100 100 100 100  99  99  99  99
   100 100 100  99  99 100  99 100  99];
mesh(x,y,z)
xlabel('X-axis(km)')
ylabel('Y-axis(km)')
zlabel('海底深度(m)')
title('海底深度图')
pause
xi=linspace(0,4,30);
yi=linspace(0,6,40);
[xxi,yyi]=meshgrid(xi,yi);
zzi=interp2(x,y,z,xxi,yyi,'linear');
mesh(xxi,yyi,zzi)
title('海底深度图(linear)')
hold on
[xx,yy]=meshgrid(x,y);
plot3(xx,yy,z+0.1,'ob')
hold off
pause
zzi=interp2(x,y,z,xxi,yyi,'bilinear');
mesh(xxi,yyi,zzi)
title('海底深度图(bilinear)')
hold on
[xx,yy]=meshgrid(x,y);
plot3(xx,yy,z+0.1,'ob')
hold off
pause
zzi=interp2(x,y,z,xxi,yyi,'cubic');
mesh(xxi,yyi,zzi)
title('海底深度图(cubic)')
hold on
[xx,yy]=meshgrid(x,y);
plot3(xx,yy,z+0.1,'ob')
hold off
pause
zzi=interp2(x,y,z,xxi,yyi,'bicubic');
mesh(xxi,yyi,zzi)
title('海底深度图(bicubic)')
hold on
[xx,yy]=meshgrid(x,y);
plot3(xx,yy,z+0.1,'ob')
hold off
pause
zzi=interp2(x,y,z,xxi,yyi,'nearest');
mesh(xxi,yyi,zzi)
title('海底深度图(nearest)')
hold on
[xx,yy]=meshgrid(x,y);
plot3(xx,yy,z+0.1,'ob')
hold off
pause
zzi=interp2(x,y,z,xxi,yyi,'cubic');
pcolor(xxi,yyi,zzi)
shading interp
hold on
contour(xxi,yyi,zzi,15,'k')
colormap(cool)
colorbar('vert')
hold off
title('海底深度等值线图')