ex1607.m

数字信号处理的MATLAB实现

%例16-7 griddata在二维插值中的应用
x = rand(100,1)*4-2; 	%产生100个[-2 2]的均匀分布数据
 y = rand(100,1)*4-2;
z = x.*exp(-x.^2-y.^2);
ti = -2:.25:2; 
[XI,YI] = meshgrid(ti,ti);	%产生精细的插值栅格点
ZI = griddata(x,y,z,XI,YI,'cubic');	%基于非栅格数据点和cubic方法估计插值点函数值
mesh(XI,YI,ZI)
hold on	%打开附加绘图模式
plot3(x,y,z,'o')	%画出原始非栅格的离散数据点
hold off	%关闭附加绘图模式