直接三角分解LU法

silence123

%直接三角分解法LU;
- v; b9 n7 g. P+ u# P1 v" Sformat rational
disp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
+ f" p9 U5 Q5 z% y    for j=1:N
" j' N  u& O/ }3 a# ^2 f        A(i,j)=input('A(i,j)=');
    end
/ L" @: I( i2 Z- i6 s( Qend
) ^3 V  X: B& q4 F7 a$ x3 X, eL=eye(N);
' [& E  F/ D6 ]' v# R/ O. `" Ifor i=1:N
4 U1 [1 T4 N, G+ i    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
7 i- u8 b/ C+ G: N4 a, rend
* I5 e! J% a. `; n* e4 a$ Kfor r=2:N
    for i=r:N
1 W* k# x/ y; r( f* Z        U(r,i)=A(r,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);
        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
- k3 N9 `3 [) }+ u        L(i,r)=L(i,r)/U(r,r);
    end
" t% P$ y% E: v0 Fend
disp('A=');
! X' L. C) ~' g6 xdisp(A);
' C3 Q  M+ d3 V6 Z' I  d) F$ x9 {. ldisp('L=');
4 E; r/ I+ x# {4 H8 y8 [* Hdisp(L);
. P4 O  U. p7 M" s) N6 }5 Z/ Wdisp('U=')
1 C& I4 {7 h% F+ B/ W5 J' L: Bdisp(U);

0.9清1.8清2.7清