直接三角分解LU法

silence123

%直接三角分解法LU;
format rational
disp('请输入矩阵的阶数');
N=input('N=');
' G  }! z  a$ A! J0 D: ?for i=1:N
# d; x( X, d. Q" ^5 S( M- w    for j=1:N
        A(i,j)=input('A(i,j)=');
    end
end
% E3 n6 f0 `& v9 {. D" d; rL=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
; M" `2 [+ v4 w. aend
9 M( e5 \! N7 J- C* ^6 Ofor r=2:N
    for i=r:N
        U(r,i)=A(r,i);
        for k=1:r-1
; @4 ^* `& @& S- @# {3 E            U(r,i)=U(r,i)-L(r,k)*U(k,i);
- Z: ^: L7 K" u+ {6 t, Z# @! @% @        end
2 H7 J1 ]9 m8 Q' A/ K. p        L(i,r)=A(i,r);
        for k=1:r-1
  {) s% T$ \) j/ F( E, T            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
        L(i,r)=L(i,r)/U(r,r);
4 H! k8 {3 m. B7 \    end
5 E5 f7 v9 l$ w2 e. nend
9 [) P0 F9 X+ `  V+ ]; d9 qdisp('A=');
disp(A);
disp('L=');
disp(L);
7 T$ [; R  U+ o" ~9 q8 S" Qdisp('U=')
disp(U);

0.9清1.8清2.7清