数据结构Ch3解非线性方程组的直接解法su教程文件.ppt

第三章 解线性方程组的直接法;§1 高斯消去法;解 ② +①×12/18 ,③ +①×1/18, 得 ; ③+②×1/1.167得 ;(1) (2) (n);回代公式;? 运算量 /* Amount of Computation */;? 运算量 /* Amount of Computation */;Step k：设 ，计算因子 且计算;消元乘除次数： 回代过程：; 前述顺序消去法是按序通过用a11,a22(1),…,an-1,n-1(n-2) (akk(k-1) ≠0)作为除数来达到消元目的的，因此应特别注 意到元素 (称之为主元素).因为如果主元素为0,消元 过程将无法进行下去.如果主元的绝对值较小,作为除数也 不合适.所以这种方法是不恰当的. ;例如用顺序消去法求解下列方程组(用四位有效数字计算):;②-×105得 (1.000-1.000×105)x2=1.000-0.6000×105 (0.0000-0.1000×106)x2=0.0000-0.0600×106 化简可得 x2=0.6000 回代入①求得 x1=105(0.6-0.6000)=0 而方程组的解应为 x1=0.4000 x2=0.6000 ; 显然用上述方法求出的解x1与方程组的实际解相差很大。若改变两个方程的顺序,即 x1+x2=1 ① 10-5 x1+x2=0.6 ② ②-①×10-5得 (1.000-1.000×10-5)x2=0.6-1.000×10-5 化简得 x2=0.6000 回代求得 x1=(1-0.6000)=0.4000; 高斯主元素消去法是顺序消去法的一种改进。它的基本思想是在逐次消元时总是选绝对值最大的元素(称之为主元)做除数,按顺序消去法的步骤消元。 这里主要介绍求解线性方程组最常用的列主元素消去法和全主元素消去法。 1.列主元素消去法 所谓列主元素消去法就是在每一步消元过程中取系数子矩阵的第一列元素中绝对值最大者作主元。对线性方程组(3.1)进行n-1次消元后,可得到上三角形方程组 ; 必须指出的是方程组(3.5)中的aij（k）(i≤j) 并非与顺序消去法得到的（3.4）相同。 然后就可对方程组(3.5)回代求解。;解 ②中-18为主元, 交换②和①得 ;②+①×12/18, ③+①×1/18得 ;③+②×1/1167得 ;列主元素消去法的计算过程为: (1)消元过程 对于k=1,2,…,n-1进行下述运算: ①选主元,确定r,使得 若ark（k-1）=0,说明系数矩阵为奇异,则停止计算否则进行下一步。 ②交换A的r、k两行，即互换arj（k-1） 、 akj（k-1）的值 ③计算; 2. 全主元素消去法 所谓全主元素消去法,就是每步消元时选取系数子矩阵中绝对值最大的元素作主元。经过n-1次消元后,方程组(3.1)可化为上三角形方程组，后回代求解。 ; 2. 全主元素消去法 所谓全主元素消去法,就是每步消元时选取系数子矩阵中绝对值最大的元素作主元。经过n-1次消元后,方程组(3.1)可化为上三角形方程组，后回代求解。 ;例 用全主元素消去法 求解方程组 ;① ② ③ ;③-②×0.944/2.333得 ;在计算机上求解时,我们仍把增广矩阵(A, b) 记为A,具体计算过程为 (1)消元过程 对k=1,2,…,n-1进行下列运算:; 图 3.2 ;? Complete Pivoting:;证明：;用矩阵的分块理论可设:;定理3.1 若方阵A的各阶主子式都不为零,则存在唯一的单位下三角矩阵L和上三角矩阵U,使得LU=A。证明：;;则;下证此分解时唯一的：设A还有另外的三角分解,即有;此定理是矩阵分解的重要定理,称作矩阵的LU分解,也叫杜利特尔分解.它是用直接法解线性方程组的理论基础.此外,可利用LU分解来计算A的行列式的值.;3.2.2