exm050402.m

来自「matlab电子书籍」· M 代码 · 共 17 行

clf
x=3*pi*(-1:1/15:1);y=x;[X,Y]=meshgrid(x,y);
R=sqrt(X.^2+Y.^2)+eps;Z=sin(R)./R;
[dzdx,dzdy]=gradient(Z);
dzdr=sqrt(dzdx.^2+dzdy.^2);						
dz2=del2(Z);										 			
subplot(1,2,1),surf(X,Y,Z,abs(dzdr))
shading faceted;
colorbar('SouthOutside')				
brighten(0.6);							
colormap hsv
title('No. 1     surf(X,Y,Z,abs(dzdr))')
subplot(1,2,2);surf(X,Y,Z,abs(dz2))	
shading faceted	
colorbar('NorthOutside')				
title('No. 2     surf(X,Y,Z,abs(dz2))')