exm044_4.m

来自「matlab的初级教程」· M 代码 · 共 16 行

%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