Question

大神门教教怎么用MATLAB验证LU分解的正确性啊

Answer 1

A = magic(3)
[L,U] = lu(A)
A2 = L*U
R = isequaln(A,A2);
if R
    disp('A和A2相等。');
else
    disp('A和A2不相等。');
end

%% 运行结果

A =

     8     1     6

     3     5     7

     4     9     2

L =

    1.0000         0         0

    0.3750    0.5441    1.0000

    0.5000    1.0000         0

U =

    8.0000    1.0000    6.0000

         0    8.5000   -1.0000

         0         0    5.2941

A2 =

     8     1     6

     3     5     7

     4     9     2

A和A2相等。