直接三角分解LU法

silence123

%直接三角分解法LU;
format rational
* j9 o/ A( u$ r3 kdisp('请输入矩阵的阶数');
N=input('N=');
* m' P8 Y: ?. y- e. cfor i=1:N
    for j=1:N
+ B0 {" j. h# d" S, c; U        A(i,j)=input('A(i,j)=');
/ M5 M* c* v, I* U4 d/ V: L. |    end
end
L=eye(N);
% i: T, t; x. [, U' lfor i=1:N
2 G0 v8 \  s9 q# m    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
/ w( F0 R% t# b& g3 e. \end
for r=2:N
    for i=r:N
9 Z+ F: E" \# E1 W# W* y& }        U(r,i)=A(r,i);
. U; Y- ?& E8 F        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
! g  ?! e- @! L. }% O        L(i,r)=A(i,r);
, v( N* a/ w2 e/ F  O7 r8 `6 S" b        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
        L(i,r)=L(i,r)/U(r,r);
    end
end
disp('A=');
disp(A);
9 @9 T( r/ r1 c7 I  v/ Z$ Pdisp('L=');
7 R) U3 [8 T* N8 Ydisp(L);
' O: a0 W6 E' K  y7 o5 ndisp('U=')
  c# T: _' R. k- ~5 q- Mdisp(U);

0.9清1.8清2.7清