直接三角分解LU法

silence123

%直接三角分解法LU;
5 Z, ]- A0 C8 V# @format rational
disp('请输入矩阵的阶数');
8 q4 u" r" s" t& e7 r$ Q8 ?N=input('N=');
for i=1:N
/ T5 C# a  p9 o1 ?, o0 K    for j=1:N
! w( S) f4 ], O6 z        A(i,j)=input('A(i,j)=');
    end
' Q3 `1 ]8 T9 a6 A+ `. rend
L=eye(N);
9 \: E0 M; e7 y% ?for i=1:N
$ [, z7 q- d# j" D& i) v: L    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
, P. G! m# e2 e  ]5 iend
for r=2:N
8 ]4 a3 U9 U3 s( p" @    for i=r:N
        U(r,i)=A(r,i);
7 c/ I7 T; r" \/ n4 C        for k=1:r-1
* a0 M6 v5 R2 X) T( L0 s& Q" ]) U+ k            U(r,i)=U(r,i)-L(r,k)*U(k,i);
( M  z" E" T( {- X, p& E        end
        L(i,r)=A(i,r);
# f5 |) ]; M# [        for k=1:r-1
! _7 H) Y+ ~% u1 `3 c3 Q- y, k; J            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
        L(i,r)=L(i,r)/U(r,r);
    end
end
) p8 m) ^- n# t4 Idisp('A=');
) C9 @, a5 n# H) j& adisp(A);
7 ~8 a- J& ~7 d# Pdisp('L=');
disp(L);
) I, P, d; I9 S+ jdisp('U=')
disp(U);

0.9清1.8清2.7清