第5章 绘图2.pptVIP

第5章 MATLAB绘图（二） 5.2 三维绘图 5.3 三维图形的精细处理 5.4 图像与动画 5.2 三维绘图5.2.1 三维曲线图 plot3(x,y,z,线型) %绘制三维曲线 说明：x，y，z必须是相同尺寸的数组，当是向量时则绘制一条三维曲线，当是矩阵时绘制多条曲线，三维曲线的条数等于矩阵的列数。 5.2.2三维曲面图 1.平面网格坐标矩阵的生成 绘制z=f(x,y)所代表的三维曲面图，先要在x-y平面选定一矩形区域，假定矩形区域D=[a,b]×[c,d],然后将[a,b]在x方向分成m份，将[c,d]在y方向分成n份,由各划分点分别作平行于两坐标轴的直线，将区域D分成m×n个小矩形，生成代表每一个小矩形顶点坐标的平面网格坐标矩阵，最后利用有关函数绘图。 (1)利用矩阵运算生成。 x=a:dx:b; y=(c:dy:d); X=ones(size(y))*x; Y=y*ones(size(x)); 矩阵X的每一行都是向量x，行数等于向量y的元素的个数，矩阵Y的每一列都是向量y,列数等于向量x的元素的个数。于是，X和Y相同位置上的元素(X(i,j),Y(i,j))恰好是区域D的(i,j)网格点的坐标。若根据每一个网格点上的x,y坐标求函数值z，则得到函数值矩阵Z.显然，X,Y,Z各列或各行所对应坐标，对应于一条空间曲线，空间曲线的集合组成空间曲面。 (2)利用meshgrid函数生成。 [X,Y]＝meshgrid(x,y) %产生XY矩形网格 说明：x和y分别是有n个和m个元素的一维数组，X和Y都是n×m的矩阵，每个（X,Y）对应一个网格点；如果y省略，则X和Y都是n×n的矩阵。 例13 已知6x30，15y36，求不定方程2x+5y=126的整数解。 程序如下： x=5:29; y=14:35; [x,y]=meshgrid(x,y); %在[5,29]×[14,35]区域生成网格坐标 z=2*x+5*y; k=find(z==126); %找出解的位置 x(k),y(k) %输出对应位置的x,y即方程的解 2. 绘制三维曲面的函数 1)三维网格图 三维网格图就是将平面上的网格点（X,Y）对应z值的顶点画出，并将各顶点用线连接起来。 mesh(X,Y,Z,C) %绘制网格点数据对应的三维网格，C是指定在不同高度下的颜色范围。 例14 绘制z=x2+y2的三维网格图 x=0:10; [X,Y]=meshgrid(x) %y省略则表示x=y Z=X.^2+Y.^2; mesh(X,Y,Z) 2) 三维曲面图 三维曲面图与网格图相似，但不同的是网格图中网格范围内的区域为空白，而三维曲面图则用颜色来填充。 surf(X,Y,Z,C)%绘制网格点数据对应的三维曲面图 3. 标准三维曲面 1)sphere函数绘制三维球面 调用格式为： [x,y,z]=sphere(n) 该函数将产生(n+1)×(n+1)矩阵x,y,z，采用这3个矩阵可以绘制出圆心位于原点、半径为1的单位球体。若在调用该函数时不带输出参数，则直接绘制所需球面。n决定了球面的圆滑程度，其缺省值为20。若n值较小，则将绘制出多面体表面图。 subplot(1,3,1); [x,y,z]=sphere;surf(x,y,z); title(‘sphere’); subplot(1,3,2); [x,y,z]=sphere(100);mesh(x,y,z); title(‘n=100’); subplot(1,3,3); [x,y,z]=sphere(5);mesh(x,y,z); title(‘n=5’); 2)cylinder函数绘制三维柱面 调用格式为： [x,y,z]= cylinder(R,n) 其中R是一个向量,存放柱面各个层次上的半径。 如：cylinder(3)生成一个圆柱； cylinder([10,1])生成一个圆锥； t=0:pi/100:4*pi;R=sin(t); cylinder(R,30)生成一个正弦型柱面。 另外，生成矩阵的大小与R向量的长度及n有关。其余与sphere函数相同。 3)peaks 函数，称为多峰函数，常用于三维曲面的演示。该函数可以用来生成绘图数据矩阵，矩阵元素由函数在矩形区域[-3,3]×[-3,3]的等分网格点上的函数值确定。 例如peaks(30)将生成一个30×30矩阵z，即分别沿x和y方向将区间[-3,3]等分成29份，并计算这些网格点上的函数值。缺省的等分数是48，即p=peaks将生成一个49×49矩阵p。也可以根据网格坐标矩阵x