3-38.m

《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计算

axis([-3 3 -3 3 -5 20])
% 确定坐标轴的范围
[X,Y] = meshgrid(-3:.5:3);
% 生成网格
Z = peaks(X,Y);
% 计算peaks函数值
figure(1)
mesh(X,Y,Z)
% 绘制原始数据曲面图，如图3-9
 

[XI,YI] = meshgrid(-3:.125:3);
% 确定插值点
ZI1 = interp2(X,Y,Z,XI,YI,'*nearest');
% 计算用nearest方法所得二维插值
figure(2)
mesh(XI,YI,ZI1)
% 绘制最近点插值法曲面图，如图3-10
  
ZI2 = interp2(X,Y,Z,XI,YI,'*linear');
% 计算用linear法所得二维插值
figure(3)
mesh(XI,YI,ZI2)
% 绘制双线性插值法曲面图，如图3-11
ZI3 = interp2(X,Y,Z,XI,YI,'*spline');
% 计算用spline法所得二维插值
figure(4)
mesh(XI,YI,ZI3)
% 绘制双样条插值法曲面图，如图3-12
ZI4 = interp2(X,Y,Z,XI,YI,'*cubic');
% 计算用cubic法所得二维插值
figure(5)
mesh(XI,YI,ZI4)
% 绘制双立方插值法曲面图，如图3-13