%interp2_example
%采用二次插值對三維高斯型分布函數進行插值
[X,Y] = meshgrid(-3:0.8:3);                %原始數據
Z = peaks(X,Y);
[XI,YI] = meshgrid(-3:.125:3);             %插值點
ZI = interp2(X,Y,Z,XI,YI,'cubic');         %插值
hold on;
mesh(X,Y,Z,'Linestyle','-.');              %畫原始數據
mesh(XI,YI,ZI+30);                         %畫插值后的點
xlabel('x');
ylabel('y');
zlabel('z');
legend('原始數據','三次多項式插值');        %添加圖例
hold off
axis([-3 3 -3 3 -5 40]);                   %定義坐標軸范圍