直接三角分解LU法

silence123

%直接三角分解法LU;
$ q0 k9 X& B! T' A7 H' kformat rational
disp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
    for j=1:N
, t9 r) _+ v" W2 W        A(i,j)=input('A(i,j)=');
    end
: P' E* Q6 m4 lend
L=eye(N);
for i=1:N
    U(1,i)=A(1,i);
* ?9 p" v* T* T1 R  |    L(i,1)=A(i,1)/U(1,1);
2 \' }) \8 L* d6 q7 b" y2 Y# k9 bend
  \6 B! D* M- i1 P) jfor r=2:N
6 ^; J+ G1 z3 ]' w) Q    for i=r:N
- W0 H5 |0 k, S- m! ^! c7 S1 A        U(r,i)=A(r,i);
3 U# _1 _8 G* U/ k2 ^/ ^: A! W' G        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
$ e; I. g- L3 ?% ?2 {# j" E        L(i,r)=A(i,r);
        for k=1:r-1
% d7 ?% U1 c6 G& X            L(i,r)=L(i,r)-L(i,k)*U(k,r);
. A. K6 u0 M* [5 M, V+ L# o        end
8 P. C4 t8 x- c        L(i,r)=L(i,r)/U(r,r);
7 N7 `) @: ]1 k& k  f, i" e    end
end
/ y4 Z( h0 M' s1 q6 }8 \disp('A=');
disp(A);
disp('L=');
disp(L);
disp('U=')
, C+ q- g9 t' ~( qdisp(U);

0.9清1.8清2.7清