exm044_4.m

北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压缩包

%exm044_4.m
clear,
dd=0.2;
x=-1:dd:1;y=x;
[X,Y]=meshgrid(x,y);
Z=(X.^2)+(Y.^2);
[DZx,DZy]=gradient(Z,dd,dd);
DZ2=4*del2(Z,dd,dd);
%与理论计算比较
DDZx0=DZx-2*X;DDZy0=DZy-2*Y;DDZ20=DZ2-4;
subplot(1,3,1),stem3(X,Y,DDZx0)
subplot(1,3,2),stem3(X,Y,DDZy0)
subplot(1,3,3),stem3(X,Y,DDZ20)
axis([-1,1,-1,1,-0.4,0.4])
xlabel('x'),ylabel('y')  
shg