1、数值计算软件实践.1第一讲 绪论数值计算软件的产生原因和应用领域。典型数值计算软件的基础知识。我们要学习的数值计算软件是工程科学工程科学界界的周润发。.2现实中的一些问题S矩阵运算z加减z乘z除z逆z1122abba 11111111.mnnnmmmnaabbaabb 11282.3现实中的一些问题S函数图形的绘制zy=sin(t)；y=cos(t)zS数据处理z插值、拟合z数值积分、微分zS微分方程的数值解.4现实中的一些问题S物理现象的数值模拟zBP墨西哥湾漏油事件S控制系统的设计.5数值计算软件实践目标：了解MATLAB，利用该软件进行矩阵运算，简单编程，简单的数据处理，基本图形绘制，了
2、解动态系统仿真工具。MATLAB是学习和科研的工具工具。（一）概述（一）概述（二）（二）MATLAB基础知识基础知识（三）矩阵与数组（三）矩阵与数组（四）绘图功能（四）绘图功能（五）数值与符号运算（五）数值与符号运算（六）（六）M语言程序与图形用户界面设计语言程序与图形用户界面设计（七）（七）Simulink仿真环境仿真环境.61 MATLAB基础知识SMATLAB是MathWorks公司开发的科学与工程计算软件，它以矩阵运算为基础，把计算、绘图及动态系统仿真等功能有机地融合在一起，同时，它又具有程序设计语言（M语言）的基本特征。S它将一个优秀软件的易用性与可靠性、通用性与专业性 、一般目的的
3、应用一般目的的应用与高深的科学技高深的科学技术应用术应用有机的结合。.71.1 MATLAB简介S1980年，美国新墨西哥大学计算机科学系主任Cleve 着手编写供学生使用的子程序接口程序S设计者的初衷是为解决线性代数线性代数课程的矩阵运算问题S取名MATLAB即MATrix LABoratory 矩阵实验室S1984年，推出了MATLAB第一个商业版本.81.1 MATLAB简介S1992年MATLAB 4.0版；S1997年MATLAB 5.0版；S2000年MATLAB 6.0版；S2004年MATLAB 7.0版；S2005年MATLAB 7.1版；S随后每年发布两次以年份命名的版本：
4、z如R2006a、R2006b、R2007a、R2007b和最新的版本R2010a。.91.1 MATLAB简介SMATLAB 的应用：z数值和符号计算z工程与科学绘图z控制系统的设计与仿真z动态系统建模、分析和仿真z通讯系统设计与仿真z数字图像处理z数字信号处理z财务与金融工程z.101.2 MATLAB安装S普通PC基本都能满足MATLAB的安装环境要求和安装普通软件没有任何区别！.111.3 MATLAB集成开发环境.121.3.1 工作界面S工作界面z菜单z工具栏z命令窗口z历史命令窗口z当前工作目录窗口z工作空间窗口.131.3.2 命令窗口SCommand WindowS用于输入命
5、令输入命令并显示除图形以外的所有执行结果。SMATLAB用“”作为运算提示符，表明系统处于准备状态。S当在提示符后输入一段程序或一段运算式后回车，MATLAB会给出计算结果，并再次进入准备状态。.141.3.2 命令窗口S命令以回车结束。命令之间以逗号逗号或者分号分号分隔。z例如：x=720, y=86x = 720y = 86x=720; y=86y = 86S若一个命令行很长，可以使用续行符，即三个三个小黑点。z例如：z=1+1/(1*2)+1/(1*2*3)+. 1/(1*2*3*4*5)命令输入采用英文字符英文字符.151.3.2 命令窗口S单击命令窗口右上角的“ ”按钮，可以使命令窗
6、口脱离主窗口而成为一个独立的窗口。.161.3.3 历史命令窗口SCommand HistoryS用于记录所有执行过的命令（可以单个或全部删除）S标明使用时间（启动软件时的时刻，只标一次启动软件时的时刻，只标一次）S可以通过用鼠标双击双击某一历史命令来重新执行该命令S可以成为一个独立的窗口通过小键盘上的“”和“”键可以在命令窗口命令窗口中依次依次显示一行一行历史命令.171.3.4 工作空间窗口SWorkspaceS工作空间窗口将显示当前当前内存内存中所有的MATLAB变量的变量名、数据结构、字节数以及类型等信息。.181.3.4 当前工作目录SCurrent DirectoryS当前工作目录
7、是指MATLAB运行文件时的目录。只有在当前工作目录或搜索路径下的文件、函数才可以被运行或调用。S如果没有特殊指明，数据文件也将存放在当前目录下。 .191.3.5 Start按钮SStart按钮.201.4 MATLAB通用命令S通用命令是MATLAB中经常使用的一组命令，这些命令可以用来管理目录、函数、变量、工作空间、文件和窗口等。z常用命令z输入内容的编辑z操作符和运算符z搜索路径与扩展.211.4.1 常用命令helphelp在线帮助文件在线帮助文件doc装入超文本说明whatM、MAT、MEX文件的目录列表type列出M文件demo运行演示程序path控制MATLAB的搜索路径who
8、列出当前变量whos列出当前变量（长表）load从磁盘文件中恢复变量save保存工作空间变量pack整理工作空间内存clearclear从内存中清除变量和函数从内存中清除变量和函数cd改变当前工作目录dir目录列表!执行DOS操作系统命令clcclc清命令窗口清命令窗口home光标置左上角format设置输出格式quit退出MATLABinfo显示系统和公司信息subscribe成为MATLAB的订购用户hostidMATLAB主服务程序的识别代号whatsnew在说明书中未包含的新信息ver版本信息.221.4.2 输入内容的编辑S命令窗口中，MATLAB提供了控制光标位置和进行简单编辑的键
9、盘按键，部分常用的键盘按键及其功能如下表所示。.231.4.3 操作符和运算符+加减*矩阵乘法.*数组乘法数组乘法矩阵幂矩阵幂.数组幂数组幂左除或反斜杠左除或反斜杠/右除或斜杠右除或斜杠./数组除数组除KronKronecker张量积张量积:冒号冒号( )圆括号，指定运算的优先级 方括号，定义矩阵.小数点小数点.父目录父目录续行符号续行符号,逗号逗号;分号，区分行及取消运行结果显示%注释标记!感叹号感叹号转置或引用转置或引用=赋值= =相等关系操作符关系操作符&逻辑与逻辑与|逻辑或逻辑或逻辑非逻辑非xor逻辑异或逻辑异或.241.4.3 操作符和运算符S常用运算的MATLAB表示 数学表示数学
10、表示 MATLAB表示表示.251.4.4 搜索路径与扩展S当用户在命令窗口输入一条命令后，MATLAB将按照下列顺序搜索所输入的命令：z检查该命令是不是一个变量。z检查该命令是不是一个内部函数。z检查该命令是否为当前目录下的M文件。z检查该命令是否为搜索路径中其他目录下的M文件。z都找不到时系统将报错！.261.4.4 搜索路径与扩展S查看搜索路径z选择MATLAB主窗口中的【File】|【Set Path】菜单，进入【设置搜索路径】对话框。通过该对话框可为MATLAB添加或删除搜索路径。S自己书写的函数有可能并没有保存在搜索路径下。要解决这个问题，只需把程序所在的目录扩展成搜索路径即可。.
11、271.4.4 搜索路径与扩展S扩展扩展搜索路径z选择主窗口中的【File】|【Set Path】|【Add Folder】。z在命令窗口执行“path(path,D:Study )”。z在命令窗口执行“addpath D:Study- end”，将新的目录加到整个搜索路径的末尾。如果将end改为begin，可以将新的目录加到整个搜索路径的开始。.281.5 帮助系统SMATLAB为用户提供了非常完善非常完善的帮助系统。z联机帮助系统z命令窗口查询帮助系统z联机演示系统 (Demo)Demo给了大家一个模仿提高的机会一定好好利用！.291.5.1 联机帮助系统S进入联机帮助系统有四种方法z按F
12、1键z选择Help菜单中的“MATLAB Help”选项。z单击主窗口工具栏中的Help按钮（小问号）。z在命令窗口中输入“helpwin”、“helpdesk”或“doc”命令。.301.5.2 命令窗口查询帮助系统S帮助命令zhelp 命令z查询函数语法的最基本方法z结果直接显示在命令窗口。z通过help加函数名来显示该函数的帮助说明。z例如，为了显示sin函数的使用方法与功能help sinzlookfor命令z该命令可以帮助用户了解函数的用法。例如，输入命令lookfor sin.311.5.3 联机演示系统S进入联机演示系统的方法z选择MATLAB主窗口菜单的【Help】|【Demo
13、s】选项；z在命令窗口输入demos；z直接在帮助页面上选择Demos页。.321.6 MATLAB 示例S下面以一个简单的例子展示如何使用MATLAB进行简单的应用。.331.6.1 数值计算S求解线性方程组z在MATLAB命令窗口输入命令：za=2,3,-1;3,-5,3;6,3,-8; zb=7;8;9; zx=inv(a)*bS也可以通过符号计算来解此方程zsyms x1 x2 x3zx1,x2,x3=solve(2*x1+3*x2-x3-7,3*x1-5*x2+3*x3-8,6*x1+3*x2-8*x3-9).341.6.2 绘图功能S绘制函数y=300sin(x)/x和y=x2的曲
14、线S只需要在命令窗口中输入：zx=-20:0.1:20;zplot(x,300*sin(x)./x,:,x,x.2);.351.6.3 编程计算S汽车性能计算软件.36小结S初步了解MATLAB。S熟悉MATLAB界面，进行简单操作。数值计算软件实践.37第二讲 基础知识介绍MATLAB的变量、数据类型以及矩阵基本操作。.382.1 MATLAB的变量S2.1.1 变量命名与赋值S2.1.2 预定义变量S2.1.3 变量的管理.392.1.1 变量命名与赋值S变量命名z在MATLAB中，变量名是以字母开头，后接字母、数字或下画线的字符序列，最多63个字符。例如：x、x_1、xy1z在MATLA
15、B中，变量名区分字母的大小写。不支持中中文文字字符符语句.402.1.1 变量命名与赋值S赋值zMATLAB赋值语句有两种格式：变量名 = 表达式；在该语句形式下，MATLAB将右表达式的值赋给左边的变量。表达式表达式；在该语句形式下， MATLAB将表达式的值赋给的预定义变量ans。z注释以%开头。.412.1.2 预定义变量S预定义变量可以被重新赋值。.422.1.3 变量的管理S显示显示内存变量z显示在工作空间工作空间中已经驻留的变量清单。z主要命令：qwho % %显示变量名显示变量名qwhoszwhos在给出变量名的同时，还给出它们的维数、所占字节数以及变量的类型。.432.1.3
16、变量的管理S删除删除MATLAB工作空间中的变量z全部删除命令：clearS工作空间窗口专门用于内存变量的管理。S变量编辑器Array Editor.442.1.3 变量的管理S保存保存内存变量z利用MAT文件可以把当前工作空间中的一些有用变量长久地保留下来。zMAT是MATLAB保存数据的一种标准格式二进制文件，扩展名是.mat。 兼容性不好.452.1.3 变量的管理S点击工作空间窗口的save按钮则保存所有变量S命令命令生成与导入MAT文件zsave 文件名 变量名表load 文件名 变量名表z例如：save mydata a bload mydataz实现变量定点定点保存.462.2
17、MATLAB的数据类型S2.2.1 数值类型S2.2.2 逻辑类型S2.2.3 字符类型和字符串S2.2.4 结构体类型.472.2.1 数值类型S数值的类型z整数z浮点数zInfzNaNS数值的输出格式.482.2.1 数值类型S整数zMATLAB支持1、2、4和8字节的有符号整数和无符号整数。z数据类型的名称、表示范围和转换函数转换函数如下表所示：.492.2.1 数值类型S浮点数zMATLAB有单精度和双精度两种浮点数。z浮点数的存储空间、表示范围和转换函数如下表所示： .502.2.1 数值类型SInf ( inf ) & NaNzInf和-Inf分别表示正无穷大和负无穷大。除法运算中
18、除数为0或者运算结果溢出都会导致Inf或-Inf的运行结果。zNaN在MATLAB中表示一个既不是实数也不是复数的数值。MATLAB 对除数为0的情况会发出警告但运算仍可进行.512.2.1 数值类型S数值的输出格式S用“format”命令设置数值输出格式。zformat 格式符格式符短格式短格式(short)短格式短格式e方式方式(short e)长格式长格式(long)长格式长格式e方式方式(long e)银行格式银行格式(bank) 十六进制格式十六进制格式(hex) S“format”命令只影响数值的输出格式输出格式，而不影响数据的计算和存储。.522.2.2 逻辑类型SMATLAB中
19、的逻辑类型（logical）包含true和false，分别由1和0表示。在MATLAB中用函数logical()将任何非零的数值非零的数值转换为true（即1），将数值数值0转换为false（即0）。.532.2.3 字符类型和字符串S在MATLAB中，数据类型（char）表示一个字符。S一个char类型的1n数组称为字符串string。.542.2.4 结构体类型S结构体类型是一种由若干属性（field）组成的MATLAB数组，其中的每个属性可以是任意数据类型。.552.3 MATLAB矩阵及其基本操作S在MATLAB中，所有的数据均以二维、三维或高维矩阵的形式存储，每个矩阵的单元可以是数值
20、类型、逻辑类型、字符类型或者其他任何数据类型。 S2.3.1 矩阵的构造S2.3.2 矩阵大小的改变S2.3.3 矩阵元素的操作S2.3.4 矩阵信息的获取S2.3.5 矩阵结构的改变.562.3.1 矩阵的构造S简单矩阵构造S特殊矩阵构造S行向量构造.572.3.1 矩阵的构造S简单矩阵构造zMATLAB矩阵构造符为“”。同一行元素依次放入矩阵构造符内，并且以空格空格或者逗逗号号分隔；行与行之间用分号分号分隔。z在MATLAB中，生成矩阵时无须对矩阵的维数和类型进行说明，MATLAB会根据用户所输入的内容自动进行配置。z矩阵元素也可是表达式(MATLAB将自动计算结果).582.3.1 矩阵
21、的构造S简单矩阵构造z已建好的小矩阵小矩阵可拼接成大矩阵大矩阵，从而利用已建好的矩阵建立更大的矩阵。.592.3.1 矩阵的构造S特殊矩阵构造z在MATLAB中还提供一些函数函数用来构造特殊矩阵。.602.3.1 矩阵的构造S特殊矩阵构造（续表）.612.3.1 矩阵的构造S特殊矩阵构造示例.622.3.1 矩阵的构造S行向量构造zMATLAB提供了向量构造符“:”。z常用的用法如下：a:b以a为起点，且所有取值在a与b之间的向量。 a:s:b 返回以a为起点，以s为步长，且所有取值在a与b之间的向量。 默认步长是默认步长是1 1.632.3.2 矩阵大小的改变S矩阵的合并合并S矩阵行列的删除
22、删除.642.3.2 矩阵大小的改变S矩阵的合并z使用构造符“”z矩阵合并函数.652.3.2 矩阵大小的改变S矩阵行列的删除z将要删除矩阵的行或者列赋予一个空矩阵 。z命令为：删除行q矩阵名 (n,:) = 删除列q矩阵名 (:,n) = .662.3.3 矩阵元素的操作S本小节将介绍存取矩阵元素元素值的方法。 z访问单个元素z线性引用元素z访问多个元素.672.3.3 矩阵元素的操作S访问单个元素z对于二维矩阵A，访问其第m行第n列元素的命令：A(m,n)z对元素进行重新赋值的命令：A(m,n) = 55.682.3.3 矩阵元素的操作S线性引用元素z对于矩阵A，线性引用元素格式为A(k)
23、。通常用于行向量或列向量，也可用于二维矩阵。zMATLAB按列优先排列的一个长列向量格式（线性引用元素）来存储存储矩阵元素。.692.3.3 矩阵元素的操作S访问多个元素z操作符“:”z若A是二维矩阵，其主要用法如下：A(:,:)返回矩阵A的所有元素A(i,:)返回矩阵A第i行的所有元素A(i,k1:k2)返回矩阵A第i行的自k1到k2列的所有元素A(:,j)返回矩阵A第j列的所有元素A(k1:k2,j)返回矩阵A第j列的自k1到k2行的所有元素z若A是多维矩阵，也可以通过类似的方法实现对其访问。.702.3.4 矩阵信息的获取S本小节介绍如何获取矩阵的信息。 z矩阵尺寸信息z元素的数据类型z
24、矩阵的数据结构.712.3.4 矩阵信息的获取S矩阵尺寸信息z矩阵尺寸函数矩阵尺寸函数.722.3.4 矩阵信息的获取S元素的数据类型z元素数据类型函数.732.3.4 矩阵信息的获取S矩阵的数据结构z数据结构函数.742.3.5 矩阵结构的改变S改变矩阵结构的函数.75小结SMATLAB的变量SMATLAB的数据类型SMATLAB的矩阵基本操作数值计算软件实践.76第三讲 矩阵与数组介绍MATLAB中的矩阵和数组常规运算法则，多项式的相关运算，矩阵的特征分析和分解。.77主要内容S3.1 矩阵运算S3.2 数组运算S3.3 多项式运算S3.4 线性方程组求解S3.5 特殊矩阵生成S3.6 矩
25、阵分析S3.7 矩阵分解S3.8 矩阵的特征值和特征向量S3.9 矩阵相似变换.783.1 矩阵运算S3.1.1 矩阵加、减（，）运算S3.1.2 矩阵乘（*）运算S3.1.3 矩阵除（/，）运算S3.1.4 矩阵乘方 An，apS3.1.5 矩阵关系运算符S3.1.6 矩阵逻辑运算符.793.1.1 矩阵加、减（，）运算S两矩阵对应元素对应元素相加减。S标量可以和任意矩阵相加减，且与所有元素分别进行加减操作。.803.1.2 矩阵乘（*）运算SA*B：A矩阵的列数列数必须等于B矩阵的行数行数S标量则与任意矩阵的所有元素相乘。.813.1.3 矩阵除（/，）运算S矩阵右除法右除法（A / B）
26、结果等效于(BA)S矩阵左除法左除法（A B）如果除法中有一如果除法中有一个是个是标量标量 a ，A / aa Aa / AA a是否等价？是否等价？A B的结果的结果是是AX = B解解Error .823.1.4 矩阵乘方（An，ap）SAn：表示A自乘n次幂；S如果p是矩阵，a是标量，ap 使用特征值特征值和特特征向量征向量自乘到p各元素次幂。乘方运算乘方运算比较费时间比较费时间.833.1.5 矩阵关系运算符S矩阵关系运算符S矩阵与矩阵关系：两个矩阵长度相同长度相同，对应元素比较，返回具有相同长度的矩阵；S矩阵与标量关系：标量与矩阵的所有元素所有元素进行运算。.843.1.5 矩阵关系
27、运算符S关系运算的运算结果只有0和1两种情况S0 ：不满足条件S1 ：满足条件S不能进行运算时MATLAB将报错.853.1.6 矩阵逻辑运算符SMATLAB提供元素方式元素方式和比特方式比特方式逻辑运算符。S元素方式元素方式逻辑运算符：.863.1.6 矩阵逻辑运算符S比特方式比特方式逻辑运算符只接受逻辑和非负逻辑和非负整数整数类型的输入变量，它是针对输入变量的二进制二进制表示进行逻辑运算。.873.2 数组运算S数组运算指矩阵元素和元素元素和元素之间的算术运算 S3.2.1 数组加减（+，-）S3.2.2 数组乘除（，/，）S3.2.3 数组乘方（）S3.2.4 运算优先级.883.2.1
28、 数组加减（+，-）S对应元素相加减（与矩阵加减等效）.893.2.2 数组乘除（*，./，.）S数组乘（.*）点乘z数组对应元素相乘对应元素相乘z数组必须有相同的行和列相同的行和列.903.2.2 数组乘除（，./，.）SA./B =B.A A的元素被B的对应元素除SB./A =A.B B的元素被A的对应元素除给出给出A A，B B对应元素间的商对应元素间的商.913.2.3 数组乘方（）S元素对元素元素对元素的幂.923.2.4 运算优先级S运算符的优先级决定表达式求值顺序.933.2.4 运算优先级S优先级不同先进行优先高的运算S优先级相同从左到右依次进行运算S可以用括号括号来改变默认的
29、运算顺序推荐利用括号来保推荐利用括号来保证运算顺序正确证运算顺序正确.943.3 多项式运算S3.3.1 多项式表示法S3.3.2 多项式求值S3.3.3 多项式求根S3.3.4 多项式乘法和除法S3.3.4 多项式微积分.953.3.1 多项式表示法S多项式表示S行向量行向量，其元素中按多项式降幂排列。S例如： f(x)=anxn+an-1xn-1+a0行向量行向量： an an-1 a0S函数poly2str()将多项式系数向量转换为完整形式。zf,len=poly2str(A,x); .963.3.2 多项式求值S函数polyval()用来计算多项式的值Sy = polyval(p,x)
30、 %p为多项式系数行向量，x代入多项式的值。SY = polyvalm(p,X) %把矩阵X代入多项式p中进行计算。.973.3.3 多项式求根S多项式求根Sr=roots(a) %结果以列向量列向量形式显示所有根。S函数poly()实现由根创建多项式Sp=poly(r) %输入r是多项式所有根，返回值为多项式的系数向量。.983.3.4 多项式乘法和除法S函数命令 conv() 和 deconv() 进行多项式乘法和除法Sw = conv(u,v) %返回结果多项式相乘后的系数系数行向量。Sq,r = deconv(u, v) %q为u除v后的整数，r为余数。.993.3.4 多项式微积分S
31、多项式的微分S函数polyder()zk = polyder(p) %返回多项式 p 微分的系数向量zk = polyder(a,b) %返回多项式 a*b 微分的系数向量zq,d = polyder(a,b) %返回多项式 a/b 微分的系数向量，q为分子，d为分母.1003.3.4 多项式微积分S多项式积分S函数polyint()计算多项式的不定积分Ss=polyint(p,k k) %返回多项式 p 不定积分的系数向量系数向量。.1013.4 线性方程组求解S3.4.1 矩阵求逆与线性代数方程组求解S3.4.2 利用左除运算符求解线性方程组.1023.4.1 矩阵求逆及线性代数方程组求解
32、S矩阵求逆矩阵求逆Sinv(A) %计算方阵A的逆矩阵S利用矩阵求逆方法解线性方程组Ax=b123123123231352253xxxxxxxxx.1033.4.2 利用左除运算符求解线性方程组S对于线性方程组Ax = b，可以利用左除运算符“”求解： x = A b。左除法要比求逆法耗时略多左除法要比求逆法耗时略多.1043.5 特殊矩阵生成S3.5.1 通用的特殊矩阵S3.5.2 面向特定应用的特殊矩阵.1053.5.1 通用的特殊矩阵S函数zeros全0矩阵。S函数ones全1矩阵。S函数eye单位矩阵。S函数rand产生01均匀分布的随机矩阵。S函数randn产生均值为0、方差为1的标
33、准正态分布随机矩阵。.1063.5.2 面向特定应用的特殊矩阵S魔方矩阵魔方矩阵，即每行、每列及两条对角线上的元素和都相等。S函数magic(n)，其功能是生成一个n阶魔方阵。S范德蒙范德蒙（Vandermonde）矩阵的最后一列全为1，倒数第倒数第二列为一个指定的向量二列为一个指定的向量，其他各列是其后列与倒数第二列的点乘积。S函数vander(V)生成以向量向量V V为基础向量的范得蒙矩阵。 .1073.5.2 面向特定应用的特殊矩阵S托普利兹托普利兹（Toeplitz）矩阵除第一行第一列外，其他每个元素都与左上角的元素相同。S函数是toeplitz(x,y)生成以x为第一列、y为第一行的
34、托普利兹矩阵。这里x、y均为向量，两者不必等长。toeplitz(x)用向量x生成一个对称的托普利兹矩阵。S希尔伯特希尔伯特（Hilbert）矩阵是一种数学变换矩阵，它的每个元素hij = 1/(i + j 1)。S函数是hilb(n)。 invhilb(n)为专门求希尔伯特矩阵的逆的函数。.1083.5.2 面向特定应用的特殊矩阵S帕斯卡帕斯卡（Pascal）矩阵是由杨辉三角形表组成的矩阵。函数pascal(n)生成一个n阶帕斯卡矩阵。.1093.6 矩阵分析S3.6.1 向量间的距离S3.6.2 矩阵的秩S3.6.3 矩阵的行列式S3.6.4 矩阵的迹S3.6.5 矩阵的化零矩阵S3.6.
35、6 矩阵的正交空间.1103.6.1 向量间的距离SMATLAB 函数norm求向量间距离 norm(a) norm(a,2).1113.6.2 矩阵的秩S矩阵 A中线性无关的列向量个数称为列秩，线性无关的行向量个数称为行秩。列秩与行秩是相等的。S函数rank( )计算矩阵的秩。.1123.6.3 矩阵的行列式SMATLAB提供了det( )函数来求矩阵的行列式。.1133.6.4 矩阵的迹S矩阵的迹指矩阵对角元素之和。S在MATLAB中用函数trace( )来计算矩阵的迹。 .1143.6.5 矩阵的化零矩阵S对于非满秩矩阵A，若存在矩阵Z使得AZ = 0 且 ZTZ = I，则称矩阵Z为矩
36、阵A的化零矩阵化零矩阵。S在MATLAB中用函数null( )来计算矩阵的化零矩阵。 .1153.6.6 矩阵的正交空间S矩阵A的正交空间Q满足QTQ = I，且矩阵Q与A具有相同的列基底。S在MATLAB中用函数orth()来计算矩阵的正交空间。BTB.1163.7 矩阵分解S矩阵分解是指根据一定的原理用某种算法将一个矩阵分解成若干个矩阵的乘积。S3.7.1 LU分解S3.7.2 QR分解S3.7.3 Cholesky分解.1173.7.1 LU分解S矩阵的LU分解又称Gauss消去分解或三角分解，就是将一个方阵表示为一个行交换下三角矩阵和一个上三角矩阵的乘积形式。S函数luzL,U = l
37、u(X) %X = LU，上三角阵U和一个变换形式的下三角阵L(行交换) zL,U,P = lu(X) %PX = LU，上三角阵U和一个下三角阵L以及一个置换矩阵P.1183.7.2 QR分解S对矩阵X进行QR分解，就是把X分解为一个正交矩阵Q和一个上三角矩阵R的乘积形式。S函数qrzQ,R=qr(X) %X=QR，正交矩阵Q和一个上三角矩阵RzQ,R,E=qr(X) %XE=QR，正交矩阵Q、一个上三角矩阵R以及一个置换矩阵E.1193.7.3 Cholesky分解SCholesky分解是把对称正定矩阵A表示为上三角矩阵R的转置与其本身的乘积，即A = RTR。S函数cholzR = ch
38、ol(X) % RR = X。zR,p = chol(X) %当X为对称正定的，则p = 0，R与上述格式得到的结果相同；否则 p为一个正整数。如果X为满秩矩阵，则R为一个阶数为p 1的上三角阵。 .1203.8 矩阵的特征值和特征向量S方阵A的特征值和其对应的特征向量X满足式：A*X=*X 。S函数eig() %计算特征值和其对应的特征向量zd = eig(A) %返回A的特征值特征值。zV,D = eig(A) %返回A的特征值特征值D和特征向量和特征向量V。 .1213.9 矩阵相似变换S矩阵相似变换是指，对于方阵 A和非奇异矩阵 B可得到相似矩阵X=B-1*A*B。S3.9.1 对角阵
39、变换S3.9.2 Jordan变换.1223.9.1 对角阵变换S对于方阵A，若V D=eig(A)得到的矩阵V非奇异，则A可经过相似变换得到对角阵，即D=V-1*A*V，也称矩阵A可对角化。 V-1*A*V.1233.9.2 Jordan变换S对于方阵A，若V D=eig(A)得到的矩阵 V 奇异，则A经过相似变换将不能得到对角阵，只能得到其对应的Jordan标准型。S函数jordan() ： zV,D = jordan(A) %满足 D=V-1*A*V。zD = jordan(A) %返回矩阵A对应的Jordan标准型。 .124小结S本节介绍了MATLAB矩阵与数值运算功能，通过学习应该
40、做到：z掌握矩阵及数组运算（+、-、*）z了解多项式运算z掌握线性方程组求解方法z了解特殊分析与分解命令z了解矩阵的特征值和特征向量的求法z了解矩阵的变换命令数值计算软件实践.125第四讲 绘图功能MATLAB有着强大的函数绘图功能，不需要过多地考虑绘图细节，只需给出一些基本参数就能得到所需图形。.126主要内容S4.1 二维绘图S4.2 三维绘图S4.3 图形窗口*.1274.1 二维绘图S二维图形的绘制是MATLAB语言图形处理的基础。S4.1.1 二维图绘制与修饰S4.1.2 特殊二维绘图函数S4.1.3 函数fplotS4.1.4 函数ezplot.1284.1.1 二维图绘制与修饰S
41、plot(x) %x为向量，x的元素值为纵坐元素值为纵坐标标。Splot自动打开一个图形窗口Figure。S如果图形窗口已经存在，plot则清除清除当前图形，绘制新图形。.1294.1.1 二维图绘制与修饰Splot(x,y) %基本格式，以y(x)的函数关系作出直角坐标图，如果y为nm的矩阵，则以x为自变量，作出m条曲线。.1304.1.1 二维图绘制与修饰Splot(x,y,s) %字符串字符串s用来定义曲线特征。.1314.1.1 二维图绘制与修饰Splot(x,y,s) %s可以是各种选项的组合组合。m-rok.1324.1.1 二维图绘制与修饰S单窗口多曲线绘图Splot(x1,y1
42、,x2,y2,) %多曲线绘图格式 plot(t,y,y1,y2).1334.1.1 二维图绘制与修饰S单窗口多曲线分图分图绘图Ssubplot(m,n,p) % 子图子图分割分割命令命令S从左至右，从上至下从左至右，从上至下排列 .1344.1.1 二维图绘制与修饰S单窗口多曲线分图绘图subplot(2,2,3,4).1354.1.1 二维图绘制与修饰S多窗口绘图Sfigure(n) %创建窗口函数，n为窗口序号序号.1364.1.1 二维图绘制与修饰S图形加注功能S将标题、坐标轴标记、网格线及文字注释加注到图形上的函数函数为：ztitle 给图形加标题zxlabel 给x轴加标注zyla
43、bel 给y轴加标注ztext 在图形指定位置加标注zgtext 将标注加到图形任意位置zgrid on/off 打开、关闭坐标网格线zlegend 添加图例zaxis 控制坐标轴的刻度.1374.1.1 二维图绘制与修饰S图形加注功能.1384.1.1 二维图绘制与修饰Sfill绘制二维多边形并填充颜色.1394.1.2 特殊二维绘图函数Sbar 绘制直方图Spolar 绘制极坐标图Shist 绘制统计直方图Sstairs 绘制阶梯图Sstem 绘制火柴杆图Srose 绘制统计扇形图Scomet 绘制彗星曲线Sscatter 离散点图Scompass 复数向量图(罗盘图)Sfeather 复
44、数向量投影图(羽毛图)Squiver 向量场图Sarea 区域图Spie 饼图Sconvhull 凸壳图 Serrorbar 绘制误差棒图.1404.1.2 特殊二维绘图函数S绘制极坐标图极坐标以弧度为单位极坐标以弧度为单位.1414.1.2 特殊二维绘图函数S绘制饼图.1424.1.3 函数fplot ()Splot将外部的离散数据转化为图形。数据数据点点越多绘制的曲线越精确，越能反映数据的规律。Sfplot构建所需要绘制函数的数据，同时也可以绘图，精确程度由误差误差决定。 Sx,y = fplot(fun,lims) %返回函数fun在x区间lims = xmin xmax的绘图点值。Sf
45、un用字符串字符串表示。.1434.1.3 函数fplot ().1444.1.4 函数ezplot()S函数ezplot用于绘制函数在某一自变量区域内的图形。zezplot(fun) %绘制默认范围内的图形zezplot(fun,min,max) %绘制指定范围内的图形zezplot(fun,xmin,xmax,ymin,ymax) %绘制函数f(x,y)指定范围内的图形zezplot(funx,funy) %绘制默认范围的参数方程组x=x(t),y=y(t)的图形zezplot(funx,funy,tmin,tmax) %绘制指定范围的参数方程组x=x(t),y=y(t)的图形.1454.
46、1.4 函数ezplot().1464.2 三维绘图S4.2.1 三维图绘制基本函数S4.2.2 三维网格图绘制S4.2.3 三维曲面图绘制S4.2.4 图形修饰命令.1474.2.1 三维图绘制基本函数Splot3 %基本的三维图形指令。Splot3(x,y,z) %x,y,z是长度相同长度相同的向量或维数相同维数相同的矩阵S二维图形的所有基本特性对三维图形全都适用。Splot3(x,y,z,s) Splot3(x1,y1,z1,s1, x2,y2,z2,s2, ) %多线条绘图Sgrid on/off %绘制三维网格Stext(x,y,z,string)%三维图形标注S子图和多窗口。.14
47、84.2.1 三维图绘制基本函数.1494.2.1 三维图绘制基本函数S三维饼图pie3.1504.2.2 三维网格图绘制S三维网格图函数mesh可绘出在某一区间内的完整曲面完整曲面Smesh(z) %z为nm的矩阵，x与y坐标为元素的下标Smesh(x,y,z) %x,y,z分别为三维空间的坐标位置.1514.2.2 三维网格图绘制S三维网格图绘制步骤z生成坐标： X,Y=meshgrid(x,y)z参数x和y为向量向量z表达式点运算点运算：Z=X.2+Y.2zX，Y是nm的矩阵，维数由向量尺寸决定z默认37.5，俯角30。.1524.2.3 三维曲面图绘制S曲面绘制函数surf绘制着色的三
48、维曲面Ssurf(X,Y,Z) %X,Y,Z为图形坐标向量曲面图线条是黑色空挡有颜色网格图线条有颜色空挡是白色.1534.2.3 三维曲面图绘制S带等高线的曲面图绘制函数surfc(X,Y,Z)S被光照射带阴影的曲面图surfl(X,Y,Z) .1544.2.4 图形修饰命令Scolormap(MAP) %色图设定函数Scolorbar() %标注Z Z图例图例Sshading faceted / flat / interp %网格修饰Shidden on/off %透视与消隐（网格图）Sview(az,el) %视角修饰，az 为方位角，el 为俯视角，缺省值：az=-37.5；el=30S
49、waterfall / meshz %水线修饰Scontour(Z,n) / contour3(Z,n) %绘制n条二维 /三维等高线；.1554.2.4 图形修饰命令Scolormap(MAP)zhsv %饱和值色图zgray %线性灰度色图zhot %暖色色图zcool %冷色色图zbone %兰色调灰色图zcopper %铜色色图zpink %粉红色图zprism %光谱色图zjet %饱和值色图IIzflag %红,白,蓝交替色图.1564.2.4 图形修饰命令Sshading faceted / flat / interp Shidden on / offSview(az,el) S
50、meshzScontour3.1574.3 图形窗口菜单工具栏绘图区.1584.3.1 菜单操作新建文件M文件、图形窗口、工作空间的变量或用户界面生成M-函数文件用于导入数据用于将图形窗口中的图形数据存储在二进制mat文件中，它们可以供其他的编程语言（如C语言等）调用。用于定义图形窗口的各种设置，包括字体、颜色等。用于打开【图形输出】话框，可以把图形以emf、ai、bmp、eps、jpg、pdf等格式保存，并设置有关图形窗口的显示等方面的参数。.1594.3.1 菜单操作S【Copy Figure】复制图形。S【Copy Option】打开【复制设置】对话框，设置图形复制的格式、图形背景颜色和
