使用Java BigDecimal实现Sin函数遇单项式过大迭代无法终止问题求助

原代码的核心问题 #

你最初的代码出现“单项式始终为负、绝对值不断增大、循环无法终止”的问题，根源在分子项的更新逻辑错误：

x = x.abs().multiply(x.abs()).multiply(x).multiply(sign);

这里你错误地修改了输入参数x本身，而且每次都对x取绝对值，导致：

• 符号逻辑完全混乱：虽然乘了sign，但x.abs()会抵消符号，后续的分子项符号无法正确交替。
• 分子增长失控：每次迭代的分子是|x|³ * |x|（相当于x的4次方）再乘原x，实际应该是在前一项分子的基础上乘以-x²，而不是反复放大原x的幂次。这直接导致分子的增长速度远超过分母（阶乘）的增长，项的绝对值越来越大，完全违背泰勒级数收敛的特性，循环自然无法终止。

定位这类问题的实用方法 #

1. 输出关键变量日志：在迭代中打印每一步的dividend（分子）、divisor（分母）、num（当前项），对比泰勒级数的理论项：
   • 泰勒正弦级数的第k项（从0开始）是：(-1)^k * x^(2k+1)/(2k+1)!
   • 比如x=0.5时，前几项应该是0.5、-0.5³/6≈-0.0208、0.5^5/120≈0.00026，如果代码输出的项不符合这个规律，就能快速定位计算逻辑错误。
2. 手动验算前两项：拿小的x值（比如0.5）手动计算前两项，和代码结果对比，很快就能发现分子或分母的计算错误。

更优的实现思路（结合你的修正代码） #

你的修正代码已经解决了核心问题，这里再补充几个关键优化点：

1. 先做角度归约，加速收敛

泰勒级数在x接近0时收敛最快，所以先利用正弦的数学性质把输入归约到[0, π/2]区间：

• 周期性：sin(x) = sin(x mod 2π)
• 奇偶性：sin(x) = -sin(x-π) 当x在(π,2π)时
• 对称性：sin(x) = sin(π-x) 当x在(π/2, π)时
  这一步能大幅减少迭代次数，避免大角度下级数收敛极慢的问题。

2. 迭代更新分子分母，避免重复计算

不要每次重新计算x^(2k+1)和(2k+1)!，而是基于前一项递推：

• 下一项的分子 = 当前分子 * (-x²)
• 下一项的分母 = 当前分母 * (2k+2) * (2k+3)
  这种递推方式计算效率更高，也能减少大幂次、大阶乘计算带来的精度损失。

3. 精准控制迭代终止条件

用当前项的绝对值小于10^-(scale+cutOff)来终止迭代，确保最终结果的精度满足要求，同时避免不必要的迭代。

完整修正代码 #

package taylorSeries;
import java.math.BigDecimal;
import java.math.MathContext;
public class Sin {
    private static final int cutOff = 20;
    private static final BigDecimal PI = Pi.getPi(100);
    public static void main(String[] args) {
        System.out.println(getSin(Pi.getPi(100).multiply(new BigDecimal("1.5")), 100)); // Should be -1
    }
    public static BigDecimal getSin(final BigDecimal x, int scale) {
        // 角度归约到[0, π/2]区间，加速级数收敛
        if (x.compareTo(PI.multiply(new BigDecimal(2))) > 0)
            return getSin(x.remainder(PI.multiply(new BigDecimal(2)), new MathContext(x.precision())), scale);
        if (x.compareTo(PI) > 0)
            return getSin(x.subtract(PI), scale).multiply(new BigDecimal("-1"));
        if (x.compareTo(PI.divide(new BigDecimal(2))) > 0)
            return getSin(PI.subtract(x), scale);
        
        BigDecimal sign = new BigDecimal("-1");
        BigDecimal divisor = BigDecimal.ONE;
        BigDecimal i = BigDecimal.ONE;
        BigDecimal num = null;
        BigDecimal dividend = x;
        BigDecimal result = dividend;
        
        do {
            // 递推更新分子：前一项分子 * (-x²)
            dividend = dividend.multiply(x).multiply(x).multiply(sign);
            // 递推更新分母：前一项分母 * (i+1) * (i+2)
            i = i.add(BigDecimal.ONE);
            divisor = divisor.multiply(i);
            i = i.add(BigDecimal.ONE);
            divisor = divisor.multiply(i);
            
            num = dividend.divide(divisor, scale + cutOff, BigDecimal.ROUND_HALF_UP);
            result = result.add(num);
        } while(num.abs().compareTo(new BigDecimal("0.1").pow(scale + cutOff)) > 0);
        
        return result.setScale(scale, BigDecimal.ROUND_HALF_UP);
    }
}

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