直接三角分解LU法

silence123

%直接三角分解法LU;
format rational
8 J3 g& F/ O8 H+ H) N# [disp('请输入矩阵的阶数');
7 K+ J% J. Y7 P8 HN=input('N=');
* i; V5 ~. k8 X8 y. N0 a3 ]6 Xfor i=1:N
    for j=1:N
        A(i,j)=input('A(i,j)=');
7 c* g3 _$ o4 t8 n    end
1 a+ _: g/ ], V0 ^& D' k  kend
3 W0 k' K  X4 ?L=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
end
for r=2:N
    for i=r:N
3 q9 l5 a( H$ ]" V4 C! u# O: e$ D        U(r,i)=A(r,i);
% s  K8 m1 c$ L* ]        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
4 @  Q4 i5 h6 I% Y& M- z, \        L(i,r)=A(i,r);
! \* U* C% O4 l        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
$ P8 W, F- k% Q9 \  ?) _3 Tend
disp('A=');
: q- n+ W& ^7 e4 idisp(A);
disp('L=');
disp(L);
9 a: c8 p. H* \4 T. sdisp('U=')
, F1 J) o; `7 S( A  k( {* {/ p& wdisp(U);

0.9清1.8清2.7清