矩阵的LU分解

矩阵的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