C语言实现三元一次方程组求解：联合消元法详解

叶初

发布时间：2026-10-13| 3人浏览过

原创

本文演示如何利用C语言结合初中代数的“联合消元法”，在不引入复杂线性代数库的情况下，手动实现三元一次方程组的求解过程。

解题思路与目标

对于形如以下方程组的三元一次方程：

A₁x + B₁y + C₁z = D₁
A₂x + B₂y + C₂z = D₂
A₃x + B₃y + C₃z = D₃

我们的目标是通过编程模拟高斯消元的简化版，逐步消去变量 $z$ 和 $y$，先求出 $x$，再回代得到 $y$ 和 $z$。这种方法仅涉及基础的加减乘除运算，非常适合作为算法入门练习。

准备条件与环境设置

编译器环境：建议使用支持标准C99的编译器（如 Visual Studio 2010 及以上版本），需包含 stdio.h、stdlib.h 和 math.h 头文件。
数据类型选择：所有系数和解均使用 double 类型，以保证计算精度。
数据结构设计：
- A[3][4]：存储原始方程组的增广矩阵（前3列为系数，第4列为常数项）。
- B[2][4]：用于存储消去 $z$ 后得到的关于 $x, y$ 的二元一次方程组。
- C[1][4]：用于存储进一步消去 $y$ 后得到的只含 $x$ 的一元一次方程。
- x[3]：数组分别存储解得的 $x, y, z$ 值。

操作步骤与代码实现

第1步：初始化方程组数据

在程序中定义并输入待求解的三个方程的系数。这里以一组测试数据为例：

#include "stdafx.h"
#include "stdio.h"
#include "stdlib.h"
#include "math.h"

double x[3]; // 存放解 x, y, z
// 原始方程组 A:
// 1x + 2y + 3z = -42
// 2x - 1y + 5z = -49
// -1x + 3y - 3z = 17
double A[3][4] = {
	{1, 2, 3, -42},
	{2, -1, 5, -49},
	{-1, 3, -3, 17}
};

double B[2][4]; // 中间结果：消去z后的二元方程
double C[1][4]; // 最终结果：消去y后的一元方程

第2步：联合消元法——消除变量 z

为了从方程组中消去 $z$，我们需要将第一个方程分别与第二、第三个方程进行处理。核心策略是：调整两个方程使得它们对应变量的系数绝对值相等或互为相反数，然后相加或相减。

函数 CombineB(int i) 负责将 A[0] 行与 A[i] 行合并，生成新的二元方程存入 B 数组：

void CombineB(int i) {
    double m0 = abs(A[i][2]);        // 取第i个方程z的系数绝对值
    int b = A[0][2]*A[i][2]>0 ? -1 : 1; // 判断符号是否相同，决定加减号
    double mi = b * abs(A[0][2]);     // 取第一个方程z的系数调整值
    
    for(int j=0; j<4; j++) {
        // 执行乘法并相加，实现消元
        B[i-1][j] = A[0][j]*m0 + A[i][j]*mi;
    }
}

在主函数中调用该函数两次，分别处理第二个和第三个方程：

CombineB(1); // A[0]与A[1]消去z，存入B[0]
CombineB(2); // A[0]与A[2]消去z，存入B[1]

第3步：再次联合消元——消除变量 y

现在得到了两个关于 $x$ 和 $y$ 的方程（存放在 B 中）。按照同样逻辑，利用第一个二元方程（B[0]）与第二个二元方程（B[1]）消去变量 $y$，从而得到一个只包含 $x$ 的方程存入 C：

void CombineC() {
    double m0 = abs(B[1][1]);       // 取B[1]中y的系数绝对值
    int b = B[0][1]*B[1][1]>0 ? -1 : 1;
    double m1 = b * abs(B[0][1]);   // 调整系数
    
    for(int j=0; j<4; j++) {
        C[0][j] = B[0][j]*m0 + B[1][j]*m1;
    }
}

此时，C[0][0] 即为 $x$ 的系数，C[0][3] 为常数项。

第4步：回代求解 x, y, z

根据消元后的结果反向计算各个变量的值：

// 1. 求解 x
x[0] = (0 - C[0][3]) / C[0][0];

// 2. 代入B[0]求解 y
x[1] = (0 - B[0][3] - B[0][0]*x[0]) / B[0][1];

// 3. 代入A[0]求解 z
x[2] = (0 - A[0][3] - A[0][0]*x[0] - A[0][1]*x[1]) / A[0][2];

for(int i=0; i<3; i++)
    printf("x[%d]=%lf\t", i, x[i]);

完成判断与排错

如何验证程序运行正确：检查终端输出的变量值是否符合原方程组的解。例如上述代码设定的测试数据，解应为 x=-3.857..., y=-2.285..., z=...（视具体计算而定）。

常见注意事项：

除以零错误：如果某个需要消去的变量系数为0，当前硬编码逻辑会报错（除零）。在实际应用中，应先判断系数是否为0，若为0则无需消元；若非0但极小，需注意浮点数精度问题。
无解或无穷多解：本段代码假设方程组有唯一解。若消元后出现 0 = 非0 的情况，说明无解；若出现 0 = 0，说明可能有无穷多解或方程不独立。