matlab追赶法求解三对角方程组,【良心出品】MATLAB 追赶法求解三对角方程组的算法原理例题与程序...

3)三对角形线性方程组

123456789104100000000141000000001410000000014100000000141000000001410000000014100000000141000000001410000000014x x x x x x x x x x -⎡⎤⎡⎤⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎢⎥--⎢⎢⎥⎢⎢⎥-⎣⎦⎣⎦

7513261214455⎡⎤⎢⎥⎢⎥⎢⎥-⎢⎥

⎢⎥

⎢⎥=⎢⎥-⎢⎥

⎢⎥⎢⎥

-⎢⎥

⎥⎢⎥⎥⎢⎥

⎥⎢⎥-⎣⎦

*(2,1,3,0,1,2,3,0,1,1)T

x =---

二、数学原理

设系数矩阵为三对角矩阵

1

122233111000000000

000000

n n n n

n b c a b c a b A a b c a b ---⎛⎫ ⎪ ⎪ ⎪=

⎪ ⎪ ⎪

⎪ ⎪⎝

⎭

则方程组Ax=f 称为三对角方程组。

设矩阵A 非奇异，A 有Crout 分解A=LU ，其中L 为下三角矩阵，U 为单位上三角矩阵，记

1

122

233

1

10

00010

000

0001000

000100,00000000

00

0001n n n

n b L U γαβγββγβ--⎛⎫⎛⎫ ⎪

⎪ ⎪ ⎪ ⎪ ⎪∂==

⎪

⎪ ⎪

⎪ ⎪ ⎪

⎪ ⎪ ⎪ ⎪∂⎝

⎭

⎝

⎭

可先依次求出L ，U 中的元素后，令Ux=y ，先求解下三角方程组Ly=f 得出y ，再求解上三角方程组Ux=y 。