直接三角分解LU法

silence123

%直接三角分解法LU;
) H* b5 E4 \5 aformat rational
disp('请输入矩阵的阶数');
N=input('N=');
for i=1:N
' i" a7 w. h; J- K    for j=1:N
        A(i,j)=input('A(i,j)=');
  y; N3 X& N2 W* Q+ }/ z    end
7 ?, J) F6 Z9 _% n4 bend
: D8 f. F5 h4 D  K! y2 s" nL=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
    for i=r:N
  B7 s$ M& ^) s0 Q        U(r,i)=A(r,i);
0 ~' y% U; D4 U" C, J        for k=1:r-1
            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
        L(i,r)=A(i,r);
        for k=1:r-1
            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
        L(i,r)=L(i,r)/U(r,r);
  v# p4 {1 V/ r  K8 i    end
end
disp('A=');
disp(A);
3 l  g* T) j1 b+ O; edisp('L=');
disp(L);
. d; H4 c! ?! m8 i0 Zdisp('U=')
* K/ N  Z  A9 fdisp(U);

0.9清1.8清2.7清