第6章==高斯消元法.ppt

例1,解方程组 2x1+x2+x3=7 4x1+5x2-x3=11 x1-2x2+x3=0 2x1+x2+x3=7 3x2-3x3=-3 2x3=-6 * 第6章 线性方程组的数值解法 6.1 问题背景——电阻网络 已知某电阻网络如图6.1－1所示： 图6.1－1电阻网络 由图6.1－1可建立以下电流方程组： （6.1－1） 由此可知，解上述方程组可以求出电流流经各电阻的电流 或简写为 （6.1－3） 形如（6.1－2）式的n阶线性代数方程组中 ， 为未知系数； 为常数。 或者式子(6.1-2)写成矩阵形式 （6.1－4） 对于此类型问题的线性方程组，其一般形式为： （6.1－2） 式中， 是由方程组的系数组成n阶矩阵，称为系数矩阵； 是由右端各系数组成的n维列向量。 由线性代数知识知，当且仅当行列式 即矩阵A非奇异时， 根据克莱姆法则，线性代数方程组（6.1－4）存在唯一解，并可表示为两个行列式之比： 其中， Di是将行列式D的第i列用右端向量b来代替所构成的行列 式. 克莱姆法则仅适用于低阶方程组，而对于高阶方程组却不适用。 （6.1－4） 求解线性方程组的数值方法大致可分为两大类：一类是直接法（也称精确解法）；另一类是迭代法。 直接法 直接法主要有高斯（Gauss）消去法和矩阵分解法等。 消去法的基本思想是：通过将一个方程乘或除以某个常数，以及将两个方程相加、减这两种手续，逐步减少方程中变元的数目，最终使每个方程仅含一个变元，从而得出所求的解。 2.迭代法 迭代法就是通过某种极限过程去逐次逼近方程组精确解的数值计算方法，特别适合于求解系数矩阵为大型稀疏矩阵的线性方程组。 6.2 高斯（Gauss）消元法 1.高斯消去法的计算过程 用高斯消去法求下列方程组 （1）消去过程 增广矩阵 （2）回代过程 从上三角形方程组逐次回代计算出线性方程组的解。 以上求解线性方程组的方法称为Gauss消去法 Gauss消去法步骤: 消元、回代 同解 即 根据Cramer(克莱姆)法则,若 解：第一步，消元。此线性方程组的增广矩阵为 对增广矩阵进行初等行变换，使系数矩阵化为上三角形。 1 1 7 4 5 -1 11 1 -2 1 0 1 1 7 4 5 -1 11 1 -2 1 0 (2)-(1)*4/2 (3)-(1)*1/2 1 1 7 0 3 -3 -3 0 -2.5 0.5 -3.5 (3)-(2)*2.5/3 1 1 7 0 3 -3 -3 0 0 -2 -6 于是得同解方程组 第2步，回代求解，得： x3=3, x2=2,x1=1 2.高斯消去法应注意的问题 （1） 全不为0 在高斯消去过程中，第k步要用 做除法运算，如果 计算将无法继续下去， 因此一定要保证 全不为0。 （2） 的绝对值不能太小 即使 不为0，但是如果其绝对值很小，舍入误差的影响也会严重地损失精度，因此在实际计算时，必须避免此类情况发生。 例6.1 考察方程组 （6.2－7） 用高斯消去法求解时，先用10-5除（6.2－7）中第一式，然后从（6.2－7）中第二式中消去x1，得 （6.2－8） 计算机按四舍五入计算，有 因而方程组（6.2－8）的实际形式是 由此回代得出 这个结果是严重错误的。 原因：a11=10-5太小，使得方程组（6.2－8）中第一式在消元过程中“吃掉”了方程组（6.2－8）中第一式。 避免此类错误的一种有效办法是，在消元前先调整方程的秩序（行交换），即将方程组（6.2－7）改写为 再进行消元则得 因而上述方程组的实际形式是 由此回代得到x1=x2=1，这个结果才是正确的。 这种行交换的方法叫列主元素高斯消元法。 (1-?)和(1-2?)都被算为1，上述方程组变为： 得到 x1=x2=1，这个结果才是正确的。 采用列主元素高斯消元法：（1）行交换；（2）消元。 (1)*(-?)+(2) 这表明：列主元素高斯消元法应用广。但也有缺点，如下列 列主元素高斯消元法或高斯消元法，上述方程组变为：