第2章1-3节消去法-分解-条件数.ppt

例 用Doolittle分 解法求解 解: 所以 §2.3 线性方程组的误差分析 考虑方程组Ax=b的系数矩阵或右端项有误差时,误差对解将产生怎样的影响. 例2: 用直接法求解下列方程组: 显然此方程组的解为 定义7:设A是非奇异矩阵,若Cond(A)远大于1,则称方程组AX=b是病态的,否则称方程组是良态的. 对于先前的例子 §2.3 Error Analysis for . ? 近似解的误差估计及改善： 设　　 的近似解为　 ，则一般有 cond (A) 误差上限 ? 改善方法： Step 1:　　　　　近似解 Step 2: Step 3: Step 4: 若　 可被精确解出，则有 就是精确解了。 经验表明：若 A 不是非常病态（例如： ），则如此迭代可达到机器精度；但若 A 病态，则此算法也不能改进。 AX = b （2.1） 第二章 线性方程组解法 例1: 用顺序消去法解方程组 解: 解此上三角方程组，得 称为上三角方程组 §1 高斯消去法 1．三角形方程组的解法 2．高斯消去法 消元 记 Step 1：设 ，计算因子 将增广矩阵的第 i 行 ? mi1 ? 第1行，得到 其中 Step k：设 ，计算因子 且计算 共进行 ？ 步 n ? 1 回代 What if ? No unique solution exists. What if ? Then we must find the smallest integer k ? i with , and interchange the k-th row with the i-th row. What if we can’t find such k ? No unique solution exists. 定理 若A的所有顺序主子式 /* determinant of leading principal submatrices */ 均不为0，则高斯消元无需换行即可进行到底，得到唯一解。 注：事实上，只要 A 非奇异，即 A?1 存在，则可通过逐次消元及行交换，将方程组化为三角形方程组，求出唯一解。 ? 选主元消去法 /* Pivoting Strategies */ 例：单精度解方程组 /* 精确解为 和 */ 8个 8个 用Gaussian Elimination计算： 8个 小主元 /* Small pivot element */ 可能导致计算失败。 ? 全主元消去法 /* Complete Pivoting */ 每一步选绝对值最大的元素为主元素，保证 。 Step k: ① 选取 ② If ik ? k then 交换第 k 行与第 ik 行; If jk ? k then 交换第 k 列与第 jk 列; ③ 消元 注：列交换改变了 xi 的顺序，须记录交换次序，解完后再换回来。 ? 列主元消去法 /* Partial Pivoting, or maximal column pivoting */ 省去换列的步骤，每次仅选一列中最大的元。 例： ? 注：列主元法没有全主元法稳定。 例： 注意：这两个方程组在数学上严格等价。 ? §2 三角分解法 /* Matrix Factorization */ ? 高斯消元法的矩阵形式 /* Matrix Form of G.E. */： Step 1: 记 L1 = ，则 Step n ? 1: 其中 Lk = §2 Matrix Factorization – Matrix Form of G.E. 记为 L 单位下三角阵 /* unitary lower-triangular matrix */ 记 U = A 的 LU 分解 /* LU factorization */ 为神马要考虑A的矩阵分解形式??! When you have to solve the system for different with a fixed A. Could you be more specific, please? Factorize A first, then for every