matlab编写下三角形方程组前推法的函数

下面是使用 MATLAB 编写的下三角形方程组前推法的函数：

function x = forwardSubstitution(L, b)
    % 获取方程组的维度
    n = size(L, 1);
    
    % 初始化解向量
    x = zeros(n, 1);
    
    % 前推法求解下三角形方程组
    for i = 1:n
        x(i) = b(i);
        for j = 1:i-1
            x(i) = x(i) - L(i, j) * x(j);
        end
        x(i) = x(i) / L(i, i);
    end
end

该函数接受两个输入参数：下三角形矩阵 L 和向量 b。函数使用前推法，通过求解下三角形方程组 Lx = b，得到解向量 x。

注意：这里假设输入的下三角形矩阵 L 是非奇异的，并且方程组有唯一解。