《科学计算引论》实验课：第四章 线性方程组数值解法.pptVIP

算法 4.2 追赶法 * * 《科学计算引论》实验 第四章 线性方程组数值解法 本章内容简介 § 4.1 Gauss 消元法 § 4.2 特殊线性方程组的解法及敏度分析 § 4.3 经典迭代方法 § 4.4 Krylov 子空间方法 A*X=b 矩阵除法 矩阵的除法：/、\ 右除和左除 若 A 为可逆方阵，则 左除：A\B == A 的逆左乘 B == inv(A)*B 右除：B/A == A 的逆右乘 B == B*inv(A) 右除：X=A\B == A*X=B 左除：X=B/A == X*A=B 通常，矩阵除法可以理解为 当 A 和 B 行数相等时即可进行左除 当 A 和 B 列数相等时即可进行右除 doc inv；%(See \；see Algorithm for details) Matlab中， X=A\B ≠ X=inv(A)*B， X=A\B更好 LU分解 利用矩阵乘法直接计算 LU 分解 L ? U = A 其中：L --- 单位下三角矩阵，U --- 上三角矩阵 LU分解Matlab命令 LU分解(列主元Gauss消元)： [L,U]=lu(A) %产生一个上三角阵U和一个变换形式的下三角阵L(行交换)，矩阵A必 须是方阵。 [L,U,P]=lu(X)： 产生一个上三角阵U和一个下三角