数值计算方法ch2—2.1.pptVIP

第2章 解线性代数方程组的直接法 本章介绍求解n阶线性代数议程组 2.1 高斯消去法 2.1.1 高斯消去法的基本步骤 高斯（Gauss）消去法其实就是一种消去法（消元法）①，只是步骤规范，便于编写计算机程序。例如要解方程组 ①中国东汉初年诞生的《九章算术》，记载了求解线性代数组的消去法，比国外早1 500年。 2.1.2 高斯消去法的运算量 高斯消去法比起克莱姆法则和约当消去法，主要长处就是运算量较少。 用克莱姆法则求解n阶线性方程组，需计算n+1个行列式；每个行列式是n!个乘积之和，而每个乘积是n个元素相乘，因此共需乘法次数 ( n + 1）· n!( n － 1)= ( n + 1) ! ( n － 1) 当n=20时( n + 1) ! ( n － 1)≈9.7×1020。要完成这么多次乘法，在每秒做一亿次乘法运算的计算机上，也需30.8万年。因此克莱姆法则在实际计算中不适用。 高斯消去法消去过程分n-1步，第k步变换n－k行：对这些行先求倍数，再从n+1－k元素减去第k行对应列的倍数；因此共需乘、除次数 回代过程求xn需1次除法，求xn-1需1次乘法、1次除法，……，求需n –1次乘法1次除法，因此共需乘除次数 两过程共需乘除次数 。 当n=20时，N=3 060，显然比克莱姆法则少得多。 线性代数教材讲授线性方程组解法时，往往提倡用行初等变换把增广矩阵化为最简形。如上例方程组，多用如下变换： 2、计算机上为什么不用克莱姆法与约当消去法？ 3、何谓主元消去法？有何优点？ * * （2—1） 的直接法，并研究舍入误差对解的影响。这里假定系数行列式不为零，方程组的解存在且唯一。所谓直接法，就是不计舍入误差时，通过有限次算术运算能求得准确解的方法。克莱姆法就是一种直接法。只是方程组阶数较高时它的运算量太大，实际无法使用。本章介绍计算机上常用、有效的直接法——高斯消去法与三角分解法。 线性方程级（2-1）常简写成矩阵形式 Ax=b 其中 A= ，x = ， b = 其步骤如下：将第一个方程乘-5，-7，分别加于第二、第三方程，水去未知量x1，得同解方程组 将所得方程组的第二方程乘 ，加到第三方程，消去未知量x2，得同解方程组 这是上三角形方程组，它很容易求解；由第三方程得x3=-2，代入第二方程可得x2=5，再 代入第一方程得x1=-3。 一般高斯消去法包括两过程：先把方程组化为同解的上三角形方程组，再按相反顺序求解上三角形方程组。前者称消去或消元过程，后者称回代过程。 消去过程实际上是对增广矩阵 作行初等变换。如上例可表示为 对一般的n阶方程组，消去过程分n-1步：第一步消去a11下方元素。第二步消去a22下方元 素，……，第n-1步消去an-1，n-1下方元素。第k步将第k行的适当倍数加于其后名称，出可说是从k+1~n行减去第k行的适当倍数，使它们的第k列元素变为零，而列元素减去第k行对应列元素的倍数。因此，如把增广矩阵 变换前后都在计算机上用同一数组A存储，则消去过程可写为： 对k=1~n-1（步）做 对i= k+1~n（行）做 令 （“适当倍数”） 对j= k+1~n+1（列）令 回代过程是解同解的上三角形方程组 注意此处增广矩阵不是原始增广矩阵，而是变换后的增广矩阵，b是变换后增广矩阵的第n+1列。回代过程可写为： 令 对k = n – 1 ~ 1令 变换结束后时增广矩阵最后一列就是解，不需回代。这种解法称为约当消去法。它容易学习掌握，也容易编写计算机程序。其计算过程可写为： 对k=1~n做 对j= k+1~n+1令 对I=1~n但i≠k做（“适当倍数”） 对j= k+1~n+1令 不过这种解法需要乘除次数 比高斯消去法多