矩阵因式分解（LU矩阵分解）与GSL实现

   矩阵的因式分解是把一个矩阵A表示为两个或更多个矩阵的乘积，是将复杂的数据进行分解，其中有多种方法，例如：LU分解，秩分解，QR分解，奇异值分解，谱分解等。这里主要介绍对LU分解的认识。

根据参考的书籍，这里的LU分解只限于一系列具有相同系数矩阵的线性方程：

Ax=b₁,   Ax=b₂, Ax=b_{p                  （1）}

当A为可逆矩阵时，可计算A^-1，然后计算A^-1 b₁，A^-1 b₂，等等。但是，真正在社会实践的运用中，又是如何计算并使用的呢？ 实际而言，（1）中的第一个方程是由行变换解出的，并同时得出矩阵A的LU分解。

设A为m×n阶矩阵，则A_m×n可进行化简为阶梯形，此时不必行对换，那么A可写成形式A=LU，L是m×m下三角矩阵，主对角线元素全是1，U是A的一个等价的m×n阶梯形矩阵。如下：

这样的一个分解称为LU分解，矩阵L是可逆的，我们称L为单位下三角矩阵。

由上，我们可知，当A=LU时，方程Ax=b可写成L（Ux）=b，把Ux写成y，可以有解下面一对方程来求解x：

                  Ly=b

                  Ux=y