C++ 代码 如何实现 二元一次方程组组求解

在C++中，我们可以使用高斯消元法（Gaussian Elimination）或者LU分解、QR分解等数值计算库来求解二元一次方程组。这里简单展示一下基于增广矩阵（包含系数矩阵和常数项）的高斯消元基本步骤：

#include <iostream>
#include <vector>

// 高斯消元函数
std::vector<double> solveLinearEquation(const std::vector<std::vector<double>>& augmentedMatrix) {
    int n = augmentedMatrix.size() - 1; // 行数
    for (int i = 0; i < n; ++i) {
        if (augmentedMatrix[i][i] == 0) { // 检查主元素是否为零
            std::cout << "Error: Division by zero in row " << i + 1 << std::endl;
            return {}; // 或者返回一个错误处理的向量
        }
        
        // 对于当前行i，除以其主元素
        for (int j = i; j <= n; ++j)
            augmentedMatrix[i][j] /= augmentedMatrix[i][i];

        // 将行i应用于下方的行
        for (int k = i + 1; k <= n; ++k) {
            double multiplier = augmentedMatrix[k][i];
            for (int l = i; l <= n; ++l)
                augmentedMatrix[k][l] -= multiplier * augmentedMatrix[i][l];
        }
    }

    // 解出最后一个非零元素对应的x和y
    double x = augmentedMatrix[n][n];
    double y = augmentedMatrix[n][n - 1] / augmentedMatrix[n][n]; // 注意这里可能会有除以零的风险

    // 返回结果作为向量 [x, y]
    std::vector<double> solution({x, y});
    return solution;
}

int main() {
    std::vector<std::vector<double>> augmentedMatrix = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}};
    std::vector<double> result = solveLinearEquation(augmentedMatrix);
    if (!result.empty()) {
        std::cout << "Solution is: x = " << result[0] << ", y = " << result[1] << std::endl;
    }
    return 0;
}

这个示例展示了如何使用高斯消元法求解基础的三元一次方程组。如果是二元一次方程组，可以直接去除最后一列并相应地调整行数。注意，实际项目中建议使用数值计算库如Eigen或Boost.Multiprecision来进行更稳定的计算。