当前位置: 首页 > article >正文

巴塞尔问题详解:计算所有正整数平方的倒数之和

1 相关历史背景

巴塞尔问题(Basel Problem)是数学史上一个著名的问题,由意大利数学家皮埃特罗·门戈利(Pietro Mengoli)在1644年首次提出。 但他未能解决,只能给出小数点后六位的近似解是1.644934,但是未能给出精确值。

Pietro Mengoli:意大利数学家。生于波伦亚,卒于同地。1650年在波伦亚大学获哲学学位,毕业后留校工作,直到去逝。门戈利在数学上的主要贡献是在他的《算术求积新法》(Novae quadraturae arithmeticae1650)中证明了调和级数 ∑ n = 1 ∞ 1 n = 1 + 1 2 + 1 3 + ⋯ \sum_{n=1}^{\infty} \frac{1}{n} = 1 + \frac{1}{2} + \frac{1}{3} + \cdots n=1n1=1+21+31+ 的发散性。他将这一级数划分为一个无穷数组,使每一组中各项的和均大于1,从而得到该级数发散的结论。这项成就比雅各布·伯努利的第一个证明(1689)早了近40年。此外,他还证明了三角形级数 1 3 + 1 6 + 1 10 + ⋯ + 1 n ( n + 1 ) / 2 + ⋯ \frac{1}{3} + \frac{1}{6} + \frac{1}{10} + \cdots + \frac{1}{n(n+1)/2} + \cdots 31+61+101++n(n+1)/21+ 的和为1,进一步又求得 ∑ n 1 n ( n + r ) = 1 n ( 1 + 1 2 + 1 3 + ⋯ + 1 r ) \sum_{n} \frac{1}{n(n+r)} = \frac{1}{n} (1 + \frac{1}{2} + \frac{1}{3} + \cdots + \frac{1}{r}) nn(n+r)1=n1(1+21+31++r1) ∑ n 1 n ( n + 1 ) = 1 4 \sum_{n} \frac{1}{n(n+1)} = \frac{1}{4} nn(n+1)1=41 等结果,并利用这些结果计算了形如 ∫ 0 t x m / 2 ( 1 − x ) m / 2 d x \int_{0}^{t} x^{m/2} (1-x)^{m/2} dx 0txm/2(1x)m/2dx 类的积分 (1659)。他在无穷级数方面的工作领先于惠更斯和莱布尼茨等人。他的工作是从卡瓦列里不可分原理向牛顿流技术和莱布尼茨微分法的过渡。(摘自数学百科:http://shuxueji.com/)

问题的核心是计算所有正整数平方的倒数之和,即:

ζ ( 2 ) = ∑ n = 1 ∞ 1 n 2 \zeta(2) = \sum_{n=1}^{\infty} \frac{1}{n^2} ζ(2)=n=1n21

这个问题之所以被称为“巴塞尔问题”,是因为它引起了瑞士巴塞尔(Basel)的数学家们的广泛关注,尤其是莱昂哈德·欧拉(Leonhard Euler)的杰出贡献。


2 欧拉的突破

欧拉(Leonhard Euler)在1734年首次解决了巴塞尔问题,当时他只有28岁。
欧拉的解法基于对函数 sin ⁡ x x \frac{\sin x}{x} xsinx 的分析。他通过将 sin ⁡ x x \frac{\sin x}{x} xsinx 展开为无穷乘积,并与泰勒展开式进行比较,得到了 ζ ( 2 ) \zeta(2) ζ(2) 的精确值。
欧拉最初的结果是:

∑ n = 1 ∞ 1 n 2 = π 2 6 \sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6} n=1n21=6π2

这个结果震惊了当时的数学界,因为 π \pi π 的出现,简直布诗人。而如果对上述答案展开,小数点后六位刚好就是1.644934,很难相信这是巧合。


3 欧拉的进一步研究

在此基础上,欧拉进一步研究了更一般的级数 ζ ( s ) = ∑ n = 1 ∞ 1 n s \zeta(s) = \sum_{n=1}^{\infty} \frac{1}{n^s} ζ(s)=n=1ns1,并发现了著名的黎曼ζ函数(Riemann Zeta Function)的雏形。还证明了 ζ ( 2 ) \zeta(2) ζ(2) 的值的正确性,并通过多种方法验证了他的结果。例如,他通过计算部分和并观察其收敛性,进一步确认了 π 2 6 \frac{\pi^2}{6} 6π2 的正确性。


4. 巴塞尔问题的意义

  • 巴塞尔问题的解决展示了数学中无穷级数的深刻性质,为后来的数学家提供了研究类似问题的工具和方法。
  • 欧拉的解法启发了许多后续研究,例如傅里叶级数复分析
  • 巴塞尔问题还与物理学中的许多问题密切相关。例如,在量子力学和统计力学中, ζ ( 2 ) \zeta(2) ζ(2) 的值经常出现在各种公式中。
  • 在现代数学中,巴塞尔问题是黎曼ζ函数研究的一个特例。
  • 巴塞尔问题的解决也为后来的数学家提供了研究无穷级数的新思路。例如,拉马努金(Ramanujan)在研究无穷级数时,受到了欧拉的启发。

5 公式推导

前面讲了一堆历史,接下来就上主菜吧。

  1. 函数选择
    欧拉选择了函数 sin ⁡ x x \frac{\sin x}{x} xsinx,因为它在 x = 0 x = 0 x=0 处有定义,并且其零点为 x = ± n π x = \pm n\pi x=± n = 1 , 2 , 3 , … n = 1, 2, 3, \dots n=1,2,3,)。

  2. 无穷乘积展开
    欧拉假设 sin ⁡ x x \frac{\sin x}{x} xsinx 可以表示为无穷乘积:

    sin ⁡ x x = ∏ n = 1 ∞ ( 1 − x 2 n 2 π 2 ) \frac{\sin x}{x} = \prod_{n=1}^{\infty} \left(1 - \frac{x^2}{n^2 \pi^2}\right) xsinx=n=1(1n2π2x2)

    这个假设基于函数在 x = ± n π x = \pm n\pi x=± 处的零点。

  3. 泰勒展开
    另一方面, sin ⁡ x x \frac{\sin x}{x} xsinx 的泰勒展开式为:

    sin ⁡ x x = 1 − x 2 3 ! + x 4 5 ! − x 6 7 ! + … \frac{\sin x}{x} = 1 - \frac{x^2}{3!} + \frac{x^4}{5!} - \frac{x^6}{7!} + \dots xsinx=13!x2+5!x47!x6+


  1. 比较系数
    将无穷乘积展开与泰勒展开式进行比较。首先,展开无穷乘积的前几项:

    ∏ n = 1 ∞ ( 1 − x 2 n 2 π 2 ) = ( 1 − x 2 π 2 ) ( 1 − x 2 4 π 2 ) ( 1 − x 2 9 π 2 ) … \prod_{n=1}^{\infty} \left(1 - \frac{x^2}{n^2 \pi^2}\right) = \left(1 - \frac{x^2}{\pi^2}\right)\left(1 - \frac{x^2}{4\pi^2}\right)\left(1 - \frac{x^2}{9\pi^2}\right) \dots n=1(1n2π2x2)=(1π2x2)(14π2x2)(19π2x2)

    展开后, x 2 x^2 x2 项的系数为:

    − ( 1 π 2 + 1 4 π 2 + 1 9 π 2 + …   ) = − 1 π 2 ∑ n = 1 ∞ 1 n 2 -\left(\frac{1}{\pi^2} + \frac{1}{4\pi^2} + \frac{1}{9\pi^2} + \dots\right) = -\frac{1}{\pi^2} \sum_{n=1}^{\infty} \frac{1}{n^2} (π21+4π21+9π21+)=π21n=1n21

    泰勒展开式中 x 2 x^2 x2 项的系数为 − 1 6 -\frac{1}{6} 61

    所以将两边的 x 2 x^2 x2 项系数相等:

    − 1 π 2 ∑ n = 1 ∞ 1 n 2 = − 1 6 -\frac{1}{\pi^2} \sum_{n=1}^{\infty} \frac{1}{n^2} = -\frac{1}{6} π21n=1n21=61

  2. 求解 ζ ( 2 ) \zeta(2) ζ(2) 得到结果

    ∑ n = 1 ∞ 1 n 2 = π 2 6 \sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6} n=1n21=6π2


另外一种思路就是由于与 π \pi π 相关,所以想到和三角函数相关,两边同除 x x x ,右边有平方差公式,但是发现发散。所以右边式子应该加入常数 C C C 使得等号成立 ,接着考虑凑平方和、极限、二次项系数。

∑ n = 1 ∞ 1 n 2 = ? \sum_{n=1}^{\infty} \frac{1}{n^2} = ? n=1n21=?

泰勒: sin ⁡ x = x − x 3 3 ! + x 5 5 ! − x 7 7 ! \sin x = x - \frac{x^3}{3!} + \frac{x^5}{5!} - \frac{x^7}{7!} sinx=x3!x3+5!x57!x7

无穷和可以转为无穷积: () × () ⋯ ()× ()\cdots ()×()

由于解是 n π n\pi ,所以应该是如下:

⋯ ( x − π ) × ( x ) × ( x + π ) × ( x + 2 π ) ⋯ \cdots (x-\pi) \times (x) \times (x+\pi) \times (x+2\pi) \cdots (xπ)×(x)×(x+π)×(x+2π)

⇒ ⋯ x ( x 2 − π 2 ) ( x 2 − 4 π 2 ) ( x 2 − 9 π 2 ) ⋯ \Rightarrow \cdots x(x^2-\pi^2)(x^2-4\pi^2)(x^2-9\pi^2) \cdots x(x2π2)(x24π2)(x29π2)

但是上述式子是发散的,所以还是有问题的,应该再乘一个常数 C C C 使得其收敛。

∴ sin ⁡ x = C ⋅ x ( x 2 − π 2 ) ⋯ \therefore \sin x = C \cdot x(x^2-\pi^2) \cdots sinx=Cx(x2π2)

两边同除 x x x ,考虑极限

sin ⁡ x x = C ⋅ x ( x 2 − π 2 ) ⋯ x \frac{\sin x}{x} = \frac{C \cdot x(x^2-\pi^2) \cdots}{x} xsinx=xCx(x2π2)

∴ lim ⁡ x → 0 sin ⁡ x x = C ⋅ x ( x 2 − π 2 ) ⋯ x \therefore \lim_{x \to 0} \frac{\sin x}{x} = \frac{C \cdot x(x^2-\pi^2) \cdots}{x} x0limxsinx=xCx(x2π2)

∴ 1 = C ⋅ ( x 2 − π 2 ) ⋯ \therefore 1 = C \cdot (x^2-\pi^2) \cdots 1=C(x2π2)

∴ C = 1 ( − π ) 2 ( − 4 π ) 2 ⋯ = ∏ n = 1 ∞ ( 1 − x 2 n 2 π 2 ) \therefore C = \frac{1}{(-\pi)^2(-4\pi)^2} \cdots = \prod_{n=1}^{\infty}(1 - \frac{x^2}{n^2\pi^2}) C=(π)2(4π)21=n=1(1n2π2x2)

现在有sinx的无穷和和无穷积:
∴ sin ⁡ x = x − x 3 3 ! + x 5 5 ! − x 7 7 ! + ⋯ \therefore \sin x = x - \frac{x^3}{3!} + \frac{x^5}{5!} - \frac{x^7}{7!} + \cdots sinx=x3!x3+5!x57!x7+

sin ⁡ x = x ∏ n = 1 ∞ ( 1 − x 2 n 2 π 2 ) \sin x = x \prod_{n=1}^{\infty}(1 - \frac{x^2}{n^2\pi^2}) sinx=xn=1(1n2π2x2)

把无穷积展开,找二次项系数
x ( 1 − ( 1 π 2 + 1 2 2 π 2 + 1 3 2 π 2 ⋯   ) x 2 − ⋯   ) x\left(1 - \left(\frac{1}{\pi^2} + \frac{1}{2^2\pi^2} + \frac{1}{3^2\pi^2} \cdots\right)x^2 - \cdots\right) x(1(π21+22π21+32π21)x2)

泰勒无穷和中的二项式系数为
1 3 ! \frac{1}{3!} 3!1

∴ ∏ n = 1 ∞ 1 n 2 π 2 = 1 6 \therefore \prod_{n=1}^{\infty} \frac{1}{n^2\pi^2} = \frac{1}{6} n=1n2π21=61
∴ ∏ n = 1 ∞ 1 n 2 = π 2 6 \therefore \prod_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6} n=1n21=6π2

在这里插入图片描述


http://www.kler.cn/a/526301.html

相关文章:

  • MATLAB的数据类型和各类数据类型转化示例
  • DeepSeek-R1 模型及GRPO算法学习
  • 模型I/O功能之模型包装器
  • 一文了解性能优化的方法
  • SQL注入漏洞之绕过[前端 服务端 waf]限制 以及 防御手法 一篇文章给你搞定
  • .NET MAUI 入门学习指南
  • DeepSeek R1本地部署详细指南
  • Java 性能优化与新特性
  • [OO ALV] OO ALV 基础显示
  • allegro修改封闭图形线宽
  • 独立成分分析 (ICA):用于信号分离或降维
  • wordpress外贸独立站常用询盘软件
  • Rust语言进阶之enumerate用法实例(九十六)
  • 第33篇:Python开发进阶:自然语言处理与文本分析
  • Java继承中的静态方法隐藏与实例变量隐藏:深入解析与最佳实践
  • 年化19.3%策略集|ctpbee_api替换成openctp整合backtrader实盘方案(代码+数据)
  • 大厂面试题备份20250129
  • dify实现原理分析-rag-检索(Retrieval)服务的实现
  • 信号处理以及队列
  • 一文讲解Java中的异常处理机制
  • 变量和简单数据类型(字符串)
  • doris:导入时实现数据转换
  • Java 分布式与微服务架构:现代企业应用开发的新范式
  • JAVASE入门十二脚-IO流charArrayReader,bufferedReader,输入与输出,采集百度网页,分块操作
  • Golang 并发机制-1:Golang并发特性概述
  • 实战:如何快速让新网站被百度收录?