2021-01-11

矩阵LU分解

存在矩阵A=$\left[\begin{array}{ccc}1& 0& 1\\ a& a& a\\ b& b& a\end{array}\right]$
L=$\left[\begin{array}{ccc}1& 0& 0\\ i1& 1& 0\\ i2& i3& 1\end{array}\right]$为对角线为1的下三角矩阵
U=$\left[\begin{array}{ccc}u11& u12& u13\\ 0& u22& u23\\ 0& 0& u33\end{array}\right]$为上三角矩阵
矩阵A乘以一个对角线为1的下三角型矩阵就是进行初等行变换
而我们目标是把A变成上三角
$A\left[\begin{array}{ccc}1& 0& 1\\ a& a& a\\ b& b& a\end{array}\right]⟹B\left[\begin{array}{ccc}1& 0& 1\\ 0& a& 0\\ b& b& a\end{array}\right]⟹C\left[\begin{array}{ccc}1& 0& 1\\ 0& a& 0\\ 0& b& a-b\end{array}\right]⟹U\left[\begin{array}{ccc}1& 0& 1\\ 0& a& a\\ 0& 0& a-b\end{array}\right]$

$E21\left[\begin{array}{ccc}1& 0& 0\\ -a& 1& 0\\ 0& 0& 1\end{array}\right]\ast A=B$

$E21中第一行\ast A=B中第一行$
$E21中第二行\ast A=B中第二行$
$E21中第三行\ast A=B中第三行$
$B二行=\left[\begin{array}{ccc}1& 0& 0\end{array}\right]\ast -a+\left[\begin{array}{ccc}-a& 1& 0\end{array}\right]\ast 1＋\left[\begin{array}{ccc}0& 0& 1\end{array}\right]\ast 0=\left[\begin{array}{ccc}0& a& 0\end{array}\right]$
$这就是矩阵中的向量乘法$

$E31\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ -b& 0& 1\end{array}\right]\ast B=C$


$E32\left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& -b/a& 1\end{array}\right]\ast C=U$


$(E32\ast E31\ast E21)A=U$
$A=E21^{-1}\ast E31^{-1}\ast E32^{-1}\ast U=LU$
$L=\left[\begin{array}{ccc}1& 0& 0\\ a& 1& 0\\ 0& 0& 1\end{array}\right]\ast \left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ b& 0& 1\end{array}\right]\ast \left[\begin{array}{ccc}1& 0& 0\\ 0& 1& 0\\ 0& b/a& 1\end{array}\right]这就是著名LU分解$

LU分解短视频
LU分解的快速求解—矩阵的LU分解步骤-待定系数法
Typora数学公式汇总（Markdown）

                                                                                                -2021-01-011