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 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