第四章-第一节-Gauss消元法.ppt

* * 第一节 消元法 四、 消元法的矩阵形式 三、 主元消元法 二、 消元法的基本思想 一、引言 五、小结 在自然科学和工程技术问题中，涉及到许多数值计算问题，最终都要归结为解线性代数方程组 。其中 是可逆的。本章和下一章分别讨论解方程组的直接方法和迭代方法。所谓直接方法就是通过有限次的精确运算能得到真解的一类数值方法。从本质上讲，直接方法的原理是找到一个可逆矩阵 ，使得 是一个上三角阵，这个过程称为“消元”过程。消元之后再进行“回代”，即求解 。实际计算过程中，不必明显地计算出矩阵 ,而只须把 和 计算出来。这类直接方法中最基本和最简单的就是 消元法，本章首先讨论 消元法和矩阵分解法，以及 消元法在各种情况下的变形，并分析其误差。 一、引言 用矩阵和向量的记号表示，则有 其中 为可逆矩阵， 消元法分消元和迭代两个过程 ，消元过程是将(1.1)化成为如下形式的上三角方程组： 二、 消元法的基本思想 考虑 阶线性代数方程组： 迭代过程是从(1.3)最后一个方程直接解出 然后依公式 依次求出 ，称为回代求解。消元过程的实质是对增广矩阵 作一系列初等行变换，最后把 化为上三角矩阵 ，得 ，因为对 每做一次初等行变换，相当于对方程式组(1.1)进行一次同解变换，所以与 相应的上三角形方程组（1.3）是（1.1）的同解方程组。 第二步，假设 令 于是 用上述方法又可把(1.5)化为 设方程组（1.1）的增广矩阵 不妨设 并令 第一步消元是用 主对角元以下的元素全化为0，得 乘第一行然后加到第 行上去，从而把第一列 (1.5) (1.6) 第三步，假设 再用 消去 如此继续，共 步即可把方程组(1.1)化为形如（1.3）的上三角方程组： (1.7) 其中 和 分别为方程组(1.3)的系数矩阵和右端向量。 这样就完成了消元过程，最后利用公式(1.4)“回代”求解。 由以上分析可以看出，消元过程的第 步共含除法运算 次，乘法运算 次，所以消元过程共含乘除法次数为 而回代过程的乘除法运算次数为 所以 消元法总的乘除法次数为 如果我们用 法则计算(1.1)的解， 要计算 个 阶行列式并作 次除法。 而计算每个行列式， 若用子式展开的方法，则有 次乘法， 所以用 法则大约需要 次乘除法运算。例如，当 时约需 次运算，而用 消元法只需 次乘法运算。 例1 方程组 解： 方程 方程(2)得： 代入方程(1)得 由此得到的解完全失真， 如果交换两个方程的顺序，得到等价方程组 的精确解为： 在 消元法计算中取 位有效 数字。 经 消元后有 由此可以看到，在有些情况下，调换方程组的次序对方程组的解是有影响的，在消元法中抑制舍入误差的增长是十分重要的。 三、 主元消元法 选主元的一种简单办法是，第 步消元时，在 的第 列元素 中选取绝对值最大者作为主元，并将其对换到 位置上，然后再进行消元计算，这样选取的主元叫列主元。