2-Matlab绘图.ppt

1、一、二维图形的基本命令 Funtool %函数作图工具; taylortool %taylor展开工具; plot(x) % x为向量，以x元素的值为纵坐 标，x的序号为横坐标绘制曲线 plot(x,y) % 以x元素为横坐标，y元素为纵坐标绘制曲线 plot(x,y1,x,y2,) % 以公共的x元素为横坐标，y1,y2，.为纵坐标绘制多条曲线,二维平面图形,例1 画出一条正弦曲线和余弦曲线。 x=0:pi/10:2*pi; y1=sin(x); y2=cos(x); plot(x,y1, x,y2),plot(x,y,颜色+线型+点型,) 颜色代码： 黄色（y），紫色（m），青色（c），

2、红色（r），绿色（g），蓝色（b）， 白色（w），黑色（k） 点型和线型代码： . 。 + * s(正方形)，d(菱形)，h(六角星)，p(五角星) - ： -. - 实线 点线 点划线 虚线,例2 绘制不同颜色的曲线。-zhao31 x=0:pi/10:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,b-,x,y2,r:),命令 subplot(m,n,k) 实现一个区域中显示mn个子图形窗口，并指定在第k个子窗口绘图。,在同一窗口绘制多图 -zhao32,x=0:pi/10:2*pi; y1=sin(x); y2=cos(x); y3=x.2+6*x+3; y4

3、=exp(x); subplot(2,2,1); plot(x,y1,bo-) subplot(2,2,2); plot(x,y2,R*:) subplot(2,2,3); plot(x,y3,g+) subplot(2,2,4); plot(x,y4,mp),hold on 实现在同一窗口中多次绘制曲线， hold off取消。 例3 -zhao33 x=0:pi/10:2*pi; y1=sin(x); y2=cos(x); y3=log(x); plot(x,y1,bo-,x,y2,r*:) hold on plot(x,y3,md-) hold off,二、图形的标识与修饰,grid 命

4、令对图形窗口加坐标网格 axis 实现坐标轴的重新设置 axis(xmin xmax ymin ymax) linewidth和markersize实现图形中线宽和点型大小的设置，默认值是1。 plot(x,y,b*-,linewidth,4,markersize,6) 4 使用title，xlabel，ylabel，zlabel，text等命令实现对图形的文字说明.,例4 二维平面图形-zhao34,x=0:pi/10:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,bo-,x,y2,R*:) grid %加注坐标网格 xlabel(x) %标记横坐标轴 ylab

5、el(y) title(Sine and Cosine Curve) %标记图名 % gtext(y=sinx) %在选定的地方标记曲线名 % gtext(y=cosx) text(2.5,0.7,y=sin(x) %在(2.5,0.7)处标记曲线名 text(1.5,0.1,y=cos(x),三 其他图形函数,fill(x,y,颜色) %对闭合图形进行填充 polar(theta,rho) %极坐标绘图 stairs(x,y) %绘制阶梯图 pie(x) %饼图,向量中各元素所占的比重 area(x,y) %面积图; bar(x,y) %绘制二维直方图 bar3(x,y) %绘制三维直方图,

6、其他图形的显示 -zhao35,例5 x=0:pi/10:2*pi; y=sin(x); subplot(2,2,1); plot(x,y); subplot(2,2,2); fill(x,y,g); subplot(2,2,3); bar(x,y); subplot(2,2,4); stairs(x,y);,例6 某大班高等数学考试，90分以上15人，80分至89分32人，70分至79分58人，60分至 69分20人，60分以下10人，画出饼图。 -zhao36,x=15,32,58,20,10; pie(x); explode=0,0,0,0,1;%让不及格部分脱离饼图 pie(x,exp

7、lode) %向量explode和x元素相同，explode中不为零的部分会被分开,例7 画面积图。 x=-2:2; y=-1,3;2,5;3,6;8,2;2,4; %x的元素个数必须与y的行数相同。 area(x,y),例8 画出极坐标方程 R=sin(2*theta)cos(2*theta)的图形。 theta=linspace(0,2*pi); rho=sin(2*theta).*cos(2*theta); polar(theta,rho,r); title(polar plot of sin(2*theta)cos(2*theta),三维图形的绘制,三维曲线的绘制 plot3(x,y,

8、z,颜色+线型+点型) 例9 画出螺旋线 t=0:0.1:50; x=sin(t); y=cos(t); z=t; plot3(x,y,z,m*) grid %加注坐标网格,三维曲面图的绘制,mesh(x,y,z) %画网格曲面，已知点和附近的点以直线连接 surf(x,y,z) %画光滑曲面，已知点和附近点以平面连接 例10 绘制高斯矩阵的曲面.-zhao37 z=peaks(40); subplot(1,2,1); mesh(z); subplot(1,2,2); surf(z);,还有一个重要函数 X,Y=meshgrid(x,y) 其中x,y是向量，通过上述函数就可将x和y指定的区域转

9、换成为矩阵X和Y。,例11 -zhao38 x=-7.5:0.5:7.5; y=x; x,y=meshgrid(x,y); %生成矩形域上网格节点矩阵 z=sin(sqrt(x.2+y.2)./(sqrt(x.2+y.2)+eps); %生成函数值矩阵 subplot(2,1,1) %画子图1 surf(x,y,z) %画曲面 subplot(2,1,2) %画子图2 surf(x,y,z) %画光滑曲面 shading flat %对曲面平滑并除去网格线,二、设置观察点 view(方位角,仰角) 其中方位角是观察点和坐标原点连线在x-y平面内的投影和y轴负方向的夹角, 仰角是观察点和坐标原点

10、连线与x-y平面的夹角。 对于这两个角度，三维图形的默认值是-37.5和30，二维图形的默认值是0和90。,例12 -zhao39 z=peaks(40); subplot(2,2,1); mesh(z); subplot(2,2,2); surf(z); view(-37.5,-30); subplot(2,2,3); surf(z); view(180,0); subplot(2,2,4); mesh(z); view(0,90);,例13 绘图表示实验数据,-zhao310 t=1:9; d1=12.51,13.54,15.60,15.92,20.04,24.53,30.24,50.12

11、,76.98; d2=9.87,20.11,32.56,40.23,57.77,69.13,72.24,82.50,89.88; d3=10.78,8.54,14.50,45.21,39.04,60.53,70.24,50.12,36.98; plot(t,d1,R-,t,d2,gx:,t,d3,m*-) title(Time axis(0,10,0,100) text(6.5,25.5,leftarrowdata1) text(4,65,data2rightarrow) text(3.5,25,leftarrowdata3),三、 符号表达式绘图,1、fplot(fun,lims) %绘制数

12、学函数，fun表示函数名，也可以是定义的m文件名，也可以是以x为变量的字符串。 Lims=xmin xmax ymin ymax限定x、y轴上的绘图空间。 例 fplot(sin(1./x),0.001 0.1),2、ezplot(f,xmin,xmax) %直接画出函数图形， xmin,xmax缺省值为-2pi,2pi 例 syms x; y=sin(x); ezplot(y) 或者 syms x; ezplot(sin(x),0,2*pi) 3、ezsurf(f) %画符号函数f(x,y)的曲面图， ezsurf(sin(x2+y2)/(x2+y2),四、特殊图形,peaks函数 x,y,

13、z=peaks; surf(x,y,z) sphere(N) % 产生(N+1)* (N+1)矩阵，绘制单位球面 x,y,z=sphere; %默认值N=20 surf(x,y,z) cylinder(r,N) % 产生(N+1)* (N+1)矩阵，绘制半径为r的 圆柱体 x,y,z=cylinder; %默认r=1，N=20 surf(x,y,z),例14 绘制高度为10的圆柱面x2+y2=4 -zhao311 r=2; n=1000; xb,yb,zb=cylinder(r,n); mesh(xb,yb,10*zb),x=-2:0.01:2; y=x; xx,yy=meshgrid(x,y); zz=xx.2+yy.2; i=find(zz1); %找出圆域外的坐标点i zz(i)=NaN; %对圆域外面的坐标点i处函数值赋空 surf(xx,yy,zz