高斯消去法MATLAB代码

以下是高斯消元法的MATLAB代码：

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

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

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

使用方法：

1.将上述代码保存为.m文件，例如gauss_elimination.m。

2.在MATLAB命令窗口中输入以下命令：

A = [1 2 3; 4 5 6; 7 8 9];
b = [1; 2; 3];
x = gauss_elimination(A, b);

其中A为系数矩阵，b为常数向量，x为方程组的解。