randn('state',2)							%為仿真的可重復性而設
x=-5:5;y=-5:5;[X,Y]=meshgrid(x,y);	%產生“經緯”矩陣
%以下兩條指令本可寫成一條。因Notebook無法正確處理續行號，不得不如此。
zz=1.2*exp(-((X-1).^2+(Y-2).^2))-0.7*exp(-((X+2).^2+(Y+1).^2));
Z=-500+zz+randn(size(X))*0.05;		%使數據帶標準差為0.05的隨機誤差
surf(X,Y,Z);view(-25,25) 
 

xi=linspace(-5,5,50);yi=linspace(-5,5,50);[XI,YI]=meshgrid(xi,yi);
ZI=interp2(X,Y,Z,XI,YI,'*cubic');
surf(XI,YI,ZI),view(-25,25)