数学2-用MATLAB绘制二维-三维图形.pptVIP

2021/3/27 * 3)极坐标系 t=linspace(0,2*pi); r=1+0*t; 或 r=ones(1,100) %不可写成ｒ＝１,因为必须是向量才行。 polar(t,r) 2021/3/27 * 2.2.1 三维曲线的绘制（plot3命令） plot3是二维绘图命令plot的推广,其使用方法和功能基本相同。 格式: plot3（x,y,z,颜色＋线型＋点型,...) 当（x,y,z）是一个坐标时,则在空间相应位置画出一个点。用法: plot3（x,y,z,r*) 当（x,y,z）是点列时,则画出依次连接这些点的曲线。用法: plot3（x,y,z,r－－) ％红色实线 其中,曲线的颜色和线型取值同前表。 2021/3/27 * 例2-6 画出螺旋线,x＝sin(t),y=cos(t),z=t, t ,上的一段曲线。 t=0:pi/50:10*pi; x=sin(t); y=cos(t); z=t; plot3(x,y,z,g-,linewidth,3) grid 2021/3/27 * 2.2.2 三维曲面的绘制(mesh、surf命令) 对于二元函数z＝f(x,y),设其定义域为[a,b]×[c,d],则其几何图形为空间坐标系中一片曲面,该曲面的投影域即为函数的定义域,使用mesh、surf等命令可以绘制出相应曲面。 格式: mesh（x,y,z） ％绘制网格曲面 surf（x,y,z） ％绘制光滑曲面 其中,x、y分别为投影域上网格划分节点处对应的横坐标矩阵和纵坐标矩阵,可由meshgrid命令生成; z为与投影域上网格划分节点（x,y）对应的函数值（坐标z）矩阵。具体步骤: 2021/3/27 * 1)对投影域进行划分: x＝a:p1:b y＝c:p2:d 2）按上述划分生成投影域上全部网格节点的坐标矩阵。 [X,Y]=meshgrid(x,y) 3)根据函数表达式生成全部网格节点出对应的函数值矩阵z: z＝f(X,Y) 4）顺序连接已经产生的空间点（x,y,z)绘制相应曲面: mesh（X,Y,Z) surf（X,Y,Z) shading flat ％去除网格线。 2021/3/27 * 例2-7画出矩形域[-1,1]×[-1,1]旋转抛物面:z=x2+y2. x=linspace(-1,1,100); y=x; [X,Y]=meshgrid(x,y); %生成矩形区[-1,1]×[-1,1]的网格坐标矩阵 Z=X.^2+Y.^2; subplot(1,2,1) mesh(X,Y,Z); subplot(1,2,2) surf(X,Y,Z); shading flat; %对曲面z=x2+y2平滑并除去网格 2021/3/27 * 未进行 Shading flat 2021/3/27 * Shading Flat 之后 2021/3/27 * 2.8 在圆形域x2+y2≦1上绘制旋转抛物面:z= x2+y2. x=linspace(-1,1,300); y=x; [X,Y]=meshgrid(x,y); %生成矩形域[-1,1] × [-1,1] 网格节点坐标矩阵 Z=X.^2+Y.^2; i=find(Z1); %找出圆域x2+y2≦1之外的函数值坐标点i Z(i)=NaN; %对圆域x2+y2≦1之外的坐标点i处函数值进行“赋空” subplot(1,2,1) mesh(X,Y,Z); %网格曲面 subplot(1,2,2) surf(X,Y,Z); %光滑曲面 shading flat; 2021/3/27 * 2021/3/27 * 注:find的用方法 find（A）返回矩阵A中非零元素所在位置 A = [1 0 4 -3 0 0 0 8 6]; X = find(A) X = 1 3 4 8 9find（A5） 返回矩阵A中大于5的元素所在位置 find(A5) ans = 8 9 [i,j,v]=find(A) 返回矩阵A中非零元素所在的行i,列j,和元素的值v(按所在位置先后顺序输出) A=[3 2 0; -5 0 7; 0 0 1]; [i,j,v]=find(A) i = 1 2 1 2 3 j = 1 1 2 3 3 v = 3 -5 2 7 1 2021/3/27 * 2.2.3 特殊图形和