直接三角分解LU法

silence123

%直接三角分解法LU;
format rational
# c9 c# g9 ^* T& F9 `disp('请输入矩阵的阶数');
N=input('N=');
9 t2 |1 T1 x. N/ X2 T5 i% Jfor i=1:N
    for j=1:N
% U$ X& F9 Z  a  b+ ?        A(i,j)=input('A(i,j)=');
5 m/ s  I! u, v5 d9 F1 v) d: U/ s& n    end
+ {7 i6 J9 ]: \- n( L/ W7 t+ Mend
L=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
end
  E+ @5 g4 X" W& c+ afor r=2:N
    for i=r:N
        U(r,i)=A(r,i);
, h* c% `* M) h* \! i        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
        L(i,r)=A(i,r);
* }2 J" Z0 G: _+ _6 s. _        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
! h1 J6 R. u, w+ L( r        end
        L(i,r)=L(i,r)/U(r,r);
    end
- ^0 E7 g0 O1 X3 _end
# n/ E9 B6 X/ l! Y" v/ Gdisp('A=');
) {1 X4 V# Q) wdisp(A);
5 T# {. ^5 f$ Y- V" Ldisp('L=');
disp(L);
4 F$ t6 _4 S3 ?' o& q" n5 w3 ]disp('U=')
- P: h% [1 h. o& _8 x, G% O1 T; udisp(U);

0.9清1.8清2.7清