(完整word版)追赶法求解三对角方程及其算例.docxVIP

追赶法求解三对角方程组 要求：对于给定的三对角系数矩阵和右端项，可以求 解线性代数方程组 一、 追赶法的数学理论 设系数矩阵为三对角矩阵 b1 c1 0 0 0 0 a2 b2 c2 0 0 0 0 a3 b3 0 0 0 A 0 0 0 an 1 bn 1 cn 1 0 0 0 0 an bn 则方程组 Ax=f 称为三对角方程组。 设矩阵 A 非奇异， A 有 Crout 分解 A=LU ，其中 L 为下三角矩阵， U 为单位上 三角矩阵，记 b1 0 0 0 0 1 1 0 0 0 2 2 0 0 0 0 1 2 0 0 0 3 3 0 0 , 0 0 1 0 0 L U 0 0 0 n 1 0 0 0 0 0 n 1 0 0 0 n n 0 0 0 1 可先依次求出 L ，U 中的元素后， 令 Ux=y ，先求解下三角方程组 Ly=f 得出 y， 再求解上三角方程组 Ux=y 。 事实上，求解三对角方程组的 2 追赶法将矩阵三角分解的计算与求解两个三角 方程组的计算放在一起，使算法更为紧凑。其计算公式为： 1 b1, 1 c1 , y1 f1 1 1 对 i 2,3, , n ai, ibi ai i 1, ci i i i yi fi i yi 1 (*) i xn yn 对 i n 1,n 2, ,1 xi yi i xi 1 二、 追赶法的算法和