直接三角分解LU法

silence123

%直接三角分解法LU;
format rational
disp('请输入矩阵的阶数');
N=input('N=');
9 ^1 v2 F7 Q  ~2 t( A( ]for i=1:N
    for j=1:N
4 Y; H: d/ S& F        A(i,j)=input('A(i,j)=');
    end
end
L=eye(N);
5 G2 I' t" M( N" v) x. Lfor i=1:N
: j/ K/ j0 r/ r4 o+ x2 q    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
end
for r=2:N
    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);
4 ?' |1 j4 z. D+ ?' Z0 M        end
        L(i,r)=A(i,r);
        for k=1:r-1
; h) R. j+ a; H9 t& d! `            L(i,r)=L(i,r)-L(i,k)*U(k,r);
& v: b: N! v# y        end
/ H9 l$ F2 ]8 N% a        L(i,r)=L(i,r)/U(r,r);
    end
, G% u; v9 K; [$ O5 q& y- X, v4 jend
disp('A=');
disp(A);
disp('L=');
' T' _/ t, c: D+ H6 kdisp(L);
' q+ [6 O$ `, P& _disp('U=')
disp(U);

0.9清1.8清2.7清