matlab三维作图 plot3(x,y,z) 三维线条图 t=0:pi/50:15*pi;
plot3(sin(t),cos(t),t,'r*') 与plot相似 v=axis 返回各个轴的范围
text(0,0,0,'origin') 在某个坐标点加入文字
plot3 增加维数可以一次画多个图，使所个二维图形眼一个轴排列
三维网线图的绘制 mesh(x,y,z) 网格图
mesh(x,y,z,c) 四维作图,(x,y,z)代表空间三维，c代表颜色维 mesh(…,‘property name‘,property value,…) 设置曲面各属性的值 [x,y,z]=sphere(12);
mesh(x,y,z)，hidden off 曲面设置为透明 meshc(x,y,z) 画网格图和基本的等值线图 meshz(x,y,z) 画包含零平面的网格图
waterfall(x,y,z) 与mesh一样，只是在效果上它的网格线只在x轴一个方向出现，呈瀑布状水线 两个变量的标量指令meshgrid(x)或meshgrid(x,y) (p179)
将两个一维向量生成两个二维向量，以便进行z=f(x,y)运算，算出z的所有值，z为x y的标量指令 [X,Y]=meshgrid(x) meshgrid(x,x)的简略式 [X,Y]=meshgrid(x,y)
[X,Y,Z]=meshgrid(x,y,z) 用于三维图形的绘制 [x,y]=meshgrid([-2:0.1:2]); z=x.*exp(-x.^2-y.^2); plot3(x,y,z)
surf(x,y,z,c) 着色表面图 surf(x,y,z) 隐含着c=z
surf（z）隐含着x，y的值为surf指令根据z的尺寸自动生成 surfc 画出具有基本等值线的曲面图 surfl 画出一个具有亮度的曲面图
shading flat 网线图的某整条线段或曲面图的某个贴片都着一种颜色
shading interp 某一线段或贴片上各点的颜色由线或片的顶端颜色经线性插值而得 曲面图不能设成网格图那样透明，但需要时，可以在孔洞处将数据设成nun
等高线的绘制
在二维空间绘制等高线contour
contour(x,y,z,n) 绘制n条等值线（n可省略）
contour(x,y,z,v) 在向量v所指定的高度上绘制等高线（可省） c=contour(x,y,z) 计算等值线的高度值
c=contourc(x,y,z,n) 计算n条等高线的x－y坐标数据
c=contourc(x,y,z,v) 计算向量v所指定的等高线的x－y坐标数据 clabel(c) 给c阵所表示的等高线加注高度标识 clabel(c,v) 给向量v所指定的等高线加注高度标识
clabel(c,‘manual‘) 借助鼠标给点中的等高线加注高度标识 三维空间绘制等高线contour3(x,y,z) [x,y,z]=peaks(30); contour3(x,y,z,16,'g') 二元函数的伪彩图pcolor(x,y,z)
是指令surf的二维等效指令，代表伪彩色，可与contour单色等值线结合画彩色等值线图 [x,y,z]=peaks(30); pcolor(x,y,z); 伪彩色
shading interp 颜色插值，使颜色平均渐变 hold on,contour(x,y,z,20,'k')... 画等值线 colorbar('horiz') 水平颜色标尺 c=contour(x,y,z,8); clabel(c) 标注等高线 矢量场图(速度图)quiver
用于描述函数z=f(x,y)在点(x,y)的梯度大小和方向 [X,Y]=meshgrid(x,y) X,Y为Z阵元素的坐标矩阵
[U,V]=gradient(Z,dx,dy) U，V分别为Z对x对y的导数，dx dy是x y方向上的计算步长 quiver(X,Y,U,V,s,‘linespec‘,‘filled‘) U，V为必选项，决定矢量场图中各矢量的大小和方向，s为指定所画箭头的大小，缺省时取1,linespec为字符串，指定合法的线形和彩色，filled用于填充定义的绘图标识符
[x,y]=meshgrid(-2:.2:2,-1:.15:1); z=x.*exp(-y.^2);
[px,py]=gradient(z,.2,.15); contour(x,y,z);
hold on,quiver(x,y,px,py),axis image 多边形的填色fill(x,y,c)
c定义颜色字符串，可以是‘r‘,‘b‘等，也可以用RGB三色表示[r,g,b]值为0-1 图形的四维表现
网络资源：
1. USENET新闻组
? MATLAB的新闻组是comp.soft-sys.MATLAB。浏览器指
向
? http://www.dejanews.com ? http://www.google.com 2. 网络上的工具箱
? http://www.mathtools.com/
? http://www.mathworks.com/matlabcentral/fileexchange/loadCategory.do
3. BBS
? 哈尔滨工业大学bbs: telnet://bbs.hit.edu.cn ? 上海交通大学bbs: telnet://bbs.sjtu.edu.cn ? 清华大学bbs mathtools 4. www服务
? http://www.mathworks.com ? http://matlab.turbo.hit.edu.cn ? Matlab 大观园: http://www.matlab-world.com
5. 首先是MATLAB窗口中的help命令，它最常用、最快，往往也是最有效的。给出关键字就能找到相应的说明。Help命令简单易用，但是它要求准确给出关键字，如果记忆不准，就很难找到。这时就应该用第二种方案。
6. Help Desk：这是MATLAB的HTML格式的帮助文档，有着很好的组织。可以在命令行键入helpdesk，进入“Matlab Help Desk”。Help Desk的搜索可以使用不完全关键字，这样即使遇
到记忆模糊的情况也可以很快查到。
7. 前面两种方法所得到的帮助信息常常是不够细致的，需要用第三种途径来补充，这就是电子版的使用手册。MATLAB里包含了大量关于MATLAB安装、使用、编程以及各种工具箱等的电子版手册，全部为PDF格式,具有搜索功能，因此也可以在指定的手册中搜索关键字，使用非常方便。
编程学习】matlab实用程序百例
实例1：三角函数曲线（1）
function shili01 h0=figure('toolbar','none',... 'position',[198 56 350 300],... 'name','实例01'); h1=axes('parent',h0,... 'visible','off'); x=-pi:0.05:pi; y=sin(x); plot(x,y); xlabel('自变量X'); ylabel('函数值Y'); title('SIN( )函数曲
线'); grid on
实例2：三角函数曲线（2）
function shili02 h0=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 on xlabel('自变量X'); ylabel('函数值Y'); title('三角函数');
实例3：图形的叠加
function shili03 h0=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 on xlabel('自变量X'); ylabel('函数值Y');
title('三角函数');
实例5：单个轴窗口显示多个图形
function shili05 h0=figure('toolbar','none',... 'position',[200 150 450 250],...
'name','实例05'); t=0:pi/10:2*pi; [x,y]=meshgrid(t); subplot(2,2,1)
plot(sin(t),cos(t)) axis equal
subplot(2,2,2) z=sin(x)-cos(y); plot(t,z) axis([0 2*pi -2 2]) subplot(2,2,3) h=sin(x)+cos(y); plot(t,h) axis([0 2*pi -2 2]) subplot(2,2,4) g=(sin(x).^2)-(cos(y).^2); plot(t,g) axis([0 2*pi -1 1])
实例6：图形标注
function shili06 h0=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 shili07 h0=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 shili08 h0=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 on area(x,profits2,'facecolor',[0.9 0.85 0.7],... 'edgecolor','y',... 'linewidth',3) hold on area(x,profits3,'facecolor',[0.3 0.6 0.7],... 'edgecolor','r',... 'linewidth',3) hold on area(x,profits4,'facecolor',[0.6 0.5 0.9],... 'edgecolor','m',... 'linewidth',3) hold off set(gca,'xtick',[91:95])