ch04.MATLAB及数学运算.ppt

第4章 MATLAB 的数学运算 教学目标 掌握向量和矩阵的运算 掌握线性代数的基本函数和使用 掌握稀疏矩阵的操作 教学内容 向量、矩阵及其运算 矩阵和线性代数 稀疏型矩阵 向量、矩阵及其运算 向量的点乘、叉乘和混合积 矩阵的基本运算 特殊矩阵生成 向量和矩阵的范数 矩阵的条件数 向量的点乘、叉乘和混合积 向量的点乘 向量的点乘又称为内积，是两个向量的模和两个向量之间的夹角余弦三者的乘积。MATLAB 中，实现点乘的函数是dot。dot 函数的用法为 dot(x1,x2)，其中 x1 和 x2 的维数必须相同。 向量的叉乘 向量乘法除点乘之外还有叉乘。两个向量叉积的几何意义是指以两个向量模的乘积为模，方向和两个向量构成右手坐标系的向量。向量的叉乘不可交换。在 MATLAB 中函数 cross 用于实现向量的叉乘。 向量的混合积 向量的混合积的几何意义是：它的绝对值表示以三个向量为楞的平行六面体的体积，符号由右手法则确定。上面介绍了向量的点乘和叉乘，向量的混合积由点乘和叉乘逐步实现。 矩阵的基本运算 矩阵与常数之间的四则运算 矩阵与常数的运算与数组运算相同 矩阵和矩阵之间的四则运算 矩阵和矩阵之间的加减运算与数组运算相同 设 A 是一个 m×n 矩阵，B 是一个 p×q 矩阵，当 n＝p 时，两个矩阵可以相乘，乘积为 m×q 矩阵。矩阵乘法不可逆。在 MATLAB 中，矩阵乘法由“*”实现。 矩阵除法在实际中主要用于求解线性方程组 矩阵转置 符号“’”实现矩阵的转置操作。对于实数矩阵，“’”表示矩阵转置，对于复数矩阵，“’”实现共轭转置。对于复数矩阵，如果想要实现非共轭转置，可以使用符号“.’”。 矩阵乘方 特殊矩阵生成 （1/2） 特殊矩阵生成 （2/2） 1．对角矩阵的生成 对角矩阵指除对角线以外其他元素为 0 的矩阵。函数 diag 可以生成对角矩阵。该函数的用法为： A=diag(V,K)，其中 V 是一个向量，K 是一个整数。该函数返回一个矩阵，矩阵的第 K 个对角线为 V。K 在默认情况下为 0，表示矩阵的主对角线，K 大于 0 时表示主对角线的上方，小于 0 时为主对角线的下方。 V=diag(A,K)，其中A时一个矩阵。K 与上面的语句相同。该语句返回矩阵 A 第 K 个对角线上的元素组成的矩阵。 2．魔术矩阵的生成 魔术矩阵是一种经常遇到的矩阵，除了二阶方阵之外，魔术矩阵的每一行、每一列以及每条主对角线的元素之和都相同。在 MATLAB 中，magic 函数用于生成魔术矩阵。其调用方法为 magic(N)，其中 N 为正整数，并且 N≠2. 向量和矩阵的范数 向量的范数定义为 ，其中，最常用的值为 1、2 和无穷大。矩阵 的范数定义为 ，其中，最常用的值为 1、2 和无穷大。 向量和矩阵的范数可以通过函数 norm 求解。该函数的调用格式为 n = norm(A,p)，其中 p 用于指定范数的类型。p 可以为所有大于 1 的常数，最常用的为 1、2、inf 和 fro，fro 为求解矩阵 A 的 Frobenius 范数。当 p 省略时，默认值为 2。 矩阵的条件数 矩阵的条件数是在矩阵的逆和矩阵范数的基础上定义的，是用于衡量矩阵病态程度的关键量。一个矩阵的条件数越大，表明该矩阵的病态程度越严重。 MATLAB中求解矩阵条件数的函数是cond，在命令窗口输入cond(A)即可得到矩阵A的条件数。 矩阵和线性代数 线性方程组 矩阵逆、秩及行列式 矩阵分解 矩阵指数函数和幂函数 矩阵特征值、特征向量与特征行列式 矩阵的标准正交基 矩阵奇异值分解 线性方程组（1/2） 利用矩阵求逆的方法求解方程组 我们首先求系数矩阵的逆，然后利用矩阵的逆求解方程组的解。 利用矩阵的左除符号“\”或者右除符号“/”求解方程组 线性方程组（2/2） 利用左除符号和右除符号求解线性方程组，避免了矩阵求逆操作，因此系数矩阵 不必为方阵。如果系数矩阵 的维数为 ，则有三种情况： ，此时方程组为超定方程组，MATLAB 将给出最小二乘解； ，此时方程组为方阵系统，MATLAB 给出精确解； ，此时方程组为欠约束方程组，MATLAB 将给出一组基解，该解中包含最多 个非零元素。 在采用除法符号（包括左除和右除）求解线性方程组时，MATLAB 采用 因式分解法求解方程组。尽管 MATLAB 提供了两种方法，一般更倾向于采用第二种方法，该方法用到较少的浮点数运算，执行速度较快，另外，由于采用 分解法，得出的结果要精确的多。 矩阵逆、秩及行列式 矩阵求逆 对于非奇异方阵，如果存在方阵，满足 并且 ，则称为矩阵的逆