C++高斯约当消去法求解线性方程组返回极小非0值问题咨询

1. 极小非0值是否由C++浮点数精度特性导致 #

是。单精度float类型遵循IEEE 754标准，仅能提供6_{7位有效十进制精度，高斯约当消去法包含大量连续的乘除、加减运算，每一步运算都会产生微小的舍入误差，多步累计后就会在本应输出0的位置出现1e-7}1e-8量级的极小值，你遇到的2.3842e-08完全属于单精度浮点数的典型误差范围，就是精度特性导致的正常现象。

2. 如何在不丢失有效数据的前提下将这类极小值舍入为0 #

你可以先定义一个和精度匹配的误差阈值（epsilon），单精度场景下一般取1e-6即可，该阈值远大于误差量级、又远小于正常有效数值的量级，不会误覆盖有效数据。根据你的使用场景可以选两种处理方式：

• 仅输出时处理：不修改原矩阵的计算数据，只在打印环节做判断，不会影响后续使用该矩阵做其他计算的精度
• 计算后全局处理：如果后续要用到消元结果，可以在消元完成后遍历所有矩阵元素，将绝对值小于阈值的数值直接置为0

参考修改后的打印函数代码（需要额外包含<cmath>头文件调用fabs）：

#include <cmath>
const float EPS = 1e-6; // 单精度浮点数误差阈值

void print(float x[][N], int n){
    for(int i=0;i<n;i++){
        for(int j=0;j<=n;j++) {
            float val = x[i][j];
            // 绝对值小于阈值的数判定为0
            if(fabs(val) < EPS) {
                val = 0.0f;
            }
            cout << setprecision(5) << setw(15) << val;
        }
        cout << endl;
    }
    cout << endl;
}

3. 出现-0输出的原因是什么 #

IEEE 754浮点数标准中明确定义了负零：符号位为1、数值位全为0的浮点数就是负零。你代码中出现-0是消元过程中出现了0减去一个极小正数、极小负数乘负数这类运算，运算结果的数值位为0、符号位保留为负，输出时就会显示为-0。

如果要消除-0的显示，只需要在上述极小值判断逻辑中加一步处理，所有判定为0的数值统一赋值为正0即可，参考代码中的val = 0.0f就已经可以解决这个问题。

内容的提问来源于stack exchange，提问作者Circuit_Breaker0.7