[欧拉数学]伯努利级数及相关级数的总结
By 苏剑林 | 2016-03-20 | 106838位读者 |最近在算路径积分的时候,频繁地遇到了以下两种无穷级数:
∑n1n2±ω2和∏n(1±ω2n2)
当然,直接用Mathematica可以很干脆地算出结果来,但是我还是想知道为什么,至少大概地知道。
伯努利级数 #
当ω=0的时候,第一个级数变为著名的伯努利级数
∑n1n2=1+14+19+116+…
既然跟伯努利级数有关,那么很自然想到,从伯努利级数的求和入手。
伯努利级数求和的最为人称颂的方法,当属欧拉通过大胆猜测、类比而给出的“证明”了。欧拉考虑级数
sin√x√x=1−x6+x2120−x35040+…
从左端的式子,容易看出sin√x√x=0的根为n2π2,n=1,2,3,…。
接着,我们考虑一元n次方程
1+a1x+a2x2+⋯+anxn=0
假设它有n个非零根x1,x2,…,xn,那么我们可以将它看成是关于1x的方程
(1x)n+a1(1x)n−1+⋯+an=0
对应的根变成了1/x1,1/x2,…,1/xn,那么根据韦达定理得到
1x1+1x2+⋯+1xn=−a1
欧拉的大胆之处,就是它将上述的韦达定理,直接推广到无穷级数中去。他说,既然sin√x√x=0的根为n2π2,n=1,2,3,…,那么就有
1π2+122π2+132π2+⋯=−a1=16
于是他告诉我们
1+14+19+116+⋯=π26
大胆、巧妙而不失简洁,欧拉完美地体现了冒险者的精神和令人陶醉的美!无怪乎John Stellwell在他的《数学及其历史》中,称“欧拉或许是级数运算的最伟大的学者”。
谨慎,审视 #
记得第一次看到这个证明,笔者就对欧拉佩服得五体投地了。这是怎样的头脑才能想出这样的方法!在那个级数运算才刚刚发展起来的年代,他就敢于将简单的韦达定理推广到无穷。更要命的是,他得到了正确的结果!正确的结果!~
不过如今学习更加深入,现在重看他的过程,就会带着审视的眼光了,比如哪里不严格了,哪里可以推广,等等。将韦达定理推广到无穷,并非不可能,但是有几个问题需要解决,比如怎么保证列出了所有根了呢?韦达定理中的那些根,可以是复根,而n2π2,n=1,2,3,…最多也就是全部实数根,还有复根吗?不用sin√x√x,换成其他函数,有没有可能得到其它类似的有趣结果呢?
当然,在这里只是稍微一提,这不是本文的主题,因此并不多做讨论。事实上,欧拉的方法能衍生出丰富的结果,它的启发性、创造性已经远远掩盖了它的不严格性了。
连乘积的结果 #
事实上,欧拉的方法所告诉我们的,远远多于伯努利级数和这一特例。从代数理论中我们还知道,如果一元n次方程1+a1x+a2x2+⋯+anxn=0的n个非零根为x1,x2,…,xn,那么就有
1+a1x+a2x2+⋯+anxn=an(x−x1)(x−x2)…(x−xn)
将x=0代进去,得到
1=an(−x1)(−x2)…(−xn)
用这个关系,可以将上式改写成
1+a1x+a2x2+⋯+anxn=(1−xx1)(1−xx2)…(1−xxn)
欧拉又说,这个能推广到无穷!所以有
sin√x√x=(1−xπ2)(1−x22π2)(1−x32π2)…
或者换一下变量
sinxx=(1−x2π2)(1−x222π2)(1−x232π2)…
代入x=ωπ,得到
sinωπωπ=(1−ω212)(1−ω222)(1−ω232)…
这就是文章一开头提及到的连乘积之一。如果把ω换成iω,就得到
sin(iωπ)iωπ=(1+ω212)(1+ω222)(1+ω232)…
利用sin(ix)=isinhx,就得到
sinhωπωπ=(1+ω212)(1+ω222)(1+ω232)…
这是连乘积的另一种情况。
回来继续求和 #
继续我们的求和,如果把sin√x√x换为
ωsin(π√x+ω2)√x+ω2sinωπ=1+πωcotωπ−12ω2x+…
这一次,方程
ωsin(π√x+ω2)√x+ω2sinωπ=0
的所有根为x=n2−ω2,n=1,2,3,…,因此有
∑n1n2−ω2=1−πωcotωπ2ω2
将ω换为iω,得到
∑n1n2+ω2=1−iπωcotiωπ2(iω)2=πωcothωπ−12ω2
这就解决了两个求和问题了。
欧拉,真乃神人也! #
百科上说到:欧拉从19岁开始发表论文,直到76岁,共写下了886本书籍和论文,其中在世时发表了700多篇论文。彼得堡科学院为了整理他的著作,整整用了47年!不得不说,这才是当之无愧的学神呐!
转载到请包括本文地址:https://spaces.ac.cn/archives/3680
更详细的转载事宜请参考:《科学空间FAQ》
如果您还有什么疑惑或建议,欢迎在下方评论区继续讨论。
如果您觉得本文还不错,欢迎分享/打赏本文。打赏并非要从中获得收益,而是希望知道科学空间获得了多少读者的真心关注。当然,如果你无视它,也不会影响你的阅读。再次表示欢迎和感谢!
如果您需要引用本文,请参考:
苏剑林. (Mar. 20, 2016). 《[欧拉数学]伯努利级数及相关级数的总结 》[Blog post]. Retrieved from https://spaces.ac.cn/archives/3680
@online{kexuefm-3680,
title={[欧拉数学]伯努利级数及相关级数的总结},
author={苏剑林},
year={2016},
month={Mar},
url={\url{https://spaces.ac.cn/archives/3680}},
}
March 21st, 2016
欧拉的著作由瑞士科学院欧拉委员会整理,现在还没有整理完全。欧拉委员会官方网站:http://www.leonhard-euler.ch/
至今还没有整理完全?
June 24th, 2018
方程sin(√x)/(√x)=0没有复数根,这一点很容易证明。就象exp(z)不可能等于0一样。
September 20th, 2018
[图片]苏兄,最后那个换是怎么得来的啊?
icotix=cothx
谢谢苏兄的回复,我是问你
继续我们的求和,如果把sin(√x)/(√x)
换为,下面那个东西怎么想到的?
就是类比sin√x√x求1+1/4+1/9+…的过程,构造所有根为x=n2−ω2,n=1,2,3,…的一个函数咯。有了sin√x√x这个例子,应该不难理解吧。
是的,我也这么想,应该是sin(pi√x+w^2)/(pi√x+w^2)才是,但是你还乘了w,除了sin(wpi),就没想明白了
sinπ√x+ω2π√x+ω2|x=0=sinωπωπ
要除以sinωπωπ,它的泰勒展开的第一项才为1
苏兄,我想问你个积分,你能算下吗?想了好久,都没弄出来,
觉得你很强,函数ln(1+x) /[x(1+x^2)]从0到1的积分。
September 25th, 2018
苏兄非常感谢你的耐心指导!!!
September 26th, 2018
@李刚|comment-9827
我确认下,是不是这个?
∫10ln(1+x)x(1+x2)dx
是的,苏兄!
令
F(a)=∫10ln(1+ax)x(1+x2)dx
那么原问题就是求F(1),并且我们有F(0)=0。对a求导得到
F′(a)=∫101(1+ax)(1+x2)dx
这个积分可以直接积出来,并不困难,结果是
F′(a)=−aln4+4aln(a+1)+π4a2+4
从而可以继续对a求积分得到
F(1)=∫10F′(a)da=116(π2−ln24)+∫10aln(1+a)1+a2da
而
∫10aln(1+a)1+a2da=∫10[ln(1+a)a−ln(1+a)a(1+a2)]da=∫10ln(1+a)ada−F(1)
从而
F(1)=132(π2−ln24)+12∫10ln(1+a)ada
现在剩下的积分,如果将ln(1+a)展开为a的泰勒级数,然后逐项积分,会发现等于
F(1)=132(π2−ln24)+12(1−122+132−142+…)
括号内的级数求和欧拉已经算出过了,就是π/12,所以综合起来就是
F(1)=7π296−ln228
(自己乱想的,可能不规范,答案应该是对的)
September 28th, 2018
苏兄强啊,我这样算,计算能力不行啊,我再算算,苏兄现在是干嘛啊?
有网友提供了一个答案,结果几乎一样。原来我这样乱搞,还跟标准答案几乎一样了。
函数f(x)在[0,1]连续,(0,1)可导,f(0)=0,f(1)=2,证明:在(0,1)内至少存在ξ,f'(ξ)=2f(ξ)+1,苏兄,你看看这道题呗,想了好久一直构造不出来。
不用做了,这是个错题,苏兄
September 28th, 2018
尴尬,苏兄比我小,你还在研究生在读呢!
March 7th, 2020
我发现把伯努利级数看成了伯努利数。。。(尴尬)不过也挺好看的,感谢。( •̆ ᵕ •̆ )◞♡
^_^