2.mesh和surf函数

mesh函数用于绘制三维网格图，surf用于绘制三维曲面图，各线条之间的补面用颜色填充。以下介绍mesh的使用格式 (surf与mesh的使用格式相同）。

mesh(X,Y,Z,C) ： 其中，X,Y,Z 分别是由曲面上mxn个点的x坐标、y坐标和z坐标构成的mxn的矩阵，即： (X(i,j),Y(i,j),Z(i,j))是第(i,j)个网格点Pij的坐标， C(i,j)是Pij的颜色值。

mesh(X,Y,Z) ：使用C=Z（颜色值正比于高度值）

mesh(x,y,Z)或mesh(x,y,Z,C) ：使用两个向量代替两个矩阵，要求： length(x)=n, length(y)=m ,且[m,n]=size(Z)，第(i,j)个网格点Pij的坐标为(x(j),y(i),Z(i,j))

mesh(Z)：等价于mesh(1:n,1:m,Z)。

mesh(Z,C)：等价于mesh(1:n,1:m,Z,C)。

例3.绘制函数z=x*exp(-x^2-y^2)在[-2,2]x[-3,3]内的图形。下面两段代码效果是完全一样的。

%% figure1figure(1)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);mesh(X,Y,Z);title('figure 1')colorbar

%% figure2figure(2)x=-2:0.1:2;m=size(x,2);y=-3:0.1:3; n=size(y,2);Z=zeros(n,m);for i=1:m    for j=1:n        Z(j,i)=x(i)*exp(-x(i)^2-y(j)^2);    endendmesh(x,y,Z);colorbartitle('figure 2')

 结果展示：

 曲面图绘制（只需将mesh改为surf即可），这十张图片体现不同的颜色映射风格，标题为对应的colormap的值。

代码如下：

%% figure 1figure(1)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('hsv')colormap('hsv')shading flat%% figure 2figure(2)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('jet')colormap('jet')shading interp%% figure 3figure(3)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('copper')colormap('copper')shading interp%% figure 4figure(4)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('bone')colormap('bone')shading interp%% figure 5figure(5)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('flag')colormap('flag')shading interp%% figure 6figure(6)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('pink')colormap('pink')shading interp%% figure 7figure(7)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('spirng')colormap('spring')shading interp%% figure 8figure(8)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('summer')colormap('summer')shading interp%% figure 9figure(9)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('autumn')colormap('autumn')shading interp%% figure 10figure(10)x=-2:0.1:2;   y=-3:0.1:3; [X,Y]=meshgrid(x,y);Z=X.*exp(-X.^2-Y.^2);surf(X,Y,Z);title('winter')colormap('winter')shading interp

图像结果展示：

 例4.绘制三维曲面图z = sinx^2 + cosy^2，x:[0,pi]，y:[0,pi/2]。

代码：

[x,y]=meshgrid(0:pi/100:pi, 0:pi/100:pi/2);z=sin(x.^2)+cos(y.^2);surf(x,y,z);axis([0 4 0 1.8 -1.5 1.5]);colormap('prism')shading interpcolorbar

 运行结果：

这节先说这么多，下期将继续分享！

来源地址：https://blog.csdn.net/higerwy/article/details/129435292