在R语言中log1p(x)与log(x+1)结果为何不完全一致？

我原本预期log1p(x)和log(x+1)会返回相同结果，但二者存在2.775558e-17的差值。为何这两个函数结果不一致？其中哪一个精度更高？

R语言测试代码及输出 #

> y <- log({32/123}+1)
> y
#> [1] 0.2312408

> x <- log1p(32/123)
> x
#> [1] 0.2312408

> x==y
#> [1] FALSE

> y-x
#> [1] 2.775558e-17
> sprintf("%a", y-x)
#> [1] "0x1p-55"

差异原因 #

• 浮点数运算的固有精度限制：计算机采用二进制浮点数存储数值，无法精确表示所有十进制小数，运算过程中必然产生微小舍入误差。log(x+1)需要先执行x+1的加法操作，这一步就可能引入精度损失；而log1p是专门为计算log(1+x)设计的优化函数，它跳过了显式计算x+1的步骤，直接针对该数学表达式的特性设计计算逻辑，从根源上减少了中间误差。

精度对比 #

log1p的精度更高，核心原因是：

• 当x的数值较小时，x+1的浮点数计算会丢失x的部分有效位。比如当x远小于1时，x的低位二进制位会被1的高位覆盖，导致x+1的结果无法准确反映x的真实值，后续取对数后会放大这个误差。而log1p的算法专门针对这种场景做了优化，能保留更多有效精度。
• 即使像测试中x≈0.26这样不算极小的数值，x+1的加法仍会产生微小的精度损失，log运算进一步放大了这个误差；而log1p直接计算log(1+x)的近似值，避免了中间步骤的误差累积，结果更接近真实数学值。

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