MATLAB 列主元高斯消元法求解线性方程组

列主元高斯消元法 (Column pivoting Gaussian elimination) 是一种改进的高斯消元法，用于求解线性方程组。它通过在每一步消元过程中选择绝对值最大的主元来提高算法的稳定性。

下面是使用 MATLAB 实现列主元高斯消元法的代码示例：

function [x] = columnPivotingGaussian(A, b)
    % 获取方程组的维度
    n = size(A, 1);
    
    % 扩展系数矩阵
    Ab = [A, b];
    
    % 初始化列主元数组
    p = 1:n;
    
    % 遍历每一列
    for k = 1:n-1
        % 选择主元
        [~, maxIndex] = max(abs(Ab(k:n, k)));
        maxIndex = maxIndex + k - 1;
        
        % 交换行
        Ab([k, maxIndex], :) = Ab([maxIndex, k], :);
        p([k, maxIndex]) = p([maxIndex, k]);
        
        % 消元
        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);
    for k = n:-1:1
        x(k) = (Ab(k, n+1) - Ab(k, k+1:n) * x(k+1:n)) / Ab(k, k);
    end
    
    % 根据主元交换结果对 x 进行重新排序
    x(p) = x;
end

使用时，将系数矩阵 A 和常数向量 b 传递给 columnPivotingGaussian 函数，该函数将返回线性方程组的解向量 x。请注意，此代码假设输入的线性方程组具有唯一解。

以下是一个示例，演示如何使用列主元高斯消元法求解线性方程组：

% 定义系数矩阵 A 和常数向量 b
A = [1, -2, 1; 2, 1, -1; -3, 1, 2];
b = [3; 4; 1];

% 调用列主元高斯消元法求解
x = columnPivotingGaussian(A, b);

% 输出结果
disp('线性方程组的解：');
disp(x);

请注意，如果输入的线性方程组无解或有无穷多解，则该方法可能会失败。在实际应用中，建议在使用列主元高斯消元法之前，首先进行矩阵的可逆性和解的存在性检查。