2、受广大科技工作者的欢迎。qMatlab是一种广泛应用于工程计算及数值分析领域的新型高级语言,自1984年推向市场以来,历经二十多年的发展与竞争,现已成为国际公认的最优秀的工程应用开发环境。4Matlab简介qMatlab:MatrixLaboratory矩阵实验室qMatlab的发展l1980年,Moler教授用Fortran语言编写了集命令翻译、科学计算于
3、一身的一套交互式软件系统。l1984年,Moler等成立了TheMathWorks的公司,用C语言完全改写Matlab,并推出第一个商业版。l增添图形图像处理、符号运算、以及与其他流行软件的接口功能,使得Matlab的功能越来越强大。l到九十年代,在国际上30几个数学类科技应用软件中,Matlab在数值计算方面独占鳌头。5q目前,Matla
4、b已成为世界顶尖的数学应用软件,以其强大的工程计算、算法研究、工程绘图、应用程序开发、数据分析和动态仿真等功能,在航空航天、机械制造和工程建筑等领域发挥着越来越重要的作用。就影响而言,至今仍然没有一个别的计算软件可与Matlab匹敌。qMatlab的发行1984年,Matlab1.0(DOS版,182K,20多个函数)1992年,Matlab4
5、.0(93年推出Windows版,加入simulink)1994年,Matlab4.2(得到广泛重视和应用)1999年,Matlab5.3(真正实现32位运算)2002年,Matlab6.5(采用JIT加速器)2004年,Matlab7.0自2006年起,Matlab每年更新两次Matlab简介6q购买(下载)Matlab软件Matlab的安装q插入光盘,自动运行;或点击安装程序setup.exen必须安装的
6、部件MatlabSymbolicMathToolboxExtendedSymbolicMathl这里以Matlab7.0.4为例l自R2008a开始增加了激活要求7Matlab的特点与功能qMatlab具有很强的数值计算功能lMatlab以矩阵作为数据操作的基本单位,但无需预先指定矩阵维数(动态定维)l按照IEEE的数值计算标准进行计算l提供十分丰富的数值计算函数,方便计算
7、,提高效率lMatlab命令与数学中的符号、公式非常接近,可读性强,容易掌握qMatlab是一个交互式软件系统输入一条命令,立即就可以得出该命令的结果8Matlab的特点与功能qMatlab符号计算功能Matlab和著名的符号计算语言Maple相结合qMatlab的编程功能Matlab具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序
8、语言特征,而且简单易学、编程效率高。通过Matlab进行编程完成特定的任务qMatlab的绘图功能Matlab提供丰富的绘图命令,很方便实现数据的可视化9qMatlab丰富的工具箱(toolbox)Matlab的特点与功能根据专门领域中的特殊需要而设计的各种可选工具箱qMatlab的Simulink动态仿真集成环境提供建立系统模型、选择仿真参数和数值
9、算法、启动仿真程序对该系统进行仿真、设置不同的输出方式来观察仿真结果等功能SymbolicMathPDEOptimizationSignalprocessImageProcessStatisticsControlSystemSystemIdentification10Matlab的基本用法qMatlab系统的启动l使用Windows“开始”菜单l运行Matlab系统启动程序matlabl双击Matlab快捷图标qMatlab系统
10、的退出l在Matlab主窗口File菜单中选择ExitMatlabl在Matlab命令窗口输入exit或quitl单击Matlab主窗口的“关闭”按钮11Matlab的工作界面命令窗口当前工作目录当前工作空间输入命令的历史记录命令提示符12q定义矩阵:直接输入法l矩阵用方括号“”括起例:A=123;456;789l矩阵同一行中的元素之间用空格或逗号分隔
11、l矩阵行与行之间用分号分开l直接输入法中,分号可以用回车代替例:矩阵qMatlab的操作对象是l13q矩阵元素可以是任何数值表达式例:x=-1.3,sqrt(3),(1+2+3)*4/5矩阵元素赋值q矩阵元素的单独赋值例:x(5)=abs(x(1)例:x(5)=abs(x(6)自动将向量x的长度扩展到5,并将未赋值部分置零。14q大矩阵可以把小矩阵作为其元素矩阵
12、元素赋值例:A=A;111213在原矩阵的下方加一行如何在原矩阵的右边添加一列?15q单个元素的引用例:A(2,3)矩阵元素的引用q多个元素的引用:冒号的特殊用法利用小括弧和元素所在的位置(下标)x(i):向量x中的第i个元素A(i,j):矩阵A中的第i行,第j列元素a:b:c产生一个由等差序列组成的向量;a是首项,b是公差,
13、c确定最后一项;若b=1,则b可以省略。例:x=1:2:5y=1:2:6例:x=2:1:5y=2:5例:x=3:2:116例:x(1:3)A(3,1:3)矩阵元素的引用A(i:j,m:n)表示由矩阵A的第i到第j行和第m到第n列交叉线上的元素组成的子矩阵。可利用冒号提取矩阵的整行或整列。例:A(1,:)A(:,1:3)A(:,:)17Matlab帮助系统q联机帮助lhelp显示指定命
15、u以字母开头u后面可以跟字母、数字和下划线u长度不超过63个字符(6.5版本以前为19个)u变量名区分字母的大小写qMatlab语句的通常形式变量=表达式表达式是用运算符将有关运算量连接起来的式子,其结果被赋给赋值号“=”左边的变量20q分号和续行符的作用Matlab变量l若不想在屏幕上输出结果,可以在语句最后加分号l如果语
16、句很长,可用续行符“”(三个点)续行续行符的前面最好留一个空格例:21q变量的查询Matlab变量lwho显示工作空间中的所有变量lwhos查看工作空间中变量的详细属性22q系统预定义变量Matlab变量upi:圆周率,其值为imag(log(-1)uinf,Inf:无穷大unan,NaN:Not-a-Number,一个不定值,如0/0ueps:浮点运算相对精度q
17、特殊变量ansui,j:虚部单位,即1应尽量避免给系统预定义变量重新赋值!23q数与算术表达式Matlab数值运算uMatlab中的数默认是双精度实数,表示方法同C语言3,-9,0.4,1.603e-12,3.23e+20u浮点运算的相对误差为epsu浮点数表示范围为:10-30810308u复数的输入lz=3+4i(4与i之间不能有空格)lz=3+4*i复数作为矩阵元
18、素输入时,加号两边不能有空格!24q数学运算符u+加法Matlab数值运算u-减法u*乘法u/和除法(右除和左除)u幂运算q命令分隔符:逗号和分号25q输出格式Matlab的输出uMatlab以双精度执行所有的运算,运算结果可以在屏幕上输出,同时赋给指定变量;若无指定变量,则系统会自动将结果赋给变量“ans”uMatlab中数的输出格
19、式可以通过format命令指定format只改变变量的输出格式,但不会影响变量的值!26各种format格式格式解释例format短格式(缺省显示格式),同short3.1416formatshort短格式(缺省显示格式),只显示5位3.1416formatlong长格式,双精度数15位,单精度数7位3.14159265358979formatshorte短格式e方式(科学计数格式)3.141
20、6e+000formatlonge长格式e方式3.141592653589793e+000formatshortg短格式g方式3.1416formatlongg长格式g方式3.14159265358979DrWangZhengsheng-LectureNotes27变量的存储q存储当前工作空间中的变量usave将所有变量存入文件matlab.matusavemydata将所有变量存入指定文件mydata.matq存储指定的变量usav
21、emydata.mat将所有变量存入文件mydata.matsave文件名变量名列表例:savemydataAxz变量名列表中各变量之间用空格分隔28变量的读取q将数据文件中的变量载入当前工作空间uloadmydata载入数据文件中的所有变量uloadmydataAx从数据文件中提取指定变量q清除当前工作空间中的变量uclear清除当前工作空间中的所有变量
22、uclearAx清除指定的变量DrWangZhengsheng-LectureNotes29qMatlab的命令记忆功能:上下箭头键q命令补全功能:Tab键可以先输入命令的前几个字符,再按上下键缩小搜索范围q用Esc键删除命令行DrWangZhengsheng-LectureNotes30Lecture2MATLAB矩阵(数值)运算DrWangZhengsheng-LectureNotes31n利用
23、函数建立数值矩阵:MATLAB提供了许多生成和操作矩阵的函数,可以利用它们去建立矩阵。n例如:reshape函数和diag函数等。nreshape函数用于建立数值矩阵。ndiag函数用于产生对角阵。n利用M文件建立矩阵:对于比较大且比较复杂的矩阵,可以为它专门建立一个M文件。其步骤为:n第一步:使用编辑程序输入文件内容。n第二步:把输
24、入的内容以纯文本方式存盘(设文件名为mymatrix.m)。n第三步:在MATLAB命令窗口中输入mymatrix,就会自动建立一个名为AM的矩阵,可供以后显示和调用。DrWangZhengsheng-LectureNotes32n利用M文件建立矩阵:对于比较大且比较复杂的矩阵,可以为它专门建立一个M文件。其步骤为:n第一步:使用编辑程序输入文件内容。n第二步
25、:把输入的内容以纯文本方式存盘(设文件名为mymatrix.m)。n第三步:在MATLAB命令窗口中输入mymatrix,就会自动建立一个名为AM的矩阵,可供以后显示和调用。n子矩阵操作DrWangZhengsheng-LectureNotes33n2矩阵的基本运算n()矩阵转置n()矩阵加和减n()矩阵乘法n()矩阵除法Ab=inv(A)*bn()矩阵的乘方a2DrWan
26、gZhengsheng-LectureNotes34n3矩阵的函数nhelpmatfunnMatrixfunctions-numericallinearalgebra.nnMatrixanalysis.nnorm-Matrixorvectornorm.nnormest-Estimatethematrix2-norm.nrank-Matrixrank.ndet-Determinant.ntrace-Sumofdiagonalelements.nnull-Nullspace.north-O
27、rthogonalization.nrref-Reducedrowechelonform.nsubspace-Anglebetweentwosubspaces.DrWangZhengsheng-LectureNotes35n3矩阵的函数nLinearequations.nand/-Linearequationsolution;usehelpslash.ninv-Matrixinverse.nrcond-LAPACKreciprocalconditionestimatorncond-Conditio
28、nnumberwithrespecttoinversion.ncondest-1-normconditionnumberestimate.nnormest1-1-normestimate.nchol-Choleskyfactorization.ncholinc-IncompleteCholeskyfactorization.nlu-LUfactorization.nluinc-IncompleteLUfactorization.nqr-Orthogonal-triangulardecomposition.nlsqno
29、nneg-Linearleastsquareswithnonnegativityconstraints.npinv-Pseudoinverse.nlscov-Leastsquareswithknowncovariance.n36n3矩阵的函数nEigenvaluesandsingularvalues.neig-Eigenvaluesandeigenvectors.nsvd-Singularvaluedecomposition.ngsvd-Generalizedsingularvaluedecomposi
30、tion.neigs-Afeweigenvalues.nsvds-Afewsingularvalues.npoly-Characteristicpolynomial.npolyeig-Polynomialeigenvalueproblem.ncondeig-Conditionnumberwithrespecttoeigenvalues.nhess-Hessenbergform.nqz-QZfactorizationforgeneralizedeigenvalues.nschur-Schurdecomposi
31、tion.nn37n3矩阵的函数nMatrixfunctions.nexpm-Matrixexponential.nlogm-Matrixlogarithm.nsqrtm-Matrixsquareroot.nfunm-Evaluategeneralmatrixfunction.nnFactorizationutilitiesnqrdelete-DeleteacolumnorrowfromQRfactorization.nqrinsert-InsertacolumnorrowintoQRfa
32、ctorization.nrsf2csf-Realblockdiagonalformtocomplexdiagonalform.ncdf2rdf-Complexdiagonalformtorealblockdiagonalform.nbalance-Diagonalscalingtoimproveeigenvalueaccuracy.nplanerot-Givensplanerotation.ncholupdate-rank1updatetoCholeskyfactorization.nqrupdate-
33、rank1updatetoQRfactorization.DrWangZhengsheng-LectureNotes38n4建立矩阵的函数n常用函数有:neye(size(A)产生与A矩阵同阶的单位矩阵nzeros(m,n)产生0矩阵nones(m,n)产生幺矩阵nrand(m,n)产生随机元素的矩阵nSize(a)返回包含两个元素的向量。nLength(a)返回向量的长度。DrWangZhengsheng-Lecture
34、Notes39常见矩阵生成函数zeros(m,n)生成一个m行n列的零矩阵,m=n时可简写为zeros(n)ones(m,n)生成一个m行n列的元素全为1的矩阵,m=n时可写为ones(n)eye(m,n)生成一个主对角线全为1的m行n列矩阵,m=n时可简写为eye(n),即为n维单位矩阵diag(X)若X是矩阵,则diag(X)为X的主对角线向量若X是向量,diag(X)
35、产生以X为主对角线的对角矩阵tril(A)提取一个矩阵的下三角部分triu(A)提取一个矩阵的上三角部分rand(m,n)产生01间均匀分布的随机矩阵m=n时简写为rand(n)randn(m,n)产生均值为0,方差为1的标准正态分布随机矩阵m=n时简写为randn(n)其它特殊矩阵生成函数:magic、hilb、pascalDrWangZhengsheng-LectureNotes40n5数
36、组运算n(1)数组的加和减n(2)数组的乘和除n(3)数组的乘方DrWangZhengsheng-LectureNotes416应用举例n1、行列式计算;n2、求解线性方程组;n3、解特征值问题DrWangZhengsheng-LectureNotes42Lecture3MATLAB图形与可视化(Graphic)DrWangZhengsheng-LectureNotes43(1)掌握图形窗口的创建与控制,
37、以及图形窗口的基本操作;(2)熟练掌握二维和三维绘图基本的命令、线型控制;(3)初步掌握用特殊的图形来表现特殊数据的性质,如面积图、直方图、饼图等。(4)掌握坐标轴的控制和图形标注命令及其用法。DrWangZhengsheng-LectureNotes44A、二维平面图形与坐标系1.几个基本的绘图命令a.线性坐标曲线plot函数命令pl
38、ot是MATLAB二维曲线绘图中最简单、最重要、使用最广泛的一个线性绘图函数。它可以生成线段、曲线和参数方程曲线的函数图形。命令格式:plot(X,Y)plot(x1,y1,x2,y2,):综合调用方式DrWangZhengsheng-LectureNotes45用命令plot(x,y)绘制函数y=cos(x)在两个周期内的图形。x=0:0.01:2*pi;y=cos(x);plot(x,y)在同一图形窗
39、口中用命令plot(x,y)绘出正弦余弦函数的图形。x=0:0.01:2*pi;y=sin(x);cos(x);plot(x,y)DrWangZhengsheng-LectureNotes46二维函数曲线专用命令fplot用plot绘图在确定自变量的取值间隔时,一般采用平均间隔,有时会因某处间距太大,而不能反映出函数的变化情况。fplot是绘制函数y=f(x)图形的专用命令,它的数据点是自
40、适应产生的,对那些导数变化较大的函数,用fplot函数绘出的曲线比等分取点所画出的曲线更加接近真实。fplot函数命令的调用格式为:X,Y=fplot(fun,lims)fun:函数名字符串;lims:定义x的取值区间,lims=xmin,xmax;二维函数曲线专用命令ezplotDrWangZhengsheng-LectureNotes472.线型和颜色plot函数可以设置曲线的
41、线段类型、定点标记和线段颜色。常用的线段、颜色与定点标记参数DrWangZhengsheng-LectureNotes48调用格式:plot(x,y,s),s为类型说明参数,是字符串。s字符串可以是三种类型的符号之一,也可以是线型与颜色和定点标记与颜色的组合;如果没有s参数,plot将使用缺省设置(实线,前七种颜色顺序着色)绘制曲线;在当前
42、坐标系中绘图时,每调入一次绘图函数,MATLAB将擦掉坐标系中已有的图形对象。可以用holdon命令在一个坐标系中增加新的图形对象。注意MATLAB会根据新图形的大小,重新改变坐标系的比例。用不同的线型和标注来绘制两条曲线。t1=0:0.1:2*pi;t2=0:0.1:6;y1=sin(t1);y2=sqrt(t2);plot(t1,y1,:hb,t2,y2,-g)DrWangZhengsheng
43、-LectureNotes493.图形窗口的分割有时需要在一个图形窗口中显示几幅图,以便对几个函数进行直观、便捷的比较。由于每个绘图命令在绘制数据图像时都会将已有图形覆盖掉,而用hold命令不能实现同时显示几个不同坐标尺寸下的图形,用figure命令再创窗口又很难同时比较由不同的数据绘得的图像。实现在同一个窗口中同时显示多
44、个图像的命令subplot。使用格式为:subplot(m,n,i)其含义为:把图形窗口分割为m行n列子窗口,然后选定第i个窗口为当前窗口。subplot命令不仅用于二维图形,对三维图形一样适用。其本质是将figure窗口分为几个区域,再在每个区域内分别绘图。DrWangZhengsheng-LectureNotes50用subplot函数把两种不同的图形综合在
45、一个图形窗口中。subplot(2,2,1)t=0.1:0.1:2*pi;y=sin(t);semilogx(t,y)gridonsubplot(2,2,2)t=0:0.1:4*pi;y=sin(t);plot(t,y)subplot(2,2,3)x=1:0.01:5;y=exp(x);plotyy(x,y,x,y,semilogx,plot)subplot(2,2,4)x=1:0.1:10;y=sqrt(x);plot(x,y,:rd)DrWangZhengsheng-LectureNotes514.坐标系的调整实
46、现坐标系的调整的命令是axis函数。调用格式为:axis(xmin,xmax,ymin,ymax,zmin,zmax)坐标的最小值(xmin,ymin,zmin)必须小于相应的最大值(xmax,ymax,zmax),否则会出错。自动坐标系与用axis函数调整后的坐标系的比较。subplot(2,1,1)t=0:0.1:4*pi;y=sin(t);plot(t,y)subplot(2,1,2)t=0:0.1:4*pi;y=sin(t);plot(t
47、,y)axis(0,max(t),min(y),max(y)DrWangZhengsheng-LectureNotes52B、三维绘图1.三维曲线绘图命令三维函数plot3主要用来表现单参数的三维曲线,与二维绘图函数plot相比,只多了第三维数据。其调用格式为:plot3(X1,Y1,Z1,s1,X2,Y2,Z2,s2,)参数的含义如下:Xn、Yn、Zn:第一到三维数据,是尺寸相等的向量/矩阵;s、s1、
48、s2:是字符串,用来设置线型、颜色、数据点标记。DrWangZhengsheng-LectureNotes53x、y、z是向量时,plot3命令的使用t=0:0.1:8*pi;plot3(sin(t),cos(t),t)title(绘制螺旋线)%用命令title对图形主题进行标注xlabel(sin(t),FontWeight,bold,FontAngle,italic)ylabel(cos(t),FontWeight,bold,FontAngle,italic)zl
49、abel(t,FontWeight,bold,FontAngle,italic)%命令zlabel用来指定z轴的数据名称gridonx、y、z都是矩阵时,plot3命令的使用X,Y=meshgrid(-pi:0.1:pi);Z=sin(X)+cos(Y);plot3(X,Y,Z)DrWangZhengsheng-LectureNotes54为了绘制定义在平面区域D=x0,xmy0,yn上的三维曲面z=f(x,y),首先将x0,xm在x方向分成m份,将y0
50、,yn在y方向分成n份,由各划点分别作平行于坐标轴的直线,将区域D分成mn个小矩形;对于每个小矩形,计算出网格点的函数值,决定出空间中四个顶点(xi,yi,f(xi,yi)),连接四个顶点得到一个空间的四边形片;所有四边形片连在一起构成函数z=f(x,y)定义在区域D上的空间网格曲面。因此,三维曲面绘图命令可分为平面网格点的生成
51、、在平面网格基础上绘制三维网格及对三维表面进行处理三个步骤。DrWangZhengsheng-LectureNotes55a.平面网格点的生成函数命令meshgrid用来生成x-y平面上的网格点矩阵。调用形式为:X,Y=meshgrid(x,y)X,Y=meshgrid(x)等价于X,Y=meshgrid(x,x)参数含义如下:x:是区间x0,xm上分划的向量;y:是区间y0,yn上分划的向量;
52、X,Y:输出变量矩阵,矩阵X的行向量都是向量x,矩阵Y的列向量都是向量y。函数meshgrid将由两个向量决定的区域转换为对应的网格点矩阵。DrWangZhengsheng-LectureNotes56函数,定义区域为-2,2-2,2。生成网格并计算其网格点上的函数值。X,Y=meshgrid(-2:2:2,-2:2:2);X,Y%将划分结果输出至矩阵ans=-202-2
53、-2-2-202000-202222Z=X.*exp(-X.2-Y.2);%计算网格点上的函数值赋予变量ZZ=-0.000700.0007-0.036600.0366-0.000700.000722yxexzDrWangZhengsheng-LectureNotes57利用函数mesh生成网格曲面。调用格式为:mesh(X,Y,Z,C):X、Y、Z、C是同维数的矩阵,X、Y、Z对应空间上的网格点,网格
54、线颜色由C决定;mesh(X,Y,Z):相当于上面的C=Z的情况;mesh(x,y,Z,C):x和y是向量,Z和C是同维数的矩阵,网格曲面的网格顶点是(x(j),y(i),Z(i,j)),网格线的颜色由矩阵C决定;mesh(x,y,Z):相当于上面的C=Z的情况;mesh(Z,C):等价于mesh(x,y,Z,C),此时向量x=1:n,向量y=1:m;mesh(Z):相当于上面的C=
56、otes59函数surf可实现对网格曲面片进行着色,将网格曲面转化为实曲面。surf命令的调用格式与mesh相同。利用三维网格表面命令surf绘制图形。z=peaks;%绘制山峰的图像,将函数值赋予变量zsurf(z)%对山峰的图像进行着色处理shadinginterp%函数shading改变着色方式DrWangZhengsheng-LectureNotes603.等高线图形的
57、绘制4.三维视图可视效果的控制5.柱面和球面的三维表达(书)DrWangZhengsheng-LectureNotes61a.柱面的表达cylindercylinder命令中,柱面的轴线定义为z轴,只要给出母线的描述就可完成一个柱面。调用格式为:X,Y,Z=cylinder(R,N);X,Y,Z=cylinder(R):缺省值N=20;X,Y,Z=cylinder:缺省值N=20,R=1,1。R:是一描述
58、柱面母线的向量;N:是旋转柱面上的分割线条数;X,Y,Z:是返回的x,y,z坐标向量。绘制一个柱面。t=pi:0.01:3*pi;r=sin(t)+t;cylinder(r,30)shadinginterpDrWangZhengsheng-LectureNotes62b.球面的表达sphere调用格式为:X,Y,Z=sphere(N):产生一个(N+1)(N+1)的矩阵,然后用函数surf命令绘制一个单位的球面,
59、N为设置分割线的条数;X,Y,Z=sphere:缺省值N=20。画一个球面。X,Y,Z=sphere;surf(X,Y,Z)DrWangZhengsheng-LectureNotes63为了将抽象的数据表达得更形象,除了绘制二维、三维图形外,还要用到直方图、面积图、饼图等特殊图形。特殊图形指令DrWangZhengsheng-LectureNotes641.面积图命令area表现各个不同部分
60、对整体所作的贡献area(X,Y):与plot的命令的使用方法相似,将连线图到x轴的那部分填上了颜色;area(Y):缺省值X=1:SIZE(Y);area(X,Y,LEVEL)或area(Y,LEVEL):填色部分为由连线图到y=level的水平线之间的部分。绘制一面积图X=-2:2;Y=3,5,2,4,1;5,4,2,3,5;3,4,5,2,1;area(X,Y)legend(因素1,因素2,因素3)gridonDrWangZ
61、hengsheng-LectureNotes652.直方图命令bar直方图常用于统计数据的作图,有bar、bar3、barh和bar3h几种函数,其调用格式类似。以函数bar为例:bar(X,Y):X是横坐标向量,Y可以是向量或矩阵。Y是向量时,每一个元素对应一个竖条;Y是m行n列矩阵时,将画出m组竖条,每组包括n个竖条;bar(Y):横坐标使用缺省值X=1
62、:M;bar(X,Y,WIDTH)或bar(Y,WIDTH):用WIDTH指定竖条的宽度,如果WIDTH1,条与条之间将重合。缺省宽度为0.8;bar(.,grouped):产生缺省的组合直方图;bar(.,stacked):产生累积的直方图;bar(.,linespec):指定条的颜色;H=bar(.):返回条形图对象的句柄。DrWangZhengsheng-LectureNotes66用绘制直方图的几种命令绘
63、制直方图。X=-2:2;Y=3,5,2,4,1;5,4,2,3,5;3,4,5,2,1;subplot(2,2,1)bar(X,Y,r)xlabel(x)ylabel(y)colormap(cool)subplot(2,2,2)barh(X,Y,grouped)xlabel(y)ylabel(x)colormap(cool)subplot(2,2,3)bar(X,Y,stacked)xlabel(x)ylabel(Sigmay)colormap(summer)subplot(2,2,4)barh(X,Y,stacked)xlabel(y);yla
64、bel(Sigmax)colormap(summer)DrWangZhengsheng-LectureNotes67绘制三维直方图。X=-2:2;Y=3,5,2,4,1;5,4,2,3,5;3,4,5,2,1;subplot(2,2,1)bar3(X,Y,r)zlabel(y)ylabel(x)colormap(cool)subplot(2,2,2)bar3h(X,Y,grouped)ylabel(x)zlabel(y)colormap(cool)subplot(2,2,3)bar3(X,Y,stacked)ylabel(x)zla
65、bel(Sigmay)colormap(summer)subplot(2,2,4)bar3h(X,Y,stacked)zlabel(x)ylabel(Sigmay)colormap(summer)DrWangZhengsheng-LectureNotes683.饼图命令pie饼图又叫扇形图,用于显示向量中元素所占向量元素总和的百分比。pie和pie3分别用于绘制二维和三维饼图。调用格式:pie(X):向量X的饼图。把X的每一个元
66、素在所有元素总和中占的比例表达出来;pie(X,EXPLODE):向量EXPLODE(和向量X长度相等)用于指定饼图中抽出一部分的块(非零值对应的块);pie(.,LABELS):LABELS是用于标注饼图的字符串数组,其长度必须和向量X相等;H=pie(.):返回包括饼图和文本对象句柄。DrWangZhengsheng-LectureNotes69用函数pie和pie3绘制饼图。x=200,360,120,400,320;subplot(2,2,1),pie(x,00010)subplot(2,2,2),pie3(x,00010)subplot(2,2,3),pie(x(2:5)subplot(2,2,4),x=0.1,0.12,0.21,0.34,0.11;pie3(x,A,B,C,D,E)DrWangZhengsheng-Lec