5.1高斯消去法教案.ppt

线性方程组的概念 线性方程组的概念（续） 线性方程组的数值解法 例1（续） 二、Gauss消元法的计算量 Gauss法与Cramer法则的计算量比较 四、主元素法 例2（续1） 例3（续2） 例3两种解法的误差分析 4.2 列主元素法 列主元素法 4.3 全主元素法 主元素法举例 由于在计算机中作乘除运算量所需时间远大于作加减运算所需时间，故只考虑作乘除运算量。 由消元法步骤知，第k次消元需作n?k次除法，作 (n ? k)(n ? k + 1)次乘法，故消元过程中乘除法运算量为： 所以Gauss 消去法的 乘除法总运算量为： Gauss 消元法的乘 除法总运算量为： 与我们曾经介绍的Cramer法则的乘除法总运算量 (n2?1)n!+n 相比，由下表可知：当阶数越高时，Gauss消 元法所需乘除法次数比Cramer法则要少得多： 方程组阶数 3 10 20 50 Gauss消元法运算量 17 430 3060 44150 Cramer法则运算量 51 359251210 9.7×1020 7.6×1067 ? Gauss 消元法的优缺点： 但其计算过程中，要求akk(k)（称为主元素）均不为零，因而适用范围小，只适用于从1到n ?1阶顺序主子式均不为零的矩阵A，计算实践还表明，Gauss消元法的数值稳定性差，当出现小主元素时，会严重影响计算结果的精度，甚至导出错误的结果。 Gauss消元法简单易行。 三、高斯消去法实现的条件 4.1 引入主元素的必要性 对线性方程组AX = b，若其系数行列式 det(A) ? 0，则该方程组有唯一 解，但是这一条件 不能保证所有主元素都不等于零，只要某一主元 素等于零，就不能用Gauss消元法求解该方程组， 即使所有主元素不等于零，但 某一主元素的绝对 值很小时，Gauss消元法也是不适用的。如下例： 例3 解：为减小误差，计算过程中保留3位有效数字。 按Gauss消元法步骤： 第一次消元后得同解方程组： 第二次消元后得同解方程组 回代得解，x3 = 2.02，x2 = 2.40，x1 = ?5.80。 容易验证，方程组（1）的准确解为： x1 = ?2.60， x2 = 1.00，x3 = 2.0。 显然两种结果相差很大。 若在解方程组前，先交换方程的次序，如将（1）交换一行与二行改写成如下所示： 再用Gauss消元法，顺序消元后得同解方程组： 回代得解：x3 = 2.00，x2 = 1.00，x1 = ?2.60 ——与准确解相同。 在例3中，对(1)的方程进行顺序消元时， 主元a(1)11=0.50，a(2)22=0.100都比较小，以它们作 除数就增长了舍入误差,而导致计算结果不准确。 产生上述现象的原因在于舍入误差，当|a(k)kk| 很小时，进行第k次消元，要用|a(k)kk|作除数，这 样就可能增大舍入误差造成溢出停机，或者导致 错误的结果。 为了在计算过程中，抑制舍入误差的增长， 应尽量避免小主元的出现。如例3的第二种解法, 通过交换方程次序，选取绝对值大的元素作主元 基于这种思想而导出主元素法。 为简便起见，对方 程组（5-1），用 其增 广矩阵： 表示，并直接在增广矩阵上进行运算。 列主元素法的具体步骤如下： 如此经过n?1步，增广矩阵（5-3）被化成上三角形，然后由回代过程求解。 在上述过程中，主元是按列选取的，故称为列主元法。例3中的第二种解法就是按列主元法进行的。 经过n?k次消元后，得到与方程组（5-1）同解的上三角形方程组，再由回代过程求解。 例4 计算过程保留三位小数。 此例的计算结果表明，全主元素法的精度略优于列主元法，这是由于全主元法是在全体元素中选主元，故它对控制舍入误差十分有效。但是全主元法在计算过程中，需同时作行与列的互换，因而程序比较复杂,计算时间较长。列主元法的精度虽稍低于全主元法，但其计算简单，工作量大为减少，且计算经验与理论分析均表明，它与全主元法同样具有良好的数值稳定性，故列主元法是求解中小型稠密线性方程组的最好方法之一。 方法之一。 1. 用高斯消元法解方程