直接三角分解LU法

silence123

%直接三角分解法LU;
1 C% O4 b& l# j) Kformat rational
$ p; S0 N1 L; U& C5 E( Ndisp('请输入矩阵的阶数');
% k2 f* i9 T$ \* f! O5 sN=input('N=');
for i=1:N
    for j=1:N
) ^3 \: \4 r9 y5 _* i# Z' w) o6 M        A(i,j)=input('A(i,j)=');
    end
end
+ I: ]$ f9 d7 V! u% A4 I' CL=eye(N);
for i=1:N
' Z! E7 d& j3 K% N% M% j    U(1,i)=A(1,i);
* \  @2 [' b# J4 s9 x    L(i,1)=A(i,1)/U(1,1);
end
for r=2:N
& \1 ]% ^+ N. e' g0 q! J4 S    for i=r:N
        U(r,i)=A(r,i);
; `* m! ~3 {1 ?        for k=1:r-1
! g0 w1 r/ K$ a! b* J            U(r,i)=U(r,i)-L(r,k)*U(k,i);
( N  X/ }  |7 q0 G        end
+ R8 j( A% w. k2 n9 ]        L(i,r)=A(i,r);
3 U' n' w# }7 e) Y% z: h        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
  [" H% S2 L& R+ f6 q        L(i,r)=L(i,r)/U(r,r);
: H0 a' U. m! y- Z    end
& Y0 e4 K: n$ [5 M; N6 Q+ ]8 C) Fend
disp('A=');
: J4 J9 a3 H  `; odisp(A);
disp('L=');
disp(L);
; k6 {& S8 @7 b& Tdisp('U=')
disp(U);

0.9清1.8清2.7清