直接三角分解LU法

silence123

%直接三角分解法LU;
. r/ x/ n( R9 @  O$ _- Fformat rational
disp('请输入矩阵的阶数');
N=input('N=');
" k$ Q/ @, J& U; yfor i=1:N
9 X1 m" O+ t( n. X/ J/ u    for j=1:N
        A(i,j)=input('A(i,j)=');
" x9 U: q/ E* Y- _    end
end
5 |  C$ d9 @' _3 pL=eye(N);
for i=1:N
    U(1,i)=A(1,i);
    L(i,1)=A(i,1)/U(1,1);
: G" M1 R' T' i% o( \: a- ]+ Gend
for r=2:N
    for i=r:N
8 B; d( B4 Q8 O7 A- K8 S  Q. _        U(r,i)=A(r,i);
  G6 N$ K# L) |4 }        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
$ H6 i$ e/ I/ `/ ^        L(i,r)=A(i,r);
4 p% [" [' a& y: o$ T: y  }1 H2 h2 n        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
2 B% D1 d$ R( v' ?: I* V4 n        end
& K/ \' _" l5 |$ \        L(i,r)=L(i,r)/U(r,r);
    end
7 K+ a1 e7 I! `0 K, O9 s' V% o- Yend
& Z# ^' [. T: _disp('A=');
disp(A);
" C5 f# I  E( h; ]" v8 U/ [disp('L=');
disp(L);
disp('U=')
0 Y, @5 x: S  q0 ]# pdisp(U);

0.9清1.8清2.7清