exm04112_1.m

来自「同上」· M 代码 · 共 12 行

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)