MATLAB绘图.doc

实验二MATLAB绘图与处理一实验目的掌握matlab常用的二维图形和三维图形的绘制方法,并会对图形进行处理；掌握符号函数(显函数、隐函数和参数方程)的画图,掌握空间曲线和空间曲面的绘图,会对所绘图形进行加格栅,图例和标注等一些简单的处理；掌握散点数据绘图基本方法。二实验原理与方法21 二维曲线图Matlab作图是通过描点、连线来实现的，故在画一个曲线图形之前，必须先取得该图形上的一系列的点的坐标（即横坐标和纵坐标），然后将该点集的坐标传给Matlab函数画图.1使用函数（命令）为plot，其基本使用格式为：（1）plot（y）功能：绘制数组y对应于其元素序数的二维曲线图,如果y为复数向量, 则绘制虚部对于实部的二维曲线图.例: y=0 0.6 2.3 5 8.3 11.7 15 17.7 19.4 20; plot(y) 由于y矢量有10个元素,x坐标自动定义为1 2 3 4 5 6 7 8 9 10.（2）. plot(x,y)功能：绘制由x,y所确定的曲线. 1)x,y是两组向量,且它们的长度相等,则plot(x,y)可以直观地绘出以x为横坐标,y为纵坐标的图形. 例:画正弦曲线: x=0:0.1:2*pi; y=sin(x); plot(x,y) 2)当plot(x,y)中,x是向量,y是矩阵时,则绘制y矩阵中各行或列对应于向量x的曲线.如果y阵中行的长度与x向量的长度相同,则以y的行数据作为一组绘图数据;如果y阵中列的长度与x向量的 长度相同,则以y的列数据作为一组绘图数据;如果y阵中行,列均与x向量的长度相同,则以y的每列数据作为一组绘图数据.例：下面的程序可同时绘出三条曲线.MATLAB在绘制多条曲线时,会按照一定的规律自动变化每条曲线的的颜色. x=0:pi/50:2*pi;y(1,:)=sin(x); y(2,:)=0.6*sin(x); y(3,:)=0.3*sin(x);%以上语句生成矩阵y. plot(x,y)或者还可以这样用: x=0:pi/50:2*pi; y= sin(x); 0.6*sin(x); 0.3*sin(x); plot(x,y) 3) 如果x,y是同样大小的矩阵,则plot(x,y)绘出y中各列相应于x中各列的图形. 例:x(1,:)=0:pi/50:2*pi; x(2,:)=pi/4:pi/50:2*pi+pi/4; x(3,:)=pi/2:pi/50:2*pi+pi/2; y(1,:)=sin(x(1,:); y(2,:)=0.6*sin(x(2,:); y(3,:)=0.3*sin(x(3,:); plot(x,y) x=x; y=y; figure plot(x,y)在这个例子中,x-3x101,y-3x101,所以第一个plot按列画出101条曲线,每条3个点;而x-101x3,y-101x3,所以第二个plot按列画出3条曲线,每条101个点. （3）plot(x,y，line specifiers)功能：Plot等画图命令（所有的高级画图函数，除了ezplot等以外）可以添加参数line specifiers来确定以下属性:Line style 线条的类型.Color 线条的颜色Marker symbol 数据点的形状Line Style SpecifiersSpecifierLine Style-实线-虚线：点线.-点划线Marker SpecifiersSpecifierMarker symbol+加号o圆圈*星号.点x乘号square or s正方块diamond or d菱形上三角v下三角右三角左三角pentagram or p五角星hexagram or h六角星ColorSpecifieColorr红色g绿色b蓝色c蓝绿色m品红色y黄色k黑色w白色注： 1）这些选项可以连在一起用,如:-.g表示绘制绿色的点划线,g+s表示用绿色的+号绘制曲线,数据点用方块标记。2)表示属性的符号必须放在同一对单引号中，3)与先后顺序无关; 4)指定的属性中,同一种属性不能有两个以上. 例：t=0:0.5:2*pi; y1=sin(t);plot(t,y1, s-r) （4）. plot(,PropertyName,PropertyValue,.)功能：设置包括线的宽度(LineWidth)、标记点的边缘颜色（MarkerEdgeColor)、填充颜色(MarkerFaceColor)、标记点的大小(MarkerSize)等其它绘图属性. 例：t=0:pi/20:pi; y=sin(4*t).*sin(t)/2; plot(t,y,-bs,LineWidth,2,. MarkerEdgeColor,k,. MarkerFaceColor,y,. MarkerSize,10)以上画图命令设置线的宽度为2，设置标记点边缘颜色为黑色，设置标记点填充颜色为黄色，设置标记点的尺寸为10 （5）. plot(x1, y1, 选项1, x2, y2, 选项2, )功能：多条曲线画在同一个坐标轴下。例:t=0:0.1:2*pi; y1=sin(t); y2=cos(t); y3=sin(t).*cos(t); plot(t,y1, -r,t,y2, :g,t,y3, *b) 2. ezplot，其基本使用格式如下： （1）ezplot(fun,a,b) 表示在axb范围内绘制fun(x)表示的函数的图形，如不指定a,b，则默认为-2x2.fun可以是编写的函数式M文件的函数名（也为文件名），也可以是描述一元函数的字符串。例：ezplot(exp(x)*cos(x) %字符串表示的函数例: 在-1，2上画的 图形编写函数式文件myfun.mfunction y=myfun(x)y=exp(0.5*x)+sin(3*x2);然后调用ezplot画图ezplot（myfun,-1,2）（2) ezplot(f(x,y),xmin,xmax,ymin,ymax) 表示在区间xminxxmax和 yminyymax绘制隐函数f(x,y)=0的函数图，如不指定xmin,xmax,ymin,ymax，则默认画图区域为：-2x2,-2y2.例：使用ezplot画由二元方程x2-y4=0所确定的隐函数的图像ezplot(x2-y4)注：在传递给ezplot的表达式中，其中涉及到的数组或矩阵的乘除法以及指数运算可以不必顾及维数，系统自动转换。例如，字符x2 在被ezplot使用时，自动处理为x.2 例 在-2，0.5，0，2上画隐函数的图 ezplot(exp(x)+sin(x.*y),-2,0.5,0,2)注：字符串中的 .*可以改为 *，ezplot自动转成 .*（3）ezplot(x(t),y(t),tmin,tmax)表示在区间tmint=1（默认值为1），则至少画出n+1个点（即至少把范围limits分成n个小区间），最大步长不超过(xmax-xmin)/n。（6）fplot(fun,lims,) 功能：允许可选参数tol，n和LineSpec以组合方式输入，顺序可以任意。例如：fplot(sin(x),0,2,100， 2e-4,r)与fplot(sin(x),0,2,2e-4,r，100)有完全相同的效果（7）X,Y= fplot(fun,limits,)功能： 返回横坐标与纵坐标的值给变量X和Y，此时fplot不画出图形。若想画出，可用命令plot(X,Y)。（8） .fplot(fun,limits,tol,n,LineSpec,P1,P2,)功能： 允许用户直接给函数fun输入参数P1，P2等，其中函数fun的定义形式为：y=fun(x,P1,P2,)若想用缺省的tol，n或LineSpec值，只需将空矩阵传递给函数即可。例如:在区间-2,2上画函数f(x)=a*x3+b*x2+c*x其中a=2,b=3,c=1编写函数式文件myfun.mfunction y=myfun(x,a,b,c)y=a*x3+b*x2+c*x;end编写脚本式M文件（主程序）a=2;b=3;c=1;fplot(mufun,-2,2,r,a,b,c)4对数坐标图在很多工程问题中,通过对数据进行对数转换可以更清晰地看出数据的某些特征,在对数坐标系中描绘数据点的曲线,可以直接地表现对数转换.对数转换有双对数坐标转换和单轴对数坐标转换两种.用loglog函数可以实现双对数坐标转换,用semilogx和semilogy函数可以实现单轴对数坐标转换(具体用法查看帮助文件)loglog() 表示 x、y坐标都是对数坐标系semilogx() 表示 x坐标轴是对数坐标系semilogy() 表示y坐标轴是对数坐标系例： x = 0:.1:10;semilogy(x,10.x) %所画曲线已将实际纵坐标做了对数变换例： x=logspace(-1,2); % loglog(x,exp(x),-s) 5具有两个纵坐标标度的图形在MATLAB中，如果需要绘制出具有不同纵坐标标度的两个图形，可以使用plotyy绘图函数。基本调用格式为：plotyy(x1,y1,x2,y2)其中x1,y1对应一条曲线，x2,y2对应另一条曲线。横坐标的标度相同，纵坐标有两个，左纵坐标用于x1,y1数据对，右纵坐标用于x2,y2数据对。例：用不同标度在同一坐标内绘制曲线y1=0.2e-0.5*x*cos(4x) 和y2=2e-0.5xcos(x)。程序如下：x = 0:0.01:20;y1 = 200*exp(-0.05*x).*sin(x);y2 = 0.8*exp(-0.5*x).*sin(10*x);plotyy(x,y1,x,y2);6创建二维多边形：fill功能：用颜色填充二维多边形。例：t = (1/16:1/8:1)*2*pi;x = sin(t);y = cos(t);fill(x,y,r)例：t = (1/16:1/8:1)*2*pi; x = exp(t).*sin(t); y = t.*cos(t);fill(x,y,k) 7极坐标图形 用polar函数可以画出极坐标图形, 该函数有两种表达形式：polar(theta, rho) ：创建一个幅角theta相对于半径rho的极坐标图polar(theta, rho, LineSpec) ：LineSpec为绘出的图形指定线型、颜色和标识例：x=0:0.01*pi:4*pi;y=sin(x/2)+x;polar(x,y,k-) 8 二维特殊函数图下列函数可以绘制其它的二维特殊函数图形：area填充绘图；bar条形图；barh水平条形图；comet彗星图；ezpolar简单绘制极坐标图；feather矢量图； hist直方图；pie饼状图；rose极坐标系下的柱状图等等，以上各函数均有不同的调用方法，详细内容同学们可以通过MATLAB在线帮助获得。小整理：其他各种二维绘图函数 bar 长条图errorbar 图形加上误差范围 fplot 较精确的函数图形 polar 极座标图hist 累计图rose 极座标累计图stairs 阶梯图stem 针状图fill 实心图feather 羽毛图compass 罗盘图quiver 向量场图22 三维曲线绘图基本命令: plot3使用方法：plot3函数与plot函数用法十分相似，其基本调用格式为：plot3(x1,y1,z1,选项1,x2,y2,z2,选项2,xn,yn,zn,选项n) 其中每一组x,y,z组成一组曲线的坐标参数，选项的定义和plot函数相同，可以指定所画曲线的线型和颜色等；当x,y,z是同维向量时，则x,y,z 对应元素构成一条三维曲线。当x,y,z是同维矩阵时，则以x,y,z对应列元素绘制三维曲线，曲线条数等于矩阵列数。例 在区间0，10*pi画出参数曲线x=sin(t),y=cos(t),z=t. t=0:pi/50:10*pi; plot3(sin(t),cos(t),t，k,linewidth,5)例：t=0:pi/100:20*pi;x=sin(t);y=cos(t);z=t.*sin(t).*cos(t);plot3(x,y,z);例：x=0 0 0 0 0 ;y=0 0 1 1 0;z=0 1 1 0 0;plot3(x,y,z)23 三维曲面绘图命令 1 meshgrid功能： 生成二元函数z = f(x,y)中x-y平面上的矩形定义域中数据点矩阵X和Y，或者是三元函数u = f(x,y,z)中立方体定义域中的数据点矩阵X,Y和Z。将由向量x，y（可以是不同方向的）指定的区域min(x)，max(x)，min(y)，max(y)用直线x=x(i),y=y(j)（i=1,2,length(x)，j=1,2,length(y)）进行划分。这样，得到了length(x)*length(y)个点所构成的点阵。这些点的横坐标用矩阵X表示，X的每个行向量与向量x相同；这些点的纵坐标用矩阵Y表示，Y的每个列向量与向量y相同。用法 a：X,Y = meshgrid(x,y) b：X,Y = meshgrid(x) c：X,Y,Z = meshgrid(x,y,z)说明： 对于形式a，输入向量x为x-y平面上矩形定义域的矩形分割线在x轴的值，向量y为x-y平面上矩形定义域的矩形分割线在y轴的值。输出向量X为x-y平面上矩形定义域的矩形分割点的横坐标值矩阵，输出向量Y为x-y平面上矩形定义域的矩形分割点的纵坐标值矩阵。对于形式b，等价于形式a：X,Y = meshgrid(x) = meshgrid(x,x)。对于形式c，输入向量x为立方体定义域的立方体分割平面在x轴上的值，输入向量y为立方体定义域的立方体分割平面在y轴上的值，输入向量z为立方体定义域的立方体分割平面在z轴上的值。输出向量X为立方体定义域中分割点的x轴坐标值,Y为立方体定义域中分割点的y轴坐标值,Z为立方体定义域中分割点的z轴坐标值。例：考察以下语句的作用x=2:4; %在X轴上取了横坐标为2,3,4的三个点y=6:9; %在Y轴上取了纵坐标为6,7,8，9的四个点。X,Y=meshgrid(x,y);Meshgird利用数组X和Y，生成XOY平面上的一个矩形网格点阵,该点阵行数为length(y)，列数为length(x)列。X为该点阵的点的横坐标所构成的矩阵，Y为该点阵的点的纵坐标所构成的矩阵。X和Y的维数均length(y)* length(x).运行结果如下：X= 2 3 4 2 3 4 2 3 4 2 3 4Y= 6 6 6 7 7 7 8 8 8 9 9 9命令2 mesh基本格式：mesh(X,Y,Z),mesh(Z)，meshc(),meshz()功能：画出由X,Y Z确定的三维网格图说明：1若X与Y均为向量，length（X）=n，length（Y）=m，而m，n=size（Z），空间中的点 (X(j),Y(i),Z(i,j) 为所画曲面网线的交点，分别地，X对应于z的列，Y对应于z的行。2若X与Y均为矩阵，则空间中的点 (X(i,j),Y(i,j),Z(i,j)为所画曲面的网线的交点。3.mesh(Z) :由n，m = size（Z）得，X=1：n与Y=1：m，其中z为定义在矩形划分区域上的网格点阵上的单值函数。4.meshc():画出等高线的三维网格图.5.meshz():画出带有底座的三维网格图.例：x=-3:0.125:3;y=-3:0.125:3X,Y = meshgrid(x，y);Z = peaks(X,Y); %得到山峰函数的函数值矩阵Z， peaks %称为山峰数函数，是matlab函数库的二元函数mesh(X,Y,Z); %山峰曲面图或者meshc(X,Y,Z);%带有等高线的山峰函数曲面或者meshz(X,Y,Z);%带有底座的山峰函数曲面例:绘出带有等高线的马鞍面解:从命令窗口输入如下命令x=-8:8;y=-8:8;X,Y=meshgrid(x,y);Z=X.2/42-Y.2/52;Meshc(X,Y,Z)命令 2 surf基本格式：surf(X,Y,Z),surf(Z)，surfc(),surfz()功能：surf同mesh命令的用法和格式相同,不同之处在于所得图形是一个彩色曲面而非彩色网格，其中色彩的添加由高度矩阵Z（函数值矩阵）确定，颜色与高度是对应的例：使用surf绘出马鞍面x=-8:8;y=-8:8;X,Y=meshgrid(x,y);Z=X.2/42-Y.2/52;surf(X,Y,Z)命令3 contour功能:曲面的等高线图格式一: contour(x,y,z,n)说明:(x,y)是平面xoy平面上矩形网格点阵的坐标矩阵,z为相应点的高度值矩阵（函数值矩阵）,等高曲线是一个平面曲线,n为等高线条数。例: X,Y = meshgrid(-2:.2:2,-2:.2:3);Z = X.*exp(-X.2-Y.2);contour(X,Y,Z,4)格式二: contour(x,y,z,h)说明:h为一向量，向量的长度即为所画等高线的条数，向量每个元素的值为每条等高线的高度值。注意,高度值必须在曲面的高度范围以内。X,Y = meshgrid(-2:.2:2,-2:.2:3);Z = X.*exp(-X.2-Y.2);contour(X,Y,Z,0.1,0.3) %两条等高线，高度值分别为0.1,0.3格式二: clabel(C,h)说明:在从命令contour生成的二维等高线结构C的位置上添加标签h。例 X,Y = meshgrid(-2:.2:2,-2:.2:3);Z = X.*exp(-X.2-Y.2);c,h=contour(X,Y,Z,0.1 ,0.3);clabel(c,h)命令4. contour3功能：绘制三维曲面等高线。格式:contour3(x,y,z,n)说明: (x,y)是平面xoy平面上矩形网格点阵的坐标矩阵, z为相应点的高度值矩阵（函数值矩阵），该函数绘制n条等高线的三维图,用法与从头人contour相同例 X,Y = meshgrid(-2:.2:2,-2:.2:3);Z = X.*exp(-X.2-Y.2);c,h=Contour3(X,Y,Z，4)Clabel(c,h)命令4 cylinder功能 ： 生成圆柱图形。该命令生成一单位圆柱体的x，y，z轴的坐标值。用户可以用命令surf或命令mesh画出圆柱形对象，或者用没有输出参量的形式而立即画出图形。格式及说明：(1) X,Y,Z = cylinder 返回一半径为1、高度为1的圆柱体的x，y，z轴的坐标值，圆柱体的圆周有20个距离相同的点。(2) X,Y,Z = cylinder(r) 返回一半径为r、高度为1的圆柱体的x，y，z轴的坐标值，为r圆柱体的圆周有20个距离相同的点。若r为一正实数，则为一直上直下的圆柱体；若r为一数组，则所得柱体不同高度上有不同半径，此时得到一个母线绕Z轴旋转而成的旋转曲面，数组r指定母线到Z轴的距离图7-36(3) X,Y,Z = cylinder(r,n) 返回一半径为r、高度为1的圆柱体的x，y，z轴的坐标值，圆柱体的圆周有指定的n个距离相同的点。(4) cylinder() : 没有任何的输出参量，直接画出圆柱体。例：x,y,z = cylinder(3)mesh(x,y,z)例：x,y,z = cylinder(1:10);mesh(x,y,z)例：t =1:4;X,Y,Z = cylinder(2+cos(t),5);surf(X,Y,Z)命令5 sphere功能: 生成球面格式:x,y,z=sphere(n)说明:生成三维直角坐标系中的单位球体坐标,该单位球体有个面.该命令并未画图,只是返回坐标矩阵,用户可用surf(x,y,z)或mesh(x,y,z)画图.例:x,y,z=sphere(50);mesh(x,y,z);3图形处理的基本技术除了提供强大的绘图功能外，MATLAB语言还有极为强大的图形处理能力下面介绍一些图形处理技术，包括图形控制、图形标注、图形保持以及子图的绘制等(1). 图形控制MATLAB语言中较常用的图形控制函数有坐标轴控制函数axis、坐标轴缩放函数zoom和坐标网格函数grid等(i). axis函数控制坐标轴的特征在缺省情况下MATLAB自动选择图形的横、纵坐标的比例，如果你对这个比例不满意，可以用axis命令控制，常用的有：axis(xmin xmax ymin ymax)： 中分别给出了x轴和y轴的最小、最大值；axis equal 或 axis(equal)： x轴和y轴单位长度相同；axis square 或 axis(square)： 图框呈方形；axis off 或 axis(off)： 清除坐标刻度例：x=0:0.025:pi/2; plot(x,tan(x),-ko) %使用axis命令设定坐标轴之前的图axis(0 pi/2 0 5) %使用axis命令设定坐标轴之后的图形 (ii). zoom函数控制坐标轴的缩放：zoom函数可以实现对二维图形的缩放，该函数在处理局部较为密集的图形中有很大作用常用的调用格式有：Zoom： 在zoom on 和zoom off之间切换；zoom on ：允许对图形进行缩放；zoom off ：禁止对图形进行缩放；zoom xon ：允许x 轴缩放；zoom yon ：允许y轴缩放；zoom out：恢复进行的一切缩放当zoom处于on 状态时，可以通过鼠标进行图形缩放，单击鼠标左键将光标处的图形放大一倍；而单击鼠标右键将光标处的图形缩小一倍；双击鼠标左键则将会恢复缩放前的状态，即取消一切缩放操作应当注意，对图形的缩放不会影响图形的原始尺寸，也不会影响图形的横纵坐标的比例，即不会改变图形的基本结构(iii). grid函数控制平面图形的坐标网格：MATLAB提供了平面网图函数grid用于绘制坐标网格，提高图形显示效果grid函数的调用格式如下：grid on： 在图形中绘制坐标网格；grid off ： 取消坐标网格单独的函数grid将实现grid on 与grid off两种状态之间的转换例：x=0:0.1*pi:2*pi; y=sin(x);plot(x,y) grid on (2). 图形的标注MATLAB语言还提供了丰富的图形标注函数供用户自由地标注所绘制的图形(i).坐标轴标注和图形标题xlabel， ylabel： 为x, y坐标轴添加标注title ： 为图形添加标题xlabel(标注内容，属性1, 属性值1, 属性2, 属性值2,) %属性包括标注文本的属性，包括字体大小、字体名等三个函数的调用结果的区别仅在于标注所处的位置不同，title 给出的标注将置于图的顶部，而xlabel 和ylabel则分置于相应的坐标轴的边上例如：x=1:0.1*pi:2*pi; y=sin(x);plot(x,y) xlabel(x(0-2pi), fontweight, bold);title(y=sin(x),fontsize, 12, fontweight, bold) 在标注过程中经常会遇到特殊符号的输入问题，为了解决这个问题，MATLAB语言提供了相应的字符转换,如：alpha；beta；gamma；delta；epsilon；zeta；pi；omega；Omega等等alphaepsiloninftybetaetaintgamma GammapartialdeltaDeltaleftarrowthetaThetarightarrowlambdaLambdadownarrowxiXiuparrowpiPidivomegaOmegatimessigmaSigmapmphiPhileqpsiPsigeqrhotauneqmuzetaforallnuchiexists用户也可以对文本标注进行显示控制，如：bf 黑体，it 斜体，rm 标准形式.例如：x=-10:0.1:10; y=exp(-x.2/2);plot(x,y, -)title(bf y=e-x2/2) (ii). 文本标注MATLAB对图形进行文本注释所提供的函数为text和gtexttext函数的基本调用格式：text(x,y, 标注文本及控制字符串)，其中(x,y)给定标注文本在图中添加的位置.例如：.x=1:0.1*pi:2*pi; y=sin(x);.plot(x,y).text(3*pi/4,sin(3*pi/4), leftarrow sin(3pi/4)=0.707).text(5*pi/4,sin(5*pi/4),sin(5pi/4)=,num2str(sin(5*pi/4),rightarrow,HorizontalAlignment, right, Fontsize, 12). %属性HorizontalAlignment用来控制文本标注输入起点是在标注本身的左侧还是右侧图11 图12交互式文本输入函数gtext. 使用该函数，用户可以通过使用鼠标来选择文本输入的点，单击后，系统将把指定的文本输入到所选的位置上x=1:0.1*pi:2*pi; y=sin(x);plot(x,y)gtext(y=sin(x), Fontsize, 12)执行该函数时，将鼠标放在图形上会出现“+”字型交叉线供用户添加标注的点，选择添加标注的位置后，单击鼠标左键即可在该位置上添加标注(iii). 图例标注在对数值计算结果进行绘图时，经常会出现在同一张图形中绘制多条曲线的情况，这时可以使用legend命令为曲线添加图例以便于区别它们legend 函数能够为图形中的所有曲线进行自动标注，以其输入变量作为标注文本，具体调用格式如下：legend(标注1, 标注2,): 标注1，标注2等分别对应绘图过程中按绘制先后顺序所生成的曲线例：x=0:0.1*pi:2*pi; y=sin(x); z=cos(x); plot(x,y,k-o,x,z,k-h)legend(sin(x), cos(x) 可以用鼠标拖动图例框改变其在图中的位置也可以在legend函数调用时进行简单的定位设置：legend(标注1, 标注2, 定位代号)MATLAB 给出了6个定位代号，具体说明如下：0: 自动定位，使得图标与图形重复最少；1: 置于图形的右上角（默认值）；2: 置于图形的左上角；3: 置于图形的左下角；4: 置于图形的右下角；-1: 置于图形的右外侧关于标注位置，没有必要记住，可以通过在线帮助获得help legend图例标注后，也可以用鼠标来调整图例标注的位置(3). 图形保持与子图(i).图形保持在绘图过程中，经常会遇到在已存在的一张图中添加新的曲线，这就要求保持已存在的图形，MATLAB语言中实现该功能的函数是hold.hold on 启动图形保持功能，此后绘制的图形将添加到当前的图形窗口中，并自动调整坐标轴的范围；hold off 关闭图形保持功能，新绘制图形将覆盖原图形hold 在hold on 和 hold off之间切换例如：x=0:0.1*pi:2*pi; y=sin(x); z=cos(x);plot(x,y,k-*) hold onplot(x,z, k-o) plot(x,y+z,k-h)legend(sin(x), cos(x), sin(x)+cos(x), 0)hold off(ii).子图在绘图过程中，经常需要将几个图形在同一图形窗口中表示出来，但又不在同一个坐标系中绘制，此时要用到函数subplot调用格式如下：subplot(m,n,p) 将一个图形窗口分割成m*n个小窗口，可以通过参数p分别对若干子绘图区域进行操作，子绘图区域的编号为按行从左至右编号如果p是一个向量，则创建一坐标轴，包含所有罗列在p中的小窗口subplot(1,1,1) 返回非分割状态。例如：在四个子图中绘制不同的三角函数图x=0:0.1*pi:2*pi; subplot(2,2,1) %第1个绘图子域plot(x,sin(x), -*); axis(0 2*pi -1 1); title(