本文介绍了追赶法解决三对角线性方程组的过程，通过逐行和逐列更新系数来逐步求解。并通过Matlab函数实现，并与Gauss列主元消去法进行了时间消耗对比，展示了追赶法在高阶方程组中的优势。
一、问题描述
对于线性方程组
A x = b , A = ( b 1 c 1 a 2 b 2 c 2 ⋱ ⋱ ⋱ ⋱ ⋱ ⋱ a n − 1 b n − 1 c n − 1 a n b n ) , b = ( f 1 f 2 ⋮ f n ) Ax=b,\quad A=\begin{pmatrix}b_1&c_1&&&&\\a_2&b_2&c_2&&&\\&\ddots&\ddots&\ddots&& \\&&\ddots&\ddots&\ddots& \\&&&a_{n-1}&b_{n-1}&c_{n-1}\\&&&&a_n&b_n\end{pmatrix},\quad b=\begin{pmatrix} f_1\\f_2\\ \vdots\\ f_n \end{pmatrix} Ax=b,A=
b1a2c1b2⋱c2⋱⋱⋱⋱an−1⋱bn−1ancn−1bn
,b=
f1f2⋮fn
其中
{ ∣ b 1 ∣ > ∣ c 1 ∣ > 0 ∣ b n ∣ > ∣ a n ∣ > 0 ∣ b i ∣ ≥ ∣ a i ∣ + ∣ c i ∣ , a i c i ≠ 0 , i = 2 , 3 , ⋯ n − 1 \begin{cases} \begin{aligned}&\quad|b_1|>|c_1|>0\\&\quad|b_n|>|a_n|>0\\&\quad|b_i|\geq|a_i|+|c_i|,\quad a_ic_i\neq0,i=2,3,\cdots n-1\end{aligned} \end{cases} ⎩
⎨
⎧∣b1∣>∣c1∣>0∣bn∣>∣an∣>0∣bi∣≥∣ai∣+∣ci∣,aici=0,i=2,3,⋯n−1
二、追赶法解三对角方程组
A = ( b 1 c 1 a 2 b 2 c 2 ⋱ ⋱ ⋱ a n − 1 b n − 1 c n − 1 a n b n ) , L = ( 1 l 2 1 l 3 ⋱ ⋱ 1 l n 1 ) , U = ( u 1 d 1 u 2 d 2 ⋱ ⋱ u n − 1 d n − 1 u n ) \begin{equation*} A=\begin{pmatrix}b_1&c_1\\a_2&b_2&c_2\\&\ddots&\ddots&\ddots\\&&a_{n-1}&b_{n-1}&c_{n-1}\\&&&a_n&b_n\end{pmatrix}, {L}=\begin{pmatrix}1\\l_2&1\\&l_3&\ddots\\&&\ddots&1\\&&&l_n&1\end{pmatrix},{U}=\begin{pmatrix}u_1&d_1\\&u_2&d_2\\&&\ddots&\ddots\\&&&u_{n-1}&d_{n-1}\\&&&&u_n\end{pmatrix} \end{equation*} A= b1a2
最低0.47元/天 解锁文章
扫一扫
1万+
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
