高斯消元法 matlab

高斯消元法是一种用于求解线性方组的方法，它通过将方程组转化为上三角矩阵来简化求解过程。下面是一个使用Matlab实现高斯消元法的例子：

function x = gaussElimination(A, b)
    n = size(A, 1);
    Ab = [A, b];
    
    % 前向消元
    for k = 1:n-1
        for i = k+1:n
            factor = Ab(i, k) / Ab(k, k);
            Ab(i, k:n+1) = Ab(i, k:n+1) - factor * Ab(k, k:n+1);
        end
    end
    
    % 回代求解
    x = zeros(n, 1);
    x(n) = Ab(n, n+1) / Ab(n, n);
    for i = n-1:-1:1
        x(i) = (Ab(i, n+1) - Ab(i, i+1:n) * x(i+1:n)) / Ab(i, i);
    end
end

这个函数接受两个参数：系数矩阵A和常数向量b。它返回一个解向量x，满足Ax=b。