自定义对数函数与Math.log结果差异及NaN问题排查

常见问题根源 #

• 泰勒级数收敛域限制：ln(1+x)的泰勒展开仅在-1 < x ≤ 1区间收敛，直接用它计算大于1的整数n会导致级数发散，结果完全偏离真实值，甚至出现计算溢出或NaN。比如输入n>2时，变量超出收敛域，迭代会产生无意义的数值。
• 定义域处理缺失：对数函数的定义域是n > 0，如果输入n=0或者负数，会直接返回NaN，必须先对输入做合法性校验。
• 底转换错误：以2为底的对数lg(n)应该是ln(n)/ln(2)，如果exercise方法没做这个转换，或转换时用错底数，结果会和内置Math.log()（自然对数）差异巨大。
• 幂运算的精度与溢出：循环实现幂运算时，高次幂可能出现浮点数精度丢失，甚至整数溢出（如果用int存中间结果），导致迭代项完全错误。

修正方案 #

1. 调整泰勒级数的适用范围 #

先把n转换到收敛域内：利用对数性质ln(n) = ln(2^k * m) = k*ln(2) + ln(m)，其中0.5 < m ≤ 1，再通过ln(m) = -ln(1/m)将其转换为符合收敛域的形式（此时(1/m)-1 ∈ (0,1]）。

2. 输入合法性校验 #

在函数开头添加校验逻辑：

if (n <= 0) {
    throw new IllegalArgumentException("n must be positive");
}

3. 正确实现底转换 #

计算以2为底的对数时，需基于自然对数做转换：

public static double lg(int n) {
    double lnN = computeNaturalLog(n);
    return lnN / Math.log(2);
}

4. 优化泰勒迭代的幂运算 #

用累乘代替循环幂运算，减少精度丢失：

double term = x; // 第一项x
double sum = term;
double sign = -1.0;
for (int k = 2; k <= 20; k++) {
    term *= x; // 每次乘x得到x^k，避免重复计算
    sum += sign * term / k;
    sign *= -1;
}

示例修正后的核心代码片段 #

public static double computeNaturalLog(int n) {
    if (n <= 0) {
        throw new IllegalArgumentException("n must be positive");
    }
    if (n == 1) {
        return 0.0;
    }
    // 把n转换到(0.5, 1]区间，记录2的幂次k
    int k = 0;
    double m = n;
    while (m > 1) {
        m /= 2;
        k++;
    }
    // 转换为符合收敛域的x值
    double x = 1/m - 1;
    double sum = 0.0;
    double term = x;
    double sign = 1.0;
    for (int i = 1; i <= 20; i++) {
        sum += sign * term / i;
        term *= x;
        sign *= -1;
    }
    double lnM = -sum;
    // 还原得到ln(n)
    return k * Math.log(2) + lnM;
}

public static double lg(int n) {
    return computeNaturalLog(n) / Math.log(2);
}

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