现代数值分析 第2章线性方程组的直接法.pptVIP

为避免过多的开方运算，在更多情况下是将 A 分解 为： A = LDL T ，其中 L 为下三角， D 为对角阵。事实上， 由前面的讨论可知： A=LU ，其中 L 为单位下三角阵， U 为上三角阵。记 U 的元素为 u ij ， D =diag( u ii ) ，由于 A 对称正定，必有 u ii ≠0( i =1,2,…, n ) ，所以 U=DU * ，其中 * 1 * T T U D U A = A L =U A= LDL ? ? 为单位上三角阵 , 由 , 知 , 因此 。 注： 将矩阵 A 作 Doolittle 分解或 Crout 分解 , 由矩阵乘 法可得 A 的 LDL T 分解。 ? 解三对角方程组的 追赶法 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? n n n n n n n f f f x x x b a c b a c b a c b 2 1 2 1 1 1 1 2 2 2 1 1 定理 2.2.3 ： 若 A 为 对角占优 的三对角阵，且满足 则方程组有唯一的 LU 分解。 0 , 0 , 0 | | | | , 0 | | | | 1 1 ? ? ? ? ? ? i i n n c a a b c b 第一步 : 对 A 作 Doolittle 分解 追赶法公式的推导：（以四阶为例） 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 1 1 2 1 2 1 2 2 3 2 3 2 3 3 4 3 4 3 4 1 1 2 2 2 1 3 3 3 2 1 1 1 1 ( 1, 2,3); i i b c u d l a b c u d A LU l a b c u d l a b u u d l u l d u d l u l d u d l u l d u u b u b l d d c i u b l d ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 比较对应位置的元素，得递推公式： 第二步： 2 2 1 3 3 2 4 4 3 4 4 4 3 / ; / / l a u l a u l a u u b l d ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1 1 1 1 2 1 2 2 2 2 2 3 3 3 3 2 3 3 4 4 4 4 3 4 4 1 1 2 2 2 1 3 3 3 2 4 4 4 3 1 1 1 1 y f y f l y y f l y f Ly f l y f l y y f y f l l y y f y f y f l y y f l y L y y f l y y f ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 第三步：由 ： 先求 该过程称为“ 追 ”的过程。 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 4 4 3 3 3 4 3 3 3 4 3 2 2 2 3 2 1 1 1 2 1 / ( ) / ( ) / ( ) / ( ) / u d x y u d x y Ux y u d x y u x y x y u x y d x u y c x u x y c x u x y c x Ux y x u ? ?? ? ? ? ? ?? ? ? ? ? ? ? ? ?? ? ? ? ? ?? ? ? ? ? ?? ? ? ? ? ?? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 第四步：再由 回代解得 解得： 该过程称为“ 赶 ”的过程。 一般情形的三对角方程组计算公式： 1 1 1 1 , 1,2, , 1 / , 2,3, , 1 , 2,3, , i i i i i i i i i d c i n u b l a u i n u b l c i n ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? L L L 计算次序为： 1 1 2 2 3 3 (1) (2) (3) i i n n d c u b l u l u l u ?