LU 分解和前向/后向代入法来解线性方程组 Ax = b

2744557306

这段代码实现了 LU 分解和前向/后向代入法来解线性方程组 Ax = b，其中 A 是一个方阵，b 是一个列向量。以下是代码的中文解释：
+ {, L0 N- y5 Q9 P6 `a = [1,2,3;2,5,2;3,1,5];
b = [14,18,20];
# [' K" v2 o* f6 rn = size(a, 1);
) Y( m2 K$ z. @# _7 ]u = zeros(size(a));
l = zeros(size(a));
  K$ Z( N( [& \! l) \4 wu(1, = a(1, ;
2 Q& u2 ^( b. ]* ]# F8 _8 N8 S! ]
% LU 分解
6 v& z1 Y+ w' J" [( z% f4 qfor i = 2:n
3 [8 L! y+ F  P& n9 m# U6 C: i    l(i, 1) = a(i, 1) / u(1, 1);
end
for r = 2:n
    for i = r:1:n
; z* N( c2 g. K, [        sum1 = 0;
; J4 r7 z) A: ?        for k = 1:r-1
- G4 Z0 \- K. j3 Y6 e1 Z) c$ V7 O* e7 o            sum1 = sum1 + l(r, k) * u(k, i);
        end
6 R/ g1 m3 C" c0 q% H        u(r, i) = a(r, i) - sum1;
    end
: h& E( U" x6 Q. }/ G- G    for i = r+1:n
- Q& W! Y) K( L        sum2 = 0;
/ W/ e; O) i& u0 Y        for k = 1:r-1
            sum2 = sum2 + l(i, k) * u(k, r);
7 O4 C: Y& U  M- H* @        end
        l(i, r) = (a(i, r) - sum2) / u(r, r);
/ p6 k  I& J( }    end
+ y/ e# X! @$ @8 q0 }+ y) D& F* z( c0 @end
% A7 O, B4 W' |. F; K3 ?
% Z* b* l0 _9 ]( P( c4 [% 设置 L 的对角线为1
  z) U9 ~& |: C% h2 B0 y( Kfor i = 1:n
$ C7 S# `) k$ ~6 j: A  `    l(i, i) = 1;
end
  x0 {& m1 @+ C- Z
% 前向代入
) B9 F& m) b3 f3 [/ j& b% p2 sy(1) = b(1);
! d8 b/ j$ F% D7 L, Qfor i = 2:n
    sum3 = 0;
) {, _! e( u. Z3 V    for k = 1:i-1
  C* Q) a. a, j+ \" K        sum3 = sum3 + l(i, k) * y(k);
    end
7 v7 |  W- G6 l. a0 p& {6 u2 r0 w    y(i) = b(i) - sum3;
end

' F" m" n. l* ^9 n3 U, }4 g% 后向代入
$ K4 I0 K& m' ix(n) = y(n) / u(n, n);
for i = n-1:-1:1
2 g5 u1 O$ h0 X$ A    sum4 = 0;
    for k = i+1:n
6 v/ Z6 A6 R2 Y% ^' j% G( P        sum4 = sum4 + u(i, k) * x(k);
    end
    x(i) = (y(i) - sum4) / u(i, i);
end
9 g1 S: Y1 Z% ^" v; H
% 输出结果
disp('解 y：');
, R, U# D- J) e4 [/ E( c* @disp(y');
disp('解 x：');
disp(x');
3 E( z; c; F( K. ]/ }
这段代码通过 LU 分解将矩阵 A 分解为下三角矩阵 L 和上三角矩阵 U，然后使用前向代入和后向代入求解线性方程组 Ax = b。最后，输出解 y 和解 x。