2、线性系统求解方法：直接法与迭代法

线性系统求解方法：直接法与迭代法

1. 线性系统直接求解方法

1.1 求解下三角系统

在求解线性系统时，下三角系统的求解是一个基础问题。以下是求解下三角系统的 MATLAB 代码：

n = length(b);
x = zeros(n, 1);
x(1) = b(1)/A(1, 1);
for j = 2 : n
    x(j) = b(j);
    for k = 1 : j-1
        x(j) = x(j) - A(j, k)*x(k);
    end
    x(j) = x(j)/A(j, j);
end

对应的 Python 代码如下：

def SolveLowerLinearSystem(A, b):
    import numpy as np
    n = len(b)
    x = np.zeros((n, 1), 'float')
    x[0] = b[0]/A[0, 0]
    for i in range(1, n):
        x[i] = b[i]
        for j in range(i):
            x[i] -= A[i, j]*x[j]
        x[i] /= A[i, i]
    return x

1.2 高斯消元法与高斯 - 约旦消元法

高斯消元法和高斯 - 约旦消元法都通过对矩阵进行初等行变换来求解线性系统。高斯消元法在得到矩阵的行