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