MATLAB画图函数.docx

1、MATLAB画图函数目录MATLAB函数画图 2基本XY平面绘图命令 2Plot 2loglog 2semilogx 2semilogy 2axis 3xlabel 3ylabel 3title 3legend 3bar 3errorbar 3fplot 3polar 3hist 3rose 3stairs 3stem 3fill 3feather 3Compass 3quiver 3基本XYZ立体绘图命令 5mesh 5surf 5waterfall 5meshc 6surfc 6contour3 6contour 6plot3 6MATLAB的图视化功能 6MATLAB函数画图MATLAB

2、不但擅长於矩阵相关的数值运算，也适合用在各种科学目视表示（Scientific visualization）。本节将介绍MATLAB基本xy平面及xyz空间的各项绘图命令，包含一维曲线及二维曲面的绘制、列印及存档。基本XY平面绘图命令Plot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y座标。下例可画出一条正弦曲线：close all; %linspace(5,100,20)和5：5：100的区别在于前者已知元素总个数而不知道步长，后者已知步长不知元素个数，这两者的效果是一样的x=linspace(0, 2*pi, 100); % 100个点的x座标y=sin(

3、x); % 对应的y座标plot(x,y);=小整理：MATLAB基本绘图函数plot: x轴和y轴均为线性刻度（Linear scale）loglog: x轴和y轴均为对数刻度（Logarithmic scale）semilogx: x轴为对数刻度，y轴为线性刻度semilogy: x轴为线性刻度，y轴为对数刻度=若要画出多条曲线，只需将座标对依次放入plot函数即可：plot(x, sin(x), x, cos(x);若要改变颜色，在座标对后面加上相关字串即可：plot(x, sin(x), c, x, cos(x), g);若要同时改变颜色及图线型态（Line style），也是在座标对

4、后面加上相关字串即可：plot(x, sin(x), co, x, cos(x), g*);=小整理：plot绘图函数的叁数字元 颜色 字元 图线型态y 黄色 . 点k 黑色 o 圆w 白色 x xb 蓝色 + +g 绿色 * *r 红色 - 实线c 亮青色 : 点线M 锰紫色 -. 点虚线- 虚线=图形完成后，我们可用axis(xmin,xmax,ymin,ymax)函数来调整图轴的范围：axis(0, 6, -1.2, 1.2);此外，MATLAB也可对图形加上各种注解与处理：xlabel(Input Value); % x轴注解ylabel(Function Value); % y轴注解

5、title(Two Trigonometric Functions); % 图形标题legend(y = sin(x),y = cos(x); % 图形注解grid on; % 显示格线=我们可用subplot来同时画出数个小图形於同一个视窗之中：subplot(2,2,1); plot(x, sin(x);subplot(2,2,2); plot(x, cos(x);subplot(2,2,3); plot(x, sinh(x);subplot(2,2,4); plot(x, cosh(x);MATLAB还有其他各种二维绘图函数，以适合不同的应用，详见下表。=小整理：其他各种二维绘图函数ba

6、r 长条图errorbar 图形加上误差范围fplot 较精确的函数图形polar 极座标图hist 累计图rose 极座标累计图stairs 阶梯图stem 针状图fill 实心图feather 羽毛图Compass 罗盘图quiver 向量场图=以下我们针对每个函数举例。当资料点数量不多时，长条图是很适合的表示方式：close all; % 关闭所有的图形视窗x=1:10;y=rand(size(x);bar(x,y);=如果已知资料的误差量，就可用errorbar来表示：下例以单位标准差来做资料的误差量：x = linspace(0,2*pi,30);y = sin(x);e = std

7、(y)*ones(size(x);errorbar(x,y,e)=对于变化剧烈的函数，可用fplot来进行较精确的绘图，会对剧烈变化处进行较密集的取样，如下例：fplot(sin(1/x), 0.02 0.2); % 0.02 0.2是绘图范围=若要产生极座标图形，可用polar：theta=linspace(0, 2*pi);r=cos(4*theta);polar(theta, r);=对于大量的资料，我们可用hist来显示资料的分布情况和统计特性。下面几个命令可用来验证randn产生的高斯乱数分 ：x=randn(5000, 1); % 产生5000个 ?=0，?=1 的高斯乱数hist

8、(x,20); % 20代表长条的个数=rose和hist很接近，只不过是将资料大小视为角度，资料个数视为距离：x=randn(1000, 1);rose(x);stairs可画出阶梯图：x=linspace(0,10,50);y=sin(x).*exp(-x/3);stairs(x,y);=stems可产生针状图，常被用来绘制数位讯号：x=linspace(0,10,50);y=sin(x).*exp(-x/3);stem(x,y);=stairs将资料点视为多边行顶点，并将此多边行涂上颜色：x=linspace(0,10,50);y=sin(x).*exp(-x/3);fill(x,y,b

9、); % b为蓝色=feather将每一个资料点视复数，并以箭号画出：theta=linspace(0, 2*pi, 20);z = cos(theta)+i*sin(theta);feather(z);=compass和feather很接近，只是每个箭号的起点都在圆点：theta=linspace(0, 2*pi, 20);z = cos(theta)+i*sin(theta);compass(z);基本XYZ立体绘图命令在科学目视表示（Scientific visualization）中，三度空间的立体图是一个非常重要的技巧。本章将介绍MATLAB基本XYZ三度空间的各项绘图命令。mesh

10、和plot是三度空间立体绘图的基本命令，mesh可画出立体网状图，plot则可画出立体曲面图，两者产生的图形都会依高度而有不同颜色。下列命令可画出由函数 形成的立体网状图:=x=linspace(-2, 2, 25); % 在x轴上取25点y=linspace(-2, 2, 25); % 在y轴上取25点xx,yy=meshgrid(x, y); % xx和yy都是21x21的矩阵zz=xx.*exp(-xx.2-yy.2); % 计算函数值，zz也是21x21的矩阵mesh(xx, yy, zz); % 画出立体网状图=surf和mesh的用法类似：x=linspace(-2, 2, 25)

11、; % 在x轴上取25点y=linspace(-2, 2, 25); % 在y轴上取25点xx,yy=meshgrid(x, y); % xx和yy都是21x21的矩阵zz=xx.*exp(-xx.2-yy.2); % 计算函数值，zz也是21x21的矩阵surf(xx, yy, zz); % 画出立体曲面图为了方便测试立体绘图，MATLAB提供了一个peaks函数，可产生一个凹凸有致的曲面，包含了三个局部极大点及三个局部极小点，其方程式为：=要画出此函数的最快方法即是直接键入peaks：peaksz = 3*(1-x).2.*exp(-(x.2) - (y+1).2) .- 10*(x/5

12、- x.3 - y.5).*exp(-x.2-y.2) .- 1/3*exp(-(x+1).2 - y.2)=亦可对peaks函数取点，再以各种不同方法进行绘图meshz可将曲面加上围裙：x,y,z=peaks;meshz(x,y,z);axis(-inf inf -inf inf -inf inf);=waterfall可在x方向或y方向产生水流效果：x,y,z=peaks;waterfall(x,y,z);axis(-inf inf -inf inf -inf inf);=下列命令产生在y方向的水流效果：x,y,z=peaks;waterfall(x,y,z);axis(-inf inf

13、-inf inf -inf inf);=meshc同时画出网状图与等高线：x,y,z=peaks;meshc(x,y,z);axis(-inf inf -inf inf -inf inf);=surfc同时画出曲面图与等高线：x,y,z=peaks;surfc(x,y,z);axis(-inf inf -inf inf -inf inf);=contour3画出曲面在三度空间中的等高线：contour3(peaks, 20);axis(-inf inf -inf inf -inf inf);=contour画出曲面等高线在XY平面的投影：contour(peaks, 20);=plot3可画出

14、三度空间中的曲线：t=linspace(0,20*pi, 501);plot3(t.*sin(t), t.*cos(t), t);=亦可同时画出两条三度空间中的曲线：t=linspace(0, 10*pi, 501);plot3(t.*sin(t), t.*cos(t), t, t.*sin(t), t.*cos=MATLAB的图视化功能1. MATLAB的图视化概论数据图视化能使人们用视觉器官直接感受到数据的许多内在本质。因此，数据可视化是人们研究科学、认识世界所不可缺少的手段。MATLAB不仅数值计算方面是一个优秀的科技应用软件， 在数据可视化方面也具有上佳表现。MATLAB具有二维、三维

15、乃至四维的图形表现能力。可以从线型、边界面、色彩、渲染、光线、视角等方面把数据的特征表现出来。MAT LAB的图视化功能是建立在一组“图形对象”的基础之上的。“图形对象”的核心是图形的句柄（Granhics Handle）操作。MATLAB的有两个层次的绘图指令：（1） 底层（Lowleve）绘图指令：是直接对句柄进行操作。 底层绘图指令控制和表现数据图形的能力比高层绘图指令强。特点是灵活多变，较难掌握。（2） 高层（Highlevel）绘图指令：建立在底层指令上的绘图指令。 最常用的是高层绘图指令。高层绘图指令简单明了容易掌

16、握，本章介绍高层绘图指令。本章内容按“前易后难”的原则安排。 最常用的二个绘图指令是： plot ;mesh 2. 二维图形（1） plot函数以下例子用来体会plot 的基本的绘图原理。例：绘向量得折线图:hold onx=2.3,4.3,3,4,4.9,1.5,2.8,4.6,5.5;plot(x)plot(x,ro)注1：plot 绘图的基本素材是二维点组（x,y）(1=1,2,….n)。二维点组（x,y）(1=1,2,….n)的定义形式：*1） x=2.3,4.3,3,4,4.9,1.5,2.8,4.6,5.5;*2）

17、 y=0:0.1:5这种定义方法，默认横坐标是自然数（1，2，3，4….）*3)t=0:pi/100:2*pi x=sin(t)*4) x=1.5,2.3,2.8,3,4,4.3,4.6,4.9,5.5;y=x.2这种定义方法，要注意自变量保持升序。自变量与应变量的体积的一致。注2：plot 绘图的基本原理是依（x,y）(1=1,2,….n)排列顺序用直线连接。曲线光滑与否与点数相关。hold offt=0:pi/3:2*pi;x=sin(t);plot(t,x,r-)hold ont=0:pi/5:2*pi;x=sin(t);plot(t,x,b-)(1)

18、坐标系定制用于对坐标轴进行管理与控制，如刻度，外观，文字说明等*1）坐标轴定制指令（axis）axis用于对坐标轴刻度进行管理与控制。指令形式与作用说明如下： AXIS(XMIN XMAX YMIN YMAX) 设置 x- and y-axes刻度。AXIS(XMIN XMAX YMIN YMAX ZMIN ZMAX) 设置 x- and y-axes和 z-axes刻度。V = AXIS 返回当前图形行向量的刻度设置XMIN XMAX YMIN YMAX或(XMIN XMAX YMIN YMAX ZMIN ZMAX)。AXIS AUTO 返回刻度设置的系统默认值AXIS TIGHT

19、 依数据设置刻度AXIS IJ 设置坐标轴的原点在左上角AXIS XY 设置坐标轴的原点在左下角 AXIS EQUAL 设置坐标轴的比例因子相等。 AXIS IMAGE AXIS SQUARE AXIS NORMAL AXIS VIS3D AXIS OFF AXIS ON例：XMIN=1;XMAX=10;YMIN=10;YMAX=100;AXIS(XMIN XMAX YMIN YMAX)plot(1,50,3,60,5,20,3)*2)其它坐标系:polar例1：polar(THETA, RHO)t=0:0.1:2*pi;r=t;polar(t,r)例2:对数-对数t=0:0.1:2*pi;r

20、=t;semilogx(t,r)(2)(3) 图视效果强化例：加入格栅；坐标轴标志；文本说明等clf; hold offt=linspace(0,pi*3,30);x=sin(t);hold ony=cos(t);plot(t,x,r-,t,y,g-)grid% 加入格栅xlabel(x轴)ylabel(y轴)title(正弦与余弦曲线)text(1,0,正弦)%text(x,y,正弦)text(3,0,余弦)legend(sin(x),cos(x),3)%LEGEND(string,Pos) places the legend in the specified,%0 = Automatic

21、best placement (least conflict with data)%1 = Upper right-hand corner (default)%2 = Upper left-hand corner%3 = Lower left-hand corner%4 = Lower right-hand corner% -1 = To the right of the plot%按鼠表 left mouse button 拖legend到指定的位置(2) 子图clf; hold offt=linspace(0,pi*3,30);x=sin(exp(t);subplot (2,2,2) %(

22、n,m,p(0<p<m*n)plot (t,x,r-)y=exp(sin(t);subplot (2,2,3)plot (t,y,g-)(3) 特殊二维图形例：误差图（errorbar）clf;x=0:0.1:4;y=zeros(size(x);e=rand(size(x);yu=y+e;yd=y-e;errorbar(x,y,e)hold onplot(x,yu,r-);plot(x,yd,r-);(3) 绘图工具mmaxes prop value…修改绘图坐标轴的属性mmcxy（or）xymmcxy显示图上鼠标的xy坐标mmdraw prop value…

23、; 在图上画直线rnmfill（x，y，z，c，lb，ub）填充两条曲线间区域mmgetxy（N）使用鼠标获取xy坐标mmline prop value… 修改所画线条的属性mmtile 平铺多图形窗口mmtext（ optional text） 在图上放置或拖曳文本mrnzoom 用橡皮框缩放坐标轴mmzap object使用鼠标删除文本，线型或坐标轴mmfont prop value修改文本字体属性例：clf;x=0:0.1:4;y=zeros(size(x);e=rand(size(x);yu=y+e;yd=y-e;error

24、bar(x,y,e)hold onplot(x,yu,r-);plot(x,yd,r-);yu(1)=0;yu(41)=0;fill(x,yu,r);yd(1)=0;yd(41)=0;fill(x,yd,g);  3. 三维图形(1) plot3(三维直线函数)以下例子用来体会plot 3的基本的绘图原理。例：绘参数方程 x=t;y=sin(t);z=cos(t) 的空间曲线clft=0:0.05:100;x=t;y=sin(t);z=sin(2*t);plot3(x,y,z,b:)例：空间划线：clft=0:0.1:10;x=t;y=0*ones(size(x);z=sin(t);plot3(x,y,z,r)hold onz=0*ones(size(x);y=sin(t);plot3(x,y,z