如何在JavaScript中突破精度限制计算更多位数的π

你的问题核心在于JavaScript原生Number类型是IEEE 754双精度浮点数，最多只能提供约16位有效数字，要突破这个限制必须使用任意精度算术，以下是两种可行方案：

1. 自定义高精度数运算逻辑 #

高斯-勒让德算法需要的核心操作是：加法、乘法、除法、平方根、平方运算。你可以通过字符串或数组存储高精度数字的每一位，手动实现这些运算逻辑：

• 用数组拆分存储数字的整数部分和小数部分，避免浮点数精度丢失；
• 基于竖式计算逻辑实现高精度加法、乘法，处理进位和位数对齐；
• 用牛顿迭代法实现高精度平方根，依赖已完成的高精度加减乘除操作。

这种方案灵活性强，但需要编写大量基础运算代码，适合需要完全自定义实现的场景。

2. 使用现成的高精度数学库 #

更高效的方式是使用纯JavaScript实现的高精度浮点数库，这类库已封装好所有需要的高精度运算：

• decimal.js：支持任意精度的浮点数运算，包含平方根、幂运算等高斯-勒让德算法所需操作；
• big.js：轻量级高精度浮点数库，API简洁，适合基础高精度计算场景。

基于decimal.js的代码改写示例 #

// 引入decimal.js（需先安装或引入库文件）
const Decimal = require('decimal.js');
// 设置精度，例如保留100位小数
Decimal.set({ precision: 100 });

let a = new Decimal(1);
let b = new Decimal(1).div(Decimal.sqrt(2));
let t = new Decimal(1).div(4);
let p = new Decimal(1);

let i = 0;
// 设置收敛阈值：当a和b的差小于1e-100时停止循环
while (a.sub(b).abs().greaterThan(new Decimal('1e-100'))) {
    i++;

    let an = a.add(b).div(2);
    b = Decimal.sqrt(a.mul(b));
    t = t.sub(p.mul(a.sub(an).pow(2)));
    p = p.mul(2);

    a = an;

    if (i === 100) break;
    let pi = a.add(b).pow(2).div(t.mul(4));
    console.log(pi.toString());
}

关键修正说明 #

你原代码中的循环条件a - b < 3是错误的，高斯-勒让德算法需要判断a和b的差足够小（趋近于0）来停止迭代，否则无法正确收敛到高精度的π值。

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