北航数值分析B第二章课件ch2.3.pptVIP

§3 高斯选主元素消去法 例4 用高斯消去法解方程组 要求用具有舍入的10位浮点数进行计算。 精确到10位真解： 解法1（高斯消去法） 消元： 舍去或着说被“吃” 舍去或着说被“吃” 计算解： 解法2 用行变换的高斯消去法. 消元： 计算解: 该结果较好。该例子说明，在采用高斯消去法解方程组时，应 。对一般系数矩阵，最好保持乘数 ，因此在高斯消去法中引进选主元素技巧。 4.1 完全主元素消去法 第一步： （3）消元计算： 在A中选取绝对值最大的元素作为主元素，即确定 第k 步：重复进行，设已完成第1步—第k-1 的选主元，使 [A， ]约化为: 第k步的步骤： （3）消元计算： 二 回代求解： 算法8 是完全选主元素消去法。 工作量大。 经过上述过程，方程组约化为： 缺点： 优点： 改进方法： 列主元消去法， 设已完成第1步~第k-1步计算，得到与原方程组等价的方程组 方框内为第k步选主元素区域。 4.2?? 列主元素消去法 以下步骤类似完全选主元素消去法。 算法9 是列主元素消去法。 例5 用列主元素消去法解方程组 分析：由精确解看出有两位有效数字，因此，用4位浮点数 解 ：消元： 舍去或着说被“吃” 回代计算解： 高斯选主元消去法的步骤： 注：该解若取两位有效数字，则与真解完全相同。 优点： 数值稳定。 修正方法： 消元；回代。 列主元高斯-约当（Gauss-Jordam）消去法。 缺点： 既消元；又回代。 4.3 列主元高斯—约当（Gauss –Jordan)消去法 假设G--J消去法已完成第1步～第k-1步，得到与原方程组等价 第k步计算步骤： 的方程组 ，其中 （1）按列选主元： （2）换行： 消元： （3）消元计算： （4）计算主行（主元素所在行） 计算解： 说明： 因此，可以用来求逆矩阵。 如果用列主元G-J消去法将（A，I） 不用回代，将A化为单位矩阵，则解为常数项列。 定理9 （列主元高斯—约当法求逆矩阵） 化为（I，T）， 优点： 缺点： 因为计算量太大，但是在解多个方程组而它们的系数矩阵相同时， 该方法与高等代数中求逆矩阵方法的不同之处是有选主 注： 元，实际上选主元就是交换两行的位置，仍是初等变换，在一般的 求逆矩阵方法中也有交换两行元素。 例6 用列主元G-J消去法求 解： ? ?