2-Matlab绘图.ppt

二 三维图形的绘制 1、用plot3(x,y,z,颜色+线型+点型)命令 例10 画出螺旋线 t=0:0.1:50; x=sin(t); y=cos(t); z=t; plot3(x,y,z,m*) grid on; %加注坐标网格 基本步骤： 1. 生成二维网格点 2. 计算函数在网格点上的值 3. 绘制函数图形 2 用meshgrid，surf命令 1. 用meshgrid生成二维网格点 [X,Y]=meshgrid(x,y) 其中x,y是向量，通过上述函数就可将x和y指定的区域转换成为矩阵X和Y。 观察meshgrid指令的效果。 a=-0.98; b=0.98; c=-1; d=1; n=10; x=linspace(a,b,n); y=linspace(c,d,n); [X,Y]=meshgrid(x,y); plot(X,Y,+) 计算函数在网格点上的值 z=X.^2+Y.^2; 绘制函数图形 mesh(x,y,z) %画网格曲面，已知点和附近的点以直线连接 surf(x,y,z) %画光滑曲面，已知点和附近点以平面连接 例11 x=-7.5:0.5:7.5; y=x; [x,y]=meshgrid(x,y); %生成矩形域上网格节点矩阵 z=sin(sqrt(x.^2+y.^2))./(sqrt(x.^2+y.^2)+eps); %生成函数值矩阵 subplot(2,1,1) %画子图1 mesh(x,y,z) %画曲面 subplot(2,1,2) %画子图2 surf(x,y,z) %画光滑曲面 shading flat %对曲面平滑并除去网格线 2 ezsurf命令 ezsurf(f(x,y),[a,b,c,d]) %画函数f(x,y)的曲面图 ezsurf(sin(x^2+y^2)/(x^2+y^2)) 简捷绘图 peaks函数 % 产生N*N的高斯分布矩阵，对应的图形是凹凸有致的曲面 [x,y,z]=peaks; surf(x,y,z) sphere(N) % 产生(N+1)* (N+1)矩阵，绘制单位球面 [x,y,z]=sphere; %默认值N=20 surf(x,y,z) cylinder(r,N) % 产生(N+1)* (N+1)矩阵，绘制半径为r的圆柱体 [x,y,z]=cylinder; %默认r=1，N=20 surf(x,y,z) ellipsoid(xc,yc,zc,xr,yr,zr) %产生中心在(xc,yc,zc)截距分别为xr,yr,zr的椭球 [x,y,z]=ellipsoid(0,0,0,1,3,7); surf(x,y,z) 例12 绘制高度为10的圆柱面x^2+y^2=4 r=2; n=1000; [xb,yb,zb]=cylinder(r,n); mesh(xb,yb,10*zb) x=-2:0.01:2; y=x; [X,Y]=meshgrid(x,y) Z=X.^2+Y.^2; surf(X,Y,Z) shading flat 4 用find语句精细绘制曲面 再看 x=-2:0.01:2; y=x; [xx,yy]=meshgrid(x,y); zz=xx.^2+yy.^2; ii=find(zz1); %找出圆域外的坐标点i zz(ii)=NaN; %对圆域外面的坐标点i处函数值赋空 surf(xx,yy,zz); shading flat; clf %清除图形窗口已有的内容 shg %显示图形窗口。 figure(p) %打开第p个窗口绘图 figure(1) plot(sin(0:0.1:pi)) figure(2) plot(cos(0:0.1:pi)) 若干有用的命令 … % 续行号 clear、 clear x % 清除工作空间的已有变量。 例12