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

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.mzph.cn/bicheng/69691.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

android 圆形弹窗摄像头开发踩坑——源码————未来之窗跨平台操作

一、飘窗刷脸&#xff0c;拍照采用飘窗 刷脸认证安卓接口采用飘窗具有在不干扰用户主要操作的前提下以醒目方式引导用户完成认证&#xff0c;且能灵活定制样式以提升用户体验和认证效率的优点 二、踩坑只有一个扇形 <?xml version"1.0" encoding"utf-8&quo…

电子电气架构 --- 在智能座舱基础上定义人机交互

我是穿拖鞋的汉子&#xff0c;魔都中坚持长期主义的汽车电子工程师。 老规矩&#xff0c;分享一段喜欢的文字&#xff0c;避免自己成为高知识低文化的工程师&#xff1a; 简单&#xff0c;单纯&#xff0c;喜欢独处&#xff0c;独来独往&#xff0c;不易合同频过着接地气的生活…

图漾相机——Sample_V1示例程序

文章目录 1.SDK支持的平台类型1.1 Windows 平台1.2 Linux平台 2.SDK基本知识2.1 SDK目录结构2.2 设备组件简介2.3 设备组件属性2.4 设备的帧数据管理机制2.5 SDK中的坐标系变换 3.Sample_V1示例程序3.1 DeviceStorage3.2 DumpCalibInfo3.3 NetStatistic3.4 SimpleView_SaveLoad…

如何使用formlinker,重构微软表单创建的数字生产力法则?

仅需三步&#xff1a;上传文件-下载文件-导入文件到微软表单 凌晨两点的格式炼狱&#xff1a;被浪费的300万小时人类创造力 剑桥大学的实验室曾捕捉到一组震撼数据&#xff1a;全球教育工作者每年花在调整试题格式上的时间&#xff0c;足够建造3座迪拜哈利法塔。当北京某高校的…

idea对jar包内容进行反编译

1.先安装一下这个插件java Bytecode Decompiler 2.找到这个插件的路径&#xff0c;在idea的plugins下面的lib文件夹内&#xff1a;java-decompiler.jar。下面是我自己本地的插件路径&#xff0c;以作参考&#xff1a; D:\dev\utils\idea\IntelliJ IDEA 2020.1.3\plugins\java-d…

C#方法(练习)

1.定义一个函数&#xff0c;输入三个值,找出三个数中的最小值 2.定义一个函数&#xff0c;输入三个值,找出三个数中的最大值 3.定义一个函数&#xff0c;输入三个值,找出三个数中的平均值 4.定义一个函数&#xff0c;计算一个数的 N 次方 Pow(2, 3)返回8 5.传入十一…

1.五子棋对弈python解法——2024年省赛蓝桥杯真题

问题描述 原题传送门&#xff1a;1.五子棋对弈 - 蓝桥云课 "在五子棋的对弈中&#xff0c;友谊的小船说翻就翻&#xff1f;" 不&#xff01;对小蓝和小桥来说&#xff0c;五子棋不仅是棋盘上的较量&#xff0c;更是心与心之间的沟通。这两位挚友秉承着"友谊第…

基于STM32的智能停车场管理系统设计

目录 引言系统设计 硬件设计软件设计 系统功能模块 车辆识别与进出管理模块车位检测与引导模块计费与支付模块数据存储与查询模块远程监控与异常报警模块 控制算法 车牌识别与车辆进出管理算法车位检测与引导算法计费与支付处理算法数据存储与远程反馈算法 代码实现 车辆检测与…

单细胞-第五节 多样本数据分析,打分R包AUCell

文件在单细胞\5_GC_py\1_single_cell\3.AUCell.Rmd 1.基因 rm(list = ls()) load("g.Rdata")2.AUCell https://www.ncbi.nlm.nih.gov/pmc/articles/PMC9897923 IF: NA NA NA用这个文章里的方法,将单细胞亚群的marker基因与ros相关基因取交集,用作AUCell的基因集…

HBase-2.5.10 伪分布式环境搭建【Mac】

文章目录 前言一、搭建单节点Zookeeper1. 解压zookeeper2. 配置环境变量3. 修改配置文件4. 启动zk 二、搭建伪分布式Hbase1. 解压hbase2. 配置环境变量3. 修改配置4. 启动HBase 前言 搭建hbase伪分布式环境 提示&#xff1a;以下是本篇文章正文内容&#xff0c;下面案例仅供参…

蓝牙技术在物联网中的应用有哪些

蓝牙技术凭借低功耗、低成本和易于部署的特性&#xff0c;在物联网领域广泛应用&#xff0c;推动了智能家居、工业、医疗、农业等多领域发展。 智能家居&#xff1a;在智能家居系统里&#xff0c;蓝牙技术连接各类设备&#xff0c;像智能门锁、智能灯泡、智能插座、智能窗帘等。…

NLP深度学习 DAY5:Seq2Seq 模型详解

Seq2Seq&#xff08;Sequence-to-Sequence&#xff09;模型是一种用于处理输入和输出均为序列任务的深度学习模型。它最初被设计用于机器翻译&#xff0c;但后来广泛应用于其他任务&#xff0c;如文本摘要、对话系统、语音识别、问答系统等。 核心思想 Seq2Seq 模型的目标是将…

DeepSeek模型:开启人工智能的新篇章

DeepSeek模型&#xff1a;开启人工智能的新篇章 在当今快速发展的技术浪潮中&#xff0c;人工智能&#xff08;AI&#xff09;已经成为了推动社会进步和创新的核心力量之一。而DeepSeek模型&#xff0c;作为AI领域的一颗璀璨明珠&#xff0c;正以其强大的功能和灵活的用法&…

单细胞-第四节 多样本数据分析,下游画图

文件在单细胞\5_GC_py\1_single_cell\2_plots.Rmd 1.细胞数量条形图 rm(list ls()) library(Seurat) load("seu.obj.Rdata")dat as.data.frame(table(Idents(seu.obj))) dat$label paste(dat$Var1,dat$Freq,sep ":") head(dat) library(ggplot2) lib…

freeswitch在centos上编译过程

操作系统&#xff1a;centos9-last usr/local/freeswitch/bin/freeswitch -version FreeSWITCH version: 1.10.13-devgit~20250125T131725Z~3f1e4bf90a~64bit (git 3f1e4bf 2025-01-25 13:17:25Z 64bit)vi /etc/ssh/sshd_config ip a nmtui reboot ip a curl -o /etc/pki/rpm-…

NLP模型大对比:Transformer >Seq2Seq > LSTM > RNN > n-gram

结论 Transformer 大于 传统的Seq2Seq 大于 LSTM 大于 RNN 大于 传统的n-gram n-gram VS Transformer 我们可以用一个 图书馆查询 的类比来解释它们的差异&#xff1a; 一、核心差异对比 维度n-gram 模型Transformer工作方式固定窗口的"近视观察员"全局关联的&q…

Julius AI 人工智能数据分析工具介绍

Julius AI 是一款由 Casera Labs 开发的人工智能数据分析工具&#xff0c;旨在通过自然语言交互和强大的算法能力&#xff0c;帮助用户快速分析和可视化复杂数据。这款工具特别适合没有数据科学背景的用户&#xff0c;使数据分析变得简单高效。 核心功能 自然语言交互&#x…

H3CNE-31-BFD

Bidirectional Forwarding Dection&#xff0c;双向转发检查 作用&#xff1a;毫秒级故障检查&#xff0c;通常结合三层协议&#xff08;静态路由、vrrp、ospf、BGP等&#xff09;&#xff0c;实现链路故障快速检查。 BFD配置示例 没有中间的SW&#xff0c;接口down&#xff…

2025最新版MySQL安装使用指南

2025最新版MySQL安装使用指南 The Installation and Usage Guide of the Latest Version of Oracle MySQL in 2025 By JacksonML 1. 获取MySQL 打开Chrome浏览器&#xff0c;访问官网链接&#xff1a;https://www.mysql.com/ &#xff0c;随即打开MySQL官网主页面&#xff…

基于STM32的智能温控花盆设计

目录 引言系统设计 硬件设计软件设计 系统功能模块 空气质量传感器模块数据处理与分析模块实时显示与用户交互模块报警与提示模块远程监控与数据上传模块 控制算法 空气质量数据处理与分析算法异常检测与报警算法数据上传与历史数据回溯算法 代码实现 空气质量检测与数据处理代…