直接三角分解LU法

silence123

%直接三角分解法LU;
5 R2 J! F0 F6 g' r7 mformat rational
0 e0 r0 F3 D5 M9 Gdisp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
1 C& }+ F- Y) Z3 G    for j=1:N
# W, t; o" W8 N* {+ @% R1 d* R! e        A(i,j)=input('A(i,j)=');
    end
! v6 S5 E, N- o; lend
L=eye(N);
for i=1:N
3 e. S8 G( |# z6 S6 P& w6 E    U(1,i)=A(1,i);
! B- b4 b: q) x  i! T    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);
        end
: L" P+ \, g! A( K- Y        L(i,r)=A(i,r);
8 |4 X9 _3 g$ Q3 p) C) b- u        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
* p$ h1 u# b2 x/ {        end
1 s$ q; a0 N# J1 x5 y# @        L(i,r)=L(i,r)/U(r,r);
+ S4 l. X% Q3 ~1 L, U, b    end
7 T2 d# {& _& }3 C3 Q7 Oend
disp('A=');
# _  T  `% N7 v1 r8 u# {- J9 l" {disp(A);
disp('L=');
disp(L);
disp('U=')
- z/ Z; R( {" z- y1 U8 ?disp(U);

0.9清1.8清2.7清