LU 分解

由此，矩阵将被分解为上部三角矩阵和下三角矩阵。通常，它将用于提高 Gauß消除的性能和稳定性（如果通过排列完成）。

然而，这种方法通常不会或不好地工作，因为它不稳定。例如

A = [8 1 6
    3 5 7
    4 9 2];
[L,U] = lu(A);

添加置换矩阵就足以使 PA = LU：

[L,U,P]=lu(A);

在下文中，我们将绘制 lu 的运行时间，该运行时间取决于矩阵元素的平方根。 http://i.stack.imgur.com/cnhNK.jpg