直接三角分解LU法

silence123

%直接三角分解法LU;
/ q: j% m0 s2 p( U2 [, Bformat rational
disp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
( o* C8 p8 B2 r    for j=1:N
        A(i,j)=input('A(i,j)=');
    end
end
L=eye(N);
6 Y% B; t+ ]# ]( _) o) n3 Vfor i=1:N
8 J$ n2 g4 D9 V7 q1 S    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
+ r1 Z5 n. U& u) O7 p' Rend
1 @/ m8 ?& O3 {: A* k/ f/ }& ^6 Tfor r=2:N
$ Y2 _" @- s- E8 L) T$ Y/ K" T$ }, a    for i=r:N
. @) @, e2 u/ f/ Y: j1 l5 _        U(r,i)=A(r,i);
1 k, O% q% Q7 C- L        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
        L(i,r)=A(i,r);
        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
% T+ ?1 a' v: S. y        L(i,r)=L(i,r)/U(r,r);
    end
0 Q$ W7 `6 [! l" q+ Z2 Y6 p5 O9 Gend
disp('A=');
disp(A);
1 N8 \4 R) L; W! \1 s5 }  X% Ndisp('L=');
disp(L);
disp('U=')
disp(U);

0.9清1.8清2.7清