三角矩阵的lu分解.docVIP

?三角矩阵的lu分解

?在线性代数中，LU分解(LUDecomposition)是矩阵分解的一种，可以将一个矩阵分解为一个单位下三角矩阵和一个上三角矩阵的乘积（有时是它们和一个置换矩阵的乘积）。LU分解主要应用在数值分析中，用来解线性方程、求反矩阵或计算行列式。

什么是LU分解

如果有一个矩阵A，将A表示成下三角矩阵L和上三角矩阵U的乘积，称为A的LU分解。

?

更进一步，我们希望下三角矩阵的对角元素都为1：

?

一旦完成了LU分解，解线性方程组就会容易得多。

LU分解的步骤

对于满秩矩阵A来说，通过左乘一个消元矩阵，可以得到一个上三角矩阵U。

可以看到，L实际上就是消元矩阵的逆。容易知道二阶矩阵的逆：

?

现在假设A是一个3×3矩阵，在不考虑行交换的情况下，通过消元得到上三角矩阵的过程是：

?

LU分解的前提

并非所有矩阵都能进行LU分解，能够LU分解的矩阵需要满足以下三个条件：

矩阵是方阵（LU分解主要是针对方阵）；

矩阵是可逆的，也就是该矩阵是满秩矩阵，每一行都是独立向量；

消元过程中没有0主元出现，也就是消元过程中不能出现行交换的初等变换。

LU分解的意义

LU分解的意义在于求解大型方程组。一个方程组可以简化为Ax=b的形式，其中A是n阶方阵，x是未知数组成的向量，b是n×1矩阵，例如：

以往求解的方式有两种，一是高斯消元法，二