完整版matlab的一些画图技巧.docx

1、完整版matlab的一些画图技巧semilogx 和假设我们希望即向量y的下x和y是由计算matlab中如何在指定一点画一个填充颜色的小圆plot(1,1,r.,markersize,50)二维作图绘图命令plot绘制x-y坐标图；loglog命令绘制对数坐标图;semilogy命令绘制半对数坐标图；polor命令绘制极坐标图.基本形式如果y是一个向量，那么plot(y)绘制一个y中元素的线性图. 画出y=0., 0.48, 0.84, 1., 0.91, 6.14 则用命令：plot(y)它相当于命令：plot(x, y)，其中x=1,2,n或x=1;2;n,标编号,n为向量y的长度Matl

2、ab会产生一个图形窗口，显示如下图形，请注意：坐标 机自动绘出的.上面的图形没有加上x轴和y轴的标注，也没有标题.用xlabel, ylabel , title命令可以加上.如果x, y是同样长度的向量，plot(x,y)命令可画出相应的x元素与y元素的x-y坐标图.例：x=0:0.05:4*pi; y=s in (x); plot(x,y)grid on, title( y=sin( x )曲线图)xlabel( x = 0 : 0.05 : 4Pi )结果见下图.图4.1.1.2 y=sin(x)的图形title图形标题text标注数据点lege nd在右上角加解释文字grid 给图形加上

3、网格hold 保持图形窗口的图形表4.1.1.1 Matlab 图形命令多重线在一个单线图上，绘制多重线有三种办法 第一种方法是利用plot的多变量方式绘制:plot(x1,y1,x2,y2,.,x n,yn)x1,y1,x2,y2,.,xn,yn是成对的向量，每一对 x, y在图上产生如上方式的单 线多变量方式绘图是允许不同长度的向量显示在同一图形上.第二种方法也是利用plot绘制，但加上hold on/off命令的配合：plot(x1,y1)hold onplot(x2,y2)hold off第三种方法还是利用plot绘制，但代入矩阵:如果plot用于两个变量plot(x,y)，并且x,

4、y是矩阵，则有以下情况：(1) 如果y是矩阵，x是向量，plot(x,y)用不同的画线形式绘出y的行或列 及相应的x向量，y的行或列的方向与x向量元素的值选择是相同的.(2) 如果x是矩阵，y是向量，则除了 x向量的线族及相应的y向量外，以 上的规则也适用.(3) 如果x，y是同样大小的矩阵，plot(x,y)绘制x的列及y相应的列.还有其它一些情况，请参见 Matlab的帮助系统.线型和颜色的控制如果不指定划线方式和颜色，Matlab会自动为您选择点的表示方式及颜色.您也可以用不同的符号指定不同的曲线绘制方式.例如：plot(x,y,*) 用*作为点绘制的图形.plot(x1,y1,:,x2

5、,y2,+) 用:画第一条线，用+画第二条线.线型、点标记和颜色的取值有以下几种：线型点标记颜色实线.占y黄虚线o小圆圈m棕色点划线x叉子符c青色间断线+加号r红色*星号g绿色s方格b蓝色d菱形w白色A朝上三角k黑色v朝下三角朝右三角朝左三角p五角星h六角星表4.131线型和颜色控制符如果你的计算机系统不支持彩色显示，Matlab将把颜色符号解释为线型符号，用不同的线型表示不同的颜色. 颜色与线型也可以一起给出，即同时指定曲线的颜色和线型.例如： t=-3.14:0.2:3.14;x=si n( t); y=cos(t);plot(t,x, +r,t,y, -b)对数图、极坐标图及条形图log

6、log、semilogx、semilogy和polar的用法和plot相似.这些命令允许数据在不同的graph paper上绘制，例如不同的坐标系统.先介绍的 fplot是扩展来的可用于符号作图的函数.fplot(fname,lims)绘制fname指定的函数的图形.polar( theta, rho)使用相角theta为极坐标形式绘图，相应半径为 rho，其次可使用grid命令画出极坐标网格.loglog 用 log10-log10 标度绘图.semilogx用半对数坐标绘图，x轴是log10 , y是线性的.semilogy用半对数坐标绘图，y轴是log10，x是线性的.bar(x)显示x

7、向量元素的条形图，bar不接受多变量.hist绘制统计频率直方图.histfit(data ,n bi ns)绘制统计直方图与其正态分布拟合曲线.fplot函数的绘制区域为lims=xmin,xmax，也可以用Iims=xmin,xmax,ymin,ymax指定y轴的区域.函数表达式可以是一个函数名，如sin , tan等；也可以是带上参数x的函数表达式，如sin(x), diric(x,10)；也 可以是一个用方括号括起的函数组，如sin, cos.例 1: fplot(sin,0 4*pi)例2： fplot(sin(1 ./x), 0.01 0.1)例3: fplot(abs(exp(-

8、j*x*(0:9)*ones(10,1),0 2*pi,-o)例4： fplot(sin(x), cos(x) , tan(x),-2*pi 2*pi -2*pi 2*pi) %(图4.141)图4.1.4.1 sin,cos,tan 函数图形 图4.1.4.2半对数图下面介绍的是其它几个作图函数的应用. 例5:半对数坐标绘图t=0.001:0.002:20;y=5 + log(t) + t;semilogx(t,y, b)hold onsemilogx(t,t+5, r) %(图 4.142)例6:极坐标绘图t=0:0.01:2*pi;polar(t,si n( 6*t) % (图 4.14

9、3)图4.1.4.3极坐标绘图 图4.1.4.4正态分布的统计直方图与其正态分布拟合曲线例7:正态分布图我们可以用命令normrnd生成符合正态分布的随机数.n ormrnd(u,v,m, n)其中，u表示生成随机数的期望，v代表随机数的方差.运行：a=n ormrnd(10,2,10000,1);histfit(a) % (图 4.144)我们可以得到正态分布的统计直方图与其正态分布拟合曲线.例&比较正态分布(图4.145( 1)与平均分布(图4.145( 2)的分布图:yn=ra ndn(30000,1); % 正态分布x=min(yn) : 0.2 : max( yn);subplot(

10、121)hist( yn, x)yu=ra nd(30000,1); % 平均分布subplot(122)hist(yu, 25)4.145(2)4.145(1)在绘图过程中，经常要把几个图形在同一个图形窗口中表现出来，而不是简 单地叠加(例如上面的例8).这就用到函数subplot .其调用格式如下：subplot( m,n ,p)subplot函数把一个图形窗口分割成 m x n个子区域，用户可以通过参数 p调用个各子绘图区域进行操作.子绘图区域的编号为按行从左至右编号.例9:绘制子图x=0:0.1*pi:2*pi; subplot(2,2,1) plot(x,s in (x),-*);

11、title(si n(x); subplot(2,2,2) plot(x,cos(x),-o); title(cos(x); subplot(2,2,3) plot(x,s in (2*x),-.*); title(si n(2x); subplot(2,2,4); plot(x,cos(3*x),:d)title(cos(3x)得到图形如下:图4.1.5.1子图利用二维绘图函数patch，我们可绘制填充图绘制填充图的另一个函数为 fill.下面的例子绘出了函数humps （一个Matlab演示函数）在指定区域内的函 数图形.例10:用函数patch绘制填充图fplot(humps,0,2,b

12、) hold onpatch(0.5 0.5:0.02:1 1,0 humps(0.5:0.02:1) 0,r);hold offtitle(A region under an interesting function.)grid图4.1.6.1填充图我们还可以用函数fill来绘制类似的填充图. 例11：用函数fill绘制填充图x=0:pi/60:2*pi;y=si n( x);x1=0:pi/60:1;y1=s in( x1);plot(x,y,r);hold onfill(x1 1,y1 0,g)图4.162填充图三维作图mesh(Z)语句可以给出矩阵Z元素的三维消隐图，网络表面由Z坐标点

13、定义， 与前面叙述的x-y平面的线格相同，图形由邻近的点连接而成.它可用来显示用 其它方式难以输出的包含大量数据的大型矩阵，也可用来绘制 Z变量函数.显示两变量的函数Z=f(x,y)，第一步需产生特定的行和列的x-y矩阵.然后计 算函数在各网格点上的值.最后用 mesh函数输出.下面我们绘制sin(r)/r函数的图形.建立图形用以下方法：x=-8:.5:8;y=x；x=on es(size(y)*x;y=y* on es(size(y);R=sqrt(x.A2+y.A2)+eps;z=si n(R)./R;mesh(z) % 试运行 mesh(x,y,z)，看看与mesh(z)有什么不同之处？

14、(坐标问题)各语句的意义是：首先建立行向量 x，列向量y;然后按向量的长度建立1- 矩阵；用向量乘以产生的1-矩阵，生成网格矩阵，它们的值对应于x-y坐标平面; 接下来计算各网格点的半径；最后计算函数值矩阵 Z 用mesh函数即可以得到 图形.图4.2.1三维消隐图第一条语句x的赋值为定义域，在其上估计函数；第三条语句建立一个重复 行的x矩阵，第四条语句产生y的响应，第五条语句产生矩阵 R (其元素为各网 格点到原点的距离)用mesh方法结果如上.另外，上述命令系列中的前4行可用以下一条命令替代:x, y=meshgrid(-8:0.5:8)(1)meshc与函数mesh的调用方式相同，只是该

15、函数在 mesh的基础上又 增加了绘制相应等高线的功能下面来看一个 meshc的例子：x,y=meshgrid(-4:.5:4);z=sqrt(x.A2+y.A2);meshc(z) % 试运行 meshc(x,y,z)，看看与 meshc(z)有什么不同之处？我们可以得到图形：图 4.2.2.1 meshc 图地面上的圆圈就是上面图形的等高线.(2)函数meshz与mesh的调用方式也相同，不同的是该函数在 mesh函数的作用之上增加了屏蔽作用，即增加了边界面屏蔽例如：x,y=meshgrid(-4:.5:4);z=sqrt(x.A2+y.A2);meshz(z) % 试运行 meshz(x

16、,y,z)，看看与 meshz(z)有什么不同之处？我们得到图形：图 4.2.2.2 meshz 图(1)在Matlab中有一个专门绘制圆球体的函数 sphere，其调用格式如下:x,y,z=sphere( n)此函数生成三个(n+1) x (n+1)阶的矩阵，再利用函数surf(x,y,z)可生成单位球 面.x,y,z=sphere 此形式使用了默认值n=20sphere( n) 只绘制球面图，不返回值.运行下面程序：sphere(30);axis square;我们得到球体图形：图4.2.3.1 球面图若只输入sphere画图，则是默认了 n=20的情况.surf函数也是Matlab中常用

17、的三维绘图函数.其调用格式如下:surf(x,y, z,c)输入参数的设置与mesh相同，不同的是mesh函数绘制的是一网格图，而 surf绘制的是着色的三维表面.Matlab语言对表面进行着色的方法是，在得到 相应网格后，对每一网格依据该网格所代表的节点的色值(由变量 c控制)，来定义这一网格的颜色.若不输入 C,则默认为c=z.我们看下面的例子：%绘制地球表面的气温分布示意图.a,b,c=sphere(40);t=abs(c); %求绝对值surf(a,b,c,t);axis equalcolormap(hot)我们可以得到图形如下：图423.2等温线示意图(1)坐标轴的控制函数axis，

18、调用格式如下:axis(xmi n,xmax,ymi n,ymax,zmi n,zmax) 用此命令可以控制坐标轴的范围.与axis相关的几条常用命令还有：axis equalaxis squareaxis onaxis offaxis manual严格控制各坐标的分度使其相等使绘图区为正方形恢复对坐标轴的一切设置取消对坐标轴的一切设置以当前的坐标限制图形的绘制(2) grid on在图形中绘制坐标网格.grid off取消坐标网格.(3)xlabel, ylabel, zlabel分别为x轴,y轴,z轴添加标注.title为图形添加 标题.以上函数的调用格式大同小异，我们以 xlabel为例

19、进行介绍：xlabel(标注文本，属性1，属性值1，属性2，属性值2，)这里的属性是标注文本的属性，包括字体大小、字体名、字体粗细等.例如：x, y=meshgrid(-4:.2:4);R=sqrt(x.A2+y.A2);z=-cos(R);mesh(x,y,z)xlabel(xi n-4,4,fo ntweight,bold);ylabel(yi n-4,4,fo ntweight,bold);zlabel(z=-cos(sqrt(xA2+yA2),fo ntweight,bold);title(旋转曲面,fontsize,15,fontweight,bold,fontname,隶书);旅曲

20、曲面*44】 xtM4图424.1添加标注以上各种绘图方法的详细用法，请看联机信息.对平面上 n 个点. - L1 i，使得散点用微积分的在平面直线族 为实数中寻求一条直线到与散点相对应的在直线上的点之间的纵坐标的误差的平方和最小,方法可得：L為戸（艺禹）（乏。岸工：一0护所求得的这条直线称为回归直线.例：已知如下点列，求其回归直线，并计算最小误差平方和.x0.10.1.12.13.14.15.16.17.18.2.21.231y4243.4545.4547.495350555560555参考的程序如下：x=0.1 0.11 .12 .13 .14 .15 .16 .17 .18 .2 .21 .23; y=42 43.5 45 45.5 45 47.5 49 53 50 55 55 60; n=len gth(x);xb=mea n( x);yb=mea n( y);x2b=sum(x.A2)/n;xyb=x*y/n;b=(xb*yb-xyb)/(xbA2-x2b);a=yb-b*xb;y1=a+b.*x;plot(x,y,*,x,y1);serror=sum(y-y1).A2)