第5节解线性方程组的直接方法.ppt

解线性方程组的直接方法/ Direct Method for Solving Linear Systems / 邹昌文 引言与预备知识 线性方程组数值解法分类 直接法 迭代法 向量与矩阵基本知识复习P162 直接法—高斯消去法 顺序高斯消去法/ Gaussian Elimination / 例:解下列线性方程组 方程组等价于 由最后一方程 现考虑n阶线性方程组 主元对解的影响 若主对角线上某元素为0则无法进行消元 若主对角线上某元素绝对值接近0,则用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，将严重影响计算结果的精度。 * 讨论线性方程组 消元过程 回代过程 ? 高斯消去法一般思路： 思路 首先将A化为上三角阵 / upper-triangular matrix /，再回代求解 / backward substitution /。 = 消元 Step 1：设 ，计算因子 将增广矩阵/ augmented matrix / 第 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 存在，则可通过逐次消元及行交换，将方程组化为三角形方程组，求出唯一解。 例：单精度解方程组 /* 精确解为 和 */ 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 / 省去换列的步骤，每次仅选一列中最大的元。 例： ? 注：列主元法没有全主元法稳定。 注意：这两个方程组在数学上严格等价。 ? ? 标度化列主元消去法 / Scaled Partial Pivoting / 对每一行计算　　　　　 。为省时间，si 只在初始时计算一次。以后每一步考虑子列 中　 最大的 aik 为主元。 例： 注：稳定性介于列主元法和全主元法之间。 *