数值分析－线性代数方程组的直接解法.pptxVIP

会计学;引例: 设有电源及一些电阻组成的简单电路，求各环路电流。;;; 或写为矩阵形式即得到环路电流I1 ,I2 , I3所满足的线性方程组： (R1 + R4 + R5) -R4 - R5 I1 E1 -R4 (R2 + R4 + R6) - R6 I2 = E2 -R5 - R6 (R3 + R5 + R6) I3 E3 ;第3章 解线性方程组的直接法;线性方程组的数值解法一般有两类： 直接法：就是经过有限步算术运算，可求得方程组精确解的方法（若计算过程中没有舍入误差），如克莱姆法则就是一种直接法，直接法中具有代表性的算法是高斯(Gauss)消去法。 迭代法: ( 第四章介绍）就是用某种极限过程去逐步逼近线性方程组的精确解的方法。也就是从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。(一般有限步内得不到精确解); 3.2 解线性方程组的直接法（高斯消去法） ;（1）消元过程 第1步:将方程①乘上(-2)加到方程 ②上去,将方程 ①乘上 加到方程 ③上去,这样就消去了第2、3个方程的 项,于是就得到等价方程组 ;第2步：将方程 ④乘上 加到方程 ⑤上去，这样就消去了第3个方程的 项，于是就得到等价方程组 ;前述的消元过程相当于对原方程组 ; 由此看出,高斯消去法解方程组基本思想是设法消去方程组的系数矩阵A的主对角线下的元素,而将Ax=b化为等价的上三角形方程组,然后再通过回代过程便可获得方程组的解。换一种说法就是用矩阵行的初等变换将原方程组系数矩阵化为上三角形矩阵,而以上三角形矩阵为系数的方程组的求解比较简单,可以从最后一个方程开始,依次向前代入求出未知变量 这种求解上三角方程组的方法称为回代, 通过一个方程乘或除以某个常数,以及将两个方程相加减,逐步减少方程中的变元数,最终将方程组化成上三角方程组,一般将这一过程称为消元,然后再回代求解。;高斯消去法算法构造 我们知道,线性方程组(3.1)用矩阵形式表示为 ;⑴ 消元过程,高斯消去法的消元过程由n-1步组成： 第1步 设 ,把(3.3)中的第一列中元素 消为零,令 ;第k步 （k=2,3,…,n-1）继续上述消元过程，设第k-1次消??已经完成，得到与原方程组等价的方程组 ;只要 ,消元过程就可以进行下去,直到经过n-1次消元之后，消元过程结束，得到与原方程组等价的上三角形方程组,记为 ;;（3）高斯消去法的计算步骤： ① 消元过程;设 计算 ;(4) 高斯消去法流程图 ，见P42 (5) Gauss消去法计算量 ≈;第k 步;;设方程组系数矩阵 ，其顺序主子式 ;定义3.1 设矩阵 每一行对角元素的绝对值都大于同行其他元素绝对值之和 ;证:先考察消元过程的第1步,因A为严格对角占 优,故 故 ,又根据高斯消 去公式得 于是 ; 一般线性方程组使用高斯消去法求解时，在消元过程中可能会出现 的情况，这时消去法将无法进行；即使 ，但它的绝对值很小时，用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，将严重影响计算结果的精度。实际计算时必须避免这类情况的发生。主元素消去法就可弥补这一缺陷。 ;交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k) 为主元素，并称使用主元素的消元法为主元素法 根据主元素选取范围分为：列主元素法、行主元素法、全主元素法;主元素法的意义;这时 ;全主元素消去法 是通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为 称这样的 为主元素。尽管它的算法更稳定，但计算量较大，实际应用中大多数使用列主元素消去法即可满足需要。 ;全主元素法不是按列选主元素，而是在全体待选系数中选取，则得全主元素法。 例3.3 用全主元素法解下列线组;计算m21=-19/40=0.475，m31=4/40=0.1 (5)- m21(4), (6)- m31(4) 消去x2 得;保留有主元素的方程;列主元素法;(5)- m21