高精度计算大x值下的-atanh(1/(1+10**-x))问题求助

当x取值量级达到240时，直接计算1/(1+10**-x)会因浮点数精度限制被近似为1.0，而atanh(1)在数学上是无穷大，这就是你遇到inf或域错误的核心原因。我们可以通过等价公式转换避开这个问题，利用反双曲正切的对数定义推导近似表达式：

公式推导 #

反双曲正切的定义为：

atanh(y) = 0.5 * ln((1+y)/(1-y))

代入y = 1/(1+ε)（其中ε=10**-x，x很大时ε极小）：

(1+y)/(1-y) = [1 + 1/(1+ε)] / [1 - 1/(1+ε)] = (2+ε)/ε ≈ 2/ε = 2*10^x

因此：

atanh(1/(1+ε)) ≈ 0.5 * ln(2*10^x) = 0.5*(ln2 + x*ln10)

最终目标表达式的近似值为：

-atanh(1/(1+10**-x)) ≈ -0.5*(math.log(2) + x*math.log(10))

代码实现 #

用Python直接计算这个近似式即可得到高精度结果，无需依赖特殊库：

import math

def compute_target_value(x):
    return -0.5 * (math.log(2) + x * math.log(10))

# 测试x=240
print(compute_target_value(240))  # 输出约-276.65675，与预期的-276.65一致

精度说明 #

x越大，ε=10**-x越趋近于0，这个近似的误差就越小。对于x=240，误差小到可以忽略，完全满足精度需求。

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