Matlab绘制三维网线图 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 1

12345678910111213141516171819202122232425262728 %% ======四维表现========``%用色彩表现函数的特征``%当三维网线图、曲面图的第四个输入宗量取一些特殊矩阵时，色彩就能表现或加强函数的某特征，如梯度、曲率、方向导数等。``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 的全导数``dz2=del2(Z); % 计算曲率``figure;``subplot(``1``,``2``,``1``),surf(X,Y,Z),title(``'No. 1 surf(X,Y,Z)'``);`` shading faceted,colorbar( ``'horiz'``) ,brighten(``0.2``);``subplot(``1``,``2``,``2``),surf(X,Y,Z,R),title(``'No. 2 surf(X,Y,Z,R)'``);`` shading faceted;colorbar( ``'horiz'``);``%色彩分别表现函数的高度和半径特征``figure;``subplot(``1``,``2``,``1``),surf(X,Y,Z,dzdx) ;``shading faceted;brighten(``0.1`` );colorbar( ``'horiz'``);``title(``'No. 3 surf(X,Y,Z,dzdx)'``);``subplot(``1``,``2``,``2``),surf(X,Y,Z,dzdy);`` shading faceted;colorbar( ``'horiz'``);``title(``'No. 4 surf(X,Y,Z,dzdy)'``);``%色彩分别表现函数的 x 方向和 y 方向导数特征``figure;``subplot(``1``,``2``,``1``),surf(X,Y,Z,abs(dzdr)) ;``shading faceted;brighten(``0.6`` );colorbar( ``'horiz'``);``title(``'No. 5 surf(X,Y,Z,abs(dzdr))'``);``subplot(``1``,``2``,``2``),surf(X,Y,Z,abs(dz2));`` shading faceted;colorbar( ``'horiz'``);``title(``'No. 6 surf(X,Y,Z,abs(dz2))'``);