算法概述

高斯消元法是线性代数中用于求解线性方程组的经典算法，能够高效处理包含数千个方程和未知数的系统。该方法通过初等行变换将增广矩阵转化为上三角矩阵，再通过回代过程求解未知数。

数学基础

当方程组包含n个未知数时，必须存在n个线性无关的方程才能确定唯一解。若存在冗余方程（如0x + 0y = 0）或矛盾方程（如0x = 5），则系统无解或有无穷多解。增广矩阵将系数矩阵与常数项合并为n×(n+1)结构，例如方程组：

3x + 2y = 5
x - y = 1

对应的增广矩阵为：

[3  2 | 5]
[1 -1 | 1]

核心步骤

1. 主元选择：在当前列的剩余行中选取绝对值最大的元素作为主元，通过行交换提升数值稳定性
2. 归一化：将主元所在行除以主元值，使主元位置变为1
3. 消元：用当前行消去下方所有行的当前列元素，使下方元素归零
4. 回代求解：从最后一行开始依次计算各未知数的值

代码实现

#include <iostream>
#include <cmath>
#include <algorithm>
using namespace std;

const int MAX_SIZE = 1010;
double coefficientMatrix[MAX_SIZE][MAX_SIZE]; // 增广矩阵

int main() {
    int equationCount;
    cin >> equationCount;
    for (int i = 1; i <= equationCount; i++) {
        for (int j = 1; j <= equationCount + 1; j++) {
            cin >> coefficientMatrix[i][j];
        }
    }

    for (int currentRow = 1; currentRow <= equationCount; currentRow++) {
        // 寻找主元行
        int pivotRow = currentRow;
        double maxValue = fabs(coefficientMatrix[currentRow][currentRow]);
        for (int row = currentRow + 1; row <= equationCount; row++) {
            if (fabs(coefficientMatrix[row][currentRow]) > maxValue) {
                maxValue = fabs(coefficientMatrix[row][currentRow]);
                pivotRow = row;
            }
        }
        if (maxValue < 1e-8) {
            cout << "No Solution" << endl;
            return 0;
        }
        if (pivotRow != currentRow) {
            swap(coefficientMatrix[currentRow], coefficientMatrix[pivotRow]);
        }

        // 归一化当前行
        double pivot = coefficientMatrix[currentRow][currentRow];
        for (int col = currentRow; col <= equationCount + 1; col++) {
            coefficientMatrix[currentRow][col] /= pivot;
        }

        // 消去下方行的当前列元素
        for (int row = currentRow + 1; row <= equationCount; row++) {
            double factor = coefficientMatrix[row][currentRow];
            for (int col = currentRow; col <= equationCount + 1; col++) {
                coefficientMatrix[row][col] -= factor * coefficientMatrix[currentRow][col];
            }
        }
    }

    // 回代过程
    for (int i = equationCount; i >= 1; i--) {
        for (int j = i + 1; j <= equationCount; j++) {
            coefficientMatrix[i][equationCount + 1] -= coefficientMatrix[i][j] * coefficientMatrix[j][equationCount + 1];
        }
    }

    // 输出结果
    for (int i = 1; i <= equationCount; i++) {
        printf("%.2f\n", coefficientMatrix[i][equationCount + 1]);
    }
    return 0;
}

该实现采用列主元法提升数值稳定性，通过动态调整主元位置有效控制计算误差。实际应用中需注意浮点数精度问题，通常使用1e-8作为阈值判断主元有效性。