直接三角分解LU法

silence123

%直接三角分解法LU;
9 V9 ]4 r+ j3 X) W: d! [0 ^format rational
disp('请输入矩阵的阶数');
N=input('N=');
4 K" R  c$ |7 zfor i=1:N
    for j=1:N
; {" m4 s7 {# b6 v" E5 a# B        A(i,j)=input('A(i,j)=');
/ J7 [, C- g; c& ?, E    end
1 ^: A! e- z3 F* P( f2 z+ oend
; g1 r3 G; M5 ?+ \$ rL=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
end
& D8 o3 J, u! L/ O# H5 Rfor r=2:N
7 m% c( s: F8 ^% `1 p4 y/ b3 @    for i=r:N
        U(r,i)=A(r,i);
( Y* m# O2 k8 V0 p" B# M- C: A        for k=1:r-1
+ \( V- W; U) c; H  b' f            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
        L(i,r)=A(i,r);
/ s" `( f5 }1 h( F; W( |! Z        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
! \, x9 n8 m; G! @- u        L(i,r)=L(i,r)/U(r,r);
    end
end
disp('A=');
disp(A);
disp('L=');
$ C1 a: `( M( W. K6 M% b4 |disp(L);
' F0 w7 R# O8 ~8 ^$ \disp('U=')
disp(U);

0.9清1.8清2.7清