直接三角分解LU法

silence123

%直接三角分解法LU;
- g. n; `8 T! J# Q( G) n, yformat rational
disp('请输入矩阵的阶数');
' t$ B; o& X1 l! WN=input('N=');
3 p7 [3 c" e0 W  S8 `for i=1:N
    for j=1:N
        A(i,j)=input('A(i,j)=');
. v. v% B! W8 Q% D5 S8 A* y, y# r, A: M    end
8 U/ m5 y+ f) V. F  o9 j3 d6 t+ }end
L=eye(N);
for i=1:N
+ f7 _9 |1 f$ [    U(1,i)=A(1,i);
% G5 k  R6 ^) N# j# ]9 |4 N" S5 E    L(i,1)=A(i,1)/U(1,1);
end
for r=2:N
* |$ H. q4 n  j) k2 t    for i=r:N
        U(r,i)=A(r,i);
4 r1 }7 ~) r# a1 a        for k=1:r-1
& c% l+ E! A" c. j            U(r,i)=U(r,i)-L(r,k)*U(k,i);
        end
5 L$ `; D  b- t( b% ]; P        L(i,r)=A(i,r);
: k& Z4 Q) [& C        for k=1:r-1
0 @5 A3 _7 _  Y! u6 g* Y; G% F            L(i,r)=L(i,r)-L(i,k)*U(k,r);
        end
+ K" K$ s# X' Q+ r: R- p9 D7 u5 T) ]        L(i,r)=L(i,r)/U(r,r);
5 R7 D) H' e  V$ D! ]& F    end
+ ~$ g! S) X; J  C, B  i2 Q& i8 rend
disp('A=');
disp(A);
disp('L=');
8 j" T, u, r( g$ H. Xdisp(L);
disp('U=')
% T6 m. T7 S) d- \! [( v0 X& ^disp(U);

0.9清1.8清2.7清