直接三角分解LU法

silence123

%直接三角分解法LU;
+ [* x% }- R+ R; ]3 N3 Z1 ?format rational
disp('请输入矩阵的阶数');
% _# z9 h1 S) }5 o* r5 HN=input('N=');
$ @+ g! a/ }/ ^  M& S: A/ Jfor i=1:N
    for j=1:N
        A(i,j)=input('A(i,j)=');
7 h: m/ j0 A) @! u; e1 a! x1 C1 c    end
end
L=eye(N);
$ j- ^' R' L7 R$ w3 u9 p* Rfor i=1:N
    U(1,i)=A(1,i);
9 c' n, y0 _7 D) A$ ?) {* a$ z    L(i,1)=A(i,1)/U(1,1);
& Y, Z! y4 h1 @' V& g' uend
for r=2:N
    for i=r:N
        U(r,i)=A(r,i);
        for k=1:r-1
" A- q& \* D' b            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
- a1 s- `$ A! i& V3 m  g7 _, e" ?        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);
* K! p$ K0 K: c- ~$ M& c9 k8 l    end
end
disp('A=');
disp(A);
9 |$ o" x* X( s3 v2 y2 rdisp('L=');
/ e% F" R: t* u& bdisp(L);
disp('U=')
disp(U);

0.9清1.8清2.7清