直接三角分解LU法

silence123

%直接三角分解法LU;
9 ~9 U' g6 o4 x. V' H0 R4 M; jformat rational
disp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
6 h9 G8 z, W9 g) b) a5 w8 ~    for j=1:N
$ V7 H# w1 ~, }7 U9 o        A(i,j)=input('A(i,j)=');
+ W+ }, a! z' G% i! j6 I    end
) z( v+ g3 n# Q; {7 cend
" H8 ?+ U7 h/ u6 P1 e, iL=eye(N);
for i=1:N
! }; \/ m% f, y. K( |7 Y9 B    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
- F* m& d) u, lend
: l" [' b. [  L& m8 qfor r=2:N
    for i=r:N
/ M. T1 O1 j$ d* e1 \        U(r,i)=A(r,i);
        for k=1:r-1
0 B6 g( e' N. V. H7 Q9 V: }            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
        L(i,r)=A(i,r);
' M- w& \# m4 P' R9 ]( l        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
) E, X/ ]/ Z# P/ L8 J& v' U        L(i,r)=L(i,r)/U(r,r);
9 c) f* i# s) d. \. Z9 g    end
' \( C  L  k8 O& d5 x7 K" ]. jend
disp('A=');
disp(A);
disp('L=');
disp(L);
disp('U=')
disp(U);

0.9清1.8清2.7清