标签：begin end ddots 三角 矩阵 cdots 分解 pmatrix 21
文章目录
一、三对角方程组追赶法
A
x
=
f
Ax=f
Ax=f的系数矩阵呈对三角形
A
=
(
b
1
c
1
a
2
b
2
c
2
⋱
⋱
⋱
a
n
−
1
b
n
−
1
c
n
−
1
a
n
b
n
)
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}
A=⎝⎜⎜⎜⎜⎛b1a2c1b2⋱c2⋱an−1⋱bn−1ancn−1bn⎠⎟⎟⎟⎟⎞
A
=
L
U
A=LU
A=LU具有以下形式
L
=
(
1
l
2
1
l
3
⋱
⋱
1
l
n
1
)
L= \begin{pmatrix} 1 \\ l_{2} & 1 \\ &l_3& \ddots\\ && \ddots& 1\\ &&& l_n & 1\\ \end{pmatrix}
L=⎝⎜⎜⎜⎜⎛1l21l3⋱⋱1ln1⎠⎟⎟⎟⎟⎞
U = ( u 1 d 1 u 2 d 2 ⋱ ⋱ u n − 1 d n − 1 u n ) U= \begin{pmatrix} u_1&d_1 \\ & u_{2} & d_2 \\ &&\ddots& \ddots\\ &&&u_{n-1}& d_{n-1}\\ &&&& u_n\\ \end{pmatrix} U=⎝⎜⎜⎜⎜⎛u1d1u2d2⋱⋱un−1dn−1un⎠⎟⎟⎟⎟⎞
计算公式如下：
d
i
=
c
i
d_i=c_i
di=ci
u
1
=
b
1
u_1=b_1
u1=b1
l
i
=
a
i
/
u
i
−
1
l_i=a_i/u_{i-1}
li=ai/ui−1
u
i
=
b
i
−
l
i
c
i
−
1
u_i=b_i-l_ic_{i-1}
ui=bi−lici−1
计算次序是
u
1
,
l
2
,
u
2
,
l
3
,
u
3
,
⋯
,
l
n
,
u
n
u_1,l_2,u_2,l_3,u_3,\cdots,l_n,u_n
u1,l2,u2,l3,u3,⋯,ln,un
原方程组转化为
L
y
=
f
Ly=f
Ly=f
U
x
=
y
Ux=y
Ux=y
即可得解
求解三对角方程组前提条件是A是严格对角占优阵或对称正定
二、对称正定的Cholesky分解法
若A对称正定，则存在一个实的下三角阵L，使
A
=
L
L
T
A=LL^T
A=LLT
这种分解是唯一的
L
=
(
l
11
l
21
l
22
⋯
⋱
l
n
1
l
n
2
⋯
l
n
n
)
L= \begin{pmatrix} l_{11} \\ l_{21} & l_{22} \\ \cdots & & \ddots\\ l_{n1}&l_{n2}& \cdots&l_{nn}\\ \end{pmatrix}
L=⎝⎜⎜⎛l11l21⋯ln1l22ln2⋱⋯lnn⎠⎟⎟⎞
逐列计算L的元素，即
l
11
,
l
21
,
⋯
,
l
n
1
,
l
22
,
⋯
,
l
n
2
,
⋯
,
l
n
n
l_{11},l_{21},\cdots,l_{n1},l_{22},\cdots,l_{n2},\cdots,l_{nn}
l11,l21,⋯,ln1,l22,⋯,ln2,⋯,lnn
原方程组转化为
L
y
=
b
Ly=b
Ly=b
L
T
x
L^Tx
LTx=y
PS:Matlab有专用的函数来计算
L
L
L
L=chol(A)
为了避免相对耗时的平方根运算，将A分解为
A
=
L
D
L
T
A=LDL^T
A=LDLT
其中
L
=
(
1
l
21
1
⋯
⋱
l
n
1
l
n
2
⋯
1
)
L= \begin{pmatrix} 1 \\ l_{21} &1\\ \cdots & & \ddots\\ l_{n1}&l_{n2}& \cdots&1\\ \end{pmatrix}
L=⎝⎜⎜⎛1l21⋯ln11ln2⋱⋯1⎠⎟⎟⎞
D
=
(
d
1
d
2
⋱
d
n
)
D= \begin{pmatrix} d_1 \\ \ &d_2\\ & & \ddots\\ && &d_n\\ \end{pmatrix}
D=⎝⎜⎜⎛d1 d2⋱dn⎠⎟⎟⎞
称为改进的平方根法
标签：begin,end,ddots,三角,矩阵,cdots,分解,pmatrix,21 来源： https://blog.csdn.net/wang7301/article/details/120588354
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。