资源描述：
《线性方程组的直接解法ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、第五章线性方程组的直接解法§1Gauss消去法1.1顺序Gauss消去法1.2列主元Gauss消去法§2直接三角分解方法2.1Gauss消去法的矩阵运算2.2Doolittle分解法2.3平方根法2.4追赶法8/4/20211第五章线性方程组的直接解法在科学计算中，经常需要求解含有n个未知量的n个方程构成的线性方程组方程组还可以用矩阵形式表示为:Ax=b(5.1)8/4/20212第五章线性方程组的直接解法根据Gramer（克莱姆）法则，求解方程组（5.1）时，要计算大量的行列式，所需乘法次数大约为当n较大时,这个计算量是惊
2、人的。例如，当n=20时，约需乘法次数为N=9.7×1020若系数矩阵A非奇异，即det(A)≠0,则方程组有惟一解x=(x1,x2,…,xn)T.如果用每秒一亿次的计算机来计算，需要三十万年时间。可见Gramer法则不是一种实用的方法。因此，必须构造出适合于计算机使用的线性方程组的求解方法。N=(n+1)n!(n-1)=(n2-1)n!8/4/20213第五章线性方程组的直接解法直接方法的特点是:如果不考虑计算过程中的舍入误差，运用此类方法经过有限次算术运算就能求出线性方程组的精确解。求解线性方程组的数值方法可分为两大类：
3、直接方法和迭代方法。本章讨论直接方法，迭代方法将在下一章中讨论。需要指出，由于实际计算中舍入误差的存在，用直接方法一般也只能求得方程组的近似值。本章我们将给出直接解法的若干算法。8/4/20214第五章线性方程组的直接解法§1Gauss消去法Gauss（高斯）消去法是一种规则化的加减消元法基本思想通过逐次消元计算把需求解的线性方程组转化成上三角形方程组，也就是把线性方程组的系数矩阵转化为上三角矩阵，从而使一般线性方程组的求解转化为等价（同解）的上三角形方程组的求解。Gauss消去法由消元和回代两个过程组成，先讨论一个具体的线
4、性方程组的求解。8/4/20215第五章线性方程组的直接解法一、顺序Gauss消去法例1.用Gauss消去法解方程组用增广矩阵进行进算8/4/20216第五章线性方程组的直接解法这样，对于方程组(5.1)我们用增广矩阵表示，并给出gauss消去法的具体算法或者Ax=b8/4/20217第五章线性方程组的直接解法顺序Gauss消去法的消元过程可表述如下：第一步，设a11(1)≠0，将第一列中a11(1)以下的各元素消成零乘以矩阵[A（1），b（1）]的第一行再加到第i行，得到矩阵(i＝2，3，…，n)即依次用8/4/20218
5、第五章线性方程组的直接解法其中第二步，设a22(2)≠0，将第二列a22(2)以下各元素消成零，（i＝3，4，…，n）即依次用乘以矩阵[A（2），b（2）]的第二行再加到第i行，得到矩阵8/4/20219第五章线性方程组的直接解法其中如此继续消元下去直到第n－1步结束后，得到矩阵8/4/202110第五章线性方程组的直接解法增广矩阵[A（n），b（n）]对应如下上三角形方程组这是与原线性方程组（5.1）等价的方程组.8/4/202111第五章线性方程组的直接解法对于等价方程组进行回代求解，可以得到：8/4/202112第五章
6、线性方程组的直接解法首先写出增广矩阵于是，采用Gauss消去法求解方程组(5.1)8/4/202113第五章线性方程组的直接解法然后进行消元，采用公式最后进行回代得到方程组的解得到相似增广矩阵（i＝k+1，k+2，…，n）8/4/202114第五章线性方程组的直接解法在编程计算时，最后的增广矩阵存放的元素是：下面给出Gauss消去法的计算流程：8/4/202115第五章线性方程组的直接解法消元过程：j＝k＋1，…，n回代过程：对于i＝n－1，n－2，…，1，计算k＝1，2，…，n－1，i＝k＋1，k＋2，…，n对于执行计算置
7、8/4/202116第五章线性方程组的直接解法算法Gauss(A,a,b,n,x)1.消元Fork=1,2,…,n-11.1ifakk=0,stop;1.2Fori=k+1,k+2,…,n1.2.1lik=aik/akk=>aik1.2.2Forj=k+1，k+2,…,naij-aikakj=>aij1.2.3bi-aikbk=>bi2.回代2.1bn/an=>xn;2.2Fori=n-1,n-2,…,2,12.2.1bk=>S2.2.2Forj=k+1，k+2,…,nS–akjxj=>S2.2.3S/akk=>xka11a
8、12a13…a1nb1a21a22a23…a2nb2a31a32a33…a3nb3an1an2an3…annbn…………………………l21l31l41.ln1a22a23.a2nb2l32l42.ln2...a33.a3nb3l43.ln3a11a12a13.a1nb1a4nb4...an
