黎曼ζ函數








複平面中一矩形區域之黎曼ζ函數ζ(z){displaystyle zeta (z)}zeta (z);此圖用Matplotlib程式繪圖產生,使用到定義域著色方法。[1]


黎曼ζ函數ζ(s)的定義如下:
設一複數s,其實數部份> 1而且:


ζ(s)=∑n=1∞1ns{displaystyle zeta (s)=sum _{n=1}^{infty }{frac {1}{n^{s}}}}zeta(s) =<br />
sum_{n=1}^infin frac{1}{n^s}

它亦可以用积分定义:


ζ(s)=1Γ(s)∫0∞xs−1ex−1dx{displaystyle zeta (s)={frac {1}{Gamma (s)}}int _{0}^{infty }{frac {x^{s-1}}{e^{x}-1}}mathrm {d} x} zeta(s) = frac{1}{Gamma(s)} int_{0}^{infty} frac{x ^ {s-1}}{e ^ x - 1} mathrm{d}x

在区域{s : Re(s) > 1}上,此无穷级数收敛并为一全纯函数(其中Re表示复数的实部,下同)。欧拉在1740年考虑过s为正整数的情况,后来切比雪夫拓展到s>1。[2]波恩哈德·黎曼认识到:ζ函数可以通过解析开拓来扩展到一个定义在复数域(s, s≠ 1)上的全纯函数ζ(s)。这也是黎曼猜想所研究的函数。


虽然黎曼的ζ函数被数学家认为主要和“最纯”的数学领域数论相关,它也出现在应用统计学(参看齊夫定律(Zipf's Law)和齊夫-曼德爾布羅特定律英语Zipf–Mandelbrot law(Zipf-Mandelbrot Law))、物理,以及调音的数学理论中。




目录






  • 1 历史


    • 1.1 奥里斯姆


    • 1.2 欧拉


    • 1.3 黎曼


    • 1.4 阿达马与普森


    • 1.5 希尔伯特


    • 1.6 玻尔与兰道


    • 1.7 哈代与李特尔伍德


    • 1.8 塞尔伯格




  • 2 解析延拓


  • 3 和数论函数的关系


  • 4 佩龙公式


  • 5 和素数的关系


    • 5.1 欧拉乘积


    • 5.2 更进一步的联系


      • 5.2.1 黎曼阶梯素数计数函数


      • 5.2.2 切比雪夫函数






  • 6 零点


  • 7 函数值


    • 7.1 当s为正整数


    • 7.2 s趨近於1


    • 7.3 负整数


    • 7.4 复数值


    • 7.5 幅角


    • 7.6 函数值表


    • 7.7 临界线上的数值计算




  • 8 參考資料


  • 9 相關條目





历史



奥里斯姆


ζ函数最早出现于1350年左右,当时的尼克尔·奥里斯姆发现了调和级数发散,即
ζ(1)=1+12+13+14+...→{displaystyle zeta (1)=1+{frac {1}{2}}+{frac {1}{3}}+{frac {1}{4}}+...to infty }zeta (1)=1+{frac  {1}{2}}+{frac  {1}{3}}+{frac  {1}{4}}+...to infty







欧拉




第n个调和数(蓝点)与Log(n)+γ(红线)的图像


之后的一次进展来自莱昂哈德·欧拉,他给出了调和级数呈对数发散。






除此之外,他还在1735年给出了巴塞尔问题的解答,得到
ζ(2)=π26{displaystyle zeta (2)={frac {pi ^{2}}{6}}}zeta (2)={frac  {pi ^{2}}{6}}
的结果。欧拉最初的证明可以在巴塞尔问题中看到,然而那是他的第一个证明,因而广为人知。
事实上,那个证明虽有不严谨之处,但是欧拉仍然有自己的严格证明。[4]







欧拉在1737年还发现了欧拉乘积公式:
n=1∞1ns=∏p(1−1ps)−1{displaystyle sum _{n=1}^{infty }{frac {1}{n^{s}}}=prod _{p}(1-{frac {1}{p^{s}}})^{-1}}sum _{{n=1}}^{{infty }}{frac  {1}{n^{s}}}=prod _{{p}}(1-{frac  {1}{p^{s}}})^{{-1}}
这是ζ函数与素数的联系的朦胧征兆,其证明可以在证明黎曼ζ函数的欧拉乘积公式中看到。
通过这条公式,容易证明当 Re⁡(s)>1{displaystyle {begin{smallmatrix}operatorname {Re} (s)>1end{smallmatrix}}}{displaystyle {begin{smallmatrix}operatorname {Re} (s)>1end{smallmatrix}}} 时,ζ(s)>0{displaystyle {begin{smallmatrix}zeta (s)>0end{smallmatrix}}}{begin{smallmatrix}zeta (s)>0end{smallmatrix}}



1749年,欧拉通过大胆的计算發現了[5]
ζ(−1)=1+2+3+4+5+...=−112{displaystyle zeta (-1)=1+2+3+4+5+...=-{frac {1}{12}}}zeta (-1)=1+2+3+4+5+...=-{frac  {1}{12}}
ζ(−2)=12+22+32+42+52+...=0{displaystyle zeta (-2)=1^{2}+2^{2}+3^{2}+4^{2}+5^{2}+...=0}zeta (-2)=1^{2}+2^{2}+3^{2}+4^{2}+5^{2}+...=0
ζ(−3)=13+23+33+43+53+...=1120{displaystyle zeta (-3)=1^{3}+2^{3}+3^{3}+4^{3}+5^{3}+...={frac {1}{120}}}zeta (-3)=1^{3}+2^{3}+3^{3}+4^{3}+5^{3}+...={frac  {1}{120}}
发现ζ(s)与ζ(1-s)之间存在某些关系。



黎曼





波恩哈德·黎曼对ζ解析延拓,用于素数的分布理论


将欧拉所做的一切牢牢地置于坚石之上的是黎曼,他在1859年的论文论小于给定数值的素数个数英语On_the_Number_of_Primes_Less_Than_a_Given_Magnitude以及未发表的手稿中做出了多项进展:[6]



  • 第一积分表示: ζ(s)=1Γ(s)∫0∞xs−1ex−1dx{displaystyle zeta (s)={frac {1}{Gamma (s)}}int _{0}^{infty }{frac {x^{s-1}}{e^{x}-1}},mathrm {d} x}zeta (s)={frac  {1}{Gamma (s)}}int _{{0}}^{{infty }}{frac  {x^{{s-1}}}{e^{x}-1}},{mathrm  {d}}x

  • 完备化的ζ,即黎曼ξ函数: ξ(s)=πs2Γ(s2)ζ(s){displaystyle xi (s)=pi ^{-{frac {s}{2}}}Gamma left({frac {s}{2}}right)zeta (s)}xi(s) = pi^{-frac{s}{2}} Gammaleft(frac{s}{2}right)zeta(s) ,满足函数方程 ξ(s)=ξ(1−s){displaystyle xi (s)=xi (1-s)}xi (s)=xi (1-s)

  • 第二积分表示: φ(x)=∑n=1∞e−πn2x{displaystyle varphi (x)=sum _{n=1}^{infty }e^{-pi n^{2}x}}varphi (x)=sum _{{n=1}}^{{infty }}e^{{-pi n^{2}x}} ,则 ξ(s)=∫0∞φ(x)xs2−1dx{displaystyle xi (s)=int _{0}^{infty }varphi (x)x^{{frac {s}{2}}-1},mathrm {d} x}xi (s)=int _{{0}}^{{infty }}varphi (x)x^{{{frac  {s}{2}}-1}},{mathrm  {d}}x


  • 黎曼 - 冯·曼戈尔特公式英语Riemann–von Mangoldt formula:以0<N⁡(T){displaystyle {begin{smallmatrix}0<operatorname {N} (T)end{smallmatrix}}}{begin{smallmatrix}0<operatorname {N}(T)end{smallmatrix}}表示虚部介于0与T之间的非平凡零点数量,则 N(T)=T2πlog⁡T2πT2π+O(log⁡T){displaystyle N(T)={frac {T}{2pi }}log {frac {T}{2pi }}-{frac {T}{2pi }}+mathrm {O} (log T)}N(T)={frac  {T}{2pi }}log {frac  {T}{2pi }}-{frac  {T}{2pi }}+mathrm{O} (log T)


  • 黎曼猜想:ζ函数的所有非平凡零点的实部非常有可能均为12{displaystyle {begin{smallmatrix}{frac {1}{2}}end{smallmatrix}}}{begin{smallmatrix}{frac  {1}{2}}end{smallmatrix}}


  • 第三积分表示: ζ(s)=12π(1−s)∮γzs−1ez1−ezdz{displaystyle zeta (s)={frac {1}{2pi i}}Gamma (1-s)oint _{gamma }{frac {{z^{s-1}}{e^{z}}}{1-{e^{z}}}},mathrm {d} z}zeta (s)={frac  {1}{{2pi i}}}Gamma (1-s)oint _{gamma }{{frac  {{{z^{{s-1}}}{e^{z}}}}{{1-{e^{z}}}}}},{mathrm  {d}}z ,其中围道γ逆时针环绕负实轴



第三积分表示的围道γ




  • 黎曼-西格尔公式:给出计算ξ函数的数值的方法

  • 零点的计算:计算了虚部介于0与100的所有零点的数值

  • 素数的分布公式:引入黎曼素数计数函数,给出了它与ζ函数的关系



阿达马与普森




ζ(1+it)的图像,蓝色为实部,黄色为虚部


1896年,雅克·阿达马与普森几乎同时地证明了ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}的所有非平凡零点的实部均小于1,即Re⁡(s)=1{displaystyle {begin{smallmatrix}operatorname {Re} (s)=1end{smallmatrix}}}{displaystyle {begin{smallmatrix}operatorname {Re} (s)=1end{smallmatrix}}}上无非平凡零点,从而完成了素数定理的证明。



希尔伯特


1900年,希尔伯特在巴黎的第二届国际数学家大会上作了题为《数学问题》的演讲,提出了23道最重要的数学问题,黎曼假设在其中作为第8题出现。
之后,希尔伯特提出了希尔伯特-波利亚猜想,具体时间及场合未知。



玻尔与兰道




虚部介于0与T的零点数量(蓝点)与黎曼-冯·曼格尔特公式(红线)的图像


1914年,哈那德·玻爾和愛德蒙·蘭道证明了玻爾-蘭道定理:含有临界线的任意带状区域都几乎包含了ζ的所有非平凡零点,表明了临界线为零点汇聚的“中心位置”。



哈代与李特尔伍德


1921年,哈代和李特尔伍德证明了存在常数T,使临界线上虚部位于0与T之间的非平凡零点的数量至少为KT{displaystyle {begin{smallmatrix}KTend{smallmatrix}}}{begin{smallmatrix}KTend{smallmatrix}}



塞尔伯格


1942年,阿特勒·塞尔伯格更进一步,证明了存在常数T,使临界线上虚部位于0与T之间的非平凡零点的数量至少为KTlog⁡T{displaystyle {begin{smallmatrix}KTlog Tend{smallmatrix}}}{begin{smallmatrix}KTlog Tend{smallmatrix}},这意味着ζ函数在临界线上的非平凡零点在所有零点中占有一个正密度,而临界线Re⁡(s)=12{displaystyle {begin{smallmatrix}operatorname {Re} (s)={frac {1}{2}}end{smallmatrix}}}{begin{smallmatrix}operatorname {Re}(s)={frac  {1}{2}}end{smallmatrix}}对于临界带0<Re⁡(s)<1{displaystyle {begin{smallmatrix}0<operatorname {Re} (s)<1end{smallmatrix}}}{begin{smallmatrix}0<operatorname {Re}(s)<1end{smallmatrix}}的测度为0。



解析延拓




对ζ函数解析延拓时使用的围道


ζ函数原本定义在右半平面Re⁡s>1{displaystyle {begin{smallmatrix}operatorname {Re} s>1end{smallmatrix}}}{begin{smallmatrix}operatorname {Re}s>1end{smallmatrix}}上,并且在此区域内为全纯函数



ζ(s)=∑n=1∞1ns=1Γ(s)∫0∞xs−1ex−1dx{displaystyle zeta (s)=sum _{n=1}^{infty }{frac {1}{n^{s}}}={frac {1}{Gamma (s)}}int _{0}^{infty }{frac {x^{s-1}}{e^{x}-1}},mathrm {d} x}zeta (s)=sum _{{n=1}}^{infty }{frac  {1}{n^{s}}}={frac  {1}{Gamma (s)}}int _{{0}}^{{infty }}{frac  {x^{{s-1}}}{e^{x}-1}},{mathrm  {d}}x (Re⁡s>1){displaystyle (operatorname {Re} s>1)}(operatorname {Re}s>1)

解析延拓后在全局具有积分表达式



ζ(s)=12π(1−s)∮γzs−1ez1−ezdz{displaystyle zeta (s)={frac {1}{2pi i}}Gamma (1-s)oint _{gamma }{frac {{z^{s-1}}{e^{z}}}{1-{e^{z}}}},mathrm {d} z}zeta (s)={frac  {1}{{2pi i}}}Gamma (1-s)oint _{gamma }{{frac  {{{z^{{s-1}}}{e^{z}}}}{{1-{e^{z}}}}}},{mathrm  {d}}z


满足函数方程



ζ(1−s)=2(2π)−(s)cos⁡s2)ζ(s){displaystyle zeta (1-s)=2(2pi )^{-s}Gamma (s)cos left({frac {pi s}{2}}right)zeta (s)} zeta(1-s) = 2 (2 pi )^{-s} Gamma(s) cosleft(frac{pi s}{2}right) zeta(s)


特别地,如果考虑正规化的ζ,即黎曼ξ函数



ξ(s)=πs2Γ(s2)ζ(s){displaystyle xi (s)=pi ^{-{frac {s}{2}}}Gamma left({frac {s}{2}}right)zeta (s)}xi(s) = pi^{-frac{s}{2}} Gammaleft(frac{s}{2}right)zeta(s)


那么它满足函数方程


ξ(s)=ξ(1−s){displaystyle xi (s)=xi (1-s)}xi (s)=xi (1-s)


和数论函数的关系


黎曼ζ函数可看做是具有如下形式的级数的一个特例:


F⁡(s)=∑n=1∞f(n)ns{displaystyle operatorname {F} (s)=sum _{n=1}^{infty }{frac {f(n)}{n^{s}}}}operatorname {F}(s)=sum _{{n=1}}^{{infty }}{frac  {f(n)}{n^{s}}}

这种类型的级数被称作狄利克雷级数。当f为狄利克雷特征时,又称作狄利克雷L函数,也有与黎曼猜想相应的广义黎曼猜想


为了方便对数论函数作讨论,此处引入狄利克雷卷积 f∗g{displaystyle {begin{smallmatrix}f*gend{smallmatrix}}}{begin{smallmatrix}f*gend{smallmatrix}}
(f∗g)(n)=∑pq=nf(p)g(q){displaystyle (f*g)(n)=sum _{{text{pq}}=n}f(p)g(q)}(f*g)(n)=sum _{{{text{pq}}=n}}f(p)g(q)


F⁡(s)=∑n=1∞f(n)ns{displaystyle operatorname {F} (s)=sum _{n=1}^{infty }{frac {f(n)}{n^{s}}}}operatorname {F}(s)=sum _{{n=1}}^{{infty }}{frac  {f(n)}{n^{s}}}G⁡(s)=∑n=1∞g(n)ns{displaystyle operatorname {G} (s)=sum _{n=1}^{infty }{frac {g(n)}{n^{s}}}}operatorname {G}(s)=sum _{{n=1}}^{{infty }}{frac  {g(n)}{n^{s}}}
于是显然 F⁡(s)G⁡(s)=∑n=1∞(f∗g)(n)ns{displaystyle operatorname {F} (s)operatorname {G} (s)=sum _{n=1}^{infty }{frac {(f*g)(n)}{n^{s}}}}operatorname {F}(s)operatorname {G}(s)=sum _{{n=1}}^{{infty }}{frac  {(f*g)(n)}{n^{s}}}






于是,如果数论函数h=1∗g{displaystyle {begin{smallmatrix}h=1*gend{smallmatrix}}}{begin{smallmatrix}h=1*gend{smallmatrix}},亦即 h(n)=∑d‖ng(d){displaystyle h(n)=sum _{d|n}g(d)}h(n)=sum _{{d|n}}g(d) (此时,h(n){displaystyle {begin{smallmatrix}h(n)end{smallmatrix}}}{begin{smallmatrix}h(n)end{smallmatrix}}g(d){displaystyle {begin{smallmatrix}g(d)end{smallmatrix}}}{begin{smallmatrix}g(d)end{smallmatrix}}可通过默比乌斯反演公式相互转换)
那么 H⁡(s)=∑n=1∞h(n)ns=ζ(s)∑n=1∞g(n)ns{displaystyle operatorname {H} (s)=sum _{n=1}^{infty }{frac {h(n)}{n^{s}}}=zeta (s)sum _{n=1}^{infty }{frac {g(n)}{n^{s}}}}operatorname {H}(s)=sum _{{n=1}}^{{infty }}{frac  {h(n)}{n^{s}}}=zeta (s)sum _{{n=1}}^{{infty }}{frac  {g(n)}{n^{s}}}
通常两侧的求和有一个是相对简单的函数,或是和ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}直接相关的函数
如果对g(n){displaystyle {begin{smallmatrix}g(n)end{smallmatrix}}}{begin{smallmatrix}g(n)end{smallmatrix}}的求和较简单,可以将h(n){displaystyle {begin{smallmatrix}h(n)end{smallmatrix}}}{begin{smallmatrix}h(n)end{smallmatrix}}ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}相联系,反之可以将g(n){displaystyle {begin{smallmatrix}g(n)end{smallmatrix}}}{begin{smallmatrix}g(n)end{smallmatrix}}ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}相联系
n=1∞g(n)ns=H⁡(s)ζ(s){displaystyle sum _{n=1}^{infty }{frac {g(n)}{n^{s}}}={frac {operatorname {H} (s)}{zeta (s)}}}sum _{{n=1}}^{{infty }}{frac  {g(n)}{n^{s}}}={frac  {operatorname {H}(s)}{zeta (s)}}
如下表所示:













































目标函数名 g(n) h(n) G(s)或H(s) g(n)或h(n)与ζ函数的联系
莫比乌斯函数 μ(n)=μ(p1a1p2a2...pkak)={(−1)ka1=a2=...=ak0otherwise{displaystyle mu (n)=mu (p_{1}^{a_{1}}p_{2}^{a_{2}}...p_{k}^{a_{k}})={begin{cases}(-1)^{k}&a_{1}=a_{2}=...=a_{k}\0&mathrm {otherwise} end{cases}}}mu (n)=mu (p_{1}^{{a_{1}}}p_{2}^{{a_{2}}}...p_{k}^{{a_{k}}})={begin{cases}(-1)^{k}&a_{1}=a_{2}=...=a_{k}\0&{mathrm  {otherwise}}end{cases}} ⌊1n⌋{displaystyle leftlfloor {frac {1}{n}}rightrfloor }leftlfloor {frac  {1}{n}}rightrfloor H⁡(s)=1{displaystyle operatorname {H} (s)=1}operatorname {H}(s)=1
G⁡(s)=∑n=1∞μ(n)ns=1ζ(s){displaystyle operatorname {G} (s)=sum _{n=1}^{infty }{frac {mu (n)}{n^{s}}}={frac {1}{zeta (s)}}}operatorname {G}(s)=sum _{{n=1}}^{{infty }}{frac  {mu (n)}{n^{s}}}={frac  {1}{zeta (s)}}
欧拉函数 φ(n)=Card⁡{k|k<n,(k,n)=1}{displaystyle varphi (n)=operatorname {Card} {k,,|k<n,,(k,n)=1}quad }varphi (n)=operatorname {Card}{k,,|k<n,,(k,n)=1}quad n{displaystyle n}n H⁡(s)=ζ(s−1){displaystyle operatorname {H} (s)=zeta (s-1)}operatorname {H}(s)=zeta (s-1)
G⁡(s)=∑n=1∞φ(n)ns=ζ(s−1)ζ(s){displaystyle operatorname {G} (s)=sum _{n=1}^{infty }{frac {varphi (n)}{n^{s}}}={frac {zeta (s-1)}{zeta (s)}}}operatorname {G}(s)=sum _{{n=1}}^{{infty }}{frac  {varphi (n)}{n^{s}}}={frac  {zeta (s-1)}{zeta (s)}}
除数函数 {displaystyle n^{alpha }}n^{{alpha }} σα(n)=∑d|ndα{displaystyle sigma _{alpha }(n)=sum _{d|n}d^{alpha }}sigma _{alpha }(n)=sum _{{d|n}}d^{{alpha }} G⁡(s)=ζ(s−α){displaystyle operatorname {G} (s)=zeta (s-alpha )}operatorname {G}(s)=zeta (s-alpha )
H⁡(s)=∑n=1∞σα(n)ns=ζ(s−α(s){displaystyle operatorname {H} (s)=sum _{n=1}^{infty }{frac {sigma _{alpha }(n)}{n^{s}}}=zeta (s-alpha )zeta (s)}operatorname {H}(s)=sum _{{n=1}}^{{infty }}{frac  {sigma _{alpha }(n)}{n^{s}}}=zeta (s-alpha )zeta (s)
刘维尔函数 μ(n)=μ(p1a1p2a2...pkak)=a1+a2+...+ak{displaystyle mu (n)=mu (p_{1}^{a_{1}}p_{2}^{a_{2}}...p_{k}^{a_{k}})=a_{1}+a_{2}+...+a_{k}}mu (n)=mu (p_{1}^{{a_{1}}}p_{2}^{{a_{2}}}...p_{k}^{{a_{k}}})=a_{1}+a_{2}+...+a_{k} {1n=m20otherwise{displaystyle {begin{cases}1&n=m^{2}\0&mathrm {otherwise} end{cases}}}{begin{cases}1&n=m^{2}\0&{mathrm  {otherwise}}end{cases}} H⁡(s)=ζ(2s){displaystyle operatorname {H} (s)=zeta (2s)}operatorname {H}(s)=zeta (2s)
G⁡(s)=∑n=1∞λ(n)ns=ζ(2s)ζ(s){displaystyle operatorname {G} (s)=sum _{n=1}^{infty }{frac {lambda (n)}{n^{s}}}={frac {zeta (2s)}{zeta (s)}}}operatorname {G}(s)=sum _{{n=1}}^{{infty }}{frac  {lambda (n)}{n^{s}}}={frac  {zeta (2s)}{zeta (s)}}
冯·曼戈尔特函数 Λ(n)={log⁡(p)n=pk0otherwise{displaystyle Lambda (n)={begin{cases}log(p)&n=p^{k}\0&mathrm {otherwise} end{cases}}}Lambda (n)={begin{cases}log(p)&n=p^{k}\0&{mathrm  {otherwise}}end{cases}} log⁡n{displaystyle log n}log n H⁡(s)=ζ′(s){displaystyle operatorname {H} (s)=zeta '(s)}operatorname {H}(s)=zeta '(s)
G⁡(s)=∑n=1∞Λ(n)ns=−ζ′(s)ζ(s){displaystyle operatorname {G} (s)=sum _{n=1}^{infty }{frac {Lambda (n)}{n^{s}}}=-{frac {zeta '(s)}{zeta (s)}}}operatorname {G}(s)=sum _{{n=1}}^{{infty }}{frac  {Lambda (n)}{n^{s}}}=-{frac  {zeta '(s)}{zeta (s)}}


佩龙公式


ζ函数与数论函数存在的联系可以通过佩龙公式转化为它和数论函数的求和的关系:设



G(s)=∑n=1∞g(n){displaystyle G(s)={sum _{n=1}^{infty }}g(n)}G(s)={sum _{{n=1}}^{infty }}g(n)


则由佩龙公式,



A(x)=∑n≤x′g(n)=12πi∫c−i∞c+i∞G(z)xzzdz{displaystyle A(x)={sum _{nleq x}}'g(n)={frac {1}{2pi i}}int _{c-iinfty }^{c+iinfty }G(z){frac {x^{z}}{z}},mathrm {d} z}A(x)={sum _{{nleq x}}}'g(n)={frac  {1}{2pi i}}int _{{c-iinfty }}^{{c+iinfty }}G(z){frac  {x^{{z}}}{z}},{mathrm  {d}}z


其中右上角的'表示如果x是整数,那么求和的最后一项要乘以12{displaystyle {begin{smallmatrix}{frac {1}{2}}end{smallmatrix}}}{begin{smallmatrix}{frac  {1}{2}}end{smallmatrix}}
这样做的其中一个结果就是ζ函数和素数分布的关系。



和素数的关系



欧拉乘积



此函数和素数的关系已由欧拉所揭示:


ζ(s)=∏p11−p−s{displaystyle zeta (s)=prod _{p}{frac {1}{1-p^{-s}}}}zeta(s) = prod_{p} frac{1}{1-p^{-s}}

这是一个延展到所有的质数p的无穷乘积,被称为欧拉乘积。这是几何级数的公式和算术基本定理的一个结果。
如果对上式取对数,则可得到


log⁡ζ(s)=∑n=1∞1n∑pp−ns{displaystyle log zeta (s)=sum _{n=1}^{infty }{frac {1}{n}}sum _{p}p^{-ns}}log zeta (s)=sum _{{n=1}}^{{infty }}{frac  {1}{n}}sum _{p}p^{{-ns}}






更进一步的联系



黎曼阶梯素数计数函数




黎曼素数计数函数(蓝色)J(x)与对数积分(金色)Li(x)的图像,x<300




黎曼素数计数函数(蓝点)J(x)与对数积分(红线)Li(x)的图像,x<1 000 000


可以使用黎曼素数计数函数J⁡(x){displaystyle {begin{smallmatrix}operatorname {J} (x)end{smallmatrix}}}{begin{smallmatrix}operatorname {J}(x)end{smallmatrix}}建立ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}与素数分布的进一步联系,这也是黎曼在他的论文论小于给定数值的素数个数英语On_the_Number_of_Primes_Less_Than_a_Given_Magnitude中使用的函数,定义如下:


J⁡(x)=∑n≤(n){displaystyle operatorname {J} (x)=sum _{nleq x}kappa (n)}operatorname {J}(x)=sum _{{nleq x}}kappa (n)

其中κ(n)={1kn=pk0otherwise{displaystyle kappa (n)={begin{cases}{frac {1}{k}}&n=p^{k}\0&mathrm {otherwise} end{cases}}}kappa (n)={begin{cases}{frac  {1}{k}}&n=p^{k}\0&{mathrm  {otherwise}}end{cases}}
那么可以建立J⁡(x){displaystyle {begin{smallmatrix}operatorname {J} (x)end{smallmatrix}}}{begin{smallmatrix}operatorname {J}(x)end{smallmatrix}}ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}的零点ρ的联系,称为黎曼显式公式英语Explicit formulae (L-function)


J⁡(x)=12πi∫c−i∞c+i∞log⁡ζ(s)xssds=Li⁡(x)−ρLi⁡(xρ)+∫x∞1t(t2−1)log⁡(t)dx−log⁡2{displaystyle {begin{aligned}operatorname {J} (x)&={frac {1}{2pi i}}int _{c-iinfty }^{c+iinfty }log zeta (s){frac {x^{s}}{s}},mathrm {d} s\&=operatorname {Li} (x)-sum _{rho }operatorname {Li} (x^{rho })+int _{x}^{infty }{frac {1}{t(t^{2}-1)log(t)}},mathrm {d} x-log 2\end{aligned}}}{begin{aligned}operatorname {J}(x)&={frac  {1}{{2pi i}}}int _{{c-iinfty }}^{{c+iinfty }}log zeta (s){frac  {x^{s}}{s}},{mathrm  {d}}s\&=operatorname {Li}(x)-sum _{{rho }}operatorname {Li}(x^{rho })+int _{{x}}^{{infty }}{frac  {1}{t(t^{2}-1)log(t)}},{mathrm  {d}}x-log 2\end{aligned}}





J⁡(x){displaystyle {begin{smallmatrix}operatorname {J} (x)end{smallmatrix}}}{begin{smallmatrix}operatorname {J}(x)end{smallmatrix}}π(x){displaystyle {begin{smallmatrix}pi (x)end{smallmatrix}}}{begin{smallmatrix}pi (x)end{smallmatrix}}的联系可以通过莫比乌斯反演公式完成。
π(x)=∑n=1∞μ(n)nJ⁡(x)=J⁡(x)+O(xlog⁡log⁡x){displaystyle pi (x)=sum _{n=1}^{infty }{frac {mu (n)}{n}}operatorname {J} (x)=operatorname {J} (x)+mathrm {O} ({sqrt {x}}log log x)}pi (x)=sum _{{n=1}}^{{infty }}{frac  {mu (n)}{n}}operatorname {J}(x)=operatorname {J}(x)+mathrm{O} ({sqrt  x}log log x)
然而J⁡(x){displaystyle {begin{smallmatrix}operatorname {J} (x)end{smallmatrix}}}{begin{smallmatrix}operatorname {J}(x)end{smallmatrix}}的表达式过于复杂,如下的切比雪夫函数英语Chebyshev function更为常用。



切比雪夫函数




第二切比雪夫函数(蓝线)ψ(x)与y=x(金线)的图像,x<300




第二切比雪夫函数(蓝点)ψ(x)与y=x(红线)的图像,x<1 000 000


第一切比雪夫函数ϑ(x){displaystyle {begin{smallmatrix}vartheta (x)end{smallmatrix}}}{begin{smallmatrix}vartheta (x)end{smallmatrix}}定义为


ϑ(x)=∑p≤xlog⁡p{displaystyle vartheta (x)=sum _{pleq x}log p}vartheta (x)=sum _{{pleq x}}log p

而更常用的第二切比雪夫函数ψ(x){displaystyle {begin{smallmatrix}psi (x)end{smallmatrix}}}{begin{smallmatrix}psi (x)end{smallmatrix}}定义为


ψ(x)=∑pk≤xlog⁡p=∑n≤(n)=∑p≤x⌊logp⁡x⌋log⁡p{displaystyle psi (x)=sum _{p^{k}leq x}log p=sum _{nleq x}Lambda (n)=sum _{pleq x}lfloor log _{p}xrfloor log p}psi (x)=sum _{{p^{k}leq x}}log p=sum _{{nleq x}}Lambda (n)=sum _{{pleq x}}lfloor log _{p}xrfloor log p

其中,如前文定义的 Λ(n)={log⁡(p)n=pk0otherwise{displaystyle Lambda (n)={begin{cases}log(p)&n=p^{k}\0&mathrm {otherwise} end{cases}}}Lambda (n)={begin{cases}log(p)&n=p^{k}\0&{mathrm  {otherwise}}end{cases}}
第二切比雪夫函数与第一切比雪夫函数的关系,可看做“等同于”黎曼素数计数函数与素数计数函数的关系。
第二切比雪夫函数ψ(x){displaystyle {begin{smallmatrix}psi (x)end{smallmatrix}}}{begin{smallmatrix}psi (x)end{smallmatrix}}ζ(s){displaystyle {begin{smallmatrix}zeta (s)end{smallmatrix}}}{begin{smallmatrix}zeta (s)end{smallmatrix}}的零点ρ有如下的联系


ψ(x)=12πi∫c−i∞c+i∞ζ′(s)ζ(s)xssds=∑n≤(n)=x−ρρ12log⁡(1−1x2)−log⁡(2π){displaystyle {begin{aligned}psi (x)&={frac {1}{2pi i}}int _{c-iinfty }^{c+iinfty }-{frac {zeta '(s)}{zeta (s)}}{frac {x^{s}}{s}},mathrm {d} s\&=sum _{nleq x}Lambda (n)=x-sum _{rho }{frac {x^{rho }}{rho }}-{frac {1}{2}}log(1-{frac {1}{x^{2}}})-log(2pi )\end{aligned}}}{begin{aligned}psi (x)&={frac  {1}{{2pi i}}}int _{{c-iinfty }}^{{c+iinfty }}-{frac  {zeta '(s)}{zeta (s)}}{frac  {x^{s}}{s}},{mathrm  {d}}s\&=sum _{{nleq x}}Lambda (n)=x-sum _{{rho }}{frac  {x^{rho }}{rho }}-{frac  {1}{2}}log(1-{frac  {1}{x^{2}}})-log(2pi )\end{aligned}}

ψ(x){displaystyle {begin{smallmatrix}psi (x)end{smallmatrix}}}{begin{smallmatrix}psi (x)end{smallmatrix}}J⁡(x){displaystyle {begin{smallmatrix}operatorname {J} (x)end{smallmatrix}}}{begin{smallmatrix}operatorname {J}(x)end{smallmatrix}}的联系可以通过阿贝尔求和公式:



ψ(x)=∑n=pk≤xlog⁡p=ψ(x)=∑n=pk≤x1klog⁡n=∑n≤(n)log⁡n{displaystyle psi (x)=sum _{n=p^{k}leq x}log p=psi (x)=sum _{n=p^{k}leq x}{frac {1}{k}}log n=sum _{nleq x}{frac {kappa (n)}{log n}}}psi (x)=sum _{{n=p^{k}leq x}}log p=psi (x)=sum _{{n=p^{k}leq x}}{frac  {1}{k}}log n=sum _{{nleq x}}{frac  {kappa (n)}{log n}}


其中κ如前文所定义,则由阿贝尔求和公式


J⁡(x)=∑n≤(n)=ψ(x)log⁡x+∫2xψ(t)t2log⁡tdt=ψ(x)log⁡x+O(xlog2⁡x){displaystyle operatorname {J} (x)=sum _{nleq x}kappa (n)={frac {psi (x)}{log x}}+int _{2}^{x}{frac {psi (t)}{t^{2}log t}},mathrm {d} t={frac {psi (x)}{log x}}+mathrm {O} ({frac {x}{log ^{2}x}})}operatorname {J}(x)=sum _{{nleq x}}kappa (n)={frac  {psi (x)}{log x}}+int _{{2}}^{{x}}{frac  {psi (t)}{t^{2}log t}},{mathrm  {d}}t={frac  {psi (x)}{log x}}+mathrm{O} ({frac  {x}{log ^{2}x}})


零点


解析延拓之后的ζ函数具有零点,他们分别是分布有序的平凡零点(所有负偶数),以及临界带0<Re⁡s<1{displaystyle {begin{smallmatrix}0<operatorname {Re} s<1end{smallmatrix}}}{begin{smallmatrix}0<operatorname {Re}s<1end{smallmatrix}}内的非平凡零点。
N⁡(T){displaystyle {begin{smallmatrix}operatorname {N} (T)end{smallmatrix}}}{begin{smallmatrix}operatorname {N}(T)end{smallmatrix}}表示虚部介于0与T之间的非平凡零点数量,则0<N⁡(T){displaystyle {begin{smallmatrix}0<operatorname {N} (T)end{smallmatrix}}}{begin{smallmatrix}0<operatorname {N}(T)end{smallmatrix}}遵循黎曼 - 冯·曼戈尔特公式英语Riemann–von Mangoldt formulaN(T)=T2πlog⁡T2πT2π+O(log⁡T){displaystyle N(T)={frac {T}{2pi }}log {frac {T}{2pi }}-{frac {T}{2pi }}+mathrm {O} (log T)}N(T)={frac  {T}{2pi }}log {frac  {T}{2pi }}-{frac  {T}{2pi }}+mathrm{O} (log T)



函数值




黎曼函数在s > 1的情况


ζ函数满足如下函数方程:


ζ(s)=2sπs−1sin⁡s2)Γ(1−s)ζ(1−s){displaystyle zeta (s)=2^{s}pi ^{s-1}sin left({frac {pi s}{2}}right)Gamma (1-s)zeta (1-s)}zeta(s) = 2^spi^{s-1}sinleft(frac{pi s}{2}right)Gamma(1-s)zeta(1-s)

对于所有C{0,1}中的s成立。这裡,Γ表示Γ函数。这个公式原来用来构造解析连续性。在s = 1,ζ函数有一个简单极点其留数为1。上述方程中有sin函數,sin⁡s2){displaystyle sin left({frac {pi s}{2}}right)}sinleft(frac{pi s}{2}right)的零點為偶數s = 2n,這些位置是可能的零點,但s為正偶數時,sin⁡s2)Γ(1−s){displaystyle sin left({frac {pi s}{2}}right)Gamma (1-s)}sinleft(frac{pi s}{2}right)Gamma(1-s)為不為零的規則函數英语Regular function,只有s為負偶數時,ζ函数才有零點,稱為平凡零點。



当s为正整数


欧拉计算出ζ(2k),对于偶整数2k,使用公式


ζ(2k)=B2k(−1)k+1(2π)2k2(2k)!{displaystyle zeta (2k)={frac {B_{2k}(-1)^{k+1}(2pi )^{2k}}{2(2k)!}}}zeta(2k) = frac{B_{2k}(-1)^{k+1}(2pi)^{2k}}{2(2k)!}

其中B2k是伯努利数。从这个,我们可以看到ζ(2) = π2/6, ζ(4) = π4/90, ζ(6) = π6/945等等。(OEIS中的序列A046988/A002432)。这些给出了著名的π的无穷级数。奇整数的情况没有这么简单。拉马努金在这上面做了很多了不起的工作。
s{displaystyle s,}s,为正偶数时的函数值公式已经由欧拉计算出。但当s{displaystyle s,}s,为正奇数时,尚未找到封闭式。



ζ(1)=1+12+13+⋯=∞;{displaystyle zeta (1)=1+{frac {1}{2}}+{frac {1}{3}}+cdots =infty ;!}zeta(1) = 1 + frac{1}{2} + frac{1}{3} + cdots = infty;!
这是调和级数。



ζ(32)≈2.612;{displaystyle zeta left({frac {3}{2}}right)approx 2.612;!}zetaleft(frac{3}{2}right) approx 2.612;!   OEIS A078434
该值用于计算具有周期性边界条件的玻色-爱因斯坦凝聚的临界温度以及磁系统的自旋波物理。



ζ(2)=1+122+132+⋯26≈1.645;{displaystyle zeta (2)=1+{frac {1}{2^{2}}}+{frac {1}{3^{2}}}+cdots ={frac {pi ^{2}}{6}}approx 1.645;!}zeta(2) = 1 + frac{1}{2^2} + frac{1}{3^2} + cdots = frac{pi^2}{6} approx 1.645;!   OEIS A013661
即巴塞尔问题。这个结果的倒数回答了这个问题:随机选取两个数字而互质的概率是多少?[7]



ζ(3)=1+123+133+⋯1.202;{displaystyle zeta (3)=1+{frac {1}{2^{3}}}+{frac {1}{3^{3}}}+cdots approx 1.202;!}zeta(3) = 1 + frac{1}{2^3} + frac{1}{3^3} + cdots approx 1.202;!   OEIS A002117
称为阿培里常數。



ζ(4)=1+124+134+⋯490≈1.0823;{displaystyle zeta (4)=1+{frac {1}{2^{4}}}+{frac {1}{3^{4}}}+cdots ={frac {pi ^{4}}{90}}approx 1.0823;!}zeta(4) = 1 + frac{1}{2^4} + frac{1}{3^4} + cdots = frac{pi^4}{90} approx 1.0823;!   OEIS A0013662

黑體輻射裡的斯特藩-玻尔兹曼定律和维恩近似。



s趨近於1


limε(1+ε)+ζ(1−ε)2=γ{displaystyle lim _{varepsilon to 0}{frac {zeta (1+varepsilon )+zeta (1-varepsilon )}{2}}=gamma }lim _{{varepsilon to 0}}{frac  {zeta (1+varepsilon )+zeta (1-varepsilon )}{2}}=gamma

其中γ是歐拉-馬歇羅尼常數=0.577215...{displaystyle 0.577215...}{displaystyle 0.577215...}



负整数


同样由欧拉发现,ζ函数在负整数点的值是有理数,这在模形式中发挥着重要作用,而且ζ函数在负偶整数点的值為零。


事實上


ζ(−n)=−Bn+1n+1{displaystyle zeta (-n)=-{frac {B_{n+1}}{n+1}}}zeta (-n)=-{frac  {B_{{n+1}}}{n+1}}

Bn是白努利數。


因為 B2n+1 =0,故ζ函数在负偶整数点的值為零。



复数值



ζ(x+iy)=∑k=1∞cos⁡(yln⁡k)−isin⁡(yln⁡k)kx,y∈R{displaystyle zeta (x+{rm {i}}y)=sum _{k=1}^{infty }{frac {cos(yln k)-{rm {i}}sin(yln k)}{k^{x}}},yin {mathbb {R} }}zeta(x+{rm{i}}y)=sum_{k=1}^{infty}frac{cos(yln k) -{rm{i}}sin(yln k)}{k^x},yin{mathbb{R}},x>1。


幅角



arg⁡(x+iy)]=−arctan⁡k=1∞sin⁡(yln⁡k)kx∑k=1∞cos⁡(yln⁡k)kx{displaystyle arg[zeta (x+{rm {i}}y)]=-arctan {frac {sum _{k=1}^{infty }{frac {sin(yln k)}{k^{x}}}}{sum _{k=1}^{infty }{frac {cos(yln k)}{k^{x}}}}}}arg[zeta(x+{rm{i}}y)]=-arctanfrac{sum_{k=1}^{infty}frac{sin(yln k)}{k^x}}{sum_{k=1}^{infty}frac{cos(yln k)}{k^x}}


函数值表




ζ(−2nn∈Z+)=0{displaystyle zeta (-2n_{nin mathbb {Z} ^{+}})=0}zeta(-2n_{nin mathbb{Z}^+})=0


ζ(−9)=−1132{displaystyle zeta (-9)=-{frac {1}{132}}}zeta(-9)=-frac{1}{132}


ζ(−7)=1240{displaystyle zeta (-7)={frac {1}{240}}}zeta(-7)=frac{1}{240}


ζ(−5)=−1252{displaystyle zeta (-5)=-{frac {1}{252}}}zeta(-5)=-frac{1}{252}


ζ(−3)=1120{displaystyle zeta (-3)={frac {1}{120}}}zeta(-3)=frac{1}{120}


ζ(−1)=−112{displaystyle zeta (-1)=-{frac {1}{12}}}zeta(-1)=-frac{1}{12}


ζ(0)=−12{displaystyle zeta (0)=-{frac {1}{2}}}zeta(0)=-frac{1}{2}



ζ(1−)=−{displaystyle zeta (1^{-})=-infty }{displaystyle zeta (1^{-})=-infty },



ζ(1+)=∞{displaystyle zeta (1^{+})=infty }{displaystyle zeta (1^{+})=infty },


ζ(2)=π26{displaystyle zeta (2)={frac {pi ^{2}}{6}}}zeta(2)=frac{pi^2}{6}


ζ(4)=π490{displaystyle zeta (4)={frac {pi ^{4}}{90}}}zeta(4)=frac{pi^4}{90}


ζ(6)=π6945{displaystyle zeta (6)={frac {pi ^{6}}{945}}}zeta(6)=frac{pi^6}{945}


ζ(8)=π89450{displaystyle zeta (8)={frac {pi ^{8}}{9450}}}zeta(8)=frac{pi^8}{9450}


ζ(10)=π1093555{displaystyle zeta (10)={frac {pi ^{10}}{93555}}}zeta(10)=frac{pi^{10}}{93555}



临界线上的数值计算


临界线上的数值计算可以通过黎曼-西格尔公式完成。
与之相关的,林德勒夫猜想英语Lindelöf hypothesis:对于任意给定的實数ϵ>0{displaystyle {begin{smallmatrix}epsilon >0end{smallmatrix}}}{begin{smallmatrix}epsilon >0end{smallmatrix}}


ζ(12+it)=O(tϵ){displaystyle zeta left({frac {1}{2}}+itright)=mathrm {O} (t^{epsilon })}zetaleft(frac{1}{2}+itright)=Omicron(t^{epsilon})


參考資料





  1. ^ http://nbviewer.ipython.org/github/empet/Math/blob/master/DomainColoring.ipynb


  2. ^ Devlin, Keith. The Millennium Problems: The Seven Greatest Unsolved Mathematical Puzzles of Our Time. New York: Barnes & Noble. 2002: 43–47. ISBN 978-0760786598. 


  3. ^ Tom M.Apostol. Introduction to Analytic Number Theory. 北京: 世界图书出版公司北京公司. 2012年1月: 55 – 56. ISBN 978-7-5100-4062-7 (英语). Theorem 3.2 If x ≥ 1 we have :  使用|accessdate=需要含有|url= (帮助)


  4. ^ 御坂01034. 巴塞尔问题(Basel problem)的多种解法. 


  5. ^ 加藤和也 黑川信重 斎藤毅. 数论I. 北京: 高等教育出版社. 2009年6月: 197–199. ISBN 978-7-04-026360-2 (中文). 短暂的沉默被打破了...  使用|accessdate=需要含有|url= (帮助)


  6. ^ 加藤和也 黑川信重 斎藤毅. 数论I. 北京: 高等教育出版社. 2009年6月: 209–210. ISBN 978-7-04-026360-2 (中文). Riamann对ζ研究的全部内容...  使用|accessdate=需要含有|url= (帮助)


  7. ^ C. S. Ogilvy & J. T. Anderson Excursions in Number Theory, pp. 29–35, Dover Publications Inc., 1988 ISBN 978-0-486-25778-5




相關條目



  • 黎曼猜想

  • 狄利克雷级数

  • 狄利克雷卷积

  • 黎曼ξ函數





Popular posts from this blog

How to pass form data using jquery Ajax to insert data in database?

Guess what letter conforming each word

Run scheduled task as local user group (not BUILTIN)