sea.m

假设有一组海底深度测量数据

randn('state',2)
x=-5:5;y=-5:5;[X,Y]=meshgrid(x,y);
zz=1.2*exp(-((X-1).^2+(Y-2).^2))-0.7*exp(-((X+2).^2+(Y+1).^2));
Z=-500+zz+randn(size(X))*0.05;
surf(X,Y,Z);view(-25,25)
xi=linspace(-5,5,50);yi=linspace(-5,5,50);[XI,YI]=meshgrid(xi,yi);
ZI=interp2(X,Y,Z,XI,YI,'*cubic');
surf(XI,YI,ZI),view(-25,25)