www.gusucode.com > matlab通信工程仿真源码（张德丰等编著）程序书籍 > matlab_code/matlab通信工程仿真源码（张德丰等编著）/第4章/li4_43.m

    figure(close)
x=3*pi*(-1:1/30: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);
surf(X,Y,Z,abs(dzdr))
shading interp
colormap(spring)
alphamap('rampup')
alpha('color')
alpha('interp')