《MATLAB实用教程学习.pptx》由会员分享,可在线阅读,更多相关《MATLAB实用教程学习.pptx(1146页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、科学计算与仿真推荐参考书推荐参考书第1页/共1146页科学计算与仿真推荐参考书推荐参考书第2页/共1146页matlabmatlab是什么一个可视化的计算程序,广泛使用于从个人计算机到超级计算机范围内的各种计算机上包括命令控制、可编程,上百个预先定义命令和函数有许多强有力的命令,能完成大量的高级矩阵处理强有力的二维、三维图形工具能与其他程序一起使用25个(不断增加中)不同的工具箱应用于特殊的应用领域工业研究与开发的有力工具数学教学,尤其线代,数值分析,科学计算方面的教研工具电子学,控制理论,物理学等工程科学方面的教研工具经济学,化学和生物学等有计算问题的所有领域中的教学与研究名字取自矩阵实验室
2、(matrix laboratory)第3页/共1146页matlabmatlab不是什么不是万能的解决工具不是最高性能的编程语言受计算条件限制,不能解决超大型实际问题不能解决工具箱之外的问题种类-需要编写接口、算法甚至工具箱第4页/共1146页MATLAB 的名称源自 Matrix Laboratory,它是一种科学计算软件,专门以矩阵的形式处理数据。第5页/共1146页MATLAB 将高性能的数值计算和可视化集成在一起,并提供了大量的内置函数,从而被广泛地应用于科学计算丶控制系统丶信息处理等领域的分析丶仿真和设计工作,而且利用 MATLAB 产品的开放式结构,可以非常容易地对 MATLAB
3、 的功能进行扩充,从而在不断深化对问题认识的同时,不断完善 MATLAB 产品以提高产品自身的竞争能力。第6页/共1146页MATLAB Compiler 是一种编译工具,它能够将那些利用 MATLAB 提供的编程语言 M 语言编写的函数文件编译生成为函数库丶可执行文件 COM 组件等等。这样就可以扩展 MATLAB 功能,使 MATLAB 能够同其他高级编程语言例如 C/C+语言进行混合应用,取长补短,以提高程序的运行效率,丰富程序开发的手段。第7页/共1146页目前 MATLAB 产品族可以用来进行:数值分析 数值和符号计算 工程与科学绘图 控制系统的设计与方针 数字图像处理 数字信号处理
4、 通讯系统设计与仿真 财务与金融工程第8页/共1146页MATLAB 是 MATLAB 产品家族的基础,它提供了基本的数学算法,例如矩阵运算丶数值分析算法,MATLAB 集成了 2D 和 3D 图形功能,以完成相应数值可视化的工作,并且提供了一种交互式的高级编程语言 M 语言,利用 M 语言可以通过编写脚本或者函数文件实现用户自己的算法。第9页/共1146页利用 M 语言还开发了相应的 MATLAB 专业工具箱函数供用户直接使用。这些工具箱应用的算法是开放的可扩展的,用户不仅可以查看其中的算法,还可以针对一些算法进行修改,甚至允许开发自己的算法扩充工具箱的功能。目前 MATLAB 产品的工具箱
5、有四十多个,分别涵盖了数据获取丶科学计算丶控制系统设计与分析丶数字信号处理丶数字图像处理丶金融财务分析以及生物遗传工程等专业领域。第10页/共1146页Simulink 是基于 MATLAB 的框图设计环境,可以用来对各种动态系统进行建模丶分析和仿真,它的建模范围广泛,可以针对任何能够用数学来描述的系统进行建模,例如航空航天动力学系统丶卫星控制制导系统丶通讯系统丶船舶及汽车等等,其中了包括连续丶离散,条件执行,事件驱动,单速率丶多速率和混杂系统等等。Simulink 提供了利用鼠标拖放的方法建立系统框图模型的图形界面,而且 Simulink 还提供了丰富的功能块以及不同的专业模块集合,利用 S
6、imulink 几乎可以做到不书写一行代码完成整个动态系统的建模工作。第11页/共1146页Stateflow 是一个交互式的设计工具,它基于有限状态机的理论,可以用来对复杂的事件驱动系统进行建模和仿真。Stateflow 与 Simulink 和 MATLAB 紧密集成,可以将 Stateflow 创建的复杂控制逻辑有效地结合到 Simulink 的模型中。第12页/共1146页在 MATLAB 产品族中,自动化的代码生成工具主要有 Real-Time Workshop(RTW)和 Stateflow Coder,这两种代码生成工具可以直接将 Simulink 的模型框图和 Stateflo
7、w 的状态图转换成高效优化的程序代码。利用 RTW 生成的代码简洁丶可靠丶易读。目前 RTW 支持生成标准的 C 语言代码,并且具备了生成其他语言代码的能力。整个代码的生成丶编译以及相应的目标下载过程都是自动完成的,用户需要做得仅仅使用鼠标点击几个按钮即可。MathWorks 公司针对不同的实时或非实时操作系统平台,开发了相应的目标选项,配合不同的软硬件系统,可以完成快速控制原型(Rapid Control Prototype)开发丶硬件在回路的实时仿真(Hardware-in-Loop)丶产品代码生成等工作。第13页/共1146页另外,MATLAB 开放性的可扩充体系允许用户开发自定义的系统
8、目标,利用Real-Time Workshop Embedded Coder 能够直接将 Simulink 的模型转变成效率优化的产品级代码。代码不仅可以是浮点的,还可以是定点的。MATLAB 开放的产品体系使 MATLAB 成为了诸多领域的开发首选软件,并且,MATLAB 还具有 300 余家第三方合作伙伴,分布在科学计算丶机械动力丶化工丶计算机通讯丶汽车丶金融等领域。接口方式包括了联合建模丶数据共享丶开发流程衔接等等。第14页/共1146页MATLAB 结合第三方软硬件产品组成了在不同领域内的完整解决方案,实现了从算法开发到实时仿真再到代码生成与最终产品实现的完整过程。主要的典型应用包括:
9、控制系统的应用与开发快速控制原型与硬件在回路仿真的统一平台 dSPACE 信号处理系统的设计与开发全系统仿真与快速原型验证,TI DSP 丶 Lyrtech 等信号处理产品软硬件平台 通信系统设计与开发结合 RadioLab 3G 和 Candence 等产品 机电一体化设计与开发全系统的联合仿真,结合 Easy 5 丶 Adams 等第15页/共1146页Matlab课程的组织结构绪论基础知识数学运算基本编程数据显示 及存取符号计算数值计算图形用户界面GUISimulink仿真外部接口第16页/共1146页第一章第一章 概述概述 第17页/共1146页 MATLAB将计算、可视化和编程等功能
10、集于一个易于使用的环境,并遇有如下特点:功能强大;简单易学;编程效率高。第18页/共1146页1.1 MATLAB简介及安简介及安装装1MATLAB的发展历程 2MATLAB的特点 3MATLAB的安装第19页/共1146页1MATLAB的发展历程 1980年,美国新墨西哥大学计算机科学 系主任Cleve 着手编写供学生使用的子程 序接口程序,取名为MATLAB;1984年,推出了MATLAB第一个商业版 本;1992年,推出MATLAB 4.0版;第20页/共1146页 1997年,推出MATLAB 5.0版;2000年,推出MATLAB 6.0版;2004年,推出MATLAB 7.0版;2
11、008年,推出MATLAB 7.6版;2012年,推出MATLAB R2012A版第21页/共1146页2MATLAB的特点 开发环境 编程 数值处理 图形化 图形用户界面 文件I/O和外部应用程序接口第22页/共1146页3MATLAB的安装MATLAB对PC系统的要求 第23页/共1146页1.2 MATLAB的目录结的目录结构构 第24页/共1146页 安装MATLAB后,在安装目录下将包含如下文件夹。第25页/共1146页 续表第26页/共1146页1.3 MATLAB的工作环的工作环境境1菜单和工具栏 2命令窗口 3历史命令窗口 4当前工作目录窗口 5工作空间窗口第27页/共1146
12、页 本节介绍MATLAB的工作界面和基本的操作方法。MATLAB的工作界面主要包括:菜单;工具栏;命令窗口;历史命令窗口;当前工作目录窗口;工作空间窗口。第28页/共1146页MATLAB工作环境 第29页/共1146页1菜单和工具栏 这里只简单介绍默认情况下的菜单和工具栏。【File】菜单主要用于对文件的处理。【Edit】菜单主要用于复制、粘贴等操 作,与一般Windows程序的类似,在此 不作详细介绍。第30页/共1146页【Debug】菜单用于调试程序。【Desktop】菜单用于设置主窗口中需要 打开的窗口。【Window】菜单列出当前所有打开的窗 口。【Help】菜单用于选择打开不同的
13、帮助系 统。第31页/共1146页 当用户单击“Current Directory”窗口时,使得该窗口成为活动窗口,同时增加一个如下图所示的菜单【View】,用于设置如何显示当前目录下的文件。第32页/共1146页 当用户单击“Workspace”窗口时,使得该窗口成为活动窗口,同时增加如下图所示的菜单【View】和【Graphics】。第33页/共1146页2命令窗口 命令窗口是MATLAB的主要交互窗口,用于输入命令并显示除图形以外的所有执行结果。MATLAB命令窗口中的“”为运算提示 符,表MATLAB处于准备状态。第34页/共1146页 当在提示符后输入一段程序或一段运算式 后【Ent
14、er】键,MATLAB会给出计算结 果,并再次进入准备状态(所得结果将被 保存在工作空间窗口中)。单击命令窗口右上角的“”按钮,可以使 命令窗口脱离主窗口而成为一个独立的窗 口。第35页/共1146页3历史命令窗口 主要用于记录所有执行过的命令;保留自安装后所有使用过命令的历史记 录,并标明使用时间;可以通过用鼠标双击某一历史命令来重新 执行该命令;可以成为一个独立的窗口。第36页/共1146页 通过上下文菜单,可删除或粘贴历史记录;可为选中的表达式或命令创建一个M文件;可为表达式或命令创建快捷按钮。第37页/共1146页4当前工作目录窗口 当前工作目录是指MATLAB运行文件时的目录。只有在
15、当前工作目录或搜索路径下的文件、函数可以被运行或调用。如下图所示。第38页/共1146页5工作空间窗口 工作空间窗口将显示目前内存中所有的MATLAB变量的变量名、数据结构、字节数以及类型等信息,如下图所示。第39页/共1146页1.4 MATLAB的通用命的通用命令令 1常用命令 2输入内容的编辑 3标点 4搜索路径与扩展第40页/共1146页 通用命令是MATLAB中经常使用的一组命令,这些命令可以用来管理目录、函数、变量、工作空间、文件和窗口等。下面对这些命令进行介绍。第41页/共1146页1常用命令 常用命令的功能第42页/共1146页2输入内容的编辑 在命令窗口中,MATLAB提供了
16、控制光标位置和进行简单编辑的键盘按键,部分常用的键盘按键及其功能如下表所示。第43页/共1146页3标点 在MATLAB中,一些标点符号被赋予了特殊的功能,如下表所示。第44页/共1146页4搜索路径与扩展 用户自己书写的函数有可能并没有保存在搜索路径下。要解决这个问题,只需把程序所在的目录扩展成MATLAB的搜索路径即可。第45页/共1146页(1)查看MATLAB的搜索路径 选择MATLAB主窗口中的【File】|【Set Path】菜单,进入【设置搜索路径】对话框。通过该对话框可为MATLAB添加或删除搜索路径。第46页/共1146页 在命令窗口中输入path或genpath可得到MAT
17、LAB的所有搜索路径,具体代码如下:第47页/共1146页(2)设置MATLAB的搜索路径 方法一:在MATLAB命令窗口中输入 editpath或pathtool命令或通过【File】|【Set Path】菜单,进入“设置搜索路 径”对话框,通过该对话框编辑搜索路 径。第48页/共1146页 方法二:在命令窗口执行“path(path,D:Study)”,然后通过“设置搜索路 径”对话框查看“D:Study”是否在搜索路 径中。方法三:在命令窗口执行“addpath D:Study-end”,将新的目录加到整个 搜索路径的末尾。如果将end改为 begin,可以将新的目录加到整个搜索路 径的
18、开始。第49页/共1146页1.5 MATLAB的帮助系的帮助系统统 1命令窗口查询帮助系统 2联机帮助系统 3联机演示系统第50页/共1146页 MATLAB为用户提供了非常完善的帮助系统。下面分别介绍MATLAB的三类帮助系统。第51页/共1146页1命令窗口查询帮助系统常见的帮助命令第52页/共1146页第53页/共1146页第54页/共1146页2联机帮助系统 用户可以通过下述方法进入MATLAB的联机帮助系统。直接单击MATLAB主窗口中的“”按钮;选中Help菜单的前4项中的任意一项;在命令窗口中执行helpwin、helpdesk 或doc。第55页/共1146页3联机演示系统
19、可以通过以下方式打开联机演示统。选择MATLAB主窗口菜单的【Help】|【Demos】选项;在命令窗口输入demos;直接在帮助页面上选择Demos页。第56页/共1146页第57页/共1146页第58页/共1146页第59页/共1146页1.6 MATLAB示例示例 下面以一个简单的例子展示如何使用MATLAB进行简单的数值计算。第60页/共1146页 第61页/共1146页第62页/共1146页习习 题题1简述MATLAB的主要功能。2在命令窗口输入“w=3+2”,然后依次使用clear和clc命 令,分别观察命令窗口、工作空间窗口和历史命令窗口 的变化。3将硬盘上一已有目录,加入到搜索
20、路径,并将其设置为 当前工作目录。4通过命令窗口,查询函数sin()的用法。5通过联机帮助系统,查询函数inv()的用法。第63页/共1146页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第64页/共1146页第二章第二章 基础知识基础知识第65页/共1146页 本章着重介绍的MATLAB基础知识包括:
21、数据类型;基本矩阵操作;运算符;字符串处理函数。第66页/共1146页目录目录 2.1 数 据 类 型 2.2 基本矩阵操作 2.3 运算符和特殊符号 2.4 字符串处理函数 习 题第67页/共1146页2.1 数数 据据 类类 型型 1数值类型 2逻辑类型 3字符和字符串类型 4结构体类型第68页/共1146页 MATLAB中定义了很多种数据类型。本节讨论MATLAB中主要的数据类型及其使用方法。在MATLAB中有15种基本数据类型,每种基本数据类型均以数组/矩阵的形式出现。第69页/共1146页1数值类型 数值类型包含 整数;浮点数;复数;Inf;NaN第70页/共1146页(1)整数类型
22、 MATLAB支持1、2、4和8字节的有符号整数和无符号整数。数据类型的名称、表示范围和转换函数如下表所示。第71页/共1146页(2)浮点数类型 MATLAB有单精度和双精度两种浮点数。其名称、存储空间、表示范围和转换函数如下表所示。第72页/共1146页(3)复数类型 复数包含实部和虚部。在MATLAB中可以用i或者j来表示虚部。第73页/共1146页(4)Inf和NaN Inf和-Inf分别表示正无穷大和负无穷 大。除法运算中除数为0或者运算结果溢 出都会导致inf或-inf的运行结果。在MATLAB中用NaN(Not a Number)来表示一个既不是实数也不是复数的数 值。第74页/
23、共1146页2逻辑类型 在 MATLAB中 逻 辑 类 型 包 含 true和 false,分 别 由 1和 0表 示。在MATLAB中用函数logical()将任何非零的数值转换为true(即1),将数值0转换为false(即0)。第75页/共1146页3字符和字符串类型 在MATLAB中,数据类型(char)表示一个字符。一个char类型的1n数组称为字符串string。第76页/共1146页4结构体类型 结构体类型是一种由若干属性(field)组成的MATLAB数组,其中的每个属性可以是任意数据类型。第77页/共1146页 下图表示了一个结构体(Personel),它包括3个属性(Nam
24、e、Score和Salary),其中Name是一个字符串,Score是一个数值,Salary是一个15的向量。第78页/共1146页(1)结构体数组的构造 构造一个结构体(数组)有两种方法。利用赋值语句 通过赋值语句为结构体中的每个指定属性赋值,从而构造结构体。第79页/共1146页 利用函数struct()在MATLAB中,函数struct()的具体用法如下:其中,strArray、field和val分别表示结构体名、属性名和相应的属性值。第80页/共1146页(2)结构体数组的访问 通过结构体数组的下标引用,可以访问任意元素的所有属性,同时可以对属性进行赋值。第81页/共1146页2.2
25、基本矩阵操作基本矩阵操作 2.2.1 矩阵的构造 2.2.2 矩阵大小的改变 2.2.3 矩阵下标引用 2.2.4 矩阵信息的获取 2.2.5 矩阵结构的改变 2.2.6 稀疏矩阵第82页/共1146页 在MATLAB中,所有的数据均以二维、三维或高维矩阵的形式存储,每个矩阵的单元可以是数值类型、逻辑类型、字符类型或者其他任何数据类型。对于标量,可以用11矩阵来表示;对于一组n个数据,可以用1n矩阵来表 示;对于多维数组,可以用多维矩阵来表示。在MATLAB中,用命令whos来显示数据的类型、存储空间等信息。第83页/共1146页2.2.1 矩阵的构造 1简单矩阵构造 2特殊矩阵构造 3向量构
26、造第84页/共1146页1简单矩阵构造 最简单的方法是采用矩阵构造符“”。构造1n矩阵(行向量)时,可以将各元素依次放入矩阵构造符内,并且以空格或者逗号分隔;构造mn矩阵时,每行如上处理,并且行与行之间用分号分隔。第85页/共1146页第86页/共1146页2特殊矩阵构造 在MATLAB中还提供一些函数用来构造特殊矩阵,这些函数如下表所示。第87页/共1146页续表第88页/共1146页第89页/共1146页3向量构造 最简单的方法是采用向量构造符“:”,其常用的用法如下。(1)a:b返回以a为起点,以1为步长,且所有取值在a与b之间的向量。第90页/共1146页(2)a:s:b 返回以a为起
27、点,以s为步长,且所有取值在a与b之间的向量。第91页/共1146页2.2.2 矩阵大小的改变 1矩阵的合并 2矩阵行列的删除第92页/共1146页1矩阵的合并 矩阵的合并就是把两个或者两个以上的矩阵连接成一个新矩阵。矩阵构造符 可用于构造矩阵,并可以作为一个矩阵合并操作符。表达式C=A B在水平方向合并矩阵A和B;表达式C=A;B在竖直方向合并矩阵A和B。第93页/共1146页 具有相同行数的两个矩阵,合并为一个新矩阵 不具有相同行数的两个矩阵,不允许合并为一个新矩阵 第94页/共1146页 除了矩阵合并符“”外,还可以使用矩阵合并函数。矩阵合并函数的描述和基本调用格式如下表所示。第95页/
28、共1146页第96页/共1146页2矩阵行列的删除 要删除矩阵的某一行或者是某一列,只需将该行或者该列赋予一个空矩阵即可。第97页/共1146页2.2.3 矩阵下标引用1访问单个元素2线性引用元素3访问多个元素第98页/共1146页 本小节将介绍通过矩阵下标来存取元素值的方法,包括访问单个元素、线性引用元素和访问多个元素等。第99页/共1146页1访问单个元素第100页/共1146页第101页/共1146页2线性引用元素 对于矩阵A,线性引用元素的格式为 A(k)。通常这样的引用用于行向量或列 向量,但也可用于二维矩阵。MATLAB按列优先排列的一个长列向量格 式(线性引用元素)来存储矩阵元素
29、。第102页/共1146页第103页/共1146页3访问多个元素 操作符“:”可以用来表示矩阵的多个元素。若A是二维矩阵,其主要用法如下:A(:,:)返回矩阵A的所有元素。A(i,:)返回矩阵A第i行的所有元素。第104页/共1146页 A(i,k1:k2)返回矩阵A第i行的自k1到k2 列的所有元素。A(:,j)返回矩阵A第j列的所有元素。A(k1:k2,j)返回矩阵A第j列的自k1到k2 行的 所有元素。若A是多维矩阵,也可以通过类似的方法实现对其访问。第105页/共1146页第106页/共1146页第107页/共1146页2.2.4 矩阵信息的获取 1矩阵尺寸信息 2元素的数据类型 3矩
30、阵的数据结构第108页/共1146页 本小节介绍如何获取矩阵的信息,包括矩阵的尺寸、元素的数据类型和矩阵的数据结构等。第109页/共1146页1矩阵尺寸信息 矩阵尺寸函数可以得到矩阵的形状和大小信息,这些函数如下表所示。第110页/共1146页第111页/共1146页2元素的数据类型 查询元素数据类型信息的部分函数如下表所示。第112页/共1146页第113页/共1146页3矩阵的数据结构 判断矩阵是否为某种指定数据结构的函数如下表所示。第114页/共1146页2.2.5 矩阵结构的改变改变矩阵结构的函数表 第115页/共1146页第116页/共1146页第117页/共1146页2.2.6 稀
31、疏矩阵 1稀疏矩阵的创建 2查看稀疏矩阵 3稀疏矩阵的运算规则第118页/共1146页 在MATLAB中,可以用满矩阵存储方 式和稀疏矩阵存储方式来存储矩阵。若一个矩阵只有少数的元素非零,称为稀 疏矩阵。稀疏矩阵非零元素及其对应的下 标来表示。用户可以创建双精度、复数和逻辑等类型 的稀疏矩阵。第119页/共1146页1稀疏矩阵的创建 在MATLAB中,用函数sparse()实现满矩阵到稀疏矩阵的转换。第120页/共1146页 在MATLAB中用函数full()实现稀疏矩阵 到满矩阵的转换。在MATLAB中,还可以用函数sparse()直接创建稀疏矩阵,其具体用法如下。S=sparse(i,j,
32、s,m,n),其中,i和j分别是稀疏矩阵非零元素的行和列下标,s为相应的非零元素的值,m和n分别是矩阵的行数和列数。第121页/共1146页 MATLAB还提供一些函数用于创建特殊稀疏矩阵,这些函数如下表所示。第122页/共1146页2查看稀疏矩阵 MATLAB提供一些函数用于查看稀疏矩阵的信息,如下表所示。第123页/共1146页 下面的例子都是基于MATLAB自带的稀疏矩阵west0479。第124页/共1146页第125页/共1146页第126页/共1146页3稀疏矩阵的运算规则 在MATLAB中的各种命令和函数都可以用于稀疏矩阵的运算,并且遵循如下的一些约定。(1)把矩阵变为标量或者定
33、长向量的函数总 是给出满矩阵;(2)对于标量或者定长向量变换到矩阵的函 数,如函数zeros()、ones()、eye()、rand()等总是给出满矩阵;第127页/共1146页(3)从矩阵到矩阵的变换函数将以原矩阵的 形式出现;(4)在参与矩阵扩展的子矩阵(如 A B;C D)中,只要有一个是稀疏矩阵,那 么所得的结果也是稀疏矩阵;(5)在矩阵引用中,将仍以原矩阵形式给出 结果。第128页/共1146页2.3 运算符和特殊符号运算符和特殊符号2.3.1 算数运算符2.3.2 关系运算符2.3.3 逻辑运算符2.3.4 运算优先级第129页/共1146页 在MATLAB中提供了丰富的运算符,包
34、括算数、关系和逻辑等3种运算符。第130页/共1146页2.3.1 算数运算符 在MATLAB中,算数运算符的用法和功能如下表所示。第131页/共1146页续表第132页/共1146页补充说明AB的用法如下:当A和B都为矩阵时,此运算无定义;当A和B都是标量时,表示标量A的B次 幂;当A是标量且B为矩阵时,表示标量A的 B中各元素次幂;当A为方阵且B为正整数时,表示矩阵A 的B次乘积;第133页/共1146页 当A为方阵且B为负整数时,表示矩阵A逆 的负B次乘积;当A为可对角化的方阵且B为非整数时,有 如下表达式:第134页/共1146页第135页/共1146页第136页/共1146页2.3.
35、2 关系运算符 MATLAB中关系运算符的用法和功能如下表所示。第137页/共1146页 值得注意的是,关系运算符只针对两个相同长度的矩阵,或其中之一是标量的情况进行运算。对于前者,是指两个矩阵的对应元素进行比 较,返回具有相同长度的矩阵;对于后者,是指这个标量与另一个矩阵的每元 素进行运算。关系运算C=f(A,B)的运算结果只有0和1两种情况,其中,函数f()表示关系运算符,0表示不满足条件,1表示满足条件。第138页/共1146页第139页/共1146页2.3.3 逻辑运算符 MATLAB提供元素方式和比特方式等逻辑运算符。元素方式逻辑运算符的用法和功能如下表所示,其中例子采用如下矩阵:第
36、140页/共1146页 元素方式逻辑运算符&、|和 与函数and()、or()和not()是等价的。比特方式逻辑运算符只接受逻辑和非负整数类型的输入变量,它是针对输入变量的二进制进行逻辑运算。第141页/共1146页 比特方式逻辑运算符的用法和功能如下表所示,表中例子采用A=28和B=200,其对应的二进制分别为11100和11001000。第142页/共1146页2.3.4 运算优先级 运算符的优先级决定表达式求值顺序;具有相同优先级的运算符从左到右依次进 行运算;不同优先级的运算符采用先进行优先高的 运算。第143页/共1146页运算符的优先等级表 由表中可以看到,括号的优先级别最高,因此
37、可以用括号来改变默认的优先等级。第144页/共1146页第145页/共1146页2.4 字符串处理函数字符串处理函数 2.4.1 字符串的构造 2.4.2 字符串的比较 2.4.3 字符串的查找和替换 2.4.4 字符串与数值间的转换第146页/共1146页 MATLAB提供了丰富的字符串操作,包括字符串的创建、合并、比较、查找以及与数值之间的转换。第147页/共1146页2.4.1 字符串的构造 在MATLAB中,可以用一对单引号来表示字符串。第148页/共1146页第149页/共1146页第150页/共1146页第151页/共1146页第152页/共1146页2.4.2 字符串的比较 1字
38、符串比较函数 2用关系运算符比较字符串第153页/共1146页 在MATLAB中提供了对字符串、字符串数组和字符子串的比较功能。第154页/共1146页1字符串比较函数 在MATLAB中,字符串比较函数如下表所示。第155页/共1146页第156页/共1146页2用关系运算符比较字符串 在MATLAB中,可以对字符串运用关系运算符,但要求两个字符串具有相同的长度,或者其中一个是标量。第157页/共1146页第158页/共1146页2.4.3 字符串的查找和替换 MATLAB提供的一些字符串查找和替换函数如下表所示。第159页/共1146页第160页/共1146页2.4.4 字符串与数值间的转换
39、 MATLAB提供的一些数值转换为字符串函数如下表所示。第161页/共1146页续表第162页/共1146页 MATLAB提供的一些字符串转换为数值函数如下表所示。第163页/共1146页第164页/共1146页习习 题题1计算复数3+4i与56i的乘积。2 构 建 结 构 体 Students,属 性 包 含 Name、Age 和Email,数据包括Zhang,18,、Wang,21,和Li,,构建后读取 所有Name属性值,并且修改Zhang的Age属 性值为19。第165页/共1146页3用满矩阵和稀疏矩阵存储方式分别构造下述矩 阵:4采用向量构造符得到向量1,5,9,41。5按水平和竖
40、直方向分别合并下述两个矩阵:第166页/共1146页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。第167页/共1146页14分别用函数strcat()和矩阵合并符合并如下字 符串:The picture is 和 very good。15 创 建 字 符 串 数 组,其 中 元 素 分
41、 别 为Picture 和 Pitch 。16在第14题结果中查找字符串e。17在第15题结果中匹配字符串Pi。18将字符串very good转换为等值的整数。19将十进制的50转换为二进制的字符串。20将十六进制的字符串50转换为三进制的整 数。第168页/共1146页第三章第三章 数学运算数学运算第169页/共1146页 本章将着重介绍MATLAB中与数学运算有关的函数和概念。在MATLAB中一切数据均能以矩阵的形式表示:针对矩阵整体的数学运算,称之为矩阵运 算;针对矩阵元素的数学运算,称之为矩阵元 素运算。第170页/共1146页目录目录 3.1 矩 阵 运 算 3.2 矩阵元素运算 习
42、 题第171页/共1146页3.1 矩矩 阵阵 运运 算算 3.1.1 矩阵分析 3.1.2 线性方程组 3.1.3 矩阵分解 3.1.4 矩阵的特征值和特征向量 3.1.5 矩阵相似变换 3.1.6 非线性运算第172页/共1146页 矩阵运算是线性代数中极其重要的部分,MATLAB具有强大的矩阵运算能力。第173页/共1146页3.1.1 矩阵分析 1向量间的距离 2矩阵的秩 3矩阵的行列式 4矩阵的迹 5矩阵的化零矩阵 6矩阵的正交空间 7矩阵的简化梯形形式 8矩阵空间之间的角度第174页/共1146页 MATLAB提供的部分矩阵分析函数如下表所示。第175页/共1146页1向量间的距离
43、第176页/共1146页第177页/共1146页2矩阵的秩 矩阵A中线性无关的列向量个数称为列秩,线性无关的行向量个数称为行秩。可以证明列秩与行秩是相等的。第178页/共1146页第179页/共1146页3矩阵的行列式第180页/共1146页第181页/共1146页4矩阵的迹 矩阵的迹定义为矩阵对角元素之和。在MATLAB中用函数trace()来计算矩阵的迹。第182页/共1146页5矩阵的化零矩阵 对于非满秩矩阵A,若存在矩阵Z使得AZ=0且ZTZ=I,则称矩阵Z为矩阵A的化零矩阵。在MATLAB中用函数null()来计算矩阵的化零矩阵。第183页/共1146页6矩阵的正交空间 矩阵A的正交
44、空间Q满足QTQ=I,且矩阵Q与A具有相同的列基底。第184页/共1146页第185页/共1146页7矩阵的简化梯形形式 矩阵A的简化梯形形式为 ,其中Ir为r阶单位矩阵。第186页/共1146页第187页/共1146页8矩阵空间之间的角度 矩阵空间之间的角度代表具有相同行数的两个矩阵线性相关程度,夹角越小代表线性相关度越高。第188页/共1146页第189页/共1146页3.1.2 线性方程组 线性方程组求解问题,可以表述为给定两个矩阵A和B,求解X使得AX=B或XA=B。XA=B可以表示为AY=B,且X=Y。下面仅讨论AX=B的情况。第190页/共1146页第191页/共1146页第192
45、页/共1146页3.1.3 矩阵分解 1Cholesky分解 2LU分解 3QR分解 4奇异值分解 5Schur分解第193页/共1146页 矩阵分解是把一个矩阵分解成比较简单或者对它性质比较熟悉的若干矩阵的乘积的形式。本小节将介绍几种矩阵分解的方法。第194页/共1146页矩阵分解函数表第195页/共1146页1Cholesky分解 Cholesky分解是把对称正定矩阵A表示为上三角矩阵R的转置与其本身的乘积,即A=RTR。第196页/共1146页第197页/共1146页 对于稀疏矩阵,MATLAB中用函数cholinc()计算不完全Cholesky分解,具体用法如下:R=full(chol
46、inc(sparse (X),DROPTOL),其中DROPTOL为不 完全Cholesky分解的丢失容限;R=full(cholinc(sparse(X),0),完 全Cholesky分解。第198页/共1146页2LU分解 高斯消去法又称LU分解,将任意一个方阵A分解为一个交换下三角 矩阵L和一个上三角矩阵U的乘积,即 A=LU。交换下三角矩阵为下三角矩阵经行变换的 结果。第199页/共1146页 LU分解在MATLAB中用函数lu()来实现,具体用法如下:L,U=lu(X),X为一个方阵,L为交换 下三角矩阵,U为上三角矩阵,满足关系 X=L*U;L,U,P=lu(X),X为一个方阵,L
47、为下 三角矩阵,U为上三角矩阵,P为置换矩 阵,满足关系P*X=L*U或X=P-1 *L*U。第200页/共1146页 考虑线性方程组AX=B和矩阵A的LU分解,线性方程组可改写成L*U*X=B,由于左除算符可以快速处理三角矩阵,因此:X=U(LB)矩阵的行列式和逆也可以利用LU分解来计算,如det(A)=det(L)*det(U)inv(A)=inv(U)*inv(L)第201页/共1146页第202页/共1146页 对于稀疏矩阵,在MATLAB中提供了函数luinc()来做不完全LU分解,其具体用法如下:L U=luinc(X,DROPTOL),其中X、L和U 的含义与函数lu()中的变量
48、相同,DROPTOL为 不完全LU分解的丢失容限。当DROPTOL设为0 时,退化为完全LU分解。L,U=luinc(X,0),0级不完全LU分解。L,U,P=luinc(X,0),0级不完全LU分解。第203页/共1146页3QR分解 QR分解就是将mn的矩阵A分解为mn的矩阵Q和nn的上三角矩阵R的乘积,且Q*Q=I,即A=Q*R。第204页/共1146页 在MATLAB中QR分解是由函数qr()来实现,其具体用法如下:Q,R=qr(A)满足A=Q*R。R=qr(A),返回上三角矩阵R。第205页/共1146页第206页/共1146页4奇异值分解 奇异值分解就是将mn的矩阵A分解为A=U*
49、S*V,其中U为mm的酉矩阵,V为nn的酉矩阵,S为mn的矩阵,并可如下表示:,其中 ,第207页/共1146页 在MATLAB中奇异值分解是由函数svd()来实现,其具体用法如下:第208页/共1146页第209页/共1146页5Schur分解 Schur分解就是将复方阵A分解为A=U*L*U,其中U为酉矩阵,L为上(下)三角矩阵,其对角线元素为A的特征值。第210页/共1146页 在MATLAB中Schur分解是由函数schur()来实现,其具体用法如下:U,L=schur(A),满足A=U*L*U,其中L为上三角矩阵。L=schur(A),返回上三角矩阵L。第211页/共1146页第212页/共1146页3.1.4 矩阵的特征值和特征向量 方阵A的特征值和其对应的特征向量满足下式:A*=*第213页/共1146页 在MATLAB中用函数eig()来计算特征值和其对应的特征向量,其具体用法如下:d=eig(A),返回矩阵A的所有特征值。V,D=eig(A),返回矩阵A的特征值和 特征向量。第214页/共1146页第215页/共1146页第216页/共1146页3.1.5 矩阵相似变换 1对角阵变换 2Jordan变换第217页/共1146页 矩阵相似变换是指,对于方阵A和非奇异矩阵B可得到相似矩阵X=B-1*A*B。第218页/共1146页1对角阵变换 对于方阵A,若V D=
限制150内