ex1605.m

来自「数字信号处理的MATLAB实现」· M 代码 · 共 10 行

%例16-5 二维插值
[X,Y] = meshgrid(-3:.25:3);	%产生已知数据的栅格点
Z = peaks(X,Y);		%计算已知数据点上的函数值，peaks是MATLAB的一个二元函数
[XI,YI] = meshgrid(-3:.125:3);		%产生更精细的插值栅格点
ZI = interp2(X,Y,Z,XI,YI);
mesh(X,Y,Z)	%绘制(X,Y,Z)已知数据的栅格图
hold on	%打开附加绘图模式，以后的绘图结果叠加在现有图像上
mesh(XI,YI,ZI+15)	%绘制(XI,YI,ZI+15)插值数据的栅格图（沿Z轴上移15个单位）
hold off	%关闭附加绘图模式
axis([-3 3 -3 3 -5 20])		%调整坐标轴范围，最终结果如图16-7