82 Matrix Factorization-Tridiagonal System >追赶法解三对角方程组 Crout Reduction for Tridiagonal Linear System*/ 与GE类似,一旦 sep:对A作(a=0则算法中断,故并非任何 三对角阵都可以用此方法分解。 比较等式两边 的元素,可得到计 ∵n1算公式(p52)。 Yu a Step2:追—即解Ly=∫:y (f-rv-1) sep3:赶—即解Ux=j:xn=p,x=y-月x(=n-1,…,1§2 Matrix Factorization – Tridiagonal System ➢ 追赶法解三对角方程组 /* Crout Reduction for Tridiagonal Linear System */                 =                                 − − − n n n n n n n f f f x x x a b a b c a b c b c 2 1 2 1 1 1 1 2 2 2 1 1 Step 1: 对 A 作Crout 分解                                 = − 1 1 1 1 2 1 n n n A       直接比较等式两边 的元素，可得到计 算公式（p.52）。 Step 2: 追——即解 L y f ：   = , 1 1 1  f y = ( 2, ... , ) ( ) 1 i n f r y y i i i i i = − = −  Step 3: 赶——即解U x y ：   = , ( 1, ... , 1) xn = yn xi = yi −  i xi+1 i = n− 与G.E.类似，一旦 i = 0 则算法中断，故并非任何 三对角阵都可以用此方法分解