一.三对角
1.1问题引入
问题引入:
A={
613861...........................86186},A= \begin{Bmatrix}6& 1 & 3 &\\ 8 & 6 & 1& \\ & ... & ...&...\\ &&...&...&...\\ &&&...&...&...\\ &&&&8&6&1\\&&&&&8&6\end{Bmatrix} ,A=⎩
⎨
⎧6816...31.....................8...6816⎭
⎬
⎫, b={
715...1514}b=\begin{Bmatrix} 7\\15\\ ...\\ 15\\14\end{Bmatrix}b=⎩
⎨
⎧715...1514⎭
⎬
⎫,则方程有解,则方程有解,则方程有解XXX使AX=bAX=bAX=b,取n=10求XXX
1.2三对角算法原理
首先将矩阵进行LU分解,设方程组AX=FAX=FAX=F,即
[b1c1a2b2c2.........aibici.........an−1bn−1cn−1anbn][x1x2............xn]=[f1f2............fn]\begin{bmatrix}b_1& c_1 & \\&a_2 & b_2 & c_2& \\&&...& ... &...\\&&& a_i &b_i&c_i\\&&&&...& ... &...\\&&&&&a_{n-1} & b_{n-1} & c_{n-1}& \\&&&&&&a_n &b_n\\\end{bmatrix} \begin{bmatrix}x_1\\x_2\\...\\...\\...\\...\\x_n\\\end{bmatrix} = \begin{bmatrix}f_1\\f_2\\...\\...\\...\\...\\f_n\\\end{bmatrix}
b1c1a2b2...c2...ai...bi...ci...an−1...bn−1ancn−1bn
x1x2............xn
=
f1f2............fn
则有则有A=LUA=LUA=LU,即
L=[1l21l31.........1ln1]L=\begin{bmatrix}1& \\l_2 & 1 \\&l_3& 1\\&&&...& ... \\&&&&...& 1 & \\&&&&&l_n &1\\\end{bmatrix}L=

最低0.47元/天 解锁文章
3150





