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')