高斯消元法 matlab实现

高斯消元法是一种求解线性方程组的方法，可以用 MATLAB 实现。以下是 MATLAB 代码：

function x = gauss_elimination(A, b)
% 高斯消元法求解线性方程组 Ax=b
% 输入参数：系数矩阵 A，常数向量 b
% 输出参数：解向量 x

n = size(A, 1);
Ab = [A, b]; % 增广矩阵

% 前向消元
for k = 1:n-1
    for i = k+1:n
        Ab(i,:) = Ab(i,:) - Ab(k,:) * Ab(i,k) / Ab(k,k);
    end
end

% 回带求解
x = zeros(n, 1);
x(n) = Ab(n,n+1) / Ab(n,n);
for k = n-1:-1:1
    x(k) = (Ab(k,n+1) - Ab(k,k+1:n) * x(k+1:n)) / Ab(k,k);
end
end