π 是如何计算的？
π 是如何计算的？
π（圆周率）是通过数学公式来近似计算圆周长与直径之比得到的。由于 π 是一个无理数（无限不循环小数），其计算依赖于算法和无穷级数。以下是几种主要的古今计算方法：
1. 几何法（古代）
- 阿基米德（约公元前 250 年） 在圆的内外分别作正多边形。通过增加多边形的边数，他将 π 的范围限定在 3.1408 到 3.1429 之间。
- 祖冲之（公元 5 世纪） 使用 12288 边的多边形，得到 π ≈ 3.1415926（精确到小数点后 6 位）。
2. 无穷级数法（文艺复兴至近代）
这些公式将无穷多项相加，逐步逼近 π。
莱布尼茨级数（1674 年）：
收敛极慢（需要数百万项才能达到较好的精度）。
马钦类公式（1706 年）：
收敛速度更快；18–19 世纪曾用手算得到上百位 π 值。
拉马努金级数（1914 年）：
收敛极快（每项增加约 8 位精度）。早期计算机计算中曾使用此公式。
3. 现代计算机算法（快速收敛）
这些算法具有二次收敛性（每步迭代精度翻倍）。
高斯–勒让德算法（1975 年）：
仅需 25 次迭代，就能得到 π 至 4500 万位。
楚德诺夫斯基算法（1988 年）：
基于拉马努金的工作但更快。近年 π 计算世界纪录（如 2022 年的 100 万亿位）多采用此算法。
4. 蒙特卡洛方法（随机法）
在一个正方形内随机投点，并内切一个四分之一圆。落在圆内的点与总点数的比值近似于 (\pi/4)。
精度低、收敛慢，但概念简单易理解。
各种方法的精度对比
| 方法 | 可达精度（示例） |
|---|---|
| 多边形法（阿基米德） | 2–3 位 |
| 莱布尼茨级数 | ~5 位（100 万项） |
| 马钦类公式 | 上百位（手算） |
| 高斯–勒让德算法 | 25 步 → 4500 万位 |
| 楚德诺夫斯基算法 | 万亿位以上 |
如今，计算 π 更多用于测试超级计算机性能和数学研究，而非实际应用——工程上使用几十位 π 就已足够。