矩阵三角分解法.pdfVIP

第三节 矩阵三角分解法 1 直接三角分解法 将高斯消去法改写为紧凑形式,可以直接从矩阵A的元素得到计算L,U元素的递推公式,而不需任何中间 步骤,这就是所谓直接三角分解法。一旦实现了矩阵LU分解，那么求解Ax=b 的问题就等价于要求解两个三 角形方程组 ⑴ Ly=b,求y ; ⑵ Ux=y,求x. 设A为非奇异矩阵,且有分解式A=LU, L为单位下三角阵,U为上三角阵,即： 由(3.1)有: a1i =U1i (i=1,2…,n),得U的第1行元素; ai1 =Li1 U11 ，Li1 =ai1 /U11 (i=2, …,n),得L的第1列元素. 设已经定出U的第1行到第r-1行元素与L的第1列到第r-1列元素。由(3.1),利用矩阵乘法(注意 当rk时,Lrk =0),有 故 又由 (3.1)有 总结上述讨论，得到用直接三角分解法解Ax=b （要求A的所有顺序主子式都不为零）的计算公式。 ⑴ u1i =a1i (i=1,2,…,n), li1 =ai1 /u11 (i=2,3,…,n), 计算U的第r行，L的第r列元素(i=2,3,…,n)。 (2) (3) 求解Ly=b，Ux=y的计算公式； (4) (5) 例5 用直接三角分解法解 解:用分解公式 求解 矩阵A的分解公式(3.2),(3.3)又称为杜利特尔 (Doolittle)分解。 2. 追赶法 在一些实际问题中，例如解常微分方程边值问题，解热传导方程等，都会要求解系数矩阵为对角占优的三 对角方程组： 简记为Ax=f 。其中，当 |i-j|1时，a =0,且： ij (a) |b ||c |〉0 1 1 ( ) | |≥| | | | ( ) (b) |b |≥|a |+ |c |,a ,c ≠0 (i=2,3,…,n-1) i i i i i (c) |b ||a |〉0 n n 下面利用矩阵的直接分解发来推到接单对角线方程组 (3.12)的计算公式。 设A=LU ,其 L为下三角矩阵，U为单位上三角矩阵。即： 其 l ,u ,d 为待定系数。比较 (3.13)两边得： i i i 由(3.14)得到： 这样就确定了 {l }，{u }，{d }, 实现了A的LU分解。 i i i 求解Ax=f等价于解两个三角方程组： （1）L f 求 （2）U 求 （1）Ly=f,求y ; （2）Ux=y,求x 从而得到解三对角线方程组的追赶法公式： 1. 计算 {L },{u }的递推公式： i i 2. 解 Ly=f 3. 解Ux=y 计算系数L →L →…→L 及 y →y →…→y 的过程称为追的过程,计算方程组的解x →x 1 2 n-1 1 2 n n n-1 →…→x 的过程称为赶的过程。 1