矩阵的LU分解是线性方程组求解方法中高斯消去法的基础,在MATLAB中由函数lu来实现。lu函数常用的语法格式如下。
[L,U] = lu(X):U矩阵为上三角矩阵,满足X = L*U。
[L,U,P] = lu(X):返回的矩阵P是置换矩阵,矩阵U是上三角矩阵,矩阵L是满秩矩阵。上述矩阵满足关系L*U = P*X。
在命令窗口中输入
code.matlab
>> a=[1 2 0;2 5 -1;4 10 -1];
>> [b,c,p]=lu(a)
b =
1.0000 0 0
0.2500 1.0000 0
0.5000 0 1.0000
c =
4.0000 10.0000 -1.0000
0 -0.5000 0.2500
0 0 -0.5000
p =
0 0 1
1 0 0
0 1 0