exm05112_1.m

来自「这是本人收集的一些程序源代码」· M 代码 · 共 14 行

%exm05112_1.m
randn('state',2)							%为仿真的可重复性而设
x=-5:5;y=-5:5;[X,Y]=meshgrid(x,y);		%产生"经纬"矩阵
%
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的随机误差
clf,shg
surf(X,Y,Z);view(-25,25)
disp('按任意键，继续！'),pause
%（2）
xi=linspace(-5,5,50);yi=linspace(-5,5,50);[XI,YI]=meshgrid(xi,yi);
ZI=interp2(X,Y,Z,XI,YI,'*cubic');
shg,surf(XI,YI,ZI),view(-25,25)