Matlab图形绘制经典案例 （2）

24、绘制函数的梯度场矢量图。

>> [x,y]=meshgrid([-2:0.1:2]); %建立栅格点数据向量

>> z=3.*x.*y*exp(-x.^2-y.^2)-1; %计算函数值向量

>> [u,v]=gradient(z,0.2,0.2); %计算梯度值向量

>> quiver(x,y,u,v,2) %绘制梯度场矢量图

 

25、给定向量x,y生成网格矩阵。

>> x=[1 2 3 4];

>> y=[10 11 12 13 14];

>> [U,V]=meshgrid(x,y)

U =

     1     2     3     4

     1     2     3     4

     1     2     3     4

     1     2     3     4

     1     2     3     4

V =

    10    10    10    10

    11    11    11    11

    12    12    12    12

    13    13    13    13

    14    14    14    14

26、生成一个5阶高斯分布矩阵，并给出相应的x,y向量矩阵。

>> [X,Y,Z]=peaks(5)

X =

   -3.0000   -1.5000         0    1.5000    3.0000

   -3.0000   -1.5000         0    1.5000    3.0000

   -3.0000   -1.5000         0    1.5000    3.0000

   -3.0000   -1.5000         0    1.5000    3.0000

   -3.0000   -1.5000         0    1.5000    3.0000

Y =

   -3.0000   -3.0000   -3.0000   -3.0000   -3.0000

   -1.5000   -1.5000   -1.5000   -1.5000   -1.5000

         0         0         0         0         0

    1.5000    1.5000    1.5000    1.5000    1.5000

    3.0000    3.0000    3.0000    3.0000    3.0000

Z =

    0.0001    0.0042   -0.2450   -0.0298   -0.0000

   -0.0005    0.3265   -5.6803   -0.4405    0.0036

   -0.0365   -2.7736    0.9810    3.2695    0.0331

   -0.0031    0.4784    7.9966    1.1853    0.0044

    0.0000    0.0312    0.2999    0.0320    0.0000

 27、在-4<=x<=4,-4<=y<=4区域上绘制z=x^2+y^2的三维网格图。

>> [x,y]=meshgrid(-4:0.125:4);

>> z=x.^2+y.^2;

>> meshc(x,y,z)  

 28、绘制高斯分布函数的网格图。

>> [x,y]=meshgrid(-3:0.125:3);

>> z=peaks(x,y);

>> meshz(x,y,z)

 
 

29、用surf绘制高斯分布函数的曲面图。

>> [x,y]=meshgrid(-3:0.125:3);

>>z=peaks(x,y);

>>surf(x,y,z)

30、绘制曲线图。

>> t=[0:pi/200:10*pi];

>> x=2*cos(t);

>> y=3*sin(t);

>> z=t.^2;

>> plot3(x,y,z)

 

31、利用peaks函数产生的数据绘制其带形图。

>> [x,y]=meshgrid([-2*pi:pi/5:2*pi],[-2:1/5:2]);

>> z=peaks(x,y);

>> ribbon(y,z)

32、绘制三维饼图。

>> A=[1 2 3;4 5 6;7 8 9];

>> ex=[1 0 0;4 0 0;0 8 0];

>> pie3(A,ex)

33、在各种style参数的条件下绘制矩阵的三维条形图。

>> z=[1 2 3;4 5 6;7 8 9];

>>bar3(z,'detached')

>>title('bar3函数以detached参数绘制A=[1 2 3;4 5 6;7 8 9]的条形图')

>> bar3(z,'grouped')

>> title('bar3函数以grouped参数绘制A=[1 2 3;4 5 6;7 8 9]的条形图')

>> bar3(z,'stacked')

>> title('bar3函数以stacked参数绘制A=[1 2 3;4 5 6;7 8 9]的条形图')

34、绘制柱形图。

>> t=[0:pi/50:2*pi];

>> [x,y,z]=cylinder(t.*sin(t));

>> surf(x,y,z)

>> cylinder(t.^2)

>> title('cylinder(t^2)绘制的柱形图')

35、绘制三维散点图。

>> x=rand(500,1);

>> y=randn(500,1);

>> z=randn(500,1);

>> scatter3(x,y,z,'p','r')

36、绘制三维火柴杆图。

>> x=[1:0.5:20];

>> y=sqrt(x);

>> z=sqrt(x.^2+y.^2);

>> stem3(x,y,z,'filled')

37、绘制高斯分布函数的三维瀑布图。

>> [x,y]=meshgrid(-4:0.05:4);

>> z=peaks(x,y);

>> waterfall(x,y,z)

38、绘制等值线图。

>> [x,y]=meshgrid(-3:0.1:3);

>> z=2-x.^2-y.^2;

>> contour3(z,20)

39、绘制一个球面。

>> [x,y,z]=sphere(40);

>> surf(x,y,z)

40、绘制三角形网格图和三角形表面图。

>> [x,y]=meshgrid(-3:0.5:3);

>> z=x.*exp(-x.^2-y.^2);

>> tri=delaunay(x,y); 

%建立三角形网格

>> trimesh(tri,x,y,z)

>> trisurf(tri,x,y,z)

 

41、绘制一个三维彗星图。

>> t=[-3*pi:pi/100:3*pi];

>> x=3.*cos(t);

>> y=2.*sin(t);

>> z=t.^2;

>> comet3(x,y,z)

42、绘制曲面z的表面法向量向量图。

>> [x,y]=meshgrid([-3:0.2:3],[-2:0.5:2]);

>> z=x.*exp(-x.^2-y.*2);

>> [u,v,w]=surfnorm(x,y,z); %计算表面法向向量

>> quiver3(x,y,z,u,v,w,1.2) %绘制三维向量图

>> hold on

>> surf(x,y,z)

>> hold off

43、绘制空间立体在-2<=x<=2, -2<=y<=2, -2<=z<=2上的切片图。

>> [x,y,z]=meshgrid(-2:0.2:2);

>> v=x.*exp(-x.^2-y.^2-z.^2);

>> xi=[-1.2 0.8 2];yi=2;zi=[-2 -0.2];

>> slice(x,y,z,v,xi,yi,zi)

44、在【-pi,pi】上制作一个不断绘制正弦曲线的动画。

>> x=[-pi:0.02:pi];

>> y=sin(x);

>> h=plot(x,y,'r-')

h =

  171.0011

>> axis([-4 4 -1 1])

>> axis square

>> grid off

>> set(h,'erasemode','xor','markersize',10)

>> while 1

drawnow

x=x+0.01;

y=sin(x)-0.01;

set(h,'xdata',x,'ydata',y)

if(x>pi)|(y<-1)

x=[-pi:0.02:pi];

y=sin(x);

end

end
 
 

45、创建一个三维曲面z=x^2+y^2的动画。

>> x=[-2:0.2:2];

>> [x1,y1]=meshgrid(x);

>> z=x1.^2+y1.^2+eps;

>> surf(z);

>> ta=axis;

>> ft=moviein(40);

>> for i=1:40

surf(sin(2*pi*i/20)*z,z)

axis(ta)

ft(:,i)=getframe;

end

>> movie(ft,20)

46、通过调整Z的数值来建立peaks函数的动画。

>> z=peaks;

>> surf(z);

>> axis tight

>>set(gca,'nextplot',

'replacechildren');

>> for i=1:20

surf(sin(2*pi*i/20)*z,z)

f(i)=getframe;

end

>> movie(f,30)

47、cool色图+faceted系统默认颜色阴影和默认色图jet+interp颜色阴影绘制peaks函数图。

>> z=peaks;

>> surf(z)

>> colormap(cool)

>> shading faceted

>>

>> z=peaks;

>> surf(z)

>> colormap(jet)

>> shading interp

48、创建一个三维表面图并设置不同的视点。

>> [x,y]=meshgrid([-3:0.2:3]);

>> z=x.*exp(-x.^2-y.^2);

>> surf(z)

>> [ax,el]=view

ax =

  -37.5000

el =

    30

>> view(30,-30)

 

49、绘图工具栏介绍。