调和级数 java,再探π^2/6
我们知道无穷级数 $sum_{ngeq1} frac{1}{n}=O(nln n)$ 发散。事实上,质数的倒数和 $sum_{pin mathbb{P}} frac{1}{p}=O(nln ln n)$ 同样发散。简单证明可以看我的博客 《调和级数相关——质数的倒数和增长速度为O(ln ln n)》。
尽管如此,平方的倒数和收敛(尽管收敛地很慢,我们将看到),而且收敛于一个有趣的数。欧拉级数
这是 1734 年 Leonhard Euler 做出的一个经典、著名且重要的结果。这个事实的一个重要解释是它导出了 Riemann zeta 函数的第一个非平凡值 $zeta(2)=frac{pi^2}{6}$ 。不仅这一结果在数学史上有显赫的地位,它的几个极其优美聪明的证明也拥有自己的历史。本文就分享几个极为精妙的证明。
Euler 当然是十分聪明的,这个证明也十分优美。首先注意到 $sin x$ 的泰勒展开
$$sin x=x-frac{x^3}{3!}+frac{x^5}{5!}-frac{x^7}{7!}+...$$
从而
$frac{sin x}{x}=1-frac{x^2}{3!}+frac{x^4}{5!}-frac{x^6}{7!}+...$
令 $frac{sin x}{x}=0$ 解得 $x=k pi (k in mathbb{Z}$ 且 $k neq 0)$
故我们将 $frac{sin x}{x}$ 因式分解
$frac{sin x}{x}=(1-frac{x}{pi})(1+frac{x}{pi})(1-frac{x}{2pi})(1+frac{x}{2pi})(1-frac{x}{3pi})(1+frac{x}{3pi})...$
$=(1-frac{x^2}{pi^2})(1-frac{x^2}{4pi^2})(1-frac{x^2}{9pi^2})...$
(Euler 似乎并没有证明这个因式分解的正确性。幸运的是,一个世纪以后,Weierstrass 提出了著名的 Weierstrass 分解定理,利用复分析证明了这个结论。)
考虑这个因式分解的 $x^2$ 项前系数,对比泰勒展开式,可以知道
$-(frac{1}{pi^2}+frac{1}{4pi^2}+frac{1}{9pi^2}+...)=-frac{1}{3!}=-frac{1}{6}$
两边同乘 $-pi^2$ 即
$sum_{n geq 1} frac{1}{n^2}=frac{pi^2}{6}$
证毕
初等数学证明
在所有证明中间,这个可能是最初等的证明。只需要知道三角函数和韦达定理,就可以基本上证明了。
我们首先建立关于余切函数和多项式之间的关系
$sum_{i=1}^{m} cot^2(frac{i}{2m+1}pi) =frac{2m(2m-1)}{6}$
为了验证这个关系,我们从
$cos nx +isin nx =(cos x +isin x)^n$
开始,取虚部,得
$sin nx =C_n^1 sin xcos^{n-1} x-C_n^3sin^3 xcos^{n-3} x+...$
令 $n=2m+1$ ,两边同除 $sin^n x$ 得到
$frac{sin((2m+1)x)}{sin^{2m+1} x} =C_{2m+1}^1 cot^{2m} x-C_{2m+1}^3 xcot^{2m} x+...$
令 $x=frac{i}{2m+1}pi (i=1,2,…,m)$ ,左边 $sin((2m+1)x)=0$ ,于是都有
$0=C_{2m+1}^1 cot^{2m} x-C_{2m+1}^3 xcot^{2m} x+...$
因此方程
$C_{2m+1}^1t^m-C_{2m+1}^3t^{m-1}+...+(-1)^mC_{2m+1}^{2m+1}=0$
有 $m$ 个不同的根
$t_i=cot^2(frac{i}{2m+1}pi) (i=1,2,...,m)$
由韦达定理可知 $m$ 根之和为
$sum_{i=1}^{m} cot^2(frac{i}{2m+1}pi) =frac{C_{2m+1}^{3}}{C_{2m+1}^{1}}=frac{2m(2m-1)}{6}$
这样我们就得到了关于余切函数和多项式之间的关系。根据余割函数和余切函数之间的关系,同样的我们也可以建立余割函数和多项式之间的关系
$sum_{i=1}^{m} csc^2(frac{i}{2m+1}pi) =sum_{i=1}^{m} 1+cot^2(frac{i}{2m+1}pi) =frac{2m(2m-1)}{6}+m=frac{2m(2m+1)}{6}$
这样我们的准备工作已经全部准备好了,现在到了大展身手的时候
考虑 $yin(0,frac{pi}{2})$ 有
$0
从而
$0
这表明
$0
现在把 $t$ 的 $m$ 个不同根带入,左边用余切函数平方和,右边用余割函数平方和,可以得到
$frac{2m(2m-1)}{6}
两边同乘 $frac{(2m+1)^2}{pi^2}$ 可得
$frac{pi^2}{6}frac{2m(2m-1)}{(2m+1)^2}
当 $m$ 趋近于 $infty$ 时,左右两端都收敛到 $frac{pi^2}{6}$ 。证毕。
还有一些极为简短的证明,下面分享两个。
傅里叶分析
考虑函数 $f(x)=x^2$ ,$xin (-pi,pi) $ ,将其傅里叶展开
$f(x)=frac{pi^2}{3}+sum_{n=1}^{infty} ((-1)^nfrac{4}{n^2}cos nx)$
显而易见,带入 $f(0)$ 即可得到答案
伯努利数的证明
函数 $B(x)=frac{x}{e^x-1}$ 是伯努利数 $B_k$ 的生成函数,有 $B$ 是亚纯,且只在 $2npi i$ 处有极点,利用Mittag-Leffler定理可以展开为
$frac{x}{e^x-1}=sum_ {ninmathbb{Z}}frac{2npi i}{x-2npi i}=sum_{ninmathbb{Z}}-frac{1}{1-frac{x}{2npi i}}=-sum _{ninmathbb{Z}} (frac{x}{2npi i})^k =sum _{n geq 0} (-1)^{n+1}frac{2zeta(2n)}{(2pi)^{2n}}x^{2n}$
其中 $zeta(s)$ 为Riemann zeta 函数,对实数 $s>1$ 定义为 $zeta(s)=sum_{n geq 1} frac{1}{n^s}$
只剩下偶数项是因为几何级数正负相消了,同时我们也可以得到奇数项的伯努利数为 $0$ ,对比偶数项可知
$B_{2n}=(-1)^{n+1}frac{2zeta(2n)}{(2pi)^{2n}}$
所以 $frac{2zeta(2)}{(2pi)^{2}}=B_2=frac{1}{12}$ ,于是 $zeta(2)=frac{pi^2}{6}$ ,证毕。
应用:任意两个正整数互质的概率
为了求任意两个正整数互质的概率。记 $p_i$ 为从小到大第 $i$ 个质数。则有任意一个正整数被 $p_i$ 整除的概率为$frac{1}{p_i}$ ,故任选两个正整数都被 $p_i$ 整除的概率为 $frac{1}{p_i^2}$ 。于是可以得到任意两个正整数互质的概率为
$Pi _{i in mathbb{N}*} (1-frac{1}{p_i^2}) =frac{1}{Pi _{i in mathbb{N}*} frac{1}{1-frac{1}{p_i^2} }} =frac{1}{Pi _{i in mathbb{N*}}(1+frac{1}{p_i^2}+frac{1}{p_i^4}+frac{1}{p_i^6}+...) }=frac{1}{sum _{n in mathbb{N*}}frac{1}{n^2}}=frac{6}{pi^2}$
其中
$Pi _{i in mathbb{N*}}(1+frac{1}{p_i^2}+frac{1}{p_i^4}+frac{1}{p_i^6}+...) =sum _{n in mathbb{N*}}frac{1}{n^2}$
是由于算数基本定理,任何正整数质因数分解唯一。
故任意两个正整数互质的概率为 $frac{6}{pi^2} $ 。
调和级数 java,再探π^2/6相关推荐
- 再探C#类与结构体究竟谁快——考虑栈变量、栈分配、64位整数、密封类
上次我对C#类与结构体做了一次速度评测(http://blog.csdn.net/zyl910/article/details/6788417).经过一段时间思索,发现还可以进一步探讨-- 第一.栈变 ...
- SAP MM 工序委外场景再探 – 相关采购申请单据中的总账科目确定
SAP MM 工序委外场景再探 – 相关采购申请单据中的总账科目确定 如下的采购申请,是工序委外场景的采购申请, 它在工单1000030 Release的时候被自动触发. 问题是:这个PR上的总账科目 ...
- 【Vue.js 牛刀小试】:第九章 - 组件基础再探(data、props)
系列目录地址 一.基础知识概览 第一章 - 一些基础概念(posted at 2018-10-31) 第二章 - 常见的指令的使用(posted at 2018-11-01) 第三章 - 事件修饰符的 ...
- RxSwift学习--核心逻辑再探
前言 通过上一篇内容RxSwift学习--核心逻辑初探,对RxSwift有了些初步的认知,下面通过源码来看一下RxSwift到底有多骚 RxSwift核心逻辑再探 先把上篇中的例子代码搬过来: //第 ...
- 【CTF大赛】2021 DASCTF July cybercms 一探再探
引言 在前不久结束的 2021 DASCTF July X CBCTF 4th 比赛中,有一道名为 cybercms 的 web 题目. 预期解是从后台登录处进行 SQL 注入写入一句话木马,然而咱在 ...
- 联发科mtk手机处理器怎么样_5G手机价格将再探新低!联发科天玑700芯片发布:入门级5G处理器...
↑↑↑点击上方蓝字订阅每日最新国产手机资讯 今年5G手机的普及之所以能进展的如此迅速,联发科在5G芯片领域的的活跃功不可没.今天一早,联发科发布的全新的天玑700移动平台,标志着5G手机的价格将进一步 ...
- C++ 类的知识 | 构造函数再探、匿名对象、友元函数、内部类、类的const成员、类的static成员
文章目录 构造函数再探 以下代码共调用多少次拷贝构造函数 委托构造函数 概念 形式 匿名对象 友元 友元的声明 友元类 令成员函数作为友元 函数重载和友元 注意 内部类 特性 类的const成员 可变 ...
- 函数重载、引用再探、内联函数
文章目录 函数重载 为什么C++支持重载,C语言不支持呢? extern "C" 引用再探 引用的特性 引用的使用场景 引用和指针 引用和指针的不同点: 内联函数 什么是内联函数? ...
- 【再探backbone 02】集合-Collection
前言 昨天我们一起学习了backbone的model,我个人对backbone的熟悉程度提高了,但是也发现一个严重的问题!!! 我平时压根没有用到model这块的东西,事实上我只用到了view,所以昨 ...
- C++之再探参数绑定bind、bind1st、bind2nd、placeholders占位符
相关博文:C++头文件<functional>和bind.placeholders占位符使用简单例子 相关博文:<Essential C++>笔记之设计一个泛型算法(二) 相关 ...
最新文章
- SharePoint 2010 BI(2):使用Visio“.NET研究” Service
- QT的QColorDialog类的使用
- Flutter之运行提示Could not update files on device: Connection closed before full header was received
- TensorFlow 之快速上手详解
- 51全志R58平台Android4.4下Camera的HAL层修改
- 2016 Multi-University Training Contest 1 T4
- Go语言的goroutine
- django和celery结合应用
- matlab读取写入图像数据格式uint8,double
- 程序员常识--OJ系统及ACM测试题库大全
- 基于ssm的城市公交查询系统的设计与实现(附源码)
- AI语音外呼机器人是如何帮助电销行业获客?
- 201871010126 王亚涛 《面向对象程序设计(java)》 第四周学习总结
- 陈小龙书pHP,PHP
- MySQL数据库基础详解
- Mimikatz各模块基本用法
- 印象笔记目录导出(失败)
- java atm柜机_java课程设计-ATM柜机模拟程序.doc
- MRI脑影像分析——根据脑图谱获取感兴趣区域mask,以海马体与丘脑为例(matlab+nilearn+nibabel+REST1.8)
- anaconda安装软件包
热门文章
- Office Web Add-in的技术原理和开发常见问题剖析
- 在SCI-Hub上下载国外论文
- WPS 2019解决不能启用宏问题
- 二进制的转换(二进制、八进制、十进制、十六进制)
- Java基础 - 替罪羊树(Scapegoat Tree)
- 杭州城市交通拥堵综合治理实践
- 2021-10-21-Flask-04 模型关系的建立,增删改查的演示,一篇文章让你直接搞懂Flask模型
- 高一计算机课期中考试总结反思,信息技术期中考试试卷分析与反思
- 黑客站在 ATM 面前,机器就直接吐出钞票,他们是怎么做到的?
- matlab中fprintf整数,matlab中fprintf函数的用法