巴塞尔问题是一个著名的数论问题,这个问题首先由义大利数学家皮耶特罗·门戈利在1644年提出,瑞士数学家莱昂哈德·欧拉于1735年解决。由于这个问题难倒了以前许多的数学家,年仅二十八岁的欧拉因此一举成名。欧拉把这个问题作了一番推广,他的想法后来被德国数学家黎曼在1859年的论文《论小于给定大数的质数个数》(Ueber die Anzahl der Primzahlen unter einer gegebenen Grösse)中所采用,论文中定义了黎曼ζ函数,并证明了它的一些基本的性质。这个问题是以瑞士的第三大城市巴塞尔命名的,它是欧拉和伯努利家族的家乡。
这个问题是精确计算所有平方数的倒数的和,也就是以下级数的和:
![{\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n^{2}}}=\lim _{n\to +\infty }\left({\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+\cdots +{\frac {1}{n^{2}}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/696f778298469c92576fda48f6a23c89b75c3813)
这个级数的和大约等于1.644934(OEIS数列A013661)。巴塞尔问题是寻找这个数的准确值,并证明它是正确的。欧拉发现准确值是
,并在1735年公布;彼时他给出了一个错误的证明,真正严密的证明在1741年给出。
利用放缩法可以简便地证明级数
有界,并由单调收敛定理得到收敛性。
注意到
且
,
。故
故
。
欧拉最初推导
的方法是聪明和新颖的。他假设有限多项式的性质对于无穷级数也是成立的。然而,欧拉没有证明此一假设,且此一假设在一般情况下也是错误的。不过他计算了级数的部分和后发现,级数确实趋近
,不多不少。这给了他足够的自信心,把这个结果公诸于众。
欧拉的方法是从正弦函数的泰勒级数展开式开始:
![{\displaystyle \sin x=x-{\frac {x^{3}}{3!}}+{\frac {x^{5}}{5!}}-{\frac {x^{7}}{7!}}+\cdots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/15a9027c77b90e40215c01281b39d37730e7e537)
两边除以
,得:
![{\displaystyle {\frac {\sin x}{x}}=1-{\frac {x^{2}}{3!}}+{\frac {x^{4}}{5!}}-{\frac {x^{6}}{7!}}+\cdots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/52e2f18cc4f8100e7be16afc351e0a7c1b7adc8e)
现在,
的根出现在
,其中
我们假设可以把这个无穷级数表示为线性因子的乘积,就像把多项式因式分解一样:
![{\displaystyle {\begin{aligned}{\frac {\sin x}{x}}&{}=\left(1-{\frac {x}{\pi }}\right)\left(1+{\frac {x}{\pi }}\right)\left(1-{\frac {x}{2\pi }}\right)\left(1+{\frac {x}{2\pi }}\right)\left(1-{\frac {x}{3\pi }}\right)\left(1+{\frac {x}{3\pi }}\right)\cdots \\&{}=\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)\cdots \end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2ad0281f1661f711479e127fb13c720d89871fcf)
如果把这个乘积展开,并把所有有
的项收集在一起,我们可以看到,
的二次项系数为:
![{\displaystyle -\left({\frac {1}{\pi ^{2}}}+{\frac {1}{4\pi ^{2}}}+{\frac {1}{9\pi ^{2}}}+\cdots \right)=-{\frac {1}{\pi ^{2}}}\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7c39ab3b889b31223d9e1a3efce3eb6b95eeb1e4)
但从
原先的级数展开式中可以看出,
的系数是
。这两个系数一定是相等的;因此,
![{\displaystyle -{\frac {1}{6}}=-{\frac {1}{\pi ^{2}}}\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/137f5522a36f142a13fa24086ada72f53fd0cd70)
等式两边乘以
就可以得出所有平方数的倒数之和。
![{\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n^{2}}}={\frac {\pi ^{2}}{6}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26a915462f991c41d6da8e006e857aed0ee427cf)
黎曼ζ函数ζ(s)是数学中的一个很重要的函数,因为它与素数的分布密切相关。这个函数对于任何实数部分大于1的复数s都是有定义的,由以下公式定义:
![{\displaystyle \zeta (s)=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/31f460a787ebaf667aaf7805f2b87e542a02836b)
取s = 2,我们可以看出ζ(2)等于所有平方数的倒数之和:
![{\displaystyle \zeta (2)=\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}={\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+{\frac {1}{3^{2}}}+{\frac {1}{4^{2}}}+\cdots ={\frac {\pi ^{2}}{6}}\approx 1.644934}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd93c00613f44f50c10a0e5e6c7de96bc7ffe57e)
用以下的等式,可以证明这个级数收敛:
![{\displaystyle \sum _{n=1}^{N}{\frac {1}{n^{2}}}<1+\sum _{n=2}^{N}{\frac {1}{n(n-1)}}=1+\sum _{n=2}^{N}\left({\frac {1}{n-1}}-{\frac {1}{n}}\right)=1+1-{\frac {1}{N}}\,{\xrightarrow {N\to \infty }}\,2}](https://wikimedia.org/api/rest_v1/media/math/render/svg/432ece1ad6b27585ac27c3e2190e5590c46315ca)
因此ζ(2)的上界小于2,因为这个级数只含有正数项,它一定是收敛的。可以证明,当s是正的偶数时,ζ(s)可以用伯努利数来表示。设
,有以下公式:
![{\displaystyle \zeta (2n)={\frac {(2\pi )^{2n}(-1)^{n+1}B_{2n}}{2\cdot (2n)!}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e6b4602234b165ad103209f1d4cbdf4833068d19)
以下介绍了一个
的证明。它是目前已知最基本的证明,大部分其它的证明都需要用到傅里叶分析、复分析和多变量微积分,但这个证明连一元微积分也不需要(在证明的最后部分需要使用极限的概念)。
考虑面积,
![{\displaystyle \tan \theta >\theta >\sin \theta }](https://wikimedia.org/api/rest_v1/media/math/render/svg/009a75107ead52731644af01338a96a7800abf0e)
![{\displaystyle {\frac {1}{\tan \theta }}<{\frac {1}{\theta }}<{\frac {1}{\sin \theta }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95a9b6ef1c157cbb04dfb507bbb9b53e91e97539)
![{\displaystyle \cot ^{2}\theta <{\frac {1}{\theta ^{2}}}<\csc ^{2}\theta }](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f9b6146f0ce2b5755918e33692d67b77677fabb)
这个证明的想法是把以下的部分和固定在两个表达式之间,这两个表达式当m趋于无穷大时都趋于
。
![{\displaystyle \sum _{k=1}^{m}{\frac {1}{k^{2}}}={\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+\cdots +{\frac {1}{m^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9afbce6d74ff38101967e217387c8e873a4dfa30)
这两个表达式从余切和余割的恒等式推出。而这些恒等式则从棣莫弗定理推出。
设x为一个实数,满足0 < x < π/2,并设n为正整数。从棣莫弗定理和余切函数的定义,可得:
![{\displaystyle {\frac {\cos(nx)+i\sin(nx)}{(\sin x)^{n}}}={\frac {(\cos x+i\sin x)^{n}}{(\sin x)^{n}}}=\left({\frac {\cos x+i\sin x}{\sin x}}\right)^{n}=(\cot x+i)^{n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2dce3548c59bc3443237347783da2a0c61207986)
根据二项式定理,我们有:
![{\displaystyle (\cot x+i)^{n}={n \choose 0}\cot ^{n}x+{n \choose 1}(\cot ^{n-1}x)i+\cdots +{n \choose {n-1}}(\cot x)i^{n-1}+{n \choose n}i^{n}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c119769323602638eff639db6ce53ba08911714c)
![{\displaystyle =\left[{n \choose 0}\cot ^{n}x-{n \choose 2}\cot ^{n-2}x\pm \cdots \right]\;+\;i\left[{n \choose 1}\cot ^{n-1}x-{n \choose 3}\cot ^{n-3}x\pm \cdots \right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9d4fe91ed922b8011563ae45e34bdf126dfc3306)
把两个方程合并,由于相等的两个复数的虚数部分也一定相等,因此有:
![{\displaystyle {\frac {\sin(nx)}{(\sin x)^{n}}}=\left[{n \choose 1}\cot ^{n-1}x-{n \choose 3}\cot ^{n-3}x\pm \cdots \right]}](https://wikimedia.org/api/rest_v1/media/math/render/svg/07546ae81a253d62baee0fd3e17a5b063acdf1bd)
固定一个正整数m,设n = 2m + 1,并考虑xr = r π/(2m + 1)对于r = 1、2、……、m。那么nxr是π的倍数,因此是正弦函数的零点,所以:
![{\displaystyle 0={{2m+1} \choose 1}\cot ^{2m}x_{r}-{{2m+1} \choose 3}\cot ^{2m-2}x_{r}\pm \cdots +(-1)^{m}{{2m+1} \choose {2m+1}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/539b7f5e5f0f0ef1b0e0fa160fea00462a25ecb1)
对于所有的r = 1、2、……、m。x1、……、xm是区间(0, π/2)内不同的数。由于函数cot2 x在这个区间内是一一对应的,因此当r = 1、2、……、m时,tr = cot2 xr的值各不同。根据以上方程,这些m个"tr"是以下m次多项式的根:
![{\displaystyle p(t):={{2m+1} \choose 1}t^{m}-{{2m+1} \choose 3}t^{m-1}\pm \cdots +(-1)^{m}{{2m+1} \choose {2m+1}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc335adb7d4e1b8289d22333576795329d03eec5)
根据韦达定理,我们可以直接从这个多项式的头两项计算出所有根的和,因此:
![{\displaystyle \cot ^{2}x_{1}+\cot ^{2}x_{2}+\cdots +\cot ^{2}x_{m}={\frac {\binom {2m+1}{3}}{\binom {2m+1}{1}}}={\frac {2m(2m-1)}{6}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67363e9ae2b6f73a388e65572bd788ee4f988892)
把恒等式csc2 x = cot2 x + 1代入,可得:
![{\displaystyle \csc ^{2}x_{1}+\csc ^{2}x_{2}+\cdots +\csc ^{2}x_{m}={\frac {2m(2m-1)}{6}}+m={\frac {2m(2m+2)}{6}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d604cad116cf57eee0ded9c115f1ffc94ec869f)
现在考虑不等式cot2 x < 1/x2 < csc2 x。如果我们把对于xr = r π/(2m + 1)的所有不等式相加起来,并利用以上的两个恒等式,便可得到:
![{\displaystyle {\frac {2m(2m-1)}{6}}<\left({\frac {2m+1}{\pi }}\right)^{2}+\left({\frac {2m+1}{2\pi }}\right)^{2}+\cdots +\left({\frac {2m+1}{m\pi }}\right)^{2}<{\frac {2m(2m+2)}{6}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5fbd75e4d5f190274ece170a6c0bf5c8ac89e8c5)
把不等式乘以(π/(2m + 1))2,便得:
![{\displaystyle {\frac {\pi ^{2}}{6}}\left({\frac {2m}{2m+1}}\right)\left({\frac {2m-1}{2m+1}}\right)<{\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+\cdots +{\frac {1}{m^{2}}}<{\frac {\pi ^{2}}{6}}\left({\frac {2m}{2m+1}}\right)\left({\frac {2m+2}{2m+1}}\right)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1174cb0db421fbbae35089ce6faac41b5f6fed11)
当m趋于无穷大时,左面和右面的表达式都趋于
,因此根据夹挤定理,有:
![{\displaystyle \zeta (2)=\sum _{k=1}^{\infty }{\frac {1}{k^{2}}}=\lim _{m\to \infty }\left({\frac {1}{1^{2}}}+{\frac {1}{2^{2}}}+\cdots +{\frac {1}{m^{2}}}\right)={\frac {\pi ^{2}}{6}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cc4c3b614c578fcadc7d10a63e231f0f779ba513)
证毕。
首先考虑二重积分
的级数展开形式,其中区域D1为x∈(0,1)且y∈(0,1)的正方形区域:
![{\displaystyle {\begin{aligned}\iint \limits _{D_{1}}{\frac {\mathrm {d} x\mathrm {d} y}{1-xy}}&=\int _{0}^{1}\left(\int _{0}^{1}{\frac {\mathrm {d} x}{1-xy}}\right)\mathrm {d} y=\int _{0}^{1}\left[\int _{0}^{1}\sum _{n=1}^{\infty }(xy)^{n-1}\mathrm {d} x\right]\mathrm {d} y\\&=\int _{0}^{1}\left.\sum _{n=1}^{\infty }{\frac {x^{n}y^{n-1}}{n}}\right|_{x=0}^{x=1}\mathrm {d} y=\int _{0}^{1}\sum _{n=1}^{\infty }{\frac {y^{n-1}}{n}}\mathrm {d} y\\&=\left.\sum _{n=1}^{\infty }{\frac {y^{n}}{n^{2}}}\right|_{y=0}^{y=1}=\sum _{n=1}^{\infty }{\frac {1}{n^{2}}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f4513be7a7ae39314ebd855b68a17e19dc0bdb8)
接下来考虑做如下变换:
![{\displaystyle {\begin{cases}x=u\cos {\dfrac {\pi }{4}}-v\sin {\dfrac {\pi }{4}}={\dfrac {u-v}{\sqrt {2}}}\\y=u\sin {\dfrac {\pi }{4}}+v\cos {\dfrac {\pi }{4}}={\dfrac {u+v}{\sqrt {2}}}\\\end{cases}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4ce019d50d5ef4b6da433c427269fd1182a8eea9)
相当于将正方形区域D1旋转45°之后变成D2,但仍然保持其形状与面积,其中正方形D2的四个顶点在(u,v)下的坐标分别是
、
、
、
:
![{\displaystyle {\begin{aligned}\iint \limits _{D_{1}}{\frac {\mathrm {d} x\mathrm {d} y}{1-xy}}&=\iint \limits _{D_{2}}{\frac {\mathrm {d} u\mathrm {d} v}{1-{\frac {(u-v)(u+v)}{2}}}}=2\iint \limits _{D_{2}}{\frac {\mathrm {d} u\mathrm {d} v}{2-u^{2}+v^{2}}}\\&=2\int _{0}^{\frac {1}{\sqrt {2}}}\left(\int _{-u}^{u}{\frac {\mathrm {d} v}{2-u^{2}+v^{2}}}\right)\mathrm {d} u+2\int _{\frac {1}{\sqrt {2}}}^{\sqrt {2}}\left(\int _{u-{\sqrt {2}}}^{{\sqrt {2}}-u}{\frac {\mathrm {d} v}{2-u^{2}+v^{2}}}\right)\mathrm {d} u\\&=4\int _{0}^{\frac {1}{\sqrt {2}}}\left(\int _{0}^{u}{\frac {\mathrm {d} v}{2-u^{2}+v^{2}}}\right)\mathrm {d} u+4\int _{\frac {1}{\sqrt {2}}}^{\sqrt {2}}\left(\int _{0}^{{\sqrt {2}}-u}{\frac {\mathrm {d} v}{2-u^{2}+v^{2}}}\right)\mathrm {d} u\\&=I_{1}+I_{2}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/be62e29bd35dfdde55d4889820674697b4a2aef1)
最后一行用到的是偶函数的性质化简积分,并且将加号前后的两个积分(包括前面的系数)简记为I1与I2。先计算I1,下面变量代换设
,则有
:
![{\displaystyle {\begin{aligned}I_{1}&=4\int _{0}^{\frac {1}{\sqrt {2}}}\left.{\frac {\arctan {\frac {v}{\sqrt {2-u^{2}}}}}{\sqrt {2-u^{2}}}}\right|_{v=0}^{v=u}\mathrm {d} u=4\int _{0}^{\frac {1}{\sqrt {2}}}{\frac {\arctan {\frac {u}{\sqrt {2-u^{2}}}}}{\sqrt {2-u^{2}}}}\mathrm {d} u\\&=4\int _{0}^{\frac {\pi }{6}}{\frac {{\sqrt {2}}\cos \theta \cdot \arctan {\frac {{\sqrt {2}}\sin \theta }{\sqrt {2-2\sin ^{2}\theta }}}}{\sqrt {2-2\sin ^{2}\theta }}}\mathrm {d} \theta =4\int _{0}^{\frac {\pi }{6}}{\frac {{\sqrt {2}}\cos \theta \cdot \arctan {\frac {{\sqrt {2}}\sin \theta }{{\sqrt {2}}\cos \theta }}}{{\sqrt {2}}\cos \theta }}\mathrm {d} \theta \\&=4\int _{0}^{\frac {\pi }{6}}\arctan \tan \theta \mathrm {d} \theta =4\int _{0}^{\frac {\pi }{6}}\theta \mathrm {d} \theta =2\theta ^{2}{\bigg |}_{\theta =0}^{\theta ={\frac {\pi }{6}}}={\frac {\pi ^{2}}{18}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ca20f64775c92171a33068a65b92d0014291c5f5)
类似地,再计算I2,下面变量代换设
,则有
,并且交换积分上下限从负号变回正号:
![{\displaystyle {\begin{aligned}I_{2}&=4\int _{\frac {1}{\sqrt {2}}}^{\sqrt {2}}\left.{\frac {\arctan {\frac {v}{\sqrt {2-u^{2}}}}}{\sqrt {2-u^{2}}}}\right|_{v=0}^{v={\sqrt {2}}-u}\mathrm {d} u=\int _{\frac {1}{\sqrt {2}}}^{\sqrt {2}}{\frac {\arctan {\frac {{\sqrt {2}}-u}{\sqrt {2-u^{2}}}}}{\sqrt {2-u^{2}}}}\mathrm {d} u\\&=4\int _{0}^{\frac {\pi }{3}}{\frac {{\sqrt {2}}\sin \theta \cdot \arctan {\frac {{\sqrt {2}}-{\sqrt {2}}\cos \theta }{\sqrt {2-2\cos ^{2}\theta }}}}{\sqrt {2-2\cos ^{2}\theta }}}\mathrm {d} \theta =4\int _{0}^{\frac {\pi }{3}}{\frac {{\sqrt {2}}\sin \theta \cdot \arctan {\frac {{\sqrt {2}}(1-\cos \theta )}{{\sqrt {2}}\sin \theta }}}{{\sqrt {2}}\sin \theta }}\mathrm {d} \theta \\&=4\int _{0}^{\frac {\pi }{3}}\arctan \tan {\frac {\theta }{2}}\mathrm {d} \theta =4\int _{0}^{\frac {\pi }{3}}{\frac {\theta }{2}}\mathrm {d} \theta =\theta ^{2}{\bigg |}_{\theta =0}^{\theta ={\frac {\pi }{3}}}={\frac {\pi ^{2}}{9}}\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de2b1b9056939df05e8f5ff99bd2d33a1bc5b269)
最终可以得到:
![{\displaystyle \sum _{n=1}^{\infty }{\frac {1}{n^{2}}}=\iint \limits _{D_{1}}{\frac {\mathrm {d} x\mathrm {d} y}{1-xy}}=2\iint \limits _{D_{2}}{\frac {\mathrm {d} u\mathrm {d} v}{2-u^{2}+v^{2}}}=I_{1}+I_{2}={\frac {\pi ^{2}}{18}}+{\frac {\pi ^{2}}{9}}={\frac {\pi ^{2}}{6}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/57f32976a44ba9bb5d84bb242ee6b6d83bc04148)
证毕。
设有函数
,其定义域为
。这个函数的傅里叶级数是:
。
根据帕塞瓦尔恒等式,我们有:
![{\displaystyle {\pi ^{2} \over 3}={1 \over 2\pi }\int _{-\pi }^{\pi }f^{2}(x)\,dx=\sum _{n=1}^{\infty }{1 \over 2\pi }\int _{-\pi }^{\pi }(2{\frac {(-1)^{n+1}}{n}}\sin(nt))^{2}dt=2\sum _{n=1}^{\infty }{1 \over n^{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4de44f9857ad86374f4fa11b4d0baedcd1a55882)
因此
![{\displaystyle {\pi ^{2} \over 6}=\sum _{n=1}^{\infty }{1 \over n^{2}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e4ddfb14f950ef1dcfd623b9a92e47812aa47e88)
证毕。
- Weil, André, Number Theory: An Approach Through History, Springer-Verlag, 1983, ISBN 0-8176-3141-0 .
- Dunham, William, Euler: The Master of Us All, Mathematical Association of America, 1999, ISBN 0-88385-328-0 .
- Derbyshire, John, Prime Obsession: Bernhard Riemann and the Greatest Unsolved Problem in Mathematics, Joseph Henry Press, 2003, ISBN 0-309-08549-7 .
- Aigner, Martin; Ziegler, Günter M., Proofs from THE BOOK, Berlin, New York: Springer-Verlag, 1998
- Edwards, Harold M., Riemann's Zeta Function, Dover, 2001, ISBN 0-486-41740-9 .