直接三角分解LU法

silence123

%直接三角分解法LU;
2 P2 q) G/ M0 J# @; y, r7 R) Bformat rational
disp('请输入矩阵的阶数');
N=input('N=');
3 Z7 _# ?; m2 K% x9 M$ cfor i=1:N
8 V# l0 T0 ?$ E' x    for j=1:N
        A(i,j)=input('A(i,j)=');
) F% F7 ^& o0 ?$ N6 d    end
3 T# `: a, n6 y. C: l" }% Hend
8 ~8 }1 y# c2 z3 gL=eye(N);
for i=1:N
    U(1,i)=A(1,i);
  b5 p# Z5 ~3 {7 _" Z    L(i,1)=A(i,1)/U(1,1);
0 O! u# q1 G* X, H5 m2 ~! p- G, Mend
5 L9 F0 ]+ r/ m* ^& kfor r=2:N
    for i=r:N
+ A3 d2 j0 A, t* e5 m        U(r,i)=A(r,i);
, S9 I: E1 S! u  o$ b# ?( `        for k=1:r-1
; k, d  l. o" s4 B            U(r,i)=U(r,i)-L(r,k)*U(k,i);
8 e6 }' e( i9 z( Z, k. ^5 H7 X- |        end
- _  B+ z2 ~5 S8 n0 L        L(i,r)=A(i,r);
        for k=1:r-1
) \0 N$ L' L4 P2 X  f            L(i,r)=L(i,r)-L(i,k)*U(k,r);
0 \1 ], n( v! ]6 M        end
        L(i,r)=L(i,r)/U(r,r);
) g3 v3 t6 N6 Q$ c$ K/ Q    end
8 J3 S& h: Q" r- D3 Yend
4 C- F# `( `5 L/ E7 Edisp('A=');
, K5 s8 r; |2 E2 K5 g! B; M3 ~& Odisp(A);
disp('L=');
disp(L);
( Y4 v/ S' V8 Z: B- m- @# C- Y0 Kdisp('U=')
1 U8 }% M9 H. K2 }, X2 \8 a# }  i6 ]disp(U);

0.9清1.8清2.7清