exm064_2.m

来自「北京航空航天大学出版社张志涌、徐彦琴主编的matlab教程附带的mfile例程压」· M 代码 · 共 19 行

%exm064_2.m

clear
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);						%计算对r的全导数	<4>
dz2=del2(Z);											%计算曲率		 	<5>
subplot(1,2,1),surf(X,Y,Z,abs(dzdr))
shading faceted;
brighten(0.6);
colorbar('horiz')
title('No. 1     surf(X,Y,Z,abs(dzdr))')
subplot(1,2,2),surf(X,Y,Z,abs(dz2))	
shading faceted;colorbar('horiz')
title('No. 2     surf(X,Y,Z,abs(dz2))')