数学函数（math.h）：浮点数运算的十八般武艺

说到 math.h，这大概是 C 语言里最「实在」的头文件之一了。我入行那会儿，第一次做信号处理项目，满屏的 sin、cos、sqrt，那时候才真正体会到——数学库用得好，代码能省一半。

但说实话，很多人用数学函数只停留在「调 API」的层面。浮点数精度问题、边界情况、性能取舍……这些坑，我踩过不少。今天咱们就掰开揉碎，把 math.h 里最常用的函数讲透。

1. 绝对值与取整：fabs、ceil、floor、round

fabs 是求浮点数的绝对值。嗯，这个很简单，但有个细节：fabs(-0.0) 返回的是 0.0，不是 -0.0。IEEE 754 标准里，负零和正零在比较时是相等的，但如果你用 signbit() 检查，它们确实不同。我在做数值计算时遇到过，某些算法对符号敏感，这时候就得小心。

取整三兄弟：

  • ceil(x)：向上取整，返回 ≥ x 的最小整数。比如 ceil(3.14) = 4.0，ceil(-2.7) = -2.0。
  • floor(x)：向下取整，返回 ≤ x 的最大整数。floor(3.14) = 3.0，floor(-2.7) = -3.0。
  • round(x)：四舍五入到最近的整数。注意，round(2.5) = 3.0，round(3.5) = 4.0，这是「远离零」的舍入方式。

我个人习惯用 round 做常规舍入，但如果你需要「银行家舍入」（四舍六入五成双），那就得自己写逻辑了。我曾经在金融项目里被这个坑过——浮点数舍入方式不同，对账差了一分钱，查了半天。

2. 幂与根：pow、sqrt、exp、log

pow(x, y) 计算 x 的 y 次方。听起来简单，但要注意：pow(0, 0) 在某些实现里返回 1，有些返回 NaN。C99 标准规定返回 1，但老编译器可能不一样。我建议你写代码时主动判断一下。

sqrt(x) 求平方根。x 必须 ≥ 0，否则返回 NaN。嗯，这里有个性能小技巧：如果你只是求整数平方根，用 sqrt 没问题；但如果要开三次方，用 pow(x, 1.0/3.0) 就行。不过要注意，1.0/3.0 不是精确的 1/3，会有浮点误差。

我的经验： 在循环里频繁调用 pow 很慢。如果指数是整数，比如 x²，直接写 x * x 比 pow(x, 2) 快一个数量级。别问我怎么知道的——优化过一段图像处理代码，把 pow 全换成乘法，帧率直接翻倍。

exp(x) 计算 e^x，log(x) 计算自然对数（底数为 e）。log10(x) 是常用对数。这些函数在科学计算里很常见。我记得有一次做数据拟合，需要把指数关系转成线性，用 log 一变换，问题就简单了。

3. 三角函数：sin、cos、tan

这三个函数接受弧度值，不是角度。很多人刚学 C 时都犯过这个错——直接传角度进去，结果完全不对。转换公式很简单：弧度 = 角度 * M_PI / 180.0。注意 M_PI 在 math.h 里不一定定义，你可以自己写 #define M_PI 3.14159265358979323846。

另外，sin/cos/tan 的精度很高，但如果你传的参数非常大（比如 1e10），误差会累积。我建议做角度归一化，把值映射到 [0, 2π) 范围内再计算。

避坑指南： 我曾经在嵌入式项目里用 sin 做波形生成，发现输出有毛刺。查了半天，原来是浮点运算太慢，导致采样点间隔不均匀。后来改用查表法 + 线性插值，问题解决。所以，性能敏感的场景，别迷信数学库。

4. 取余与模运算：fmod

fmod(x, y) 返回 x 除以 y 的余数。注意，它和整数 % 不一样：fmod 的结果符号与 x 相同。比如 fmod(-5.3, 2.0) = -1.3。

如果你需要「数学模运算」（结果总是非负），可以用 fmod(x, y) + y 再取一次模。我在做周期信号处理时经常用这个技巧。

5. 浮点数精度问题：你永远算不准

这是个大话题。简单说，浮点数在计算机里是用二进制表示的，很多十进制小数（比如 0.1）无法精确表示。所以 0.1 + 0.2 != 0.3，这在 C 里是常态。

怎么办？我的建议是：

  • 不要直接用 == 比较浮点数。 用 fabs(a - b) < EPSILON，其中 EPSILON 取 1e-9 或根据场景调整。
  • 累加时注意误差累积。 比如求 10000 个 0.1 的和，结果不是 1000.0，而是 999.999... 可以用 Kahan 求和算法改善。
  • 输出时控制精度。 用 printf("%.10f", x) 而不是默认的 %f。

一个经典例子：

#include <stdio.h>
#include <math.h>

int main() {
    double a = 0.1;
    double b = 0.2;
    double c = 0.3;
    
    if (a + b == c) {
        printf("相等\n");
    } else {
        printf("不相等！差值是 %.20f\n", fabs((a + b) - c));
    }
    // 输出：不相等！差值是 0.00000000000000005551
    return 0;
}

你看，这就是浮点数的「原罪」。但别怕，理解了这个，你就能写出更健壮的代码。

知识体系总览

下面这张图总结了 math.h 里常用函数的分类和关系。我画的时候特意把「浮点数精度」放在中间，因为它贯穿所有运算。

math.h 数学函数 取整与绝对值 fabs, ceil, floor, round 幂与根 pow, sqrt, exp, log 三角函数 sin, cos, tan 取余运算 fmod 浮点数精度 比较、累加、输出 性能与替代 查表法、乘法替代 边界与异常 NaN, Inf, 负零 图：math.h 函数分类与核心关注点

总结一下

math.h 里的函数，说白了就是工具。用得好，事半功倍；用不好，bug 藏得深。我个人觉得，最重要的不是记住每个函数的参数，而是理解浮点数的本质——它不精确，但我们可以用技巧让它「足够精确」。

嗯，这一章就到这里。代码写多了，你会发现数学库其实挺有意思的。下次遇到 pow(0, 0) 或者 sin(1e10)，你至少知道该留个心眼了。

公众号:蓝海资料掘金营，微信deep3321

← 上一章 📖 返回目录 下一章 →