高斯消元法.ppt

高斯消元法.ppt

《高斯消元法.ppt》由会员分享，可在线阅读，更多相关《高斯消元法.ppt（29页珍藏版）》请在冰豆网上搜索。

第二章第二章解线性方程组的直接解法解线性方程组的直接解法高斯消元法解三对解方程组的追赶法方程组的病态问题三角分解法本章介本章介绍绍求解求解nn元元线线性代数方程性代数方程组组的直接解法，并研究舍入的直接解法，并研究舍入误误差差对对解的影响。

解的影响。

（21）这里假定系数行列式不为零，方程组的解存解存在且唯一。

在且唯一。

所谓直接法直接法，就是不计舍入误差时，通过有限次算术运算能求得准确解的方法。

克莱姆法克莱姆法就是一种直接法。

只是方程组阶数较高时它的运算量太大，实际无法使用。

线性方程级（2-1）常简写成矩阵形式Ax=b其中A=，x=，b=2.1高斯消去法高斯消去法2.1.12.1.1高斯消去法的基本步骤高斯消去法的基本步骤高斯（高斯（GaussGauss）消去法其实将一般的线消去法其实将一般的线性方程组变换为性方程组变换为三角形（上三角）方程组三角形（上三角）方程组三角形（上三角）方程组三角形（上三角）方程组求求解问题（消元法）解问题（消元法），只是步骤规范，便于，只是步骤规范，便于编写计算机程序。

编写计算机程序。

中国东汉初年诞生的中国东汉初年诞生的九章算术九章算术，记载了求解，记载了求解线性代数组的消去法，比国外早线性代数组的消去法，比国外早15001500年。

年。

例解方程组解：

步骤如下：

将第一个方程乘-5，-7，分别加于第二、第三方程，消去未知量x1，得同解方程组将所得方程组的第二方程乘，加到第三方程，消去未知量x2，得同解方程组这就是上三角方程组，它很容易求解。

由第三方程得x3=-2，代入第二方程可得x2=5，再代入第一方程得x1=-3。

一般高斯消去法包括两过程：

先把方程组化为同解的上三角形方程组，再按相反顺序求解上三角形方程组。

前者称消去消去或消元消元过程，后者称回代回代过程。

消去过程实际上是对增广矩阵作行初等变换。

如上例可表示为对一般的n阶方程组，消去过程分n-1步：

第一步消去a11下方元素。

第二步消去a22下方元素，第n-1步消去an-1，n-1下方元素。

即第k步将第k行的适当倍数加于其后各行，或可说是从k+1n行减去第k行的适当倍数，使它们的第k列元素变为零，而其余列元素减去第k行对应列元素的倍数。

对k=1n-1（步）做对i=k+1n（行）做令（“适当倍数”）对j=k+1n+1（列）令因此，如把增广矩阵变换前后都在计算机上用同一数组A存储，则消去过程可写为：

回代过程是解同解的上三角形方程组回代过程可写为：

令令对对kk=nn1111令令2.1.2高斯消去法的运算量高斯消去法的运算量高斯消去法比起克莱姆法则和约当消去法，高斯消去法比起克莱姆法则和约当消去法，主要长处就是运算量较少。

主要长处就是运算量较少。

用克莱姆法则求解用克莱姆法则求解nn阶线性方程组，需计算阶线性方程组，需计算nn+1+1个行列式；每个行列式是个行列式；每个行列式是nn!

个乘积之和，而个乘积之和，而每个乘积是每个乘积是nn个元素相乘，因此共需乘法次数个元素相乘，因此共需乘法次数（nn+1）+1）nn!

（!

（nn1）=

（1）=（nn+1）!

（+1）!

（nn1）1）当当nn=20=20时时（nn+1）!

（+1）!

（nn1）9.7101）9.7102020。

要完成这要完成这么多次乘法，在每秒做一亿次乘法运算的计算么多次乘法，在每秒做一亿次乘法运算的计算机上，也需机上，也需30.830.8万年。

因此克莱姆法则在实际万年。

因此克莱姆法则在实际计算中不适用。

计算中不适用。

高斯消去法消去过程运算量：

高斯消去法消去过程运算量：

消元过程所用乘、除法次数为消元过程所用乘、除法次数为：

对k=1n-1（步）做对i=k+1n（行）做令（“适当倍数”）对j=k+1n+1（列）令回代过程回代过程运算量：

运算量：

回代过程求回代过程求xxnn需需11次除法，求次除法，求xxnn-1-1需需11次乘法、次乘法、11次除法，次除法，求需，求需nn11次乘法次乘法11次除法，因次除法，因此共需乘除次数：

此共需乘除次数：

两过程共需乘除次数两过程共需乘除次数。

令令对对kk=nn1111令令当当n=20n=20时，时，N=3060N=3060，显然比克莱姆法则少得显然比克莱姆法则少得多。

多。

线性代数教材讲授线性方程组解法时，线性代数教材讲授线性方程组解法时，往往提倡用行初等变换把增广矩阵化为最简往往提倡用行初等变换把增广矩阵化为最简形。

如上例方程组，多用如下变换：

形。

如上例方程组，多用如下变换：

变换结束后时增广矩阵最后一列就是解，不需回代。

这种解法称为约当消去法约当消去法。

它容易学习掌握，也容易编写计算机程序。

其计算过程可写为：

对对kk=1=1nn做做对对jj=kk+1+1nn+1+1令令对对II=1=1nn但但iikk做（做（“适当倍数适当倍数”）对对jj=kk+1+1nn+1+1令令不过这种解法需要乘除次数：

比高斯消去法多相当于多高斯消去法约一半的运算量。

因此除非方程组的阶数n较小，计算机上不用约当消去法。

上面比较了三种解法的乘除次数。

一般算法乘除次数与加减次数相差不多，而且计算机上乘除运算比加减运算费时间，所以通常只比较乘除次数。

2.1.32.1.3选主元消去法选主元消去法高斯消去法的弊端：

高斯消去法的弊端：

高斯消去法消去过程中，第k步中求n-k个倍数数用到的除数，称为主元。

它若为零或接近于零，计算机将“溢出”而停止计算，或产生较大误差。

例如方程组准确到九位小数的解是x1=0.250001875，x2=0.499998749，若在5位计算机上按高斯消去法求解，则有回代解得x1=0.00000，x2=0.50000，显然严重失真。

造成这种结果的原因，就是小主元的出现。

用它用除数产生大乘数，数乘大数变大数，大数吃小数产生舍入误差（如34105=-399997-1060.4000），从而有误差。

求出x2代回第一个方程时，因10-5x1+2x2=1，10-5x1+2x2=1两式相减得10-5（x1x1）+2（x2x2）=0，可见|x1x1|=200000|x2x2|表明x1的误差被放大200000倍，x1自然失真。

列主元消去法列主元消去法为了避免出现小主元，在每次消元前进行选主元。

即每次消元前先选取所要消元的列中绝对值最大的元素作为主元，然后再消元。

具体为：

具体为：

在第k步的第k列的元素中选主元，即在其中找出绝对值最大的元素，然后交换第k行和第p行，继续进行消元过程。

交换行相当于改变方程顺序，不会影响原方程组的解。

把它应用于上例，有回代得解=0.5,=0.25，接近真解。

算法算法:

第一步：

求主元apk使得第二步：

判断apk，则输出错误信息并停机，否则转三。

第三步：

判断若，则交换增广矩阵中第p行与第k行，否则不交换u也可在第k步时，在第k行元素akk，ak,k+1，akn中选绝对值最大的元素akq，交换k、p列，然后继续消去过程。

这种消去法称行主元法行主元法u也可在第k步时，在第kn行及第kn列的所有元素akk，ak,k+1，akn，ak+1k，ak+1,k+1，ak+1n，ann中选绝对值最大的元素apq，交换k、p行和k、q列，然后继续消去过程。

这种消去法称全主元法全主元法。

不过交换两列相当于改变未知量顺序，所以需用一数组，专门记录未知量顺序及其变化，并在最后调整解的顺序，使解中n个数与未知量正确对应。

三种主元消去法选取绝对值相对最大的元素作主元，保证了倍数c的绝对值不超过1，也有利于控制误差的传播，所以稳定性较好。

三种主元法中全主元法稳定性最好，但全主元法最费时间，程序最复杂。

通常情况下稳定性彼此相差不大，所以一般通常情况下稳定性彼此相差不大，所以一般情况都只用列主元消去法。

情况都只用列主元消去法。

复习题复习题1、何谓高斯消去法？

它与一般消去法有何不同？

怎样计算行列式？

2、计算机上为什么不用克莱姆法与约当消去法？

3、何谓主元消去法？

有何优点？