直接三角分解LU法

silence123

%直接三角分解法LU;
+ A3 U; R8 l" r" j% l, X; [format rational
! B) x6 a6 [" X) f8 zdisp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
    for j=1:N
& m5 \4 [1 e, ]& p4 G$ _        A(i,j)=input('A(i,j)=');
    end
end
! n2 T; B% J2 \5 Z7 {. O0 x- t! xL=eye(N);
* a: @( [3 k6 I5 T/ C/ x! Dfor i=1:N
$ q2 P- Z0 B: m/ Q) K! `# U    U(1,i)=A(1,i);
8 x0 v; }9 D9 N0 [    L(i,1)=A(i,1)/U(1,1);
& L9 x' Z# C" X& g, s8 H) X& d4 Pend
( Z9 ?& ^& U$ P9 Wfor r=2:N
    for i=r:N
        U(r,i)=A(r,i);
$ s& e- d# k9 k        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
4 ^8 Q! i- P" f        L(i,r)=A(i,r);
# T7 Z* f: }, |        for k=1:r-1
& N, [! Q& F. ?! R( K  v' q. p            L(i,r)=L(i,r)-L(i,k)*U(k,r);
5 j/ n) q, ^. g1 o5 Y' i        end
        L(i,r)=L(i,r)/U(r,r);
    end
end
' W+ ]7 v. K3 o  q8 R, w9 a) X! ~disp('A=');
disp(A);
% l( m1 A) m& |" M+ o1 ddisp('L=');
disp(L);
: ^) `' ]+ w1 f" S$ a8 [4 U, Odisp('U=')
0 R5 k5 z5 d) z/ T+ |/ wdisp(U);

0.9清1.8清2.7清