数值分析课件-num_2.4追赶法.pdfVIP

第二章 解线性方程组的直接法 2.4 追赶法 § 2.4 追赶法 § 对角占优矩阵: 补充 若矩阵A(a ) 满足 ij nn n |a | |a | i1,2,,n ii  ij j1 ji 则称A为严格对角占优矩阵. 若矩阵A(a ) 满足 ij nn n |a | |a | i1,2,,n ii  ij j1 ji 则称A为弱对角占优矩阵. 有一类方程组,在今后要学习的插值问题和边值问题中 有着重要的作用,即三对角线方程组,其形式为: Ax f (1) 其中 b c1 1   x  f   1 1 a b c      2 2 2  x2 f2 A    x  f    a b c       n1 n1 n1      a b xn fn  n n      A称为三对角线矩阵,并且满足 (1) |b ||c |0 1 1 (2) |b ||a ||c | , a c 0 i2,,n1 i i i i i (3) |b ||a |0 n n A称为对角占优的三对角线矩阵. 显然,A非奇异,即detA0 因此A的任意k阶顺序主子式非零,即detA 0k 所以,可以将A进行LU分解 L为单位下三角阵,U为上三角阵时,为Doolittle分解 L为下三角阵,U为单位上三角阵时,为Crout分解 以下以Doolittle分解导出三对角线方程组的解法 ( 以Crout分解的三对角线方程组的解法请参考教材) 设 ALU 由 Doolittle分解得: 1  u c     1 1  l 1