第3章 解线性方程组及直接方法.ppt

数值计算方法 第三章 解线性方程组的直接方法 上述消元过程除第一个方程不变以外, 第2—第 n 个方程全消去了变量 ?1，而系数 和常数项全得到新值： 回代过程算法 高斯主元素消去法 Gauss消去法的矩阵表示 LU形式 直接计算 A 的 LU 分解(例) 一般计算公式 LU 分解求解线性方程组 §5.向量和矩阵的范数 矩阵A的谱半径 5.6 误差分析 二、误差分析 二、矩阵的范数 矩阵范数例 与前述三种向量范数相容的三种矩阵范数： 求矩阵A的各种常用范数 解: 由于 特征方程为 容易计算 计算较复杂 对矩阵元素的 变化比较敏感 使用最广泛 性质较好 Re Im ? ? ? ? ? ? ? ? ? (A) 定理1. 与假设矛盾 i+1行 i+1行 定理 若A的所有顺序主子式 /* determinant of leading principal submatrices */ 均不为0，则 A 的 LU 分解唯一（其中 L 为单位下三角阵）。 证明：由§1中定理可知，LU 分解存在。下面证明唯一性。 若不唯一，则可设 A = L1U1 = L2U2 ，推出 Upper-triangular Lower-triangular With diagonal entries 1 ? 注： L 为一般下三角阵而 U 为单位上三角阵的分解称为Crout 分解。 实际上只要考虑 A* 的 LU 分解，即 ，则 即是 A 的 Crout 分解。 计算量与 Gauss 消去法同. ? 平方根法 /* Choleski’s Method */： ——对称 /* symmetric */ 正定 /* positive definite */ 矩阵的分解法 定义 一个矩阵 A = ( aij )n?n 称为对称阵，如果 aij = aji 。 定义 一个矩阵 A 称为正定阵，如果 对任意非零向量 都成立。 因此 Diagonal:对角 为非奇异下三角阵 为非奇异上三角阵 ----------(2) --------(3) 因此 所以 综合以上分析, 则有 -------------(4) -------------(5) 定理11. (Cholesky分解) 且该分解式唯一 这种关于对称正定矩阵的分解称为Cholesky分解 -------------(6) -------------(7) -------------(8) 对称正定线性方程组的解法 线性方程组 -------------(10) -------------(11) 则线性方程组(10)可化为两个三角形方程组 -------------(12) -------------(13) ------(14) ------(15) 对称正定方程 组的平方根法 例1. 用平方根法解对称正定方程组 解: 即 所以原方程组的解为 思考 本例中出现了大量的根式运算 原因为 考虑改变分解方式 请求解例1. 平方根法的数值稳定性 用平方根法求解对称正定方程组时不需选取主元 由 可知 因此 平方根法是数值稳定的 事实上,对称正定方程组也可以用顺序Gauss消去法求解 而不必加入选主元步骤 ? 追赶法解三对角方程组 /* Crout Reduction for Tridiagonal Linear System */ Step 1: 对 A 作Crout 分解 直接比较等式两边的元素，可得到计算公式（p.194）。 Step 2: 追——即解 ： Step 3: 赶——即解 ： 与G.E.类似，一旦 ?i = 0 则算法中断，故并非任何 三对角阵都可以用此方法分解。 定理 若 A 为对角占优 /* diagonally dominant */ 的三对角阵，且满足 ，则追赶法可解以 A 为系数矩阵的方程组。 Hey, what does diagonally dominant mean??? It means that the diagonal entries of the matrix are very LARGE. Well, how large is LARGE? T