第五章线性代数方程组的直接解法2.pptVIP

* §2 Gauss消元法/*Gaussian Elimination Method*/ 一、 基本Gauss消元法 思想 通过初等变换逐步消去未知元，将 原方程组化为同解的三角方程组。 三角（系数矩阵）方程组易于求解 方程组 的增广矩阵记为： 其中 令 Step 1：如果 初等变换 下面用初等矩阵来描述上述消元过程 取 其中 记 同样记 记 上述两个计算公式可以写成统一形式： 类似地，对 中的 部分重复以上做法 其中 Step k：第k步消元过程的计算公式 整个消元过程的矩阵表示 上三角矩阵 计算 经过n-1次消元，原方程组化为等价的上三角方程组 回代过程的计算公式： 等价方程组 得到的解向量 存放在增广矩 阵的最后一列 for for for Gauss法的消元过程 回代过程 例5：用基本Gauss消元法求解下列方程组 解： 增广矩阵 ?基本Gauss消元法的工作量 消元过程： 回代过程： 加减法的次数 乘除法的次数 （基本Gauss消元法的实现条件） 全不为零的充要条件是 的顺序主子式都不等于零，即 证明： 归纳法证明(对k归纳) 设直到k-1成立,只要证明 非零时， 非零的充要条件是 即可。 在归纳假设下，基本Gauss消元法可进行到第k-1步 其中 是对角元为 的上三角矩阵 矩阵 的k阶主子式 为上三角的 均为单位下三角矩阵 其中 因此，若方程组的系数矩阵的顺序主子式均不为零， 可以采用基本Gauss消元法求解。 结论得证 小主元 可能导致计算失败 例6：在8位制计算机上解方程组 要求用Gaussian Elimination计算。 8个 解： 二、 主元素Gauss消元法/* Pivoting Strategies */ 思想 每次消元之前，在剩余元素中选择绝对值最大的 非零元素作为主元， 然后经过换行换到主元位置 ?列主元消去法/* Column Pivoting Strategies */ Step k：第k步首先选择主元 寻求 满足 然后交换矩阵 的第 行和 行，再进行消元过程 算法: Gauss列主元消去算法 求方程组Ax=b 的解. 输入:增广矩阵An?(n+1)=（A|b）. 输出: 近似解 xk=ak,n+1(k=1,2,…,n) 或失败信息. 消元过程 for k = 1,2,…,n-1 do Step 1 - Step 4 Step 1 寻找行号 ik , 使得 Step 2 如果 ，则交换第k行和ik行； 否则转Step 7 算法: Gauss列主元消去算法（续） Step 3 for i=k+1,…,n 计算 Step 4 for j=k+1,…,n+1 计算 回代过程 Step 5 Step 6 for i=n-1,…,1 计算 Step 7 Output (系数矩阵奇异); /*不成功 */ STOP.