ex1606.m

来自「MATLAB7.0课本所有程序清单,共37章,特别适合初学者,可以帮助节省大量时」· M 代码 · 共 15 行

%例16-6 二维插值方法效果比较
[x,y] = meshgrid(-3:1:3);	%产生已知数据栅格点
z = peaks(x,y);	%计算已知点上的函数值
surf(x,y,z)	%画基于已知数据点的三维表面图，如图16-8
title('graphic based on original data')		%对图形加标题
[xi,yi] = meshgrid(-3:0.25:3);	%产生更精细的插值栅格点
zi1 = interp2(x,y,z,xi,yi,'nearest');
surf(xi,yi,zi1)	%画基于最邻近法插值的三维表面图，如图16-9
title('nearest-interpolation method')
zi2 = interp2(x,y,z,xi,yi,'linear');
surf(xi,yi,zi2)		%画基于二维分段线性插值的三维表面图，如图16-10
title('linear-interpolation method')
zi3 = interp2(x,y,z,xi,yi,'cubic');
surf(xi,yi,zi3)		%画基于二维三次多项式插值的三维表面图，如图16-11
title('cubic-interpolation method')