为何C语言float类型读取浮点数值时出现异常舍入？

我在编写读取矩阵文件的程序，每行格式是行索引(int) 列索引(int) 数值(float)，但读取时发现数值出现异常舍入。比如文件里的行是：

1 1 881.41025641026

用以下代码读取：

int r, c;
float val;
for (int i=0; i<nz; i++) {
    // 其他逻辑
    fscanf(f, "%d %d %f\n", &r, &c, &val);
    // 其他逻辑
}

读取后得到的val是881.4102783203125，而非文件里的原值。试过%g等格式符，结果一样，但换成double类型存储时，数值就完全正确。

核心原因：混淆了取值范围与精度 #

你误解了float的特性：

• float的取值范围确实是1.175494351E-38至3.402823466E+38，但它的有效精度仅为6-7位十进制数字。这是因为IEEE 754单精度浮点数的尾数部分只有23位（加隐含的1位共24位），对应十进制的有效位数约为log10(2^24) ≈ 7.22位。
• 你文件中的数值881.41025641026有11位有效数字，远超过float的精度上限。float无法精确存储这么多位数的数值，只能将其舍入到最接近的可表示值，也就是你看到的881.4102783203125。
• 而double是双精度浮点数，尾数部分有52位（加隐含1位共53位），对应十进制有效位数约log10(2^53) ≈ 15.95位，完全能覆盖你这个数值的所有有效数字，所以读取后能保持原值。

解决方案 #

如果需要保留该数值的完整精度，直接改用double类型存储和读取：

int r, c;
double val;
for (int i=0; i<nz; i++) {
    // 其他逻辑
    fscanf(f, "%d %d %lf\n", &r, &c, &val);
    // 其他逻辑
}

内容的提问来源于stack exchange，提问作者lilith