Matlab图形绘制经典案例

合集下载

Matlab绘图教程(大量实例)

例4 用不同标度在同一坐标内绘制曲线 y1=0.2e-0.5xcos(4πx) 和y2=2e-0.5xcos(πx)。程序如下： x=0:pi/100:2*pi; y1=0.2*exp(-0.5*x).*cos(4*pi*x); y2=2*exp(-0.5*x).*cos(pi*x); plotyy(x,y1,x,y2);
三维图形
1 三维曲线 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对应列元素绘制三维曲线，曲线条数等于矩阵列数。
图形的可视化编辑
MATLAB 6.5版本在图形窗口中提供了可视化的图形编辑工具，利用图形窗口菜单栏或工具栏中的有关命令可以完成对窗口中各种图形对象的编辑处理。
在图形窗口上有一个菜单栏和工具栏。菜单栏包含File、Edit、View、Insert、Tools、Window和 Help共7个菜单项，工具栏包含11个命令按钮。
MATLAB绘图
二维数据曲线图
plot函数的基本调用格式为： plot(x,y) 其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。例1 在0≤x≤2区间内，绘制曲线 y=2e-0.5xcos(4πx) 程序如下： x=0:pi/100:2*pi; y=2*exp(-0.5*x).*cos(4*pi*x); plot(x,y)
例3 分析下列程序绘制的曲线。 x1=linspace(0,2*pi,100); x2=linspace(0,3*pi,100); x3=linspace(0,4*pi,100); y1=sin(x1); y2=1+sin(x2); y3=2+sin(x3); x=[x1;x2;x3]'; y=[y1;y2;y3]'; plot(x,y,x1,y1-1)

实验二MATLAB绘制图形

grid on %在所画出的图形坐标中加入栅格
绘制图形如下
50
10
1
0.8
40
10
0.6
0.4
30
10
0.2
0
1020
-0.2
-0.4
1010
-0.6
-0.8
0
10
-1
-2
0
2
-2
0
2
10
10
10
10
10
10
如果在图中不加栅格
程序如下：
clear x=logspace(-1,2);%在10^(-1)到10^2之间产生50个对数等分的行向量 subplot(121); loglog(x,10*exp(x),'-p') subplot(122); semilogx(x,cos(10.^x))
（2）plot(x,y)：基本格式，x和y可为向量或矩阵. 1. 如果x，y是同维向量，以x元素为横坐标，以y元素为纵坐标绘图. 2. 如果x是向量，y是有一维与x元素数量相等的矩阵，则以x为共同横坐标， y元素为纵坐标绘图，曲线数目为y的另一维数. 3. 如果x，y是同维矩阵，则按列以x，y对应列元素为横、纵坐标绘图，曲线数目等于矩阵列数.
y=2*exp(-0.5*x).*cos(4*pi*x);
2
plot(x,y)
1.5
1
0.5
0
-0.5
-1
-1.5
-2
0
1
2
3
4
5
6
7
例4 绘制曲线
t=(0:0.1:2*pi);
x=t.*sin(3*t);
y=t.*sin(t).*sin(t);

用matlab绘制的漂亮图形

用matlab绘制的漂亮图形1.不同坐标系下的图形对比theta=0:pi/20:4*pi;phi= theta.^2- theta;[t,p]=meshgrid(theta,phi);r=t.*p;subplot(1,2,1);mesh(t,p,r);ylabel('x');xlabel('y');zlabel('z');[x,y,z]=sph2cart(t,p,r);subplot(1,2,2);mesh(x,y,z);ylabel('x');xlabel('y');zlabel('z');2.球曲面的法线[x,y,z]=sphere;Surfnorm(x,y,z)3.x=rand(100,1)*16-8;y=rand(100,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),33);ylin=linspace(min(y),max(y),33); [X,Y]= meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y);mesh(X,Y,Z);axis tight;hold on;ylabel('x');xlabel('y');zlabel('z'); plot3(x,y,z,’r’,’MarkerSize’,15)x=rand(1000,1)*16-8;y=rand(1000,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),99);ylin=linspace(min(y),max(y),99); [X,Y]= meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y);mesh(X,Y,Z);axis tight;hold on;ylabel('x');xlabel('y');zlabel('z'); plot3(x,y,z,'r','MarkerSize',30);去掉网格后：4.圆通fnplt(fncmb(rsmak(‘cylinder’,1,2),[0 0 -1;0 1 0;1 0 0])); axis equal,axis off;shading interp5.c=8.01*10^(-6);T=9.8*0.940*10^(-7);for m=1:100a=m/10000000;for n=1:100b=n/10000000;omiga(m,n)=sqrt((10.6667*a^3*b^5+10.6667*a^5*b^3)/(2.3112*a^5*b^5*c)* T)/(2*pi);endendmesh((1:100)/10000000,(1:100)/10000000,omiga);ylabel('半短轴b');xlabel('半长轴a');zlabel('频率f');6.马鞍面x=-8:8; y=-8:8;[X,Y]= meshgrid(x,y);Z=(X.^2./4.^2-Y.^2./5.^2);meshz(X,Y,Z);ylabel('x');xlabel('y');zlabel('z');x=-8:8;y=-8:8;[X,Y]= meshgrid(x,y);Z=X.^2./4.^2-Y.^2./5.^2;surfc(X,Y,Z);ylabel('x');xlabel('y');zlabel('z');7.螺旋线t=0:pi/60:10*pi;x=sin(t);y=cos(t);plot3(x,y,t,'*-b');ylabel('x');xlabel('y');zlabel('z');t=0:pi/60:10*pi;x=sin(t);y=cos(t);fill3(x,y,t,'*-b');ylabel('x');xlabel('y');zlabel('z');8.[x,y]= meshgrid([-2:0.2:2]); f=x.^2+y.^2;[fx,fy]= gradient(f,0.2,0.2); quiver(x,y,fx,fy);hold oncontour(x,y,f);9.[x,y]= meshgrid([-4:0.1:4]); z=x.*exp(-x.^2-y.^2);plot3(x,y,z);10. [x,y]= meshgrid([-2:0.05:2]); z=x.*exp(-x.^2-y.^2);surf(x,y,z);11.牛顿环（动态的）clearR=1;N=400;lamda=200e-9;rr=0.0017; [x,y]= meshgrid(linspace(-rr,rr,N)); r=sqrt(x.^2+y.^2);phi=r.^2/R./lamda.*pi.*2+pi;I=4.*cos(phi./2).^2;H=imshow(I);t=0;k=1;while k;s=get(gcf,'currentkey');if strcmp(s,'space');clc;k=0;endt=t+0.01;pause(0.2);phi=phi+t;I=4.*cos(phi./2).^2;set(H,'CData',I);end12.衍射调制下的双孔干涉条纹clear alllambda=500e-9;D=1;d=1e-6;a=3.5*d;ymax=10*lambda*D/a;xmax=ymax;N=101;xs=linspace(-xmax,xmax,N);ys=linspace(-ymax,ymax,N);for i=1:N;for j=1:N;r1(i,j)=sqrt(D.^2+(ys(i)-a./2).^2+xs(i).^2);r2(i,j)=sqrt(D.^2+(ys(j)-a./2).^2+xs(j).^2);Dr=r2-r1;phi(i,j)=Dr(i,j)./lambda.*pi.*2;m(i,j)=pi.*d.*sin(phi(i,j))./lambda;I(i,j)=(4.*cos(phi(i,j)./2).^2).*real((2.*besselj(1,m(i,j))./m(i,j)). ^2);endendNC=225;Ir=I*NC;subplot(1,2,1);image(xs,ys,Ir);axis equalcolormap(gray(NC));subplot(1,2,2);plot(I,m),hold onlambda=600e-9;D=1;d=1e-6;a=2.5*d;ymax=10*lambda*D/a;xmax=ymax;N=101;xs=linspace(-xmax,xmax,N);ys=linspace(-ymax,ymax,N);for i=1:N;for j=1:N;r1(i,j)=sqrt(D.^2+(ys(i)-a./2).^2+xs(i).^2);r2(i,j)=sqrt(D.^2+(ys(j)-a./2).^2+xs(j).^2);Dr=r2-r1;phi(i,j)=Dr(i,j)./lambda.*pi.*2;m(i,j)=pi.*d.*sin(phi(i,j))./lambda;I(i,j)=(4.*cos(phi(i,j)./2).^2).*real((2.*besselj(1,m(i,j))./m(i,j)). ^2);endendNC=225;Ir=I*NC;colormap(gray(NC));plot(m,I),hold on13.lambda=500e-9;D=1;d=1e-6;a=3.5*d;ymax=10*lambda*D/a;xmax=ymax;N=101;xs=linspace(-xmax,xmax,N);ys=linspace(-10000*ymax,10000*ymax,N);for i=1:N;for j=1:N;r1(i,j)=sqrt(D.^2+(ys(i)-a./2).^2+xs(i).^2);r2(i,j)=sqrt(D.^2+(ys(j)-a./2).^2+xs(j).^2);Dr=r2-r1;phi(i,j)=Dr(i,j)./lambda.*pi.*2;I(i,j)=real(4.*cos(phi(i,j)./2).^2);endendNC=225;Ir=I/3.*NC;image(xs,ys,Ir);colormap(gray(NC));hold on14.沙丘figure( 'Position',[200 31 700 696],'color','w'); hold on ;axis equal ; axis offtitle('沙丘','fontsize',16,'color','k');text(180,480,'——万祎龙 ','fontsize',12,'color','k');i=1;for py=1:2:360c=py*4*pi/360;b=pi*cos(c);a=0:pi/30:6*pi;x=100*a/pi+20;y=10*sin(a+cos(a)*pi-b)*cos(c)+py;plot(x,y,'color','k');Z(i,:)=sin(a+cos(a)*pi-b)*cos(c);i=i+1;end15.%存脚本文件heart.m%%一个半径为R的圆周等分成N份，然后一等分点为圆心，以等分点到该圆垂直直径%%的距离为半径画圆。

使用matlab绘制三维图形的方法

使用matlab绘制三维图形的方法要使用MATLAB绘制三维图形，首先需要了解MATLAB中的三维绘图函数和绘图选项。

下面将介绍一些常用的绘制三维图形的方法。

1.绘制基本的三维图形要绘制基本的三维图形，可以使用以下函数：- plot3(函数：用于在三维坐标系中绘制线条。

- scatter3(函数：用于在三维坐标系中绘制散点图。

- surf(函数：用于绘制三维曲面图。

- mesh(函数：用于绘制三维网格图。

- bar3(函数：用于绘制三维条形图。

- contour3(函数：用于绘制三维等高线图。

例如，下面的代码演示了如何使用plot3(函数绘制一个三维线条图：```x = linspace(0, 2*pi, 100);y = sin(x);z = cos(x);plot3(x, y, z, 'LineWidth', 2);xlabel('X');ylabel('Y');zlabel('Z');title('3D Line Plot');```2.添加颜色和纹理在绘制三维图形时，可以使用颜色和纹理来增加图形的信息。

MATLAB 提供了一系列函数来处理颜色和纹理，如：- colormap(函数：用于设置颜色映射。

- caxis(函数：用于设置坐标轴范围。

- shading(函数：用于设置颜色插值方法。

- texturemap(函数：用于设置纹理映射方法。

例如，下面的代码展示了如何使用纹理映射来绘制一个球体：```[X, Y, Z] = sphere(50);C = colormap('jet');surface(X, Y, Z, 'FaceColor', 'texturemap', 'CData', C);axis equal;```3.绘制多个数据集要在同一张图中绘制多个数据集，可以使用hold on和hold off命令。

案例28使用matlab绘制专业图形——绘制曲线和曲面图

案例28使用MATLAB 绘制专业图形绘制曲线和曲面图案例主要信息提示 ● 案例内容：使用MATLAB 进行专业作图。

● 关键词：MATLAB ，作图，二维作图，三维作图 ● 建议课时：2课时● 适合专业：理工科各专业●光盘中的数字资源：二维曲线绘图M 文件graph_line.m三维曲面绘图M 文件graph_surface.m一、实验内容实验内容1.绘制曲线图在1个绘图窗口的4个不同子窗口中分别绘制以下曲线： ①在子窗口1中绘制两条二维曲线，分别为：0.510.2cos(4)x y e x π-=和0.522cos()x y e x π-=并显示网格线。

②在子窗口2中绘制两条二维曲线，分别为：0.512x y e -=和2cos(2)y x π=并添加标题、x/y 轴名称、图形说明和图例。

③在子窗口3中绘制一条极坐标曲线，如下：sin()cos()r t t =④在子窗口4中绘制一条极三维螺旋线。

全部4个子窗口绘制效果如图28-1所示。

图28-1曲线绘制效果图实验内容2.绘制曲面图在1个绘图窗口的4个不同子窗口中分别绘制以下曲面：①在子窗口1和子窗口2中以不同着色方式绘制两个球面。

②在子窗口3和子窗口4中以不同视角绘制两个直径相等的圆管相交的图形。

全部4个子窗口绘制效果如图28-2所示。

图28-1曲面绘制效果图二、预备知识2.1 MATLAB绘图功能简介MA TLAB软件提供了丰富的图形表达功能，包括常用的二维图形和三维图形。

其中，各种二维图形近30种，三维图形20多种。

应用MA TLAB除了能作一般的曲线、散点图、条形图等，还能绘制流线图、三维矢量图等工程实用图形。

下面我们介绍一些基本的二维和三维绘图函数。

2.2 二维图形的绘制二维图形是将平面坐标上的数据点连接起来的平面图形。

可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。

二维图形的绘制是其他绘图操作的基础。

2.2.1绘制二维曲线的plot函数在MATLAB中，plot函数是最简单、最基本而且应用最为广泛的线性绘图函数，利用它可以在二维平面上生成线段、曲线和参数方程曲线等的函数图形。

matlab编程实例100例

1-32是：图形应用篇33-66是：界面设计篇67-84是：图形处理篇85-100是：数值分析篇实例2：三角函数曲线（2）function shili02h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例02');x=-pi:0.05:pi;y=sin(x)+cos(x);plot(x,y,'-*r','linewidth',1);grid onxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例3：图形的叠加function shili03h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例03');x=-pi:0.05:pi;y1=sin(x);y2=cos(x);plot(x,y1,...'-*r',...x,y2,...'--og');grid onxlabel('自变量X');ylabel('函数值Y');title('三角函数');实例4：双y轴图形的绘制function shili04h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例04');x=0:900;a=1000;b=0.005;y1=2*x;y2=cos(b*x);[haxes,hline1,hline2]=plotyy(x,y1,x,y2,'semilogy','plot'); axes(haxes(1))ylabel('semilog plot');axes(haxes(2))ylabel('linear plot');实例6：图形标注function shili06h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例06');t=0:pi/10:2*pi;h=plot(t,sin(t));xlabel('t=0到2\pi','fontsize',16);ylabel('sin(t)','fontsize',16);title('\it{从0to2\pi 的正弦曲线}','fontsize',16)x=get(h,'xdata');y=get(h,'ydata');imin=find(min(y)==y);imax=find(max(y)==y);text(x(imin),y(imin),...['\leftarrow最小值=',num2str(y(imin))],...'fontsize',16)text(x(imax),y(imax),...['\leftarrow最大值=',num2str(y(imax))],...'fontsize',16)实例7：条形图形function shili07h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例07');tiao1=[562 548 224 545 41 445 745 512];tiao2=[47 48 57 58 54 52 65 48];t=0:7;bar(t,tiao1)xlabel('X轴');ylabel('TIAO1值');h1=gca;h2=axes('position',get(h1,'position'));plot(t,tiao2,'linewidth',3)set(h2,'yaxislocation','right','color','none','xticklabel',[]) 实例8：区域图形function shili08h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例08');x=91:95;profits1=[88 75 84 93 77];profits2=[51 64 54 56 68];profits3=[42 54 34 25 24];profits4=[26 38 18 15 4];area(x,profits1,'facecolor',[0.5 0.9 0.6],...'edgecolor','b',...'linewidth',3)hold onarea(x,profits2,'facecolor',[0.9 0.85 0.7],...'edgecolor','y',...'linewidth',3)hold onarea(x,profits3,'facecolor',[0.3 0.6 0.7],...'edgecolor','r',...'linewidth',3)hold onarea(x,profits4,'facecolor',[0.6 0.5 0.9],...'edgecolor','m',...'linewidth',3)hold offset(gca,'xtick',[91:95])set(gca,'layer','top')gtext('\leftarrow第一季度销量')gtext('\leftarrow第二季度销量')gtext('\leftarrow第三季度销量')gtext('\leftarrow第四季度销量')xlabel('年','fontsize',16);ylabel('销售量','fontsize',16);实例9：饼图的绘制function shili09h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例09');t=[54 21 35;68 54 35;45 25 12;48 68 45;68 54 69];x=sum(t);h=pie(x);textobjs=findobj(h,'type','text');str1=get(textobjs,{'string'});val1=get(textobjs,{'extent'});oldext=cat(1,val1{:});names={'商品一：';'商品二：';'商品三：'};str2=strcat(names,str1);set(textobjs,{'string'},str2)val2=get(textobjs,{'extent'});newext=cat(1,val2{:});offset=sign(oldext(:,1)).*(newext(:,3)-oldext(:,3))/2; pos=get(textobjs,{'position'});textpos=cat(1,pos{:});textpos(:,1)=textpos(:,1)+offset;set(textobjs,{'position'},num2cell(textpos,[3,2]))实例10：阶梯图function shili10h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例10');a=0.01;b=0.5;t=0:10;f=exp(-a*t).*sin(b*t);stairs(t,f)hold onplot(t,f,':*')hold offglabel='函数e^{-(\alpha*t)}sin\beta*t的阶梯图'; gtext(glabel,'fontsize',16)xlabel('t=0:10','fontsize',16)axis([0 10 -1.2 1.2])实例11：枝干图function shili11h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例11');x=0:pi/20:2*pi;y1=sin(x);y2=cos(x);h1=stem(x,y1+y2);hold onh2=plot(x,y1,'^r',x,y2,'*g');hold offh3=[h1(1);h2];legend(h3,'y1+y2','y1=sin(x)','y2=cos(x)') xlabel('自变量X');ylabel('函数值Y');title('正弦函数与余弦函数的线性组合'); 实例12：罗盘图function shili12h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例12');winddirection=[54 24 65 84256 12 235 62125 324 34 254];windpower=[2 5 5 36 8 12 76 14 10 8];rdirection=winddirection*pi/180;[x,y]=pol2cart(rdirection,windpower); compass(x,y);desc={'风向和风力','北京气象台','10月1日0：00到','10月1日12：00'};gtext(desc)实例13：轮廓图function shili13h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例13');[th,r]=meshgrid((0:10:360)*pi/180,0:0.05:1); [x,y]=pol2cart(th,r);z=x+i*y;f=(z.^4-1).^(0.25);contour(x,y,abs(f),20)axis equalxlabel('实部','fontsize',16);ylabel('虚部','fontsize',16);h=polar([0 2*pi],[0 1]);delete(h)hold oncontour(x,y,abs(f),20)实例14：交互式图形function shili14h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例14');axis([0 10 0 10]);hold onx=[];y=[];n=0;disp('单击鼠标左键点取需要的点'); disp('单击鼠标右键点取最后一个点'); but=1;while but==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('单击鼠标左键点取下一个点');x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');hold off实例14：交互式图形function shili14h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例14');axis([0 10 0 10]);hold onx=[];y=[];n=0;disp('单击鼠标左键点取需要的点'); disp('单击鼠标右键点取最后一个点'); but=1;while but==1[xi,yi,but]=ginput(1);plot(xi,yi,'bo')n=n+1;disp('单击鼠标左键点取下一个点');x(n,1)=xi;y(n,1)=yi;endt=1:n;ts=1:0.1:n;xs=spline(t,x,ts);ys=spline(t,y,ts);plot(xs,ys,'r-');hold off实例15：变换的傅立叶函数曲线function shili15h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren')h=uicontrol('style','slider','position',...[100 10 500 20],'min',1,'max',20)for j=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0 0 1 1]);movie(m,30)实例16：劳伦兹非线形方程的无序活动function shili15h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例15');axis equalm=moviein(20,gcf);set(gca,'nextplot','replacechildren') h=uicontrol('style','slider','position',...[100 10 500 20],'min',1,'max',20) for j=1:20plot(fft(eye(j+16)))set(h,'value',j)m(:,j)=getframe(gcf);endclf;axes('position',[0 0 1 1]);movie(m,30)实例17：填充图function shili17h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例17');t=(1:2:15)*pi/8;x=sin(t);y=cos(t);fill(x,y,'r')axis square offtext(0,0,'STOP',...'color',[1 1 1],...'fontsize',50,...'horizontalalignment','center') 例18：条形图和阶梯形图function shili18h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例18');subplot(2,2,1)x=-3:0.2:3;y=exp(-x.*x);bar(x,y)title('2-D Bar Chart')subplot(2,2,2)x=-3:0.2:3;y=exp(-x.*x);bar3(x,y,'r')title('3-D Bar Chart')subplot(2,2,3)x=-3:0.2:3;y=exp(-x.*x);stairs(x,y)title('Stair Chart')subplot(2,2,4)x=-3:0.2:3;y=exp(-x.*x);barh(x,y)title('Horizontal Bar Chart')实例19：三维曲线图function shili19h0=figure('toolbar','none',...'position',[200 150 450 400],...'name','实例19');subplot(2,1,1)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,y1,z1,x,y2,z2,x,y3,z3)grid onxlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure1:3-D Plot')subplot(2,1,2)x=linspace(0,2*pi);y1=sin(x);y2=cos(x);y3=sin(x)+cos(x);z1=zeros(size(x));z2=0.5*z1;z3=z1;plot3(x,z1,y1,x,z2,y2,x,z3,y3)grid onxlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure2:3-D Plot')实例20：图形的隐藏属性function shili20h0=figure('toolbar','none',...'position',[200 150 450 300],...'name','实例20');subplot(1,2,1)[x,y,z]=sphere(10);mesh(x,y,z)axis offtitle('Figure1:Opaque')hidden onsubplot(1,2,2)[x,y,z]=sphere(10);mesh(x,y,z)axis offtitle('Figure2:Transparent')hidden off实例21PEAKS函数曲线function shili21h0=figure('toolbar','none',...'position',[200 100 450 450],...'name','实例21');[x,y,z]=peaks(30);subplot(2,1,1)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfc(x,y,z)xlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure1:surfc函数形成的曲面')subplot(2,1,2)x=x(1,:);y=y(:,1);i=find(y>0.8&y<1.2);j=find(x>-0.6&x<0.5);z(i,j)=nan*z(i,j);surfl(x,y,z)xlabel('X轴');ylabel('Y轴');zlabel('Z轴');title('Figure2:surfl函数形成的曲面') 实例22：片状图function shili22h0=figure('toolbar','none',...'position',[200 150 550 350],...'name','实例22');subplot(1,2,1)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trimesh(t,x,y,z)hidden offtitle('Figure1:Triangular Surface Plot');subplot(1,2,2)x=rand(1,20);y=rand(1,20);z=peaks(x,y*pi);t=delaunay(x,y);trisurf(t,x,y,z)title('Figure1:Triangular Surface Plot'); 实例23：视角的调整function shili23h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例23');x=-5:0.5:5;[x,y]=meshgrid(x);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;subplot(2,2,1)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure1')view(-37.5,30)subplot(2,2,2)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure2')view(-37.5+90,30)subplot(2,2,3)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure3')view(-37.5,60)subplot(2,2,4)surf(x,y,z)xlabel('X-axis')ylabel('Y-axis')zlabel('Z-axis')title('Figure4')view(180,0)实例24：向量场的绘制function shili24h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例24');subplot(2,2,1)z=peaks;ribbon(z)title('Figure1')subplot(2,2,2)[x,y,z]=peaks(15);[dx,dy]=gradient(z,0.5,0.5);contour(x,y,z,10)hold onquiver(x,y,dx,dy)hold offtitle('Figure2')subplot(2,2,3)[x,y,z]=peaks(15);[nx,ny,nz]=surfnorm(x,y,z);surf(x,y,z)hold onquiver3(x,y,z,nx,ny,nz)hold offtitle('Figure3')subplot(2,2,4)x=rand(3,5);y=rand(3,5);z=rand(3,5);c=rand(3,5);fill3(x,y,z,c)grid ontitle('Figure4')实例25：灯光定位function shili25h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例25');vert=[1 1 1;1 2 1;2 2 1;2 1 1;1 1 2;12 2;2 2 2;2 1 2];fac=[1 2 3 4;2 6 7 3;4 3 7 8;15 8 4;1 2 6 5;5 6 7 8];grid offsphere(36)h=findobj('type','surface');set(h,'facelighting','phong',...'facecolor',...'interp',...'edgecolor',[0.4 0.4 0.4],...'backfacelighting',...hold onpatch('faces',fac,'vertices',vert,...'facecolor','y');light('position',[1 3 2]);light('position',[-3 -1 3]);material shinyaxis vis3d offhold off实例26：柱状图function shili26h0=figure('toolbar','none',...'position',[200 50 450 450],...'name','实例26');subplot(2,1,1)x=[5 2 18 7 39 8 65 5 54 3 2];bar(x)xlabel('X轴');ylabel('Y轴');title('第一子图');subplot(2,1,2)y=[5 2 18 7 39 8 65 5 54 3 2];barh(y)xlabel('X轴');ylabel('Y轴');title('第二子图');实例27：设置照明方式function shili27h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例27');subplot(2,2,1)spherecamlight leftcamlight rightlighting flatcolorbaraxis offtitle('Figure1')subplot(2,2,2)sphereshading flatcamlight leftcamlight rightlighting gouraudcolorbaraxis offtitle('Figure2')subplot(2,2,3)sphereshading interpcamlight rightcamlight leftlighting phongcolorbaraxis offtitle('Figure3')subplot(2,2,4)sphereshading flatcamlight leftcamlight rightlighting nonecolorbaraxis offtitle('Figure4')实例28：羽状图function shili28h0=figure('toolbar','none',...'position',[200 150 450 350],...'name','实例28');subplot(2,1,1)alpha=90:-10:0;r=ones(size(alpha));m=alpha*pi/180;n=r*10;[u,v]=pol2cart(m,n);feather(u,v)title('羽状图')axis([0 20 0 10])subplot(2,1,2)t=0:0.5:10;x=0.05+i;y=exp(-x*t);feather(y)title('复数矩阵的羽状图')实例29：立体透视（1）function shili29h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例29');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2);grid onfor i=-2:0.5:2;h1=surf(linspace(-2,2,20),...linspace(-2,2,20),...zeros(20)+i);rotate(h1,[1 -1 1],30)dx=get(h1,'xdata');dy=get(h1,'ydata');dz=get(h1,'zdata');delete(h1)slice(x,y,z,v,[-2 2],2,-2)hold onslice(x,y,z,v,dx,dy,dz)hold offaxis tightview(-5,10)drawnowend实例30：立体透视（2）function shili30h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例30');[x,y,z]=meshgrid(-2:0.1:2,...-2:0.1:2,...-2:0.1:2);v=x.*exp(-x.^2-y.^2-z.^2);[dx,dy,dz]=cylinder;slice(x,y,z,v,[-2 2],2,-2)for i=-2:0.2:2h=surface(dx+i,dy,dz);rotate(h,[1 0 0],90)xp=get(h,'xdata');yp=get(h,'ydata');zp=get(h,'zdata');delete(h)hold onhs=slice(x,y,z,v,xp,yp,zp);axis tightxlim([-3 3])view(-10,35)drawnowdelete(hs)hold offend实例31：表面图形function shili31h0=figure('toolbar','none',...'position',[200 150 550 250],...'name','实例31');subplot(1,2,1)x=rand(100,1)*16-8;y=rand(100,1)*16-8;r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;xlin=linspace(min(x),max(x),33); ylin=linspace(min(y),max(y),33); [X,Y]=meshgrid(xlin,ylin);Z=griddata(x,y,z,X,Y,'cubic');mesh(X,Y,Z)axis tighthold onplot3(x,y,z,'.','Markersize',20)subplot(1,2,2)k=5;n=2^k-1;theta=pi*(-n:2:n)/n;phi=(pi/2)*(-n:2:n)'/n;X=cos(phi)*cos(theta);Y=cos(phi)*sin(theta);Z=sin(phi)*ones(size(theta)); colormap([0 0 0;1 1 1])C=hadamard(2^k);surf(X,Y,Z,C)axis square实例32：沿曲线移动的小球h0=figure('toolbar','none',...'position',[198****8468],...'name','实例32');h1=axes('parent',h0,...'position',[0.15 0.45 0.7 0.5],...'visible','on');t=0:pi/24:4*pi;y=sin(t);plot(t,y,'b')n=length(t);h=line('color',[0 0.5 0.5],...'linestyle','.',...'markersize',25,...'erasemode','xor');k1=uicontrol('parent',h0,...'style','pushbutton',...'position',[80 100 50 30],...'string','开始',...'callback',[...'i=1;',...'k=1;,',...'m=0;,',...'while 1,',...'if k==0,',...'break,',...'end,',...'if k~=0,',...'set(h,''xdata'',t(i),''ydata'',y(i)),',...'drawnow;,',...'i=i+1;,',...'if i>n,',...'m=m+1;,',...'i=1;,',...'end,',...'end,',...'end']);k2=uicontrol('parent',h0,...'style','pushbutton',...'position',[180 100 50 30],...'string','停止',...'callback',[...'k=0;,',...'set(e1,''string'',m),',...'p=get(h,''xdata'');,',...'q=get(h,''ydata'');,',...'set(e2,''string'',p);,',...'set(e3,''string'',q)']);k3=uicontrol('parent',h0,...'style','pushbutton',...'position',[280 100 50 30],...'string','关闭',...'callback','close');e1=uicontrol('parent',h0,...'style','edit',...'position',[60 30 60 20]);t1=uicontrol('parent',h0,...'style','text',...'string','循环次数',...'position',[60 50 60 20]);e2=uicontrol('parent',h0,...'style','edit',...'position',[180 30 50 20]);t2=uicontrol('parent',h0,...'style','text',...'string','终点的X坐标值',...'position',[155 50 100 20]);e3=uicontrol('parent',h0,...'style','edit',...'position',[300 30 50 20]);t3=uicontrol('parent',h0,...'style','text',...'string','终点的Y坐标值',...'position',[275 50 100 20]);实例33：曲线转换按钮h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例33');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhuidiao=[...'if i==1,',...'i=0;,',...'y=cos(x);,',...'delete(h),',...'set(hm,''string'',''正弦函数''),',...'h=plot(x,y);,',...'grid on,',...'else if i==0,',...'i=1;,',...'y=sin(x);,',...'set(hm,''string'',''余弦函数''),',...'delete(h),',...'h=plot(x,y);,',...'grid on,',...'end,',...'end'];hm=uicontrol(gcf,'style','pushbutton',...'string','余弦函数',...'callback',huidiao);i=1;set(hm,'position',[250 20 60 20]);set(gca,'position',[0.2 0.2 0.6 0.6])title('按钮的使用')hold on实例34：栅格控制按钮h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例34');x=0:0.5:2*pi;y=sin(x);plot(x,y)huidiao1=[...'set(h_toggle2,''value'',0),',...'grid on,',...];huidiao2=[...'set(h_toggle1,''value'',0),',...'grid off,',...];h_toggle1=uicontrol(gcf,'style','togglebutton',...'string','grid on',...'value',0,...'position',[20 45 50 20],...'callback',huidiao1);h_toggle2=uicontrol(gcf,'style','togglebutton',...'string','grid off',...'value',0,...'position',[20 20 50 20],...'callback',huidiao2);set(gca,'position',[0.2 0.2 0.6 0.6])title('开关按钮的使用')实例35：编辑框的使用h0=figure('toolbar','none',...'position',[200 150 350 250],...'name','实例35');f='Please input the letter';huidiao1=[...'g=upper(f);,',...'set(h2_edit,''string'',g),',...];huidiao2=[...'g=lower(f);,',...'set(h2_edit,''string'',g),',...];h1_edit=uicontrol(gcf,'style','edit',...'position',[100 200 100 50],...'HorizontalAlignment','left',...'string','Please input the letter',...'callback','f=get(h1_edit,''string'');',...'background','w',...'max',5,...'min',1);h2_edit=uicontrol(gcf,'style','edit',...'HorizontalAlignment','left',...'position',[100 100 100 50],...'background','w',...'max',5,...'min',1);h1_button=uicontrol(gcf,'style','pushbutton',...'string','小写变大写',...'position',[100 45 100 20],...'callback',huidiao1);h2_button=uicontrol(gcf,'style','pushbutton',...'string','大写变小写',...'position',[100 20 100 20],...'callback',huidiao2);实例36：弹出式菜单h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例36');x=0:0.5:2*pi;y=sin(x);h=plot(x,y);grid onhm=uicontrol(gcf,'style','popupmenu',...'string',...'sin(x)|cos(x)|sin(x)+cos(x)|exp(-sin(x))',...'position',[250 20 50 20]);set(hm,'value',1)huidiao=[...'v=get(hm,''value'');,',...'switch v,',...'case 1,',...'delete(h),',...'y=sin(x);,',...'h=plot(x,y);,',...'grid on,',...'case 2,',...'delete(h),',...'y=cos(x);,',...'h=plot(x,y);,',...'grid on,',...'case 3,',...'delete(h),',...'y=sin(x)+cos(x);,',...'h=plot(x,y);,',...'grid on,',...'case 4,',...'delete(h),',...'y=exp(-sin(x));,',...'h=plot(x,y);,',...'grid on,',...'end'];set(hm,'callback',huidiao)set(gca,'position',[0.2 0.2 0.6 0.6]) title('弹出式菜单的使用')实例37：滑标的使用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例37');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;h0=mesh(x,y,z);h1=axes('position',...[0.2 0.2 0.5 0.5],...'visible','off');htext=uicontrol(gcf,...'units','points',...'position',[20 30 45 15],...'string','brightness',...'style','text');hslider=uicontrol(gcf,...'units','points',...'position',[10 10 300 15],...'min',-1,...'max',1,...'style','slider',...'callback',...'brighten(get(hslider,''value''))'); 实例38：多选菜单h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例38');[x,y]=meshgrid(-8:0.5:8);r=sqrt(x.^2+y.^2)+eps;z=sin(r)./r;h0=mesh(x,y,z);hlist=uicontrol(gcf,'style','listbox',...'string','default|spring|summer|autumn|winter',...'max',5,...'min',1,...'position',[20 20 80 100],...'callback',[...'k=get(hlist,''value'');,',...'switch k,',...'case 1,',...'colormap default,',...'case 2,',...'colormap spring,',...'case 3,',...'colormap summer,',...'case 4,',...'colormap autumn,',...'case 5,',...'colormap winter,',...'end']);实例39：菜单控制的使用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例39');x=0:0.5:2*pi;y=cos(x);h=plot(x,y);grid onset(gcf,'toolbar','none')hm=uimenu('label','example');huidiao1=[...'set(hm_gridon,''checked'',''on''),',...'set(hm_gridoff,''checked'',''off''),',...'grid on'];huidiao2=[...'set(hm_gridoff,''checked'',''on''),',...'set(hm_gridon,''checked'',''off''),',...'grid off'];hm_gridon=uimenu(hm,'label','grid on',...'checked','on',...'callback',huidiao1);hm_gridoff=uimenu(hm,'label','grid off',...'checked','off',...'callback',huidiao2);实例40：UIMENU菜单的应用h0=figure('toolbar','none',...'position',[200 150 450 250],...'name','实例40');h1=uimenu(gcf,'label','函数');h11=uimenu(h1,'label','轮廓图',...'callback',[...'set(h31,''checked'',''on''),',...'set(h32,''checked'',''off''),',...'[x,y,z]=peaks;,',...'contour3(x,y,z,30)']);h12=uimenu(h1,'label','高斯分布',...'callback',[...'set(h31,''checked'',''on''),',...'set(h32,''checked'',''off''),',...'mesh(peaks);,',...'axis tight']);h13=uimenu(h1,'label','Sinc函数',...'callback',[...'set(h31,''checked'',''on''),',...'set(h32,''checked'',''off''),',...'[x,y]=meshgrid(-8:0.5:8);,',...'r=sqrt(x.^2+y.^2)+eps;,',...'z=sin(r)./r;,',...'mesh(x,y,z)']);h2=uimenu(gcf,'label','色彩');hl2(1)=uimenu(h2,'label','Default',...'checked','on',...'callback',...[...'set(hl2,''checked'',''off''),',...'set(hl2(1),''checked'',''on''),',...'colormap(''default'')']);hl2(2)=uimenu(h2,'label','spring',...'callback',...[...'set(hl2,''checked'',''off''),',...'set(hl2(2),''checked'',''on''),',...'colormap(spring)']);hl2(3)=uimenu(h2,'label','Summer',...'callback',...[...'set(hl2,''checked'',''off''),',...。

Matlab绘图教程(大量实例PPT)

MATLAB绘图二维数据曲线图pplot函数的基本调用格式为：x,y) )plot(plot(x,y其中x和y为长度相同的向量，分别用于存储x坐标和y坐标数据。

数据例1 在0≤x2π区间内，绘制曲线y=2e-0.5x cos(4πx)1≤区间内绘制曲线205x(4)程序如下：x=0:pi/100:2*pi;cos(4*pi*x);0.5*x).*cos(4*pi*x);y=2*exp(--0.5*x).*y=2*exp(x,y))plot(x,yplot(x yplot(x y)例2 绘制曲线。

绘制曲线程序如下：t=0:0.1:2*pi;x=t.sin(3t);x=t*sin(3*t);y=t.*sin(t).*sin(t); plot(x,y););plot(x,y数最简单的调用格式是包含个输参数plot函数最简单的调用格式是只包含一个输入参数：p()plot(x)在这种情况下，当x是实向量时，以该向量元素的下标为横坐标，元素值为纵坐标画出条连续曲线，标为横坐标，元素值为纵坐标画出一条连续曲线，这实际上是绘制折线图。

绘制多根二维曲线1．plot函数的输入参数是矩阵形式时数的输参数是矩阵形式时(1) 当x是向量，y是有一维与x同维的矩阵时，则绘制出多根不同颜色的曲线。

曲线条数等于y矩阵的另一维数，x被作为这些曲线共同的横坐标。

(2) 当x,y是同维矩阵时，则以x,y对应列元素为横、纵坐标分别绘制曲线，曲线条数等于矩阵的列数。

纵坐标分别绘制曲线曲线条数等于矩阵的列数(3) 对只包含一个输入参数的plot函数，当输入参数对包含个输参数的数当输参数是实矩阵时，则按列绘制每列元素值相对其下标的曲线曲线条数等于输入参数矩阵的列数的曲线，曲线条数等于输入参数矩阵的列数。

当输入参数是复数矩阵时，则按列分别以元素实部和虚部为横、纵坐标绘制多条曲线。

2．含多个输入参数的plot函数含多个输参数的数调用格式为：plot(x1,y1,x2,y2,…,xn,yn)(1) 当输入参数都为向量时，x1和y1，x2和y2，…，(1)当输入参数都为向量时xn和yn分别组成一组向量对，每一组向量对的长度可以不同每一向量对可以绘制出一条曲线度可以不同。

Matlab图形绘制经典案例资料(000001)

Matlab图形绘制经典案例1、三维曲线>> t=0:pi/50:10*pi;>> plot3(sin(2*t),cos(2*t),t)>> axis square>> grid on2、一窗口多图形>> t=-2*pi:0.01:2*pi; >> subplot(3,2,1)>> plot(t,sin(t))>> subplot(3,2,2)>> plot(t,cos(t))>> subplot(3,2,3)>> plot(t,tan(t))>> axis([-pi pi -100 100]) >> subplot(3,2,4)>> plot(t,cot(t))>> axis([-pi pi -100 100]) >> subplot(3,2,5)>> plot(t,atan(t))>> subplot(3,2,6)>> plot(t,acot(t))3、图形样式、标注、题字(也可以利用菜单直接Insert) >> x=0:pi/20:2*pi;>> plot(x,sin(x),'b-.')>> hold on>> plot(x,cos(x),'r--')>> hold on>> plot(x,sin(x)-1,'g:')>> hold on>> plot(x,cos(x)-1)>> xlabel('x');>> xlabel('x轴');>> ylabel('y轴');>> title('图形样式、标注等');>> text(pi,sin(pi),'x=\pi');>> legend('sin(x)','cos(x)','sin(x)-1','cos(x)-1');>> [x1,y1]=ginput(1) %利用鼠标定位查找线上某点的值x1 =2.0893y1 =-0.5000>> gtext('x=2.5') %鼠标定位放置所需的值在线上4、>> fplot('[sin(x),cos(x),sqrt(x)-1]',[0 2*pi]) M文件：myfun.m内容如下：function y=myfun(x)y(:,1)=sin(x);y(:,2)=cos(x);y(:,3)=x^(1/2)-1;再运行：>> fplot('myfun',[0 2*pi])同样可以得到下图5、>> [x,y]=fplot('sin',[0 2*pi]); >> [x1,y1]=fplot('cos',[0 2*pi]); >> plot(x,y,'-r',x1,y1,'-.k')>> legend('y=sinx','y=cosx')6、>> x=[-2:0.2:2];>> y=exp(x)-sin(x);>> plot(x,y,'-or','linewidth',2)7、画出y1=6(sinx-cosx),y2=x2^x-1的图形>> x=[-3:0.1:3];>> y1=6*(sin(x)-cos(x));>> y2=x.*2.^x-1;>> plot(x,y1,'-r',x,y2,'-.k','linewidth',2)8、绘制心形图r=2(1-cosθ)的极坐标图形>> theta=[0:0.01:2*pi];>> polar(theta,2*(1-cos(theta)),'-k')>> polar(theta,2*(1-cos(theta)),'-or')9、用双轴对数坐标绘制y=x*3^x-30的图形>> x=logspace(-3,3);>> y=x.*3.^x-30;>> loglog(y,'-or','linewidth',2);>> grid on10、绘制数据向量的单轴对数坐标图形>> x=[1:50];>> y=[1:50];>> semilogx(x,y,'-*b')%绘制横轴为对数坐标%纵轴为线性坐标>> grid on>> semilogy(x,y,'-*b')%绘制纵轴为对数坐标%横轴为线性坐标>> grid on11、绘制矩阵的条形图，并求出句柄属性值向量。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

>> A=[1 2 3;4 5 6;7 8 9];>> h=bar(A)h =171.0031 174.0026 176.002612、绘制矩阵的水平条形图。

>> y=[3 2 -2 2 1;-1 2 3 7 1;7 2 -3 5 2];>> x=[1:3];>> barh(x,y)13、绘制矩阵的面积图。

>> y=[3 2 -2 2 1;-1 3 3 7 2;-7 5 5 9 3]; >> area(y)14、绘制矩阵的二维饼图>> x=[1 2 3;4 5 6;7 8 9];>> explode=[0 1 0 1 0 1 0 1 0];>> pie(x,explode)15、自行确定数据向量，绘制其散点图。

>> x=rand(1,100);y=randn(1,100);scatter(x,y,20)16、自行确定数据，绘制其柱形图。

>> x=[-2:0.01:4];>>y=randn(1131,1); >>hist(y,x)17、绘制y=sinx在[0,2*pi]上的误差图。

>> x=[0:pi/20:2*pi];>> y=sin(x);>> E=std(y)*ones(size(x));%条形控制>> errorbar(x,y,E)18、绘制火柴杆图。

>> x=[1 1.5 2;3 3.5 4;5 5.5 6]; >> y=[4 3 2;4 8 9;2 7 3]; >> stem(x,y,'fill')%fill意思是“实心点”19、绘制羽列图。

>> U=[-90:5:90]*pi/180; %建立等间距数据>> V=2*ones(size(U));%根据U建立数据>> [U,V]=pol2cart(U,V);转换数据为直角坐标形式>> feather(U,V)20、同一窗口绘制和在[0,30]上的图形。

>> x=[0:0.01:30];>> y1=50*exp(-0.05*x).*sin(x); >> y2=0.5*exp(-0.5*x).*cos(x);>> plotyy(x,y1,x,y2,'plot')% plotyy(x,y1,x,y2,'plot')表示：用左侧y标度绘制（x,y1）用右侧y标度绘制（x,y2）21、绘制8阶魔方矩阵的等值线图和阶梯图。

>> A=magic(8);contour(A) %绘制等值线图stairs(A) %绘制阶梯图22、绘制玫瑰花图。

>> theta=rand(1,200)*2*pi; >> rose(theta,25)23、绘制罗盘图。

>> x=rand(20,1);y=randn(20,1); >> compass(x,y)24、绘制函数的梯度场矢量图。

>> [x,y]=meshgrid([-2:0.1:2]); %建立栅格点数据向量>> z=3.*x.*y*exp(-x.^2-y.^2)-1; %计算函数值向量>> [u,v]=gradient(z,0.2,0.2); %计算梯度值向量>> quiver(x,y,u,v,2) %绘制梯度场矢量图25、给定向量x,y生成网格矩阵。

>> x=[1 2 3 4];>> y=[10 11 12 13 14];>> [U,V]=meshgrid(x,y)U =1 2 3 41 2 3 41 2 3 41 2 3 41 2 3 4V =10 10 10 1011 11 11 1112 12 12 1213 13 13 1314 14 14 1426、生成一个5阶高斯分布矩阵，并给出相应的x,y向量矩阵。

>> [X,Y,Z]=peaks(5)X =-3.0000 -1.5000 0 1.5000 3.0000-3.0000 -1.5000 0 1.5000 3.0000-3.0000 -1.5000 0 1.5000 3.0000-3.0000 -1.5000 0 1.5000 3.0000-3.0000 -1.5000 0 1.5000 3.0000Y =-3.0000 -3.0000 -3.0000 -3.0000 -3.0000-1.5000 -1.5000 -1.5000 -1.5000 -1.50000 0 0 0 01.5000 1.5000 1.5000 1.5000 1.50003.0000 3.0000 3.0000 3.0000 3.0000Z =0.0001 0.0042 -0.2450 -0.0298 -0.0000-0.0005 0.3265 -5.6803 -0.4405 0.0036-0.0365 -2.7736 0.9810 3.2695 0.0331-0.0031 0.4784 7.9966 1.1853 0.00440.0000 0.0312 0.2999 0.0320 0.000027、在-4<=x<=4,-4<=y<=4区域上绘制z=x^2+y^2的三维网格图。