直接三角分解LU法

silence123

%直接三角分解法LU;
format rational
disp('请输入矩阵的阶数');
, T, k/ D/ n0 N' B% {4 d" |, dN=input('N=');
1 A' Q5 c( \$ a* E; h  cfor i=1:N
- x6 X( @( F% p5 M8 ?# E    for j=1:N
        A(i,j)=input('A(i,j)=');
    end
4 u3 ^; Z6 Y4 [end
3 z, Y' F/ i% p1 BL=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
end
6 \2 Y# s7 Y, z; y) R* Efor r=2:N
0 _5 E* J% N( a7 T, Z3 F    for i=r:N
        U(r,i)=A(r,i);
        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
* r* i& {* O; M        end
9 ~( q9 v4 m/ ?/ r+ R        L(i,r)=A(i,r);
        for k=1:r-1
0 P; G' X( Y) i$ M0 ]0 Y            L(i,r)=L(i,r)-L(i,k)*U(k,r);
( A; O$ q& b6 M" O7 N' H        end
+ w* g  U$ _$ Z, J        L(i,r)=L(i,r)/U(r,r);
    end
end
disp('A=');
disp(A);
6 i" i# f! d/ [disp('L=');
disp(L);
: @) y$ T9 K: Kdisp('U=')
  u' A( Z! m' c5 Z$ X7 w- e/ F0 x3 ~disp(U);

0.9清1.8清2.7清