LU分解 Ax=b, A=LU

LU分解

詳細算法： https://zhuanlan.zhihu.com/p/84210687

示例：https://www.geeksforgeeks.org/l-u-decomposition-system-linear-equations/

LU分解的實現可用用SVD分解。

問題描述: 解方程Ax=b, A可以分解成 LU，L是下三角矩陣，U是上三角矩陣。

LU分解的意義：減少計算量，便於程序化。

Ax=LUx=b,令y=Ux，則Ly=b,先求出y,再根據y=Ux求出x.

通過高斯消元法(行列變換）

A --> L0A -->L1A-->L2A->U

$\\L_{2}L_{1}L_{0}A=U\\ A=L_{0}^{-1}L_{1}^{-1}L_{2}^{-1}U=LU\\ L=L_{0}^{-1}L_{1}^{-1}L_{2}^{-1}$

示例AX=b

$\\ A=\begin{pmatrix} 1 & 1 & 1\\ 4 & 3 & -1\\ 3 & 5 & 3 \end{pmatrix}, X=\begin{pmatrix} x1\\ x2\\ x3 \end{pmatrix}, b=\begin{pmatrix} 1\\ 6\\ 4\\ \end{pmatrix}\\$

做變換，L0=r2-4r1,L1=r3-3r1 得到

$\\ L_1L_0A=\begin{pmatrix} 1 & 1 & 1\\ 0 & -1 & -5\\ 0 & 2 & 0 \end{pmatrix},$

L2=r3+2r2,得到

$L_2L_1L_0A=\begin{pmatrix} 1 & 1 &1\\ 0 & -1 & -5\\ 0 & 0 & -10\\ \end{pmatrix}=U$

$L_0=\begin{pmatrix} 1 & 0 & 0\\ -4 &1 & 0\\ 0 & 0 &1\\ \end{pmatrix},L_0^{-1}=\begin{pmatrix} 1 & 0 & 0\\ 4 &1 & 0\\ 0 & 0 &1\\ \end{pmatrix}$ 逆變換即對應元素相加爲0，

$L_1=\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ -3 & 0 &1\\ \end{pmatrix},L_1^{-1}=\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ 3 & 0 &1\\ \end{pmatrix}$

$L_2=\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ 0 & 2 &1\\ \end{pmatrix},L_2^{-1}=\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ 0 & -2 &1\\ \end{pmatrix}$

$L=L_0^{-1}L_1^{-1}L_2^{-1}=\begin{pmatrix} 1 & 0 & 0\\ 4 &1 & 0\\ 0 & 0 &1\\ \end{pmatrix}\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ 3 & 0 &1\\ \end{pmatrix}\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ 0 & -2 &1\\ \end{pmatrix}=\begin{pmatrix} 1 & 0 & 0\\ 4 &1 & 0\\3 & 0 &1\\ \end{pmatrix}\begin{pmatrix} 1 & 0 & 0\\ 0 &1 & 0\\ 0 & -2 &1\\ \end{pmatrix}=\begin{pmatrix} 1 & 0 & 0\\ 4 &1 & 0\\ 3 & -2 &1\\ \end{pmatrix}$

逆變換相乘：對應元素相加

LU分解 Ax=b, A=LU

Wireshark 安裝+使用（一）

C++ 如何正確高效地定義字符串常量

stl vector push_back解析

優秀資料網站

C++全局變量的初始化過程

LU分解 Ax=b, A=LU

Mac下配置sublime實現LaTeX

https://yachay.unat.edu.pe/blog/index.php?comment_area=format_blog&comment_component=blog&comment_co

linux以太網驅動總結