直接三角分解LU法

silence123

%直接三角分解法LU;
* B- E* U/ j* |* C( l3 R. X# b* `) D( |9 cformat rational
$ h: a  T, ^0 W( g$ J" C7 X4 h4 rdisp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
    for j=1:N
        A(i,j)=input('A(i,j)=');
) f& e% z9 G6 F8 g( W8 D    end
end
L=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
end
for r=2:N
% k1 R" t6 v3 s9 A5 F" {0 X    for i=r:N
        U(r,i)=A(r,i);
3 w7 i* ^2 _: |. t7 R$ `        for k=1:r-1
) g) u; e1 g$ k, Y            U(r,i)=U(r,i)-L(r,k)*U(k,i);
3 m0 Z9 l+ H* V        end
        L(i,r)=A(i,r);
$ Q+ ~6 Y8 ?% O- A; @        for k=1:r-1
% N& V, H, ~; W; d4 T            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
        L(i,r)=L(i,r)/U(r,r);
7 c/ l% O4 O- V1 Y7 j9 N7 k    end
end
; R0 }  Y) B; Mdisp('A=');
disp(A);
disp('L=');
disp(L);
$ U2 v2 u3 U4 zdisp('U=')
! N/ }  H1 x5 e6 n2 d$ S' Bdisp(U);

0.9清1.8清2.7清