第4章MATLAB的数学运算.pptVIP

第4章 MATLAB 的数学运算 教学内容 矩阵和线性代数 稀疏型矩阵 多项式与插值 函数运算 微分方程 矩阵的基本运算 矩阵与常数之间的四则运算 矩阵与常数的运算与数组运算相同 矩阵和矩阵之间的四则运算 矩阵和矩阵之间的加减运算与数组运算相同 设 A 是一个 m×n 矩阵，B 是一个 p×q 矩阵，当 n＝p 时，两个矩阵可以相乘，乘积为 m×q 矩阵。矩阵乘法不可逆。在 MATLAB 中，矩阵乘法由“*”实现。 矩阵除法在实际中主要用于求解线性方程组 矩阵转置 符号“’”实现矩阵的转置操作。对于实数矩阵，“’”表示矩阵转置，对于复数矩阵，“’”实现共轭转置。对于复数矩阵，如果想要实现非共轭转置，可以使用符号“.’”。 矩阵乘方 特殊矩阵生成 特殊矩阵生成 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. 矩阵和线性代数 线性方程组 逆矩阵和行列式 矩阵分解 矩阵指数函数和幂函数 矩阵特征值 矩阵奇异值分解 线性方程组求解问题 工程上经常要考虑的问题是线性方程组求解。对于给定的矩阵A和B是否存在使得AX＝B或者XA＝B的X。 利用矩阵求逆的方法求解方程组 我们首先求系数矩阵的逆，然后利用矩阵的逆求解方程组的解。 利用矩阵的左除符号“\”或者右除符号“/”求解方程组 线性方程组求解 利用左除符号和右除符号求解线性方程组，避免了矩阵求逆操作，因此系数矩阵不必为方阵。如果系数矩阵的维数为m×n时 ，则有三种情况： ，此时方程组为超定方程组，MATLAB 将给出最小二乘解； ，此时方程组为方阵系统，MATLAB 给出精确解； ，此时方程组为欠约束方程组，MATLAB 将给出一组基解，该解中包含最多 个非零元素。 在采用除法符号（包括左除和右除）求解线性方程组时，MATLAB 采用 因式分解法求解方程组。尽管 MATLAB 提供了两种方法，一般更倾向于采用第二种方法，该方法用到较少的浮点数运算，执行速度较快，另外，由于采用 分解法，得出的结果要精确的多。 方程求解实例 对于这样的方程组 利用矩阵求逆inv(A)*B 利用左除和右除AX=B的解写做:x=A\B XA=B的解写做:x=A/B 逆矩阵和行列式 矩阵行列式 在 MATLAB 中，矩阵的行列式用函数 det 求解。调用格式为 det(A)，其中 A 为方阵。 矩阵求逆 对于非奇异方阵，如果存在方阵X，满足AX=E， 则称X为矩阵A的逆，记为A-1，在 MATLAB 中，用 inv(A) 来实现矩阵逆的求解。 也可以写成A^-1 矩阵分解 Cholesky 分解 LU 分解 QR 分解（正交分解） Cholesky 分解 Cholesky 分解将对称矩阵表示为一个三角矩阵与其转置的乘积的形式，即，其中为对称矩阵，为上三角矩阵。并非所有的对称矩阵都能进行 Cholesky 分解，只有正定矩阵能够进行 Cholesky 分解，如 Pascal 矩阵。在 MATLAB 中 Cholesky 分解由函数 chol 实现，该函数对输入矩阵进行 Cholesky分解，返回其对应的三角矩阵。 Cholesky 分解同样适用于复数矩阵。如果复数矩阵满足，其中表示矩阵的共轭转置。如果矩阵存在 Cholesky 分解则称其为 Hermitian 正定。 LU 分解 矩阵的 LU 分解将一个方阵表示为一个下三角置换矩阵和一个上三角矩阵乘积的形式。如，其中为下三角置换矩阵，为上三角矩阵。MATLAB 中 LU 分解可以通过函数 lu 实现。通过矩阵的 LU 分解，可以实现线性方程组的快速求解。 另外矩阵的 LU 分解可用于矩阵快速求逆和求行列式，有 det(A) = det(L)*det(U) 和 inv(A) = inv(U)*inv(L)。 QR 分