第四章 高斯消元法与选主元.pptVIP

;关键词：

高斯消去法，主元消去法;一问题的描述;其中,A是由系数按次序排列构成的一个n阶矩阵,

称为方程组的系数矩阵,x和b都是n维向量,称为方程

组的解向量和右端向量，即

使方程组(1)中每一个方程都成立的一组数x1*,x2*,…,xn*称为式(1)的解,把它记为向量的形式,称为解向量.

;我们总是希望方程组有解,且有唯一解.由线性代数的克莱姆(cramer)规则可知,如果方程组(1)的系数矩阵A的行列式(一般记为D=|A|)不等于零,那末,这个方程组有唯一解,而且它们可以表示为

xi=Di/D(i=1,2,…,n)

这里,Di是指D中第i列元素用右端b1,b2,…bn代替构成的行列式.

如果方程组(1)有唯一解,我们按上面的等式求解,就必须计算n+1个n阶行列式.由行列式的定义,n阶行列式包含有n!项,每一项含有n个因子,计算一个n阶行列式就需要做(n-1)n!次乘法.而我们一共要计算n+1个n阶行列式,共需做(n2-1)n!次乘法.此外,还要做n次除法才能算出xi(i=1,2,…,n).也就是说,用这个办法求解就要做

N=(n2-1)n!+n

次乘除法运算,这个计算量是大得惊人的.例如,当n=10(即求解一个含10个未知量的方程组),乘除法的运算次数共;;;当n=40,乘除法运算次数可达3.18??1049次.对于上百个未知量的方程组,次数运算量就更大了.因此克莱姆规则在理论上尽管是完善的,但在实际计算中却没有什么实用价值.本文我们将重点讨论求解线性方程组的一种有效的数值方法.

(二)求解线性方程组的消元法

消(元)去法是求解线性方程组

Ax=b(3)

和满秩矩阵的逆阵A-1的一种直接方法.尽管它比较古老,但它具有演算步骤,推算公式都系统化的特点(对其中选主元消去法,还可以证明是稳定的).因此,它至今仍然是求解方程组的一种有效的方法.

消去法可以引出几种计算方法,下面按三角形方程组和一般线性方程组的顺序来讨论.

;1.三角形方程组的解法

三角形方程组包括上三角形方程组和下三角形方程组，它是最简单的线性方程组之一，实际上消元法就是通过简化一般线性方程组为三角形方程组后再求解的。上三角方程组的一般形式是：

;;;;2.Gauss消元法;因为

可得第一步消元后的方程组为;第二步:???,取做(消去第i个方程组的x2,i=3,4,…n)

mi2?第二个方程+第i个方程i=3,4,…n

类似可得第二步消元后的方程组为

;第k步:设,取做(消去第i个方程组的xk,i=k+1,k+2,…,n)

mik?第k个方程+第i个方程

i=k+1,k+2,…n

类似可得第k步消元后的方程组为

;继续下去到第n-1步消元,可将线性方程组化为如下上三角方程组:;;;;;;3.选主元;例题4:讨论下面方程组的解法;列主元基本思想及程序框图

用高斯消去法求解线性方程组时,应避免小的主元.在实际计算中,进行第k步消去前,应该在第k列元素aik(i=k,…,n)中找出绝大值最大者,例如

∣a∣=max∣a∣

再把第p个方程与第k个方程组进行交换,使apk(k-1)成为主元.我们称这个过程为选主元.由于只在第k列元素中选主元,通常也称为按列选主元(或称部分选主元).

如果在第k步消去前,在第k个方程到第n个方程所有的xk到xn的系数a(i=k,…,n;j=k,…,n)中,找出绝对值最大者,例如

∣a∣=max∣a∣

再交换第k,p两个方程和第k,q两个未知量的次序,使a成为主元.称这个过程为完全选主元.

不论是哪种方式选出主元,而后再按上面介绍的计算步骤进行消去的计算,一般都称为选主元的高斯消去法.在实际计算中,常用按列选主元的高斯消去法.

;如：经第一次消元后，增广矩阵为：

则：

列主元消元法：从这一列中选取一个主元

素，作为第二次消元的开始；

全主元消元法：从

这一子矩阵中选取一个主元素，作为第二次消元的开始；

;用列主元消去法解线性方程组AX=