几何级数与π

问题提出

“各位同学早上好！“我走进教室，看到呼噜星人们依旧保持着他们那独特的三只眼睛和毛茸茸的外表，但今天我能感觉到他们脸上多了几分期待的笑容。

“今天我们要讨论一个有趣的话题：几何级数能否用于计算π？“我微笑着说，“π，这个我们熟悉的圆周率，它神秘而又无处不在。让我们看看几何级数这个工具能否帮助我们揭开它的面纱。”

呼噜星人们面面相觑，其中一只蓝毛色的学生举起了手：“老师，什么是几何级数呢？我听说过算术级数，但是几何级数是什么？”

“好问题！“我赞许地点点头，“算术级数是每一项与前一项的差相等，而几何级数是每一项与前一项的比相等。让我们从基础开始，一步步探索。“

观察与猜想

“让我们先来看看一个简单的例子，“我在黑板上写下了数列：1, 1/2, 1/4, 1/8, 1/16, …

“大家观察一下，这个数列有什么特点？“我问道。

“每一项都是前一项的一半！“绿毛色的学生兴奋地回答。

“完全正确！“我赞赏地说，“这就是几何级数，公比r = 1/2。现在让我们思考一个问题：如果这个数列无限延续下去，它的和会是多少？”

呼噜星人们开始小声讨论。我引导道：“让我们看看前n项的和：”

$S_n = 1 + \frac{1}{2} + \frac{1}{4} + \frac{1}{8} + \cdots + \frac{1}{2^{n-1}}$

“现在我们来看看当n趋近于无穷大时会发生什么，“我继续分析，“当|r| < 1时，$r^n$会趋近于0，所以：”

$S = \lim_{n \to \infty} S_n = \lim_{n \to \infty} a \cdot \frac{1-r^n}{1-r} = a \cdot \frac{1}{1-r} = \frac{a}{1-r}$

“太棒了！“紫毛色的学生激动地说，“所以我们得到了一个强大的公式：$\sum ar^n = \frac{a}{1-r}$（|r| < 1）！”

“现在让我们看看如何用这个公式来计算π，“我转换了话题，“π的值大约是3.1415926535…，我们需要找到一种方法来表示它。“

严格证明

“首先，让我们证明这个几何级数求和公式，“我在黑板上开始推导，“设：”

$S = a + ar + ar^2 + ar^3 + \cdots$

“我们将两边都乘以r：”

$rS = ar + ar^2 + ar^3 + ar^4 + \cdots$

“然后用第一个式子减去第二个式子：”

$S - rS = a + (ar - ar) + (ar^2 - ar^2) + (ar^3 - ar^3) + \cdots$

“除了第一项，其他所有项都抵消了：”

$S - rS = a$

“提取S：”

$S(1 - r) = a$

“最后得到：”

$S = \frac{a}{1-r}$

“这个证明成立的前提是什么？“我问道。

“哦！是|r| < 1！“红毛色的学生抢答道，“如果|r| ≥ 1的话，级数不会收敛！”

“完全正确！“我赞赏地说，“现在让我们看看如何用几何级数来逼近π。”

“让我们从tan的泰勒展开开始，“我解释道，“我们知道：”

$\tan^{-1}(x) = x - \frac{x^3}{3} + \frac{x^5}{5} - \frac{x^7}{7} + \cdots = \sum_{n=0}^{\infty} \frac{(-1)^n x^{2n+1}}{2n+1}$

“这个级数在|x| ≤ 1时收敛。当x = 1时：”

$\tan^{-1}(1) = 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} - \cdots$

“而$\tan^{-1}(1) = \frac{\pi}{4}$，所以：”

$\frac{\pi}{4} = 1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} - \cdots$

“因此我们得到了第一个用级数表示π的公式：”

$\pi = 4 \left(1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} - \cdots\right) = 4 \sum_{n=0}^{\infty} \frac{(-1)^n}{2n+1}$

“这个公式看起来很优美，“我说，“但让我们看看它收敛得有多快。”

“让我们计算前几项看看精度：”

$\pi \approx 4 \left(1 - \frac{1}{3}\right) = 4 \times \frac{2}{3} \approx 2.6667$

$\pi \approx 4 \left(1 - \frac{1}{3} + \frac{1}{5}\right) = 4 \times \frac{15-5+3}{15} = 4 \times \frac{13}{15} \approx 3.4667$

$\pi \approx 4 \left(1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7}\right) = 4 \times \frac{105-35+21-15}{105} = 4 \times \frac{76}{105} \approx 2.8952$

“大家可以看到，即使计算了前4项，我们的结果离π的真实值还有很大差距，“我解释道，“这个级数收敛得非常慢。”

“老师，有没有更快的收敛方法呢？“黄毛色的学生问道。

“好问题！“我鼓励地说，“让我们尝试一些其他的几何级数方法。”

“我们知道$\sin^{-1}(x)$的泰勒展开：”

$\sin^{-1}(x) = x + \frac{1}{2} \cdot \frac{x^3}{3} + \frac{1 \cdot 3}{2 \cdot 4} \cdot \frac{x^5}{5} + \frac{1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6} \cdot \frac{x^7}{7} + \cdots$

“当x = 1/2时：”

$\sin^{-1}\left(\frac{1}{2}\right) = \frac{\pi}{6} = \frac{1}{2} + \frac{1}{2} \cdot \frac{1}{2^3 \cdot 3} + \frac{1 \cdot 3}{2 \cdot 4} \cdot \frac{1}{2^5 \cdot 5} + \frac{1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6} \cdot \frac{1}{2^7 \cdot 7} + \cdots$

“虽然这个级数比莱布尼茨级数收敛得快一些，但仍然不够理想。”

“让我们尝试另一种思路，“我说，“利用反正切函数的加法公式。我们知道：”

$\tan^{-1}(a) + \tan^{-1}(b) = \tan^{-1}\left(\frac{a+b}{1-ab}\right) \quad \text{当} \quad ab < 1$

“如果我们选择合适的a和b值，可以得到更快的收敛速度。例如：”

$\frac{\pi}{4} = \tan^{-1}(1) = \tan^{-1}\left(\frac{1}{2}\right) + \tan^{-1}\left(\frac{1}{3}\right)$

“这是因为：”

$\frac{\frac{1}{2} + \frac{1}{3}}{1 - \frac{1}{2} \cdot \frac{1}{3}} = \frac{\frac{5}{6}}{1 - \frac{1}{6}} = \frac{\frac{5}{6}}{\frac{5}{6}} = 1$

“所以我们可以得到：”

$\pi = 4 \left[\tan^{-1}\left(\frac{1}{2}\right) + \tan^{-1}\left(\frac{1}{3}\right)\right]$

“现在让我们用泰勒展开来计算这两个反正切函数：”

$\tan^{-1}\left(\frac{1}{2}\right) = \frac{1}{2} - \frac{1}{2^3 \cdot 3} + \frac{1}{2^5 \cdot 5} - \frac{1}{2^7 \cdot 7} + \cdots$

$\tan^{-1}\left(\frac{1}{3}\right) = \frac{1}{3} - \frac{1}{3^3 \cdot 3} + \frac{1}{3^5 \cdot 5} - \frac{1}{3^7 \cdot 7} + \cdots$

“让我们计算前几项看看效果：”

$\tan^{-1}\left(\frac{1}{2}\right) \approx \frac{1}{2} - \frac{1}{24} = 0.5 - 0.04167 = 0.45833$

$\tan^{-1}\left(\frac{1}{3}\right) \approx \frac{1}{3} - \frac{1}{81} = 0.33333 - 0.01235 = 0.32098$

$\frac{\pi}{4} \approx 0.45833 + 0.32098 = 0.77931$

$\pi \approx 4 \times 0.77931 = 3.11724$

“这个结果比之前的计算好一些，“我评价道，“但仍然不够精确。”

“让我们继续计算更多项：”

$\tan^{-1}\left(\frac{1}{2}\right) \approx \frac{1}{2} - \frac{1}{24} + \frac{1}{160} = 0.45833 + 0.00625 = 0.46458$

$\tan^{-1}\left(\frac{1}{3}\right) \approx \frac{1}{3} - \frac{1}{81} + \frac{1}{1215} = 0.32098 + 0.00082 = 0.32180$

$\frac{\pi}{4} \approx 0.46458 + 0.32180 = 0.78638$

$\pi \approx 4 \times 0.78638 = 3.14552$

“现在我们的结果更接近π的真实值了！“

结论与应用

“通过刚才的探索，我们发现了几何级数确实可以用于计算π，“我总结道，“但也发现了它的一些局限性。”

“那么，几何级数在计算π方面是否就毫无用处了呢？“我问道，“答案是：不是这样的！”

“虽然简单的几何级数收敛缓慢，但通过巧妙地利用各种恒等式和组合，我们可以构建更高效的级数。比如数学家们已经发现了很多更快的收敛方法。”

“让我们来看一个例题，“我在黑板上写道，例题：使用几何级数计算π的近似值，要求精确到小数点后4位。

解题思路：

1. 选择合适的级数展开
2. 计算足够的项以达到所需精度
3. 分析误差范围

解答： 我们选择使用莱布尼茨级数：

$\pi = 4 \left(1 - \frac{1}{3} + \frac{1}{5} - \frac{1}{7} + \frac{1}{9} - \frac{1}{11} + \cdots\right)$

为了保证精确到小数点后4位，我们需要计算足够的项。这个级数是一个交错级数，根据交错级数的误差估计，误差不超过下一项的绝对值。

我们需要找到n使得：

$\frac{4}{2n+1} < 0.00005$

$2n+1 > \frac{4}{0.00005} = 80000$

$n > 39999.5$

这意味着我们需要计算大约40,000项才能达到所需的精度！

实际计算结果： 计算前1000项： $\pi \approx 3.1405926538$

计算前5000项： $\pi \approx 3.1415902764$

计算前10000项： $\pi \approx 3.1415925536$

可以看到，即使计算了大量项，收敛速度仍然很慢。

更优的方法： 如果我们选择使用：

$\frac{\pi}{4} = 4 \tan^{-1}\left(\frac{1}{5}\right) - \tan^{-1}\left(\frac{1}{239}\right)$

收敛速度会快很多，只需要较少的项就能达到高精度。

“通过这个例题，我们看到了几何级数在计算π时的挑战，“我总结道，“但也看到了数学家们如何通过创新的方法来克服这些挑战。”

“今天的课告诉我们，数学不仅是知识的积累，更是智慧的体现。面对困难的问题，我们需要灵活运用各种工具，寻找最优的解决方案。“

呼噜星人的收获

“老师，我今天学到了很多东西！“蓝毛色的学生兴奋地说，“我知道了几何级数的基本概念，学会了如何用它来求和，还知道了虽然它理论上可以计算π，但在实践中有很多困难。”

“最重要的是，我理解了为什么数学家们一直在寻找更高效的计算方法，“绿毛色的学生补充道，“就像开普勒说的那样，‘数学是自然界的语言’，我们需要学会用最高效的方式来理解它。”

“我觉得数学就像是一个工具箱，“紫毛色的学生说，“有时候简单的工具就能解决大问题，但有时候我们需要更精密的工具。几何级数给我们提供了一个重要的工具，但我们还需要学习如何更好地使用它。”

“总结得很好！“我欣慰地说，“今天的课程让我们明白了几何级数在计算π方面的潜力与局限。虽然简单的几何级数收敛缓慢，但通过巧妙的数学技巧，我们可以构建更高效的算法。这充分体现了数学的魅力所在——理论上的可能性与现实中的实用技巧之间的完美结合。”

“记住，数学不仅仅是关于公式的推导和计算，更是关于如何用优雅的方式解决复杂问题的艺术。希望今天的学习能让你们对数学有更深的理解和欣赏！”

下一节：46-introducing-series