2.1高斯消元法课程教学课件.pptxVIP

第二章 解线性方程组的直接解法;本章介绍求解n元线性代数方程组 的??接解法，并研究舍入误差对解的影响。 ;; 线性方程级（2-1）常简写成矩阵形式 Ax=b 其中 A= ，x = ， b = ;§ 2.1 高斯消去法 2.1.1 高斯消去法的基本步骤 高斯（Gauss）消去法其实将一般的线性方程组变换为三角形（上三角）方程组求解问题（消元法）①，只是步骤规范，便于编写计算机程序。 ①中国东汉初年诞生的《九章算术》，记载了求解线性代数组的消去法，比国外早1500年。; 例 解方程组 解：步骤如下：① 将第一个方程乘-5，-7，分别加于第二、第三方程，消去未知量x1，得同解方程组;② 将所得方程组的第二方程乘 ，加到第三方程，消去未知量x2，得同解方程组; 一般高斯消去法包括两过程：先把方程组化为同解的上三角形方程组，再按相反顺序求解上三角形方程组。前者称消去或消元过程，后者称回代过程。 消去过程实际上是对增广矩阵 作行初等变换。如上例可表示为; 对一般的n阶方程组，消去过程分n-1步： 第一步消去a11下方元素。第二步消去a22下方元;素，……，第n-1步消去an-1，n-1下方元素。即第k步将第k行的适当倍数加于其后各行，或可说是从k+1~n行减去第k行的适当倍数，使它们的第k列元素变为零，而其余列元素减去第k行对应列元素的倍数。; 对k=1~n-1（步）做 对i= k+1~n（行）做 令 （“适当倍数”） 对j= k+1~n+1（列）令;回代过程是解同解的上三角形方程组;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万年。因此克莱姆法则在实际计算中不适用。 ;高斯消去法消去过程运算量：; 回代过程运算量： 回代过程求 xn需1次除法，求xn-1需1次乘法、1次除法，……，求需n –1次乘法1次除法，因此共需乘除次数： 两过程共需乘除次数 。; 当n=20时，N=3 060，显然比克莱姆法则少得多。 线性代数教材讲授线性方程组解法时，往往提倡用行初等变换把增广矩阵化为最简形。如上例方程组，多用如下变换：; 变换结束后时增广矩阵最后一列就是解，不需回代。这种解法称为约当消去法。它容易学习掌握，也容易编写计算机程序。其计算过程可写为：; 对k=1~n做 对j= k+1~n+1令 对I=1~n但i≠k做（“适当倍数”） 对j= k+1~n+1令; 因此除非方程组的阶数n较小，计算机上不用约当消去法。 上面比较了三种解法的乘除次数。一般算法乘除次数与加减次数相差不多，而且计算机上乘除运算比加减运算费时间，所以通常只比较乘除次数。 ; 2.1.3 选主元消去法 高斯消去法的弊端： 高斯消去法消去过程中，第k步中求n-k个倍数数 用到的除数 ，称为主元。它若为零或接近于零，计算机将“溢出”而停止计算，或产生较大误差。 ;; 造成这种结果的原因，就是小主元的出现。用它用除数产生大乘数，数乘大数变大数，大数吃小数产生舍入误差（如3－4×105=-399 997→ -106×0.400 0），从而有误差。 求出x2代回第一个方程时，因 10-5x1+2x2=1， 10-5x1’+2x2’=1 两式相减得10-5(x1－ x1’)+2(x2－ x2’)=0，可见 | x1－ x1’ |=200 000| x2－ x2’| 表明 x1的误差被放大200 000倍， x1’自然失真。; 列主元