1、推荐参考书推荐参考书推荐参考书推荐参考书 一个可视化的计算程序,广泛使用于从个人计算机到超级计算机范围内的各种计算机上 包括命令控制、可编程,上百个预先定义命令和函数 有许多强有力的命令,能完成大量的高级矩阵处理 强有力的二维、三维图形工具 能与其他程序一起使用 25个(不断增加中)不同的工具箱应用于特殊的应用领域 工业研究与开发的有力工具 数学教学,尤其线代,数值分析,科学计算方面的教研工具 电子学,控制理论,物理学等工程科学方面的教研工具 经济学,化学和生物学等有计算问题的所有领域中的教学与研究 名字取自矩阵实验室(matrix laboratory)不是万能的解决工具不是最高性能的编程语
2、言受计算条件限制,不能解决超大型实际问题不能解决工具箱之外的问题种类-需要编写接口、算法甚至工具箱 数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与方针 数字图像处理 数字信号处理 通讯系统设计与仿真 财务与金融工程绪论基础知识数学运算基本编程数据显示 及存取符号计算数值计算图形用户界面GUI Simulink仿真外部接口第一章第一章 概述概述 MATLAB将计算、可视化和编程等功能集于一个易于使用的环境,并遇有如下特点:功能强大;简单易学;编程效率高。1.1 MATLAB简介及安装1MATLAB的发展历程 2MATLAB的特点 3MATLAB的安装 1980年,美国新墨西哥大学计算
3、机科学 系主任Cleve 着手编写供学生使用的子程 序接口程序,取名为MATLAB;1984年,推出了MATLAB第一个商业版 本;1992年,推出MATLAB 4.0版;1997年,推出MATLAB 5.0版;2000年,推出MATLAB 6.0版;2004年,推出MATLAB 7.0版;2008年,推出MATLAB 7.6版;2012年,推出MATLAB R2012A版 开发环境 编程 数值处理 图形化 图形用户界面 文件I/O和外部应用程序接口MATLAB对PC系统的要求 1.2 MATLAB的目录结构 续表1.3 MATLAB的工作环境1菜单和工具栏 2命令窗口 3历史命令窗口 4当前
4、工作目录窗口 5工作空间窗口 本节介绍MATLAB的工作界面和基本的操作方法。MATLAB的工作界面主要包括:菜单;工具栏;命令窗口;历史命令窗口;当前工作目录窗口;工作空间窗口。这里只简单介绍默认情况下的菜单和工具栏。【File】菜单主要用于对文件的处理。【Edit】菜单主要用于复制、粘贴等操 作,与一般Windows程序的类似,在此 不作详细介绍。【Debug】菜单用于调试程序。【Desktop】菜单用于设置主窗口中需要 打开的窗口。【Window】菜单列出当前所有打开的窗 口。【Help】菜单用于选择打开不同的帮助系 统。当用户单击“Current Directory”窗口时,使得该窗口
5、成为活动窗口,同时增加一个如下图所示的菜单【View】,用于设置如何显示当前目录下的文件。当用户单击“Workspace”窗口时,使得该窗口成为活动窗口,同时增加如下图所示的菜单【View】和【Graphics】。命令窗口是MATLAB的主要交互窗口,用于输入命令并显示除图形以外的所有执行结果。MATLAB命令窗口中的“”为运算提示 符,表MATLAB处于准备状态。当在提示符后输入一段程序或一段运算式 后【Enter】键,MATLAB会给出计算结 果,并再次进入准备状态(所得结果将被 保存在工作空间窗口中)。单击命令窗口右上角的“”按钮,可以使 命令窗口脱离主窗口而成为一个独立的窗 口。主要用
6、于记录所有执行过的命令;保留自安装后所有使用过命令的历史记 录,并标明使用时间;可以通过用鼠标双击某一历史命令来重新 执行该命令;可以成为一个独立的窗口。通过上下文菜单,可删除或粘贴历史记录;可为选中的表达式或命令创建一个M文件;可为表达式或命令创建快捷按钮。当前工作目录是指MATLAB运行文件时的目录。只有在当前工作目录或搜索路径下的文件、函数可以被运行或调用。如下图所示。工作空间窗口将显示目前内存中所有的MATLAB变量的变量名、数据结构、字节数以及类型等信息,如下图所示。1.4 MATLAB的通用命令 1常用命令 2输入内容的编辑 3标点 4搜索路径与扩展 通用命令是MATLAB中经常使
7、用的一组命令,这些命令可以用来管理目录、函数、变量、工作空间、文件和窗口等。下面对这些命令进行介绍。常用命令的功能 在命令窗口中,MATLAB提供了控制光标位置和进行简单编辑的键盘按键,部分常用的键盘按键及其功能如下表所示。在MATLAB中,一些标点符号被赋予了特殊的功能,如下表所示。用户自己书写的函数有可能并没有保存在搜索路径下。要解决这个问题,只需把程序所在的目录扩展成MATLAB的搜索路径即可。(1)查看MATLAB的搜索路径 选择MATLAB主窗口中的【File】|【Set Path】菜单,进入【设置搜索路径】对话框。通过该对话框可为MATLAB添加或删除搜索路径。在命令窗口中输入pa
8、th或genpath可得到MATLAB的所有搜索路径,具体代码如下:(2)设置MATLAB的搜索路径 方法一:在MATLAB命令窗口中输入 editpath或pathtool命令或通过【File】|【Set Path】菜单,进入“设置搜索路 径”对话框,通过该对话框编辑搜索路 径。方法二:在命令窗口执行“path(path,D:Study)”,然后通过“设置搜索路 径”对话框查看“D:Study”是否在搜索路 径中。方法三:在命令窗口执行“addpath D:Study-end”,将新的目录加到整个 搜索路径的末尾。如果将end改为 begin,可以将新的目录加到整个搜索路 径的开始。1.5
9、MATLAB的帮助系统 1命令窗口查询帮助系统 2联机帮助系统 3联机演示系统 MATLAB为用户提供了非常完善的帮助系统。下面分别介绍MATLAB的三类帮助系统。常见的帮助命令 用户可以通过下述方法进入MATLAB的联机帮助系统。直接单击MATLAB主窗口中的“”按钮;选中Help菜单的前4项中的任意一项;在命令窗口中执行helpwin、helpdesk 或doc。可以通过以下方式打开联机演示统。选择MATLAB主窗口菜单的【Help】|【Demos】选项;在命令窗口输入demos;直接在帮助页面上选择Demos页。020406080100-40-2002040Magnitude(dB)Fr
10、equency(Hertz)00.20.40.60.81-1-0.500.51WaveformTime(Seconds)Click and drag waveform to changefundamental frequency and amplitude1.6 MATLAB示例 下面以一个简单的例子展示如何使用MATLAB进行简单的数值计算。习 题1简述MATLAB的主要功能。2在命令窗口输入“w=3+2”,然后依次使用clear和clc命 令,分别观察命令窗口、工作空间窗口和历史命令窗口 的变化。3将硬盘上一已有目录,加入到搜索路径,并将其设置为 当前工作目录。4通过命令窗口,查询函数si
11、n()的用法。5通过联机帮助系统,查询函数inv()的用法。6通过联机演示系统,查询并运行“Control Systems Toolboxes”下的“Case Studies”中的“Yaw Damper for a 747 Aircraft”演示程序。7在命令窗口依次执行“w=5;”、“p=2*w”和“q=p+w”。8在命令窗口同时执行下述代码:w=5;p=2*w q=p+w第二章第二章 基础知识基础知识 本章着重介绍的MATLAB基础知识包括:数据类型;基本矩阵操作;运算符;字符串处理函数。目录 2.1 数 据 类 型 1数值类型 2逻辑类型 3字符和字符串类型 4结构体类型 MATLAB中
12、定义了很多种数据类型。本节讨论MATLAB中主要的数据类型及其使用方法。在MATLAB中有15种基本数据类型,每种基本数据类型均以数组/矩阵的形式出现。数值类型包含 整数;浮点数;复数;Inf;NaN MATLAB支持1、2、4和8字节的有符号整数和无符号整数。数据类型的名称、表示范围和转换函数如下表所示。MATLAB有单精度和双精度两种浮点数。其名称、存储空间、表示范围和转换函数如下表所示。复数包含实部和虚部。在MATLAB中可以用i或者j来表示虚部。Inf和-Inf分别表示正无穷大和负无穷 大。除法运算中除数为0或者运算结果溢 出都会导致inf或-inf的运行结果。在MATLAB中用NaN
13、(Not a Number)来表示一个既不是实数也不是复数的数 值。在MATLAB中逻辑类型包含true和false,分别由1和0表示。在MATLAB中用函数logical()将任何非零的数值转换为true(即1),将数值0转换为false(即0)。在MATLAB中,数据类型(char)表示一个字符。一个char类型的1n数组称为字符串string。结构体类型是一种由若干属性(field)组成的MATLAB数组,其中的每个属性可以是任意数据类型。下图表示了一个结构体(Personel),它包括3个属性(Name、Score和Salary),其中Name是一个字符串,Score是一个数值,Sal
14、ary是一个15的向量。构造一个结构体(数组)有两种方法。利用赋值语句 通过赋值语句为结构体中的每个指定属性赋值,从而构造结构体。利用函数struct()在MATLAB中,函数struct()的具体用法如下:其中,strArray、field和val分别表示结构体名、属性名和相应的属性值。通过结构体数组的下标引用,可以访问任意元素的所有属性,同时可以对属性进行赋值。2.2 基本矩阵操作 2.2.1 矩阵的构造 2.2.2 矩阵大小的改变 2.2.3 矩阵下标引用 2.2.4 矩阵信息的获取 2.2.5 矩阵结构的改变 2.2.6 稀疏矩阵 在MATLAB中,所有的数据均以二维、三维或高维矩阵的
15、形式存储,每个矩阵的单元可以是数值类型、逻辑类型、字符类型或者其他任何数据类型。对于标量,可以用11矩阵来表示;对于一组n个数据,可以用1n矩阵来表 示;对于多维数组,可以用多维矩阵来表示。在MATLAB中,用命令whos来显示数据的类型、存储空间等信息。1简单矩阵构造 2特殊矩阵构造 3向量构造 最简单的方法是采用矩阵构造符“”。构造1n矩阵(行向量)时,可以将各元素依次放入矩阵构造符内,并且以空格或者逗号分隔;构造mn矩阵时,每行如上处理,并且行与行之间用分号分隔。在MATLAB中还提供一些函数用来构造特殊矩阵,这些函数如下表所示。续表 最 简 单 的 方 法 是 采 用 向 量 构 造
16、符“:”,其常用的用法如下。(1)a:b返回以a为起点,以1为步长,且所有取值在a与b之间的向量。(2)a:s:b 返回以a为起点,以s为步长,且所有取值在a与b之间的向量。1矩阵的合并 2矩阵行列的删除 矩阵的合并就是把两个或者两个以上的矩阵连接成一个新矩阵。矩阵构造符 可用于构造矩阵,并可以作为一个矩阵合并操作符。表达式C=A B在水平方向合并矩阵A和B;表达式C=A;B在竖直方向合并矩阵A和B。具有相同行数的两个矩阵,合并为一个新矩阵 不具有相同行数的两个矩阵,不允许合并为一个新矩阵 除了矩阵合并符“”外,还可以使用矩阵合并函数。矩阵合并函数的描述和基本调用格式如下表所示。要删除矩阵的某
17、一行或者是某一列,只需将该行或者该列赋予一个空矩阵即可。1访问单个元素2线性引用元素3访问多个元素 本小节将介绍通过矩阵下标来存取元素值的方法,包括访问单个元素、线性引用元素和访问多个元素等。对于矩阵A,线性引用元素的格式为 A(k)。通常这样的引用用于行向量或列 向量,但也可用于二维矩阵。MATLAB按列优先排列的一个长列向量格 式(线性引用元素)来存储矩阵元素。操作符“:”可以用来表示矩阵的多个元素。若A是二维矩阵,其主要用法如下:A(:,:)返回矩阵A的所有元素。A(i,:)返回矩阵A第i行的所有元素。A(i,k1:k2)返回矩阵A第i行的自k1到k2 列的所有元素。A(:,j)返回矩阵
18、A第j列的所有元素。A(k1:k2,j)返回矩阵A第j列的自k1到k2 行的 所有元素。若A是多维矩阵,也可以通过类似的方法实现对其访问。1矩阵尺寸信息 2元素的数据类型 3矩阵的数据结构 本小节介绍如何获取矩阵的信息,包括矩阵的尺寸、元素的数据类型和矩阵的数据结构等。矩阵尺寸函数可以得到矩阵的形状和大小信息,这些函数如下表所示。查询元素数据类型信息的部分函数如下表所示。判断矩阵是否为某种指定数据结构的函数如下表所示。改变矩阵结构的函数表 1稀疏矩阵的创建 2查看稀疏矩阵 3稀疏矩阵的运算规则 在MATLAB中,可以用满矩阵存储方 式和稀疏矩阵存储方式来存储矩阵。若一个矩阵只有少数的元素非零,
19、称为稀 疏矩阵。稀疏矩阵非零元素及其对应的下 标来表示。用户可以创建双精度、复数和逻辑等类型 的稀疏矩阵。在MATLAB中,用函数sparse()实现满矩阵到稀疏矩阵的转换。在MATLAB中用函数full()实现稀疏矩阵 到满矩阵的转换。在MATLAB中,还可以用函数sparse()直接创建稀疏矩阵,其具体用法如下。S=sparse(i,j,s,m,n),其中,i和j分别是稀疏矩阵非零元素的行和列下标,s为相应的非零元素的值,m和n分别是矩阵的行数和列数。MATLAB还提供一些函数用于创建特殊稀疏矩阵,这些函数如下表所示。MATLAB提供一些函数用于查看稀疏矩阵的信息,如下表所示。下面的例子都
20、是基于MATLAB自带的稀疏矩阵west0479。在MATLAB中的各种命令和函数都可以用于稀疏矩阵的运算,并且遵循如下的一些约定。(1)把矩阵变为标量或者定长向量的函数总 是给出满矩阵;(2)对于标量或者定长向量变换到矩阵的函 数,如函数zeros()、ones()、eye()、rand()等总是给出满矩阵;(3)从矩阵到矩阵的变换函数将以原矩阵的 形式出现;(4)在参与矩阵扩展的子矩阵(如 A B;C D)中,只要有一个是稀疏矩阵,那 么所得的结果也是稀疏矩阵;(5)在矩阵引用中,将仍以原矩阵形式给出 结果。2.3 运算符和特殊符号2.3.1 算数运算符2.3.2 关系运算符2.3.3 逻
21、辑运算符2.3.4 运算优先级 在MATLAB中提供了丰富的运算符,包括算数、关系和逻辑等3种运算符。在MATLAB中,算数运算符的用法和功能如下表所示。续表补充说明AB的用法如下:当A和B都为矩阵时,此运算无定义;当A和B都是标量时,表示标量A的B次 幂;当A是标量且B为矩阵时,表示标量A的 B中各元素次幂;当A为方阵且B为正整数时,表示矩阵A 的B次乘积;当A为方阵且B为负整数时,表示矩阵A逆 的负B次乘积;当A为可对角化的方阵且B为非整数时,有 如下表达式:MATLAB中关系运算符的用法和功能如下表所示。值得注意的是,关系运算符只针对两个相同长度的矩阵,或其中之一是标量的情况进行运算。对
22、于前者,是指两个矩阵的对应元素进行比 较,返回具有相同长度的矩阵;对于后者,是指这个标量与另一个矩阵的每元 素进行运算。关系运算C=f(A,B)的运算结果只有0和1两种情况,其中,函数f()表示关系运算符,0表示不满足条件,1表示满足条件。MATLAB提供元素方式和比特方式等逻辑运算符。元素方式逻辑运算符的用法和功能如下表所示,其中例子采用如下矩阵:元素方式逻辑运算符&、|和 与函数and()、or()和not()是等价的。比特方式逻辑运算符只接受逻辑和非负整数类型的输入变量,它是针对输入变量的二进制进行逻辑运算。比特方式逻辑运算符的用法和功能如下表所示,表中例子采用A=28和B=200,其对
23、应的二进制分别为11100和11001000。运算符的优先级决定表达式求值顺序;具有相同优先级的运算符从左到右依次进 行运算;不同优先级的运算符采用先进行优先高的 运算。由表中可以看到,括号的优先级别最高,因此可以用括号来改变默认的优先等级。2.4 字符串处理函数 2.4.1 字符串的构造 2.4.2 字符串的比较 2.4.3 字符串的查找和替换 2.4.4 字符串与数值间的转换 MATLAB提供了丰富的字符串操作,包括字符串的创建、合并、比较、查找以及与数值之间的转换。在MATLAB中,可以用一对单引号来表示字符串。1字符串比较函数 2用关系运算符比较字符串 在MATLAB中提供了对字符串、
24、字符串数组和字符子串的比较功能。在MATLAB中,字符串比较函数如下表所示。在MATLAB中,可以对字符串运用关系运算符,但要求两个字符串具有相同的长度,或者其中一个是标量。MATLAB提供的一些字符串查找和替换函数如下表所示。MATLAB提供的一些数值转换为字符串函数如下表所示。续表 MATLAB提供的一些字符串转换为数值函数如下表所示。习 题1计算复数3+4i与56i的乘积。2构建结构体Students,属性包含Name、Age 和Email,数据包括Zhang,18,、Wang,21,和Li,,构建后读取 所有Name属性值,并且修改Zhang的Age属 性值为19。3用满矩阵和稀疏矩阵
25、存储方式分别构造下述矩 阵:4采用向量构造符得到向量1,5,9,41。5按水平和竖直方向分别合并下述两个矩阵:6 分别删除第5题两个结果的第2行。7 分别将第5题两个结果的第2行最后3列的数值 改为11 12 13。8 分别查看第5题两个结果的各方向长度。9 分别判断pi是否为字符串和浮点数。10分别将第5题两个结果均转换为29的矩阵。11计算第5题矩阵A的转秩。12分别计算第5题矩阵A和B的A+B、A.*B和 AB。13判断第5题矩阵A和B中哪些元素值不小于4。14分别用函数strcat()和矩阵合并符合并如下字 符串:The picture is 和 very good。15创建字符串数组
26、,其中元素分别为Picture 和 Pitch 。16在第14题结果中查找字符串e。17在第15题结果中匹配字符串Pi。18将字符串very good转换为等值的整数。19将十进制的50转换为二进制的字符串。20将十六进制的字符串50转换为三进制的整 数。第三章第三章 数学运算数学运算 本章将着重介绍MATLAB中与数学运算有关的函数和概念。在MATLAB中一切数据均能以矩阵的形式表示:针对矩阵整体的数学运算,称之为矩阵运 算;针对矩阵元素的数学运算,称之为矩阵元 素运算。目录 3.1 矩 阵 运 算 3.2 矩阵元素运算 习 题3.1 矩 阵 运 算 3.1.1 矩阵分析 3.1.2 线性方
27、程组 3.1.3 矩阵分解 3.1.4 矩阵的特征值和特征向量 3.1.5 矩阵相似变换 3.1.6 非线性运算 矩阵运算是线性代数中极其重要的部分,MATLAB具有强大的矩阵运算能力。1向量间的距离 2矩阵的秩 3矩阵的行列式 4矩阵的迹 5矩阵的化零矩阵 6矩阵的正交空间 7矩阵的简化梯形形式 8矩阵空间之间的角度 MATLAB提供的部分矩阵分析函数如下表所示。矩阵A中线性无关的列向量个数称为列秩,线性无关的行向量个数称为行秩。可以证明列秩与行秩是相等的。矩阵的迹定义为矩阵对角元素之和。在MATLAB中用函数trace()来计算矩阵的迹。对于非满秩矩阵A,若存在矩阵Z使得AZ=0且ZTZ=
28、I,则称矩阵Z为矩阵A的化零矩阵。在MATLAB中用函数null()来计算矩阵的化零矩阵。矩阵A的正交空间Q满足QTQ=I,且矩阵Q与A具有相同的列基底。矩阵A的简化梯形形式为 ,其中Ir为r阶单位矩阵。*0*rI 矩阵空间之间的角度代表具有相同行数的两个矩阵线性相关程度,夹角越小代表线性相关度越高。线性方程组求解问题,可以表述为给定两个矩阵A和B,求解X使得AX=B或XA=B。XA=B可以表示为AY=B,且X=Y。下面仅讨论AX=B的情况。1Cholesky分解 2LU分解 3QR分解 4奇异值分解 5Schur分解 矩阵分解是把一个矩阵分解成比较简单或者对它性质比较熟悉的若干矩阵的乘积的形
29、式。本小节将介绍几种矩阵分解的方法。矩阵分解函数表 Cholesky分解是把对称正定矩阵A表示为上三角矩阵R的转置与其本身的乘积,即A=RTR。对于稀疏矩阵,MATLAB中用函数cholinc()计算不完全Cholesky分解,具体用法如下:R=full(cholinc(sparse (X),DROPTOL),其中DROPTOL为不 完全Cholesky分解的丢失容限;R=full(cholinc(sparse(X),0),完 全Cholesky分解。高斯消去法又称LU分解,将任意一个方阵A分解为一个交换下三角 矩阵L和一个上三角矩阵U的乘积,即 A=LU。交换下三角矩阵为下三角矩阵经行变换的
30、 结果。LU分解在MATLAB中用函数lu()来实现,具体用法如下:L,U=lu(X),X为一个方阵,L为交换 下三角矩阵,U为上三角矩阵,满足关系 X=L*U;L,U,P=lu(X),X为一个方阵,L为下 三角矩阵,U为上三角矩阵,P为置换矩 阵,满足关系P*X=L*U或X=P-1 *L*U。考虑线性方程组AX=B和矩阵A的LU分解,线性方程组可改写成L*U*X=B,由于左除算符可以快速处理三角矩阵,因此:X=U(LB)矩阵的行列式和逆也可以利用LU分解来计算,如det(A)=det(L)*det(U)inv(A)=inv(U)*inv(L)对于稀疏矩阵,在MATLAB中提供了函数luinc
31、()来做不完全LU分解,其具体用法如下:L U=luinc(X,DROPTOL),其中X、L和U 的含义与函数lu()中的变量相同,DROPTOL为 不完全LU分解的丢失容限。当DROPTOL设为0 时,退化为完全LU分解。L,U=luinc(X,0),0级不完全LU分解。L,U,P=luinc(X,0),0级不完全LU分解。QR分解就是将mn的矩阵A分解为mn的矩阵Q和nn的上三角矩阵R的乘积,且Q*Q=I,即A=Q*R。在MATLAB中QR分解是由函数qr()来实现,其具体用法如下:Q,R=qr(A)满足A=Q*R。R=qr(A),返回上三角矩阵R。奇异值分解就是将mn的矩阵A分解为A=U
32、*S*V,其中U为mm的酉矩阵,V为nn的酉矩阵,S为mn的矩阵,并可如下表示:,其中 ,000S12,rdiag()rrank A0(1,2,)iir 在MATLAB中奇异值分解是由函数svd()来实现,其具体用法如下:Schur分解就是将复方阵A分解为A=U*L*U,其中U为酉矩阵,L为上(下)三角矩阵,其对角线元素为A的特征值。在MATLAB中Schur分解是由函数schur()来实现,其具体用法如下:U,L=schur(A),满足A=U*L*U,其中L为上三角矩阵。L=schur(A),返回上三角矩阵L。方阵A的特征值和其对应的特征向量满足下式:A*=*在MATLAB中用函数eig()
33、来计算特征值和其对应的特征向量,其具体用法如下:d=eig(A),返回矩阵A的所有特征值。V,D=eig(A),返回矩阵A的特征值和 特征向量。1对角阵变换 2Jordan变换 矩阵相似变换是指,对于方阵A和非奇异矩阵B可得到相似矩阵X=B-1*A*B。对于方阵A,若V D=eig(A)得到的矩阵V非奇异,则A可经过相似变换得到对角阵,即D=V-1*A*V,也称矩阵A可对角化。对于方阵A,若V D=eig(A)得到的矩阵V奇异,则A经过相似变换将不能得到对角阵,只能得到其对应的Jordan标准型。在MATLAB中用函数jordan()来实现Jordan变换,其具体用法如下:V,D=jordan
34、(A),满足 D=V-1*A*V。D=jordan(A),返回矩阵A对应的 Jordan标准型。1矩阵指数运算 2矩阵对数运算 3矩阵开平方运算 4通用矩阵运算 MATLAB提供一些矩阵的非线性运算函数,其功能如下表所示。矩阵对数运算是矩阵指数运算的逆运算,在MATLAB中用函数logm()来计算矩阵对数,其具体用法如下:L=logm(A),返回矩阵A的对数。对于方阵A,可以计算它的开平方得到矩阵X,即满足X*X=A。如果矩阵A是奇异的,它有可能不存在平方根X。在MATLAB中,有两种计算矩阵A平方根的方法,即A0.5和sqrtm(A)。函数sqrtm()比A0.5的运算精度更高,其具体用法如
35、下:X=sqrtm(A),返回矩阵A的平方根X。MATLAB提供通用矩阵运算的函数funm(),其具体用法如下:F=funm(A,fun),将指定函数fun作用 在方阵A上。可以使用的指定函数fun如下表所示。3.2 矩阵元素运算 3.2.1 三角函数 3.2.2 指数和对数函数 3.2.3 复数函数 3.2.4 截断和求余函数 3.2.5 特殊函数 本节将介绍矩阵元素的数学函数:包括三角函数;指数/对数函数;复数函数;截断/求余函数;特殊函数。MATLAB提供一些三角函数,其功能如下表所示。续表 MATLAB提供一些指数和对数函数。指数和对数函数表 MATLAB提供一些复数函数,其功能如下表
36、所示。MATLAB提供一些截断和求余函数,其功能如下表所示。1工程函数 2数论函数 3坐标变换函数 本小节介绍一些用途比较特殊的数学函数,包括应用于工程上的函数(工程函数)、数论函数和坐标变换函数。工程函数经常在数学、物理和工程等问题中出现,在MATLAB中提供一些工程函数。表工程函数 MATLAB提供一些数论函数,其功能如下表所示。在MATLAB中提供一些坐标变换函数,其功能如下表所示。习 题1计算矩阵A的范数、行列式、秩、化零空间和正 交空间。2求解线性方程组AX=B,其中A如第1题所示,B=1 1 1 1 1T。3对矩阵A进行LU分解和Schur分解,其中A如第 1题。17 24 1 8
37、 50 23 5 7 14 49 4 6 13 20 43 10 12 19 21 62 11 18 25 2 56A4对矩阵A的前4行进行QR分解和奇异值分解,其 中A如第1题。5计算矩阵A的特征根及对应的特征向量,判断矩 阵A是否可对角化,其中A如第1题。6计算矩阵A的指数、开平方和余弦值,其中A如 第1题。7计算矩阵A每个元素的指数、开平方和余弦值 (元素单位为度),其中A如第1题。8 计算复数矩阵C每个元素的模、相角和共轭。9 分别使用函数fix()、floor()、ceil()和 round(),计算第8题中的相角结果。10将2-i的模结果近似为有理数,并以数值形式 显示。11计算,
38、其中m=4!和n是42与35的最大公因 数。12将球坐标系中的点(1,1,1)分别转换到笛卡尔 坐标系和极坐标系。3+4i 2-i -i 2 -2 0C第四章第四章 基本编程基本编程 MATLAB提供丰富的函数库,并可以进行程序设计,编写扩展名为.m的M文件,实现各种程序设计功能。MATLAB提供大量的函数,包括内建函数和自带函数。用户也可以利用M文件来创建函数、函数库和脚本。目录 当需要创建一个M文件时,可以通过【File】|【New】|【M-File】菜单进入文本编辑器界面。当需要打开一个M文件时,可以通过【File】|【Open】菜单进入文件选择的目录界面。本章将着重介绍M文件编程的规则
39、和注意事项。4.1 M文件基础 4.1.1 函数 4.1.2 脚本 4.1.3 子函数与私有函数 4.1.4 伪代码 M文件有函数(Functions)和脚本(Scripts)两种格式。M文件可以在任意的文本编辑器中进行 创 建 和 编 辑,但 推 荐 和 默 认 使 用MATLAB自带的文本编辑器和Notebook。MATLAB中许多常用的函数都是函数式M文件。函数被调用时,通过获取外部参数进行运算,并向外部返回运算结果。函数内部的变量都是隐含的,存放在函数本身的工作空间内。脚本也是扩展名为.m的文件,可包含MATLAB的各种命令。在命令窗口中直接键入此文件的文件名,MATLAB可逐一执行文
40、件内的所有命令。将脚本所在目录设置为当前工作目录,并在文本编辑器中打开脚本后,可通过如下方法运行脚本:单击【Debug】菜单下的【Run】;单击快捷键F5;单击工具条上的图标 。脚本运行过程所产生的变量都是全局变量,都驻留在MATLAB工作空间内。结合上例,对函数和脚本进行总结如下。(1)函数名必须与文件名相同。(2)脚本没有输入参数和返回值。(3)函数可以包括零个或多个输入参数和返 回值,如函数nargin和nargout包含输 入参数和返回值的个数。(4)函数被调用时,MATLAB会为它开辟一 函数工作空间(Function workspace),用来存放中间变量,当执行完函数文件的最后一
41、条命令或 者遇到return命令时就结束该函数的 运行,返回函数的输出,同时将临时 工作空间清空。(5)在M文件中,从开头到第一个非注释行 之间的第一个注释行组(之间无空行)是帮助文本。(6)函数中的变量(除特殊声明外)都是局部变 量,而脚本中的变量都是全局变量,关于变量 的有关内容详见第4.2节。(7)在函数中调用脚本文件,等价于在函数中将 脚本文件的内容粘贴在调用的位置。一个M文件可以包含多个函数,其中,一个是主函数,其他是子函数。主函数必须出现在最上方,其后可有若干子函数,子函数的次序可随意调整。主函数和各子函数的工作空间都是彼此独立的。私有函数是主函数的一种,它只能够在一个特定的限定函
42、数群中可见。当M文件中需要调用某一个函数时,MATLAB是按照以下顺序来搜寻的:检查此函数是否是子函数;检查此函数是否为私有目录的函数;从所设定的搜寻路径搜索此函数。搜索过程中,只要找到与第一个文件名相符的函数就会立即取用而停止搜索。一个M文件首次被调用时,MATLAB将对该M文件进行语法分析,并把生成的相应伪代码(P码)存放在内存中。此后当再次调用该M文件时,将直接运行该文件在内存中的伪代码。伪代码文件和原码文件具有相同的文件名,但其扩展名为.p。伪代码文件不是只有当M文件被调用时才产生,也可使用pcode命令预先生成。4.2 变量和语句 4.2.1 变量类型 4.2.2 程序控制结构 MA
43、TLAB的主要功能虽然是数值运算,但是它也是一个完整的程序语言,包括各种语句格式和语法规则。但与C语言不同的是,MATLAB中的变量是不需要事先定义的,如前3章的例子。在MATLAB中,变量名必须以字母开头,之后可以是任意字母、数字或下划线,但之间不能有空格;区分大小写;第63个字符之后的部分将被忽略。除了上述命名规则外,MATLAB还提供一些特殊的变量,如下表所示。除命名规则外,变量命名时还需要注意以下两个方面:(1)变量名不能与已有函数名相同,否则在 内存中有该变量时将不能调用同名函 数;(2)变量名不能与MATLAB预留的关键字和 特殊变量名相同,否则系统会显示错误 信息。变量按照作用范
44、围分为局部变量和全局变量。在MATLAB中使用命令global声明全局变量,如下所示:需要使用指定全局变量的M文件,都必须在各自的代码中声明此全局变量。只要存在声明某全局变量的M文件,则全局变量存在。在使用全局变量中需要注意以下几个方面。在使用之前必须首先定义,建议将定义放在函 数体的首行位置。虽然对全局变量的名称并没有特别的限制,但 是为了提高程序的可读性,建议采用大写字符 命名全局变量。全局变量会损坏函数的独立性,使程序的书写 和维护变得困难,尤其是在大型程序中,不利 于模块化,不推荐使用。1顺序结构 2循环结构 3选择结构 4其他控制结构 对于实现任何功能的程序,均可由顺序、循环和选择3
45、种基本结构组合实现。为了更方便地设计程序,还需要一些特殊的控制结构,如跳出循环结构。下面将分别进行介绍。顺序结构就是由前至后依次执行程序的各条代码,直至最后一条代码。脚本文件就是典型的顺序结构。01234567891000.10.20.30.40.50.60.7Step ResponseTime(sec)Amplitude 循环结构是按照给定的条件,重复执行指定的代码。该结构一般用于有规律的重复运算。在MATLAB中包括for循环和while循环。选择结构是根据给定的条件来执行不同的代码。在MATLAB中有if-else-end和switch-case-otherwise两种结构。在程序设计中
46、经常遇到提前终止循环、跳出子程序、显示出错信息等情况,主要有continue、break、return、echo、error、trycatch等。下面分别对各结构进行介绍。其作用是结束本次循环,即跳过循环体中尚未执行的代码,接着进行下一次是否执行循环的判断。它的作用是终止本次循环,跳出所在层循环。此命令可使正在运行的函数正常退出,并返回调用它的代码段继续运行,它也可强制结束该函 数 的 执 行。如 M A T L A B 自 带 函 数showopcevents()。执行M文件时,通常在命令窗口是看不到执行过程的,但在特殊情况下(比如演示)要求M文件的每条命令都要显示出来。这时可以用echo命
47、令实现这样的功能。对于脚本,echo命令可以用以下方式来实现:对于函数,echo命令可以用以下方式来实现:此命令用来提示并接收用户从键盘输入数据、字符串或表达式的值,语法规则如下:此命令将停止文件的执行并将控制权交给键盘,此时命令窗口的提示符由“”变成“K”,当输入return后,控制权将交回文件。该命令对程序的调试和在程序运行中修改变量值都很方便。此命令用于暂时中止程序的运行。该命令对程序的调试和查询中间变量值时很方便。该命令的语法规则如下:4.3 程 序 调 试4.3.1 直接调试法4.3.2 工具调试法 在编译和运行程序时出现错误(警告)无法避免,因此掌握程序调试的方法和技巧对提高工作效
48、率很重要。下面针对语法错误和逻辑错误推荐两种调试方法,即直接调试法和工具调试法。对于简单的程序往往采用直接调试法,通常采取的措施如下。(1)通过分析后,将重点怀疑语句后的分号 删掉,将结果显示出来,然后与预期 值进行比较。(2)单独调试函数时,将函数声明行注释掉,并 定义输入变量的值,然后以脚本方式运行,这 样可保存中间变量,进而进行分析和找出错 误。(3)在程序中的适当位置添加输出变量值的 代码。(4)在程序中的适当位置添加keyboard命 令。复杂的程序,必须采用工具调试法,即 借 助 M AT L A B 提 供 的 工 具 调 试 器(Debugger)进行调试。1以命令行为主的程序
49、调试 2以图形界面为主的程序调试 MATLAB提供了调试程序的工具,利用这些工具可以提高编程的效率,包括命令行的调试函数和图形界面的菜单命令。以命令行为主的程序调试手段具有通用性,可以适用于各种平台,它主要是应用MATLAB提供的调试命令。在命令窗口输入help debug可以看到对于这些命令的简单描述,下面分别进行介绍。这是最重要的部分,可以利用它来指定程序代码的断点,使得程序在断点前停止执行,并进入调试模式,从而可以检查当前各个变量的值。dbstop in mfile 在文件名为mfile的M文件的第一个可执行语句前设置断点。dbstop in mfile at lineno 在文件名为mfile的M文件的第lineno行设置断点。如果第lineno行为非执行语句,则在其后的第一个可执行语句前设置断点。dbstop in mfile at subfun 在文件名为mfile的M文件的子程序subfun的第一个可执行语句前设置断点。dbstop if error 在程序运行遇到错误时,自动设置断点。这里的错误不包括trycatch之间的错误。dbstop if all error 在程序运行遇到错误时,自动设置断点。这里的错误包括trycatch之间的错误。dbstop if warning 在程序运行遇到警告时,自动设置断点。dbstop if caught error 在程