randn('state',2)							%為仿真的可重復(fù)性而設(shè)
x=-5:5;y=-5:5;[X,Y]=meshgrid(x,y);	%產(chǎn)生“經(jīng)緯”矩陣
%以下兩條指令本可寫成一條。因Notebook無(wú)法正確處理續(xù)行號(hào)，不得不如此。
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;		%使數(shù)據(jù)帶標(biāo)準(zhǔn)差為0.05的隨機(jī)誤差
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)