Matlab与相关工具箱.ppt

单击此处编辑母版标题样式,单击此处编辑母版文本样式,第二级,第三级,第四级,第五级,*,第,2,章,Matlab,与相关工具箱,1 MATLAB,操作基础,2 MATLAB,矩阵及其运算,3 MATLAB,程序设计,4,MATLAB,工,具箱,主要内容,1.1 MATLAB,概述,1.2 MATLAB,的运行环境与安装,1.3 MATLAB,集成环境,1.4 MATLAB,帮助系统,1 MATLAB,操作基础,1.1 MATLAB,概述,1.1.1 MATLAB,的发展,1984,年，,MATLAB,第,1,版,(DOS,版,)1992,年，,MATLAB4.0,版,1997,年，,MATLAB 5.0,版,2000,年，,MATLAB 6.0,版,2002,年，,MATLAB 6.5,版,2004,年，,MATLAB 7.0,版,2006,年，,MATLAB 2006,版,2007,年，,MATLAB 2007,版,1.1.2 MATLAB,的主要功能,1,数值计算和符号计算功能,MATLAB,以矩阵作为数据操作的基本单位，还提供了十分丰富的数值计算函数。,MATLAB,和著名的符号计算语言,Maple,相结合，使得,MATLAB,具有符号计算功能。,2,绘图功能,MATLAB,提供了两个层次的绘图操作：一种是对图形句柄进行的低层绘图操作，另一种是建立在低层绘图操作之上的高层绘图操作。,3,编程语言,MATLAB,具有程序结构控制、函数调用、数据结构、输入输出、面向对象等程序语言特征，而且简单易学、编程效率高。,4,MATLAB,工具箱,MATLAB,包含两部分内容：基本部分和各种可选的工具箱。,MATLAB,工具箱分为两大类：功能性工具箱和学科性工具箱。,1.1.3,初识,MATLAB,例,1-1,绘制正弦曲线和余弦曲线。,x=0:0.5:360*pi/180;plot(x,sin(x),x,cos(x);,例,1-2,求方程,3x,4,+7x,3,+9x,2,-23=0,的全部根。,p=3,7,9,0,-23;%,建立多项式系数向量,x=roots(p)%,求根,例,1-3,求积分,quad(x.*log(1+x),0,1),例,1-4,求解线性方程组。,a=2,-3,1;8,3,2;45,1,-9;b=4;2;17;x=inv(a)*b,1.2 MATLAB,的运行环境与安装,1.2.1 MATLAB,的运行环境,硬件环境：,(1)CPU(2),内存,(3),硬盘,(4)CD-ROM,驱动器和鼠标。软件环境：,(1)Windows 98/NT/2000,或,Windows XP(2),其他软件根据需要选用,1.2.2 MATLAB,的安装,安装,MATLAB 6.5,系统，需运行系统自带的安装程序,setup.exe,，可以按照安装提示依次操作。,1.3 MATLAB,集成环境,1.3.1,启动与退出,MATLAB,集成环境,1,MATLAB,系统的启动与一般的,Windows,程序一样，启动,MATLAB,系统有,3,种常见方法：,(1),使用,Windows“,开始”菜单。,(2),运行,MATLAB,系统启动程序,matlab.exe,。,(3),利用快捷方式。,启动,MATLAB,后，将进入,MATLAB 6.5,集成环境。,MATLAB 6.5,集成环境包括,MATLAB,主窗口、命令窗口,(Command Window),、工作空间窗口,(Workspace),、命令历史窗口,(Command History),、当前目录窗口,(Current Directory),和启动平台窗口,(Launch Pad),。,当,MATLAB,安装完毕并首次启动时，展现在屏幕上的界面为,MATLAB,的默认界面，如右图所示。,2,MATLAB,系统的退出要退出,MATLAB,系统，也有,3,种常见方法：,(1),在,MATLAB,主窗口,File,菜单中选择,Exit MATLAB,命令。,(2),在,MATLAB,命令窗口输入,Exit,或,Quit,命令。,(3),单击,MATLAB,主窗口的“关闭”按钮。,1.3.2,主窗口,MATLAB,主窗口是,MATLAB,的主要工作界面。主窗口除了嵌入一些子窗口外，还主要包括菜单栏和工具栏。,1,菜单栏在,MATLAB 6.5,主窗口的菜单栏，共包含,File,、,Edit,、,View,、,Web,、,Window,和,Help 6,个菜单项。,(1)File,菜单项：,File,菜单项实现有关文件的操作。,(2)Edit,菜单项：,Edit,菜单项用于命令窗口的编辑操作。,(3)View,菜单项：,View,菜单项用于设置,MATLAB,集成环境的显示方式。,(4)Web,菜单项：,Web,菜单项用于设置,MATLAB,的,Web,操作。,(5)Window,菜单项：主窗口菜单栏上的,Window,菜单，只包含一个子菜单,Close all,，用于关闭所有打开的编辑器窗口，包括,M-file,、,Figure,、,Model,和,GUI,窗口。,(6)Help,菜单项：,Help,菜单项用于提供帮助信息。,2,工具栏,MATLAB 6.5,主窗口的工具栏共提供了,10,个命令按钮。这些命令按钮均有对应的菜单命令，但比菜单命令使用起来更快捷、方便。,1.3.3,命令窗口,命令窗口是,MATLAB,的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。,MATLAB,命令窗口中的“,”,为命令提示符，表示,MATLAB,正在处于准备状态。在命令提示符后键入命令并按下回车键后，,MATLAB,就会解释执行所输入的命令，并在命令后面给出计算结果。,一般来说，一个命令行输入一条命令，命令行以回车结束。但一个命令行也可以输入若干条命令，各命令之间以逗号分隔，若前一命令后带有分号，则逗号可以省略。例如,p=15,m=35p=15;m=35,如果一个命令行很长，一个物理行之内写不下，可以在第一个物理行之后加上,3,个小黑点并按下回车键，然后接着下一个物理行继续写命令的其他部分。,3,个小黑点称为续行符，即把下面的物理行看作该行的逻辑继续。在,MATLAB,里，有很多的控制键和方向键可用于命令行的编辑。,1.3.4,工作空间窗口,工作空间是,MATLAB,用于存储各种变量和结果的内存空间。在该窗口中显示工作空间中所有变量的名称、大小、字节数和变量类型说明，可对变量进行观察、编辑、保存和删除。,1.3.5,当前目录窗口和搜索路径,1,当前目录窗口当前目录是指,MATLAB,运行文件时的工作目录，只有在当前目录或搜索路径下的文件、函数可以被运行或调用。在当前目录窗口中可以显示或改变当前目录，还可以显示当前目录下的文件并提供搜索功能。将用户目录设置成当前目录也可使用,cd,命令。例如，将用户目录,c:mydir,设置为当前目录，可在命令窗口输入命令：,cd c:mydir,2,MATLAB,的搜索路径当用户在,MATLAB,命令窗口输入一条命令后，,MATLAB,按照一定次序寻找相关的文件。基本的搜索过程是：,(1),检查该命令是不是一个变量。,(2),检查该命令是不是一个内部函数。,(3),检查该命令是否当前目录下的,M,文件。,(4),检查该命令是否,MATLAB,搜索路径中其他目录下的,M,文件。,用户可以将自己的工作目录列入,MATLAB,搜索路径，从而将用户目录纳入,MATLAB,系统统一管理。设置搜索路径的方法有：,(1),用,path,命令设置搜索路径。例如，将用户目录,c:mydir,加到搜索路径下，可在命令窗口输入命令：,path(path,c:mydir),(2),用对话框设置搜索路径 在,MATLAB,的,File,菜单中选,Set Path,命令或在命令窗口执行,pathtool,命令，将出现搜索路径设置对话框。通过,Add Folder,或,Add with Subfolder,命令按钮将指定路径添加到搜索路径列表中。在修改完搜索路径后，则需要保存搜索路径,。,1.3.6,命令历史记录窗口,在默认设置下，历史记录窗口中会自动保留自安装起所有用过的命令的历史记录，并且还标明了使用时间，从而方便用户查询。而且，通过双击命令可进行历史命令的再运行。如果要清除这些历史记录，可以选择,Edit,菜单中的,Clear Command History,命令。,1.3.7,启动平台窗口和,Start,按钮,MATLAB 6.5,的启动平台窗口可以帮助用户方便地打开和调用,MATLAB,的各种程序、函数和帮助文件。,MATLAB 6.5,主窗口左下角还有一个,Start,按钮，单击该按钮会弹出一个菜单，选择其中的命令可以执行,MATLAB,产品的各种工具，并且可以查阅,MATLAB,包含的各种资源。,1.4 MATLAB,帮助系统,1.4.1,帮助窗口,进入帮助窗口可以通过以下,3,种方法：,(1),单击,MATLAB,主窗口工具栏中的,Help,按钮。,(2),在命令窗口中输入,helpwin,、,helpdesk,或,doc,。,(3),选择,Help,菜单中的“,MATLAB Help”,选项。,1.4.2,帮助命令,MATLAB,帮助命令包括,help,、,lookfor,以及模糊查询。,1,help,命令在,MATLAB 6.5,命令窗口中直接输入,help,命令将会显示当前帮助系统中所包含的所有项目，即搜索路径中所有的目录名称。同样，可以通过,help,加函数名来显示该函数的帮助说明。,2,lookfor,命令,help,命令只搜索出那些关键字完全匹配的结果，,lookfor,命令对搜索范围内的,M,文件进行关键字搜索，条件比较宽松。,lookfor,命令只对,M,文件的第一行进行关键字搜索。若在,lookfor,命令加上,-all,选项，则可对,M,文件进行全文搜索。,3,模糊查询,MATLAB 6.0,以上的版本提供了一种类似模糊查询的命令查询方法，用户只需要输入命令的前几个字母，然后按,Tab,键，系统就会列出所有以这几个字母开头的命令。,1.4.3,演示系统,在帮助窗口中选择演示系统,(Demos),选项卡，然后在其中选择相应的演示模块，或者在命令窗口输入,Demos,，或者选择主窗口,Help,菜单中的,Demos,子菜单，打开演示系统。,1.4.4,远程帮助系统,在,MathWorks,公司的主页,(),上可以找到很多有用的信息，国内的一些网站也有丰富的信息资源。,2.1,变量和数据操作,2.2 MATLAB,矩阵,2.3 MATLAB,运算,2.4,矩阵分析,2.5,字符串,2.7,结构数据和单元数据,2 MATLAB,矩阵及其运算,2.1,变量和数据操作,2.1.1,变量与赋值,1,变量命名在,MATLAB 6.5,中，变量名是以字母开头，后接字母、数字或下划线的字符序列，最多,63,个字符。在,MATLAB,中，变量名区分字母的大小写。,2,赋值语句,(1),变量,=,表达式,(2),表达式其中表达式是用运算符将有关运算量连接起来的式子，其结果是一个矩阵。,例,2-1,计算表达式的值，并显示计算结果。在,MATLAB,命令窗口输入命令：,x=1+2i;y=3-sqrt(17);z=(cos(abs(x+y)-sin(78*pi/180)/(x+abs(y),其中,pi,和,i,都是,MATLAB,预先定义的变量，分别代表代表圆周率,和虚数单位。输出结果是：,z=-0.3488+0.3286i,2.1.2,预定义变量,在,MATLAB,工作空间中，还驻留几个由系统本身定义的变量。例如，用,pi,表示圆周率,的近似值，用,i,，,j,表示虚数单位。预定义变量有特定的含义，在使用时，应尽量避免对这些变量重新赋值。,2.1.3,内存变量的管理,1,内存变量的删除与修改,MATLAB,工作空间窗口专门用于内存变量的管理。在工作空间窗口中可以显示所有内存变量的属性。当选中某些变量后，再单击,Delete,按钮，就能删除这些变量。当选中某些变量后，再单击,Open,按钮，将进入变量编辑器。通过变量编辑器可以直接观察变量中的具体元素，也可修改变量中的具体元素。,clear,命令用于删除,MATLAB,工作空间中的变量。,who,和,whos,这两个命令用于显示在,MATLAB,工作空间中已经驻留的变量名清单。,who,命令只显示出驻留变量的名称，,whos,在给出变量名的同时，还给出它们的大小、所占字节数及数据类型等信息。,2,内存变量文件利用,MAT,文件可以把当前,MATLAB,工作空间中的一些有用变量长久地保留下来，扩展名是,.mat,。,MAT,文件的生成和装入由,save,和,load,命令来完成。常用格式为：,save,文件名,变量名表,-append-asciiload,文件名,变量名表,-ascii,其中，文件名可以带路径，但不需带扩展名,.mat,，命令隐含一定对,.mat,文件进行操作。变量名表中的变量个数不限，只要内存或文件中存在即可，变量名之间以空格分隔。当变量名表省略时，保存或装入全部变量。,-ascii,选项使文件以,ASCII,格式处理，省略该选项时文件将以二进制格式处理。,save,命令中的,-append,选项控制将变量追加到,MAT,文件中。,2.1.4 MATLAB,常用数学函数,MATLAB,提供了许多数学函数，函数的自变量规定为矩阵变量，运算法则是将函数逐项作用于矩阵的元素上，因而运算的结果是一个与自变量同维数的矩阵。函数使用说明：,(1),三角函数以弧度为单位计算。,(2)abs,函数可以求实数的绝对值、复数的模、字符串的,ASCII,码值。,(3),用于取整的函数有,fix,、,floor,、,ceil,、,round,，要注意它们的区别。,(4)rem,与,mod,函数的区别。,rem(x,y),和,mod(x,y),要求,x,y,必须为相同大小的实矩阵或为标量。,2.1.5,数据的输出格式,MATLAB,用十进制数表示一个常数，具体可采用日常记数法和科学记数法两种表示方法。在一般情况下，,MATLAB,内部每一个数据元素都是用双精度数来表示和存储的。数据输出时用户可以用,format,命令设置或改变数据输出格式。,format,命令的格式为：,format,格式符其中格式符决定数据的输出格式,2.2 MATLAB,矩阵,2.2.1,矩阵的建立,1,直接输入法 最简单的建立矩阵的方法是从键盘直接输入矩阵的元素。具体方法如下：将矩阵的元素用方括号括起来，按矩阵行的顺序输入各元素，同一行的各元素之间用空格或逗号分隔，不同行的元素之间用分号分隔。,2,利用,M,文件建立矩阵 对于比较大且比较复杂的矩阵，可以为它专门建立一个,M,文件。下面通过一个简单例子来说明如何利用,M,文件创建矩阵。,例,2-2,利用,M,文件建立,MYMAT,矩阵。,(1),启动有关编辑程序或,MATLAB,文本编辑器，并输入待建矩阵：,(2),把输入的内容以纯文本方式存盘,(,设文件名为,mymatrix.m),。,(3),在,MATLAB,命令窗口中输入,mymatrix,，即运行该,M,文件，就会自动建立一个名为,MYMAT,的矩阵，可供以后使用。,3,利用冒号表达式建立一个向量 冒号表达式可以产生一个行向量，一般格式是：,e1:e2:e3,其中,e1,为初始值，,e2,为步长，,e3,为终止值。在,MATLAB,中，还可以用,linspace,函数产生行向量。其调用格式为：,linspace(a,b,n),其中,a,和,b,是生成向量的第一个和最后一个元素，,n,是元素总数。显然，,linspace(a,b,n),与,a:(b-a)/(n-1):b,等价。,4,建立大矩阵大矩阵可由方括号中的小矩阵或向量建立起来。,2.2.2,矩阵的拆分,1,矩阵元素 通过下标引用矩阵的元素，例如,A(3,2)=200,采用矩阵元素的序号来引用矩阵元素。矩阵元素的序号就是相应元素在内存中的排列顺序。在,MATLAB,中，矩阵元素按列存储，先第一列，再第二列，依次类推。例如,A=1,2,3;4,5,6;A(3)ans=2,显然，序号,(Index),与下标,(Subscript),是一一对应的，以,mn,矩阵,A,为例，矩阵元素,A(i,j),的序号为,(j-1)*m+i,。其相互转换关系也可利用,sub2ind,和,ind2sub,函数求得。,2,矩阵拆分,(1),利用冒号表达式获得子矩阵 ,A(:,j),表示取,A,矩阵的第,j,列全部元素；,A(i,:),表示,A,矩阵第,i,行的全部元素；,A(i,j),表示取,A,矩阵第,i,行、第,j,列的元素。,A(i:i+m,:),表示取,A,矩阵第,i,i+m,行的全部元素；,A(:,k:k+m),表示取,A,矩阵第,k,k+m,列的全部元素，,A(i:i+m,k:k+m),表示取,A,矩阵第,i,i+m,行内，并在第,k,k+m,列中的所有元素。此外，还可利用一般向量和,end,运算符来表示矩阵下标，从而获得子矩阵。,end,表示某一维的末尾元素下标。,(2),利用空矩阵删除矩阵的元素 在,MATLAB,中，定义,为空矩阵。给变量,X,赋空矩阵的语句为,X=,。注意，,X=,与,clear X,不同，,clear,是将,X,从工作空间中删除，而空矩阵则存在于工作空间中，只是维数为,0,。,2.2.3,特殊矩阵,1,通用的特殊矩阵常用的产生通用特殊矩阵的函数有：,zeros,：产生全,0,矩阵,(,零矩阵,),。,ones,：产生全,1,矩阵,(,幺矩阵,),。,eye,：产生单位矩阵。,rand,：产生,0,1,间均匀分布的随机矩阵。,randn,：产生均值为,0,，方差为,1,的标准正态分布随机矩阵。,例,2-3,分别建立,33,、,32,和与矩阵,A,同样大小的零矩阵。,(1),建立一个,33,零矩阵。,zeros(3)(2),建立一个,32,零矩阵。,zeros(3,2)(3),设,A,为,23,矩阵，则可以用,zeros(size(A),建立一个与矩阵,A,同样大小零矩阵。,A=1 2 3;4 5 6;%,产生一个,23,阶矩阵,Azeros(size(A)%,产生一个与矩阵,A,同样大小的零矩阵,例,2-4,建立随机矩阵：,(1),在区间,20,50,内均匀分布的,5,阶随机矩阵。,(2),均值为,0.6,、方差为,0.1,的,5,阶正态分布随机矩阵。命令如下：,x=20+(50-20)*rand(5)y=0.6+sqrt(0.1)*randn(5),此外，常用的函数还有,reshape(A,m,n),，它在矩阵总元素保持不变的前提下，将矩阵,A,重新排成,mn,的二维矩阵。,2,用于专门学科的特殊矩阵,(1),魔方矩阵魔方矩阵有一个有趣的性质，其每行、每列及两条对角线上的元素和都相等。对于,n,阶魔方阵，其元素由,1,2,3,n2,共,n2,个整数组成。,MATLAB,提供了求魔方矩阵的函数,magic(n),，其功能是生成一个,n,阶魔方阵。,例,2-5,将,101125,等,25,个数填入一个,5,行,5,列的表格中，使其每行每列及对角线的和均为,565,。,M=100+magic(5),(2),范得蒙矩阵范得蒙,(Vandermonde),矩阵最后一列全为,1,，倒数第二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。可以用一个指定向量生成一个范得蒙矩阵。在,MATLAB,中，函数,vander(V),生成以向量,V,为基础向量的范得蒙矩阵。例如，,A=vander(1;2;3;5),即可得到上述范得蒙矩阵。,(3),希尔伯特矩阵在,MATLAB,中，生成希尔伯特矩阵的函数是,hilb(n),。使用一般方法求逆会因为原始数据的微小扰动而产生不可靠的计算结果。,MATLAB,中，有一个专门求希尔伯特矩阵的逆的函数,invhilb(n),，其功能是求,n,阶的希尔伯特矩阵的逆矩阵。,例,2-6,求,4,阶希尔伯特矩阵及其逆矩阵。命令如下：,format rat%,以有理形式输出,H=hilb(4)H=invhilb(4),(4),托普利兹矩阵托普利兹,(Toeplitz),矩阵除第一行第一列外，其他每个元素都与左上角的元素相同。生成托普利兹矩阵的函数是,toeplitz(x,y),，它生成一个以,x,为第一列，,y,为第一行的托普利兹矩阵。这里,x,y,均为向量，两者不必等长。,toeplitz(x),用向量,x,生成一个对称的托普利兹矩阵。例如,T=toeplitz(1:6),(5),伴随矩阵,MATLAB,生成伴随矩阵的函数是,compan(p),，其中,p,是一个多项式的系数向量，高次幂系数排在前，低次幂排在后。例如，为了求多项式的,x3-7x+6,的伴随矩阵，可使用命令：,p=1,0,-7,6;compan(p),(6),帕斯卡矩阵我们知道，二次项,(x+y)n,展开后的系数随,n,的增大组成一个三角形表，称为杨辉三角形。由杨辉三角形表组成的矩阵称为帕斯卡,(Pascal),矩阵。函数,pascal(n),生成一个,n,阶帕斯卡矩阵。,例,2-7,求,(x+y),5,的展开式。在,MATLAB,命令窗口，输入命令：,pascal(6),矩阵次对角线上的元素,1,5,10,10,5,1,即为展开式的系数。,2.3 MATLAB,运算,2.3.1,算术运算,1,基本算术运算,MATLAB,的基本算术运算有：,(,加,),、,(,减,),、*,(,乘,),、,/(,右除,),、,(,左除,),、,(,乘方,),。注意，运算是在矩阵意义下进行的，单个数据的算术运算只是一种特例。,(1),矩阵加减运算 假定有两个矩阵,A,和,B,，则可以由,A+B,和,A-B,实现矩阵的加减运算。运算规则是：若,A,和,B,矩阵的维数相同，则可以执行矩阵的加减运算，,A,和,B,矩阵的相应元素相加减。如果,A,与,B,的维数不相同，则,MATLAB,将给出错误信息，提示用户两个矩阵的维数不匹配。,(2),矩阵乘法 假定有两个矩阵,A,和,B,，若,A,为,mn,矩阵，,B,为,np,矩阵，则,C=A*B,为,mp,矩阵。,(3),矩阵除法在,MATLAB,中，有两种矩阵除法运算：,和,/,，分别表示左除和右除。如果,A,矩阵是非奇异方阵，则,AB,和,B/A,运算可以实现。,AB,等效于,A,的逆左乘,B,矩阵，也就是,inv(A)*B,，而,B/A,等效于,A,矩阵的逆右乘,B,矩阵，也就是,B*inv(A),。对于含有标量的运算，两种除法运算的结果相同，如,3/4,和,43,有相同的值，都等于,0.75,。又如，设,a=10.5,25,，则,a/5=5a=2.1000 5.0000,。对于矩阵来说，左除和右除表示两种不同的除数矩阵和被除数矩阵的关系。对于矩阵运算，一般,ABB/A,。,(4),矩阵的乘方 一个矩阵的乘方运算可以表示成,Ax,，要求,A,为方阵，,x,为标量。,2,点运算 在,MATLAB,中，有一种特殊的运算，因为其运算符是在有关算术运算符前面加点，所以叫点运算。点运算符有,.*,、,./,、,.,和,.,。两矩阵进行点运算是指它们的对应元素进行相关运算，要求两矩阵的维参数相同。,2.3.2,关系运算,MATLAB,提供了,6,种关系运算符：,(,小于,),、,(,大于,),、,=(,大于或等于,),、,=(,等于,),、,=(,不等于,),。它们的含义不难理解，但要注意其书写方法与数学中的不等式符号不尽相同。,关系运算符的运算法则为：,(1),当两个比较量是标量时，直接比较两数的大小。若关系成立，关系表达式结果为,1,，否则为,0,。,(2),当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由,0,或,1,组成。,(3),当参与比较的一个是标量，而另一个是矩阵时，则把标量与矩阵的每一个元素按标量关系运算规则逐个比较，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由,0,或,1,组成。,例,2-8,产生,5,阶随机方阵,A,，其元素为,10,90,区间的随机整数，然后判断,A,的元素是否能被,3,整除。,(1),生成,5,阶随机方阵,A,。,A=fix(90-10+1)*rand(5)+10)(2),判断,A,的元素是否可以被,3,整除。,P=rem(A,3)=0,其中，,rem(A,3),是矩阵,A,的每个元素除以,3,的余数矩阵。此时，,0,被扩展为与,A,同维数的零矩阵，,P,是进行等于,(=),比较的结果矩阵。,2.3.3,逻辑运算,MATLAB,提供了,3,种逻辑运算符：,&(,与,),、,|(,或,),和,(,非,),。逻辑运算的运算法则为：,(1),在逻辑运算中，确认非零元素为真，用,1,表示，零元素为假，用,0,表示。,(2),设参与逻辑运算的是两个标量,a,和,b,，那么，,a&b a,b,全为非零时，运算结果为,1,，否则为,0,。,a|b a,b,中只要有一个非零，运算结果为,1,。,a,当,a,是零时，运算结果为,1,；当,a,非零时，运算结果为,0,。,(3),若参与逻辑运算的是两个同维矩阵，那么运算将对矩阵相同位置上的元素按标量规则逐个进行。最终运算结果是一个与原矩阵同维的矩阵，其元素由,1,或,0,组成。,(4),若参与逻辑运算的一个是标量，一个是矩阵，那么运算将在标量与矩阵中的每个元素之间按标量规则逐个进行。最终运算结果是一个与矩阵同维的矩阵，其元素由,1,或,0,组成。,(5),逻辑非是单目运算符，也服从矩阵运算规则。,(6),在算术、关系、逻辑运算中，算术运算优先级最高，逻辑运算优先级最低。,例,2-9,建立矩阵,A,，然后找出大于,4,的元素的位置。,(1),建立矩阵,A,。,A=4,-65,-54,0,6;56,0,67,-45,0(2),找出大于,4,的元素的位置。,find(A4),2.4,矩阵分析,2.4.1,对角阵与三角阵,1,对角阵只有对角线上有非,0,元素的矩阵称为对角矩阵，对角线上的元素相等的对角矩阵称为数量矩阵，对角线上的元素都为,1,的对角矩阵称为单位矩阵。,(1),提取矩阵的对角线元素设,A,为,mn,矩阵，,diag(A),函数用于提取矩阵,A,主对角线元素，产生一个具有,min(m,n),个元素的列向量。,diag(A),函数还有一种形式,diag(A,k),，其功能是提取第,k,条对角线的元素。,(2),构造对角矩阵设,V,为具有,m,个元素的向量，,diag(V),将产生一个,mm,对角矩阵，其主对角线元素即为向量,V,的元素。,diag(V),函数也有另一种形式,diag(V,k),，其功能是产生一个,nn(n=m+),对角阵，其第,k,条对角线的元素即为向量,V,的元素。,例,2-10,先建立,55,矩阵,A,，然后将,A,的第一行元素乘以,1,，第二行乘以,2,，,，第五行乘以,5,。,A=17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;.11,18,25,2,19;D=diag(1:5);D*A%,用,D,左乘,A,，对,A,的每行乘以一个指定常数,2,三角阵三角阵又进一步分为上三角阵和下三角阵，所谓上三角阵，即矩阵的对角线以下的元素全为,0,的一种矩阵，而下三角阵则是对角线以上的元素全为,0,的一种矩阵。,(1),上三角矩阵求矩阵,A,的上三角阵的,MATLAB,函数是,triu(A),。,triu(A),函数也有另一种形式,triu(A,k),，其功能是求矩阵,A,的第,k,条对角线以上的元素。例如，提取矩阵,A,的第,2,条对角线以上的元素，形成新的矩阵,B,。,(2),下三角矩阵在,MATLAB,中，提取矩阵,A,的下三角矩阵的函数是,tril(A),和,tril(A,k),，其用法与提取上三角矩阵的函数,triu(A),和,triu(A,k),完全相同。,2.4.2,矩阵的转置与旋转,1,矩阵的转置转置运算符是单撇号,(),。,2,矩阵的旋转利用函数,rot90(A,k),将矩阵,A,旋转,90,的,k,倍，当,k,为,1,时可省略。,3,矩阵的左右翻转对矩阵实施左右翻转是将原矩阵的第一列和最后一列调换，第二列和倒数第二列调换，,，依次类推。,MATLAB,对矩阵,A,实施左右翻转的函数是,fliplr(A),。,4,矩阵的上下翻转,MATLAB,对矩阵,A,实施上下翻转的函数是,flipud(A),。,2.4.3,矩阵的逆与伪逆,1,矩阵的逆对于一个方阵,A,，如果存在一个与其同阶的方阵,B,，使得：,AB=BA=I(I,为单位矩阵,),则称,B,为,A,的逆矩阵，当然，,A,也是,B,的逆矩阵。求一个矩阵的逆是一件非常烦琐的工作，容易出错，但在,MATLAB,中，求一个矩阵的逆非常容易。求方阵,A,的逆矩阵可调用函数,inv(A),。例,2-11,用求逆矩阵的方法解线性方程组。,Ax=b,其解为：,x=A,-1,b,2,矩阵的伪逆如果矩阵,A,不是一个方阵，或者,A,是一个非满秩的方阵时，矩阵,A,没有逆矩阵，但可以找到一个与,A,的转置矩阵,A,同型的矩阵,B,，使得：,ABA=ABAB=B,此时称矩阵,B,为矩阵,A,的伪逆，也称为广义逆矩阵。在,MATLAB,中，求一个矩阵伪逆的函数是,pinv(A),。,2.4.4,方阵的行列式把一个方阵看作一个行列式，并对其按行列式的规则求值，这个值就称为矩阵所对应的行列式的值。在,MATLAB,中，求方阵,A,所对应的行列式的值的函数是,det(A),。,2.4.5,矩阵的秩与迹,1,矩阵的秩矩阵线性无关的行数与列数称为矩阵的秩。在,MATLAB,中，求矩阵秩的函数是,rank(A),。,2,矩阵的迹矩阵的迹等于矩阵的对角线元素之和，也等于矩阵的特征值之和。在,MATLAB,中，求矩阵的迹的函数是,trace(A),。,2.4.6,向量和矩阵的范数矩阵或向量的范数用来度量矩阵或向量在某种意义下的长度。范数有多种方法定义，其定义不同，范数值也就不同。,1,向量的,3,种常用范数及其计算函数在,MATLAB,中，求向量范数的函数为：,(1)norm(V),或,norm(V,2),：计算向量,V,的,2,范数。,(2)norm(V,1),：计算向量,V,的,1,范数。,(3)norm(V,inf),：计算向量,V,的,范数。,2,矩阵的范数及其计算函数,MATLAB,提供了求,3,种矩阵范数的函数，其函数调用格式与求向量的范数的函数完全相同。,2.4.7,矩阵的条件数在,MATLAB,中，计算矩阵,A,的,3,种条件数的函数是：,(1)cond(A,1),计算,A,的,1,范数下的条件数。,(2)cond(A),或,cond(A,2),计算,A,的,2,范数数下的条件数。,(3)cond(A,inf),计算,A,的,范数下的条件数。,2.4.8,矩阵的特征值与特征向量在,MATLAB,中，计算矩阵,A,的特征值和特征向量的函数是,eig(A),，常用的调用格式有,3,种：,(1)E=eig(A),：求矩阵,A,的全部特征值，构成向量,E,。,(2)V,D=eig(A),：求矩阵,A,的全部特征值，构成对角阵,D,，并求,A,的特征向量构成,V,的列向量。,(3)V,D=eig(A,nobalance),：与第,2,种格式类似，但第,2,种格式中先对,A,作相似变换后求矩阵,A,的特征值和特征向量，而格式,3,直接求矩阵,A,的特征值和特征向量。,2.5,字符串 在,MATLAB,中，字符串是用单撇号括起来的字符序列。,MATLAB,将字符串当作一个行向量，每个元素对应一个字符，其标识方法和数值向量相同。也可以建立多行字符串矩阵。,字符串是以,ASCII,码形式存储的。,abs,和,double,函数都可以用来获取字符串矩阵所对应的,ASCII,码数值矩阵。相反，,char,函数可以把,ASCII,码矩阵转换为字符串矩阵。,例,2-13,建立一个字符串向量，然后对该向量做如下处理：,(1),取第,1,5,个字符组成的子字符串。,(2),将字符串倒过来重新排列。,(3),将字符串中的小写字母变成相应的大写字母，其余字符不变。,(4),统计字符串中小写字母的个数。,命令如下：,ch=ABc123d4e56Fg9;subch=ch(1:5)%,取子字符串,revch=ch(end:-1:1)%,将字符串倒排,k=find(ch=a%,找小写字母的位置,ch(k)=ch(k)-(a-A);%,将小写字母变成相应的大写字母,char(ch)length(k)%,统计小写字母的个数,与字符串有关的另一个重要函数是,eval,，其调用格式为：,eval(t),其中,t,为字符串。它的作用是把字符串的内容作为对应的,MATLAB,语句来执行。,2.7,结构数据和单元数据,2.7.1,结构数据,1,结构矩阵的建立与引用结构矩阵的元素可以是不同的数据类型，它能将一组具有不同属性的数据纳入到一个统一的变量名下进行管理。建立一个结构矩阵可采用给结构成员赋值的办法。具体格式为：结构矩阵名,.,成员名,=,表达式其中表达式应理解为矩阵表达式。,2,结构成员的修改可以根据需要增加或删除结构的成员。例如要给结构矩阵,a,增加一个成员,x4,，可给,a,中任意一个元素增加成员,x4,：,a(1).x4=410075;,但其他成员均为空矩阵，可以使用赋值语句给它赋确定的值。要删除结构的成员，则可以使用,rmfield,函数来完成。例如，删除成员,x4,：,a=rmfield(a,x4);3,关于结构的函数除了一般的结构数据的操作外，,MATLAB,还提供了部分函数来进行结构矩阵的操作。,2.7.2,单元数据,1,单元矩阵的建立与引用建立单元矩阵和一般矩阵相似，只是矩阵元素用大括号括起来。可以用带有大括号下标的形式引用单元矩阵元素。例如,b3,3,。单元矩阵的元素可以是结构或单元数据。可以使用,celldisp,函数来显示整个单元矩阵，如,celldisp(b),。还可以删除单元矩阵中的某个元素。,2,关于单元的函数,MATLAB,还提供了部分函数用于单元的操作。,1 MATLAB,操作基础,2 MATLAB,矩阵及其运算,3 MATLAB,程序设计,4 MATLAB,工,具箱,3.1 M,文件,3.2,程序控制结构,3.3,函数文件,3.4,程序举例,3.5,程序调试,3 MATLAB,程序设计,3.1 M,文件,3.1.1 M,文件概述,用,MATLAB,语言编写的程序，称为,M,文件。,M,文件可以根据调用方式的不同分为两类：命令文件,(Script File),和函数文件,(Function File),。,例,3-1,分别建立命令文件和函数文件，将华氏温度,f,转换为摄氏温度,c,。,程序,1,：,首先建立命令文件并以文件名,f2c.m,存盘。,clear;%,清除工作空间中的变量,f=input(Input Fahrenheit temperature,：,);,c=5*(f-32)/9,然后在,MATLAB,的命令窗口中输入,f2c,，将会执行该命令文件，执行情况为：,Input Fahrenheit temperature,：,73,c=,22.7778,程序,2,：,首先建立函数文件,f2c.m,。,function c=f2c(f),c=5*(f-32)/9,然后在,MATLAB,的命令窗口调用该函数文件。,clear;,y=input(Input Fahrenheit temperature,：,);,x=f2c(y),输出情况为：,Input Fahrenh