MATLAB LU函数

　　高斯消元法求解线性方程，包括把增广矩阵转换为三角矩阵形式的过程，消去阶段工作

步骤是把矩阵A分解成为下三角L和上三角U的乘积。这种计算L,U的过程称为LU分解法。

 

lu实现对矩阵的分解。

[L,U] = lu(A)

%%将矩阵A分解的上三角矩阵保存在U当中，将一个“心理学上的”下三角矩阵（例如一个下三角矩阵和置换矩阵的乘积）保存在L中，满足A=L*U，注意A不必须是方阵。

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

%%返回三个矩阵，下三角矩阵L、上三角矩阵U和一个置换矩阵P，满足P*A=L*U。

[L,U,p] = lu(A,'vector')

%%将使用向量而不是矩阵的方式返回置换信息，p是一个行向量，满足A(p,:) = L*U，如果使用[L,U,P] = lu(A,'matrix')则返回的P是向量，这也是默认的形式。

 

[L1,U1]=lu(x)
[L2,U2,P]=lu(x)
[L3,U3,P,Q] = lu(X)
MATLAB中[L1,U1]=lu(x)的结果：L是下三角的置换矩阵即L1=p*L2，U是上三角阵。Matlab中LU分解采用高斯消元法，结果是不唯一的，只要[L1,U1]=lu(x)满足L1*U1=x, [L2,U2,P]=lu(x)满足L2*U2=p*x，[L3,U3,P,Q] = lu(X)满足 L3*U3= P*X*Q就行了