消元和LU分解

矩阵乘法还有一种理解方式是
A=BC可以视为$A=\sum L_i\ast R_i$$L_i$是A的列向量,$R_i$自然是C的行向量

给定一个矩阵A
$\left[\begin{array}{cccc}1& 2& 3& 4\\ 2& 5& 6& 5\\ 3& 8& 6& 9\end{array}\right]$
先进行一波消元得到
$\left[\begin{array}{cccc}0& 0& 0& 0\\ 0& 1& 0& -3\\ 0& 2& -3& -3\end{array}\right]$
这个操作用矩阵语言描述就是,原始矩阵减去下面这个矩阵
$\left[\begin{array}{c}1\\ 2\\ 3\end{array}\right]$ * $\left[\begin{array}{cccc}1& 2& 3& 4\end{array}\right]$

同理,对下面两行再消元得到
$\left[\begin{array}{cccc}0& 0& 0& 0\\ 0& 0& 0& 0\\ 0& 0& -3& 3\end{array}\right]$

所以原始矩阵用消元的方式可以表示成如下形式
$\left[\begin{array}{c}1\\ 2\\ 3\end{array}\right]$ * $\left[\begin{array}{cccc}1& 2& 3& 4\end{array}\right]$ + $\left[\begin{array}{c}0\\ 1\\ 2\end{array}\right]$ $\left[\begin{array}{cccc}0& 1& 0& -3\end{array}\right]$+$\left[\begin{array}{c}0\\ 0\\ 1\end{array}\right]$$\left[\begin{array}{cccc}0& 0& -3& 3\end{array}\right]$

回到最初对矩阵乘法理解方式,可以知道
A = $\left[\begin{array}{ccc}1& 0& 0\\ 2& 1& 0\\ 3& 2& 1\end{array}\right]$*$\left[\begin{array}{cccc}1& 2& 3& 4\\ 0& 1& 0& -3\\ 0& 0& -3& 3\end{array}\right]$
这就是LU分解,把A拆成下三角矩阵L和上三角矩阵U的积