含有高阶广义调和数的无穷级数恒等式
Infinite Series Identities Containing High-Order Generalized Harmonic Numbers
DOI: 10.12677/pm.2024.146239, PDF, HTML, XML,    科研立项经费支持
作者: 王晓元, 孙艾莹:大连交通大学理学院,辽宁 大连
关键词: 调和数广义调和数Abel分部求和法Harmonic Numbers Generalized Harmonic Numbers Abel’s Method on Summation by Parts
摘要: 在本文中,定义一类m阶广义调和数利用组合分析中的Abel分部求和法,我们将推导出一些含有2阶广义调和数和3阶广义调和数的无穷级数与交错级数恒等式。进一步对参数a和b取特殊值,获得一些新的关于π,π2,π3,ζ(3),Catalan常数和ln2的无穷级数表达式。
Abstract: In this paper, we define a class of generalized harmonic numbers of order m by  Applying Abel’s method on summation by parts in combination analysis, we shall derive several infinite series and alternating series containing generalized harmonic numbers of order 2 and 3. Furthermore, by selecting special values for parameters a and b, several new infinite series ex-pressions are obtained for π, π2, π3, ζ(3), Catalan constant and ln2 as consequences.
文章引用:王晓元, 孙艾莹. 含有高阶广义调和数的无穷级数恒等式[J]. 理论数学, 2024, 14(6): 178-192. https://doi.org/10.12677/pm.2024.146239

1. 引言

调和数的研究由来已久,它出现在各种特殊函数的表达式中,在组合数学各分支、数论、计算机科学、计算机生物学和理论物理学等领域中都有重要应用。

众所周知,著名的经典调和数 H n 定义为

H 0 =0 H n = k=1 n 1 k = 0 1 1 t n 1t dt =γ+ψ( n+1 ) n=1,2,

其中欧拉常数 γ=0.5772 ψ= Γ /Γ 表示Euler Gamma函数的对数导数,其显式形式为

ψ( x )=γ+ k=0 { 1 k+1 1 k+x } .

1775年,Euler (见[1],p. 252)发现了与Riemann Zeta函数有关的调和数恒等式:

n=1 H n ( n+1 ) 2 = 1 2 n=1 H n n 2 =ζ( 3 ) 5 4 n=1 H n n 3 =ζ( 4 )

其中Riemann Zeta函数 ζ( s ) 定义为

ζ( s )={ k=1 1 k s = 1 1 2 s k=1 1 ( 2k1 ) s ( ( s )>1 ); 1 1 2 1s k=1 ( 1 ) k1 k s ( ( s )>0;s1 ),

为调和数的研究做出了巨大贡献。这项工作也激发了越来越多的数学工作者对调和数相关恒等式的研究热情。在过去的三十年里,关于调和数恒等式的证明方法也多种多样。De Doelder [2]通过计算双伽玛函数获得了大量有关调和数的求和恒等式。Borwein [3]利用发生函数和Parseval恒等式推导出调和数与 ζ( 4 ) 之间的求和表达式。Flajolet-Salvy [4]基于积分和留数计算得到线性与非线性的欧拉和,并通过一些求和公式给出Riemann Zeta函数值的表达式。陈永川等[5]利用与差分算子相关的Abel-Zeilberger算法证明了与调和数有关的若干有限级数求和公式。魏传安[6] [7]运用导数算子方法证明一些含有调和数与广义调和数的组合恒等式。Chu-De Donno [8]基于几个著名的超几何级数求和公式获得大量关于调和数的封闭求和公式。Gencev [9]利用Beta型积分推导出含有广义调和数和中心二项式系数的无穷级数求和公式。刘红梅[10]基于超几何级数中两个Gauss求和公式,推导出一些与中心二项式系数和广义调和数有关的无穷级数恒等式。闫庆伦等[11]运用部分分式法证明了一组关于欧拉型求和的组合恒等式,并计算了二项式系数倒数的乘积与高阶shifted调的有限求和形式。

最近,Chu [12] [13]利用Abel分部求和法重新证明了一些重要的超几何级数和基本超几何级数求和公式。2012年,他成功的将Abel分部求和法应用到调和数恒等式的研究领域,证明了一些含有经典调和数及其变换形式的无穷级数恒等式(见[14]),例如:

k=2 ( 1 ) k ( 2 k 2 +2k1 ) H ¯ k 2 k( k+1 )( k1 )( k+2 ) =2ln2 17 12

k=1 ( 1 ) k ( 4 k 2 3 ) O ¯ k 2 ( 2k+1 )( 2k+3 )( 2k1 )( 2k3 ) = π 32 .

2018年,Wang-Chu [15]非常巧妙地运用Abel分部求和法推导出含有调和数平方和立方的无穷级数公式,并获得一些漂亮的 π 2 ζ( 3 ) 的级数表达式,例如:

k1 H k 2 ( k+1 )( k+2 )( k+3 ) = π 2 24 3 16 k1 O k 2 ( 2k+1 )( 2k+3 )( 2k+5 ) = π 2 256 1 96

k2 ( 2k+1 ) H k 3 k( k+1 )( k+2 )( k1 ) = 6ζ( 3 )+19 12 k1 k O k 3 ( 2k3 )( 2k1 )( 2k+1 )( 2k+3 ) = 7ζ( 3 )4 512 .

2020年,Chen-Chen [16]定义了一类高阶广义调和数

h n ( m ) ( z )= k=1 n 1 ( k+z ) m ,其中 m,n z ,但 z1,2,,n

同样利用Abel分部求和法和部分分式分解方法证明了如下含有高阶广义调和数的无穷级数恒等式

n1 H n ( p ) ( n+a )( n+a+1 ) n1 H n ( p ) ( n+k k ) n1 H n ( p ) n( n+k k )

鉴于Abel分部求和法可以非常有效地证明与调和数有关的组合恒等式,因此,我们将进一步探索该方法在高阶广义调和数恒等式中的应用。

文中定义含有两个正实数ab的一类m阶广义调和数

h 0 m ( a,b )=0 h n m ( a,b )= k=1 n 1 ( aka+b ) m n,m=1,2,

不难看出 H n = h n 1 ( 1,1 ) 。当参数 a=b=1 a=2,b=1 时,得到两个特殊情况:

H n m := h n m ( 1,1 )= k=1 n 1 k m

O n m := h n m ( 2,1 )= k=1 n 1 ( 2k1 ) m

值得一提的是,在含有高阶调和数 H n m 的无穷级数恒等式中,有三个结构非常简洁漂亮的无穷级数:

n1 H n m n m n1 H n m ( n+1 )( n+a ) n1 H n m n( n+a )

其中第一个公式可以追溯到Euler时代(见文献[17]),另外两个公式出自Sofo [18]的研究工作。这也启发了作者想从事这方面的研究,获得更多这类广义调和数求和公式。因此,本文将利用组合分析中的Abel分部求和法,推导出一些含有2阶广义调和数和3阶广义调和数的无穷级数与交错级数恒等式,同时建立一些新的关于 π π 2 π 3 ζ( 3 ) ,Catalan常数和 ln2 的无穷级数表达式。

在接下来的推导过程中,我们需要利用Catalan常数G的表达式

G= k=0 ( 1 ) k / ( 2k+1 ) 2 =0.9159

对数 ln( 1+x ) 的泰勒展开式

ln( 1+x )= k=1 ( 1 ) k1 x k k

以及一些已知的Euler和进行计算,比如:

k=0 ( 1 ) k 2k+1 = π 4 k=1 1 k 2 = π 2 6 k=0 1 ( 2k+1 ) 2 = π 2 8 k=0 ( 1 ) k ( 2k+1 ) 3 = π 3 32

2. Abel分部求和法

组合分析学中的Abel分部求和法是检验无穷级数收敛的基本方法。2006年,Chu [19]首次将经典分析中的Abel分部求和法应用到了基本超几何级数领域,证明了基本超几何级数理论中最著名的公式之一:Bailey的 ψ 6 6 -级数恒等式。2012年,Chu [14]再次扩大Abel分部求和法的研究领域,并利用该方法推导出含有调和数的无穷级数恒等式。本文我们将继续利用该方法对高阶广义调和数进行研究。为了便于后续应用,在此我们重新回顾一下Abel分部求和引理。

对于任意的复数序列 { τ k } ,分别定义向前和向后差分算子 Δ

Δ τ k = τ k τ k+1 τ k = τ k τ k1

需要指出的是,这里的 Δ 与通常的向前差分算子仅相差一个负号,则有如下公式成立。

引理1 (Abel分部求和引理)

k1 B k A k = [ AB ] + A 0 B 1 + k1 A k Δ B k (1)

其中极限 [ AB ] + = lim n A n B n+1 存在,且公式中的一个非终止型级数是收敛的。

证明 根据向后差分的定义,对任意自然数n,有

k=1 n B k A k = k=1 n B k { A k A k1 } = k=1 n A k B k k=1 n A k1 B k

将上式最后一个求和流标k替换为 k+1 ,得到

k=1 n B k A k = A n B n+1 A 0 B 1 + k=1 n A k { B k B k+1 } = A n B n+1 A 0 B 1 + k=1 n A k Δ B k

n ,便可得到Abel分部求和引理的表达式,证毕。

下面举例说明如何利用Abel分部求和引理推导含有高阶广义调和数的无穷级数恒等式。

首先选取序列 A k = ( 1 ) k ak+b B k = h k 2 ( a,b ) ,不难计算向后和向前差分为

A k = A k A k1 = ( 1 ) k ak+b ( 1 ) k1 aka+b = ( 1 ) k ( 2aka+2b ) ( ak+b )( aka+b )

Δ B k = B k B k+1 = h k 2 ( a,b ) h k+1 2 ( a,b )= 1 ( ak+b ) 2

计算极限和初值

[ AB ] + =0 A 0 B 1 = 1 b 3

然后根据公式(1),最后得到含有二阶广义调和数的无穷级数恒等式

k1 ( 1 ) k ( 2aka+2b ) h k 2 ( a,b ) ( ak+b )( aka+b ) = 1 b 3 + k1 ( 1 ) k+1 ( ak+b ) 3 = k0 ( 1 ) k+1 ( ak+b ) 3

在上式中,令参数 a=b=1 ,由Riemann Zeta函数 ζ( s ) 的定义可得

k0 ( 1 ) k+1 ( k+1 ) 3 = k1 ( 1 ) k k 3 = 3 4 ζ( 3 )

则推导出新的关于 ζ( 3 ) 的无穷级数表达式

k1 ( 1 ) k ( 2k+1 ) H k 2 k( k+1 ) = 3 4 ζ( 3 )

再令参数 a=2,b=1 ,由Euler和公式

k=0 ( 1 ) k ( 2k+1 ) 3 = π 3 32

则推导出新的关于 π 3 的无穷级数表达式

k1 ( 1 ) k k O k 2 ( 2k1 )( 2k+1 ) = π 3 128

由上述例子可以发现在利用Abel分部求和法去研究组合公式时,关建是要找到恰当的差分对,即序列 { A k } { B k } 。本文中,为了便于接下来的计算方便快捷,我们固定其中一个序列 B k = h k m ( a,b ) ,不难计算其向前差分算子为

Δ B k = h k m ( a,b ) h k+1 m ( a,b )= 1 ( ak+b ) m

同时,文中所选择的序列 { A k } 要求满足当 n 时, A n 0 ,又因为 B n 是有界的,则 [ AB ] + =0 。注意到,当 B k = h k m ( a,b ) 时, B 0 =0 。因此,我们可以得到下面计算公式

k1 h k m ( a,b )A = k0 A k ( ak+b ) m (2)

接下来在本文的第三节和第四节,我们将主要利用公式(2)推导含有二阶广义调和数 h k 2 ( a,b ) 和三阶广义调和数 h k 3 ( a,b ) 的无穷级数与交错级数求和公式,进一步获得一些新的关于 π π 2 π 3 ζ( 3 ) ,Catalan常数和 ln2 的表达式,且相关结果均已用数学软件Mathematica程序进行验证,以确保其正确性。

3. 含有 h k 2 ( a,b ) 的无穷级数恒等式

在本节中,固定序列 { B k } B k = h k 2 ( a,b ) ,利用Abel分部求和法获得两个含有二阶广义调和数的

无穷级数恒等式和两个含有二阶广义调和数的交错级数恒等式。

定义序列 A k = 1 ( ak+a+b )( ak+2a+b ) ,计算其向后差分为

A k = 2a ( ak+b )( ak+a+b )( ak+2a+b )

根据公式(2),得到如下恒等式

k1 2a h k 2 ( a,b ) ( ak+b )( ak+a+b )( ak+2a+b ) = k0 1 ( ak+b ) 2 ( ak+a+b )( ak+2a+b )

由部分分式分解

1 ( ak+b ) 2 ( ak+a+b )( ak+2a+b ) = 1 2 a 2 ( ak+b ) 2 3 4 a 3 ( ak+b ) + 1 a 3 ( ak+a+b ) 1 4 a 3 ( ak+2a+b )

可以得到

k0 1 ( ak+b ) 2 ( ak+a+b )( ak+2a+b ) = k0 { 1 2 a 2 ( ak+b ) 2 3 4 a 3 ( ak+b ) + 1 a 3 ( ak+a+b ) 1 4 a 3 ( ak+2a+b ) } = 1 2 a 2 k0 1 ( ak+b ) 2 3 4 a 3 k0 1 ak+b + 1 a 3 k1 1 ak+b 1 4 a 3 k2 1 ak+b = 1 2 a 2 k0 1 ( ak+b ) 2 3a+2b 4 a 3 b( a+b )

则推出下面求和公式。

定理1 (含有 h k 2 ( a,b ) 的无穷级数恒等式)

k1 2a h k 2 ( a,b ) ( ak+b )( ak+a+b )( ak+2a+b ) = 1 2 a 2 k0 1 ( ak+b ) 2 3a+2b 4 a 3 b( a+b )

令参数 a=b=1 a=2,b=1 ,可以推出含有二阶调和数及其变换形式的表达式。

推论2 ( π 2 的级数表达式)

k1 H k 2 ( k+1 )( k+2 )( k+3 ) = 5 16 + π 2 24

k1 O k 2 ( 2k+1 )( 2k+3 )( 2k+5 ) = 1 48 + π 2 256

定义序列 A k = 1 ( ak+b )( 2ak+a+2b ) ,计算其向后差分为

A k = 4 a 2 k a 2 +4ab ( aka+b )( ak+b )( 2aka+2b )( 2ak+a+2b )

根据公式(2),得到如下恒等式

k1 ( 4 a 2 k a 2 +4ab ) h k 2 ( a,b ) ( aka+b )( ak+b )( 2aka+2b )( 2ak+a+2b ) = k0 1 ( ak+b ) 3 ( 2ak+a+2b )

由部分分式分解

1 ( ak+b ) 3 ( 2ak+a+2b ) = 1 a ( ak+b ) 3 2 a 2 ( ak+b ) 2 + 4 a 3 ( ak+b ) 8 a 3 ( 2ak+a+2b )

可以得到

k0 1 ( ak+b ) 3 ( 2ak+a+2b ) = k0 { 1 a ( ak+b ) 3 2 a 2 ( ak+b ) 2 + 4 a 3 ( ak+b ) 8 a 3 ( 2ak+a+2b ) } = 1 a k0 1 ( ak+b ) 3 2 a 2 k0 1 ( ak+b ) 2 + 8 a 3 k0 1 2ak+2b 8 a 3 k0 1 2ak+a+2b = 1 a k0 1 ( ak+b ) 3 2 a 2 k0 1 ( ak+b ) 2 + 8 a 3 k0 ( 1 ) k ak+2b

则推出下面求和公式。

定理3 (含有 h k 2 ( a,b ) 的无穷级数恒等式)

k1 ( 4 a 2 k a 2 +4ab ) h k 2 ( a,b ) ( aka+b )( ak+b )( 2aka+2b )( 2ak+a+2b ) = 1 a k0 1 ( ak+b ) 3 2 a 2 k0 1 ( ak+b ) 2 + 8 a 3 k0 ( 1 ) k ak+2b

令参数 a=b=1 a=2,b=1 ,可以推出含有二阶调和数及其变换形式的表达式。

推论4 ( π 2 ζ( 3 ) ln2 的级数表达式)

k1 ( 4k+3 ) H k 2 k( k+1 )( 2k+1 )( 2k+3 ) =8 π 2 3 8ln2+ζ( 3 )

k1 ( 4k+1 ) O k 2 k( k+1 )( 2k1 )( 2k+1 ) = π 2 4 +2ln2+ 7 4 ζ( 3 )

接下来,将给出含有二阶调和数及其变换形式的交错级数求和公式。

定义序列 A k = ( 1 ) k ak+2a+b ,计算其向后差分为

A k = ( 1 ) k ( 2ak+3a+2b ) ( ak+a+b )( ak+2a+b )

根据公式(2),得到如下恒等式

k1 ( 1 ) k ( 2ak+3a+2b ) h k 2 ( a,b ) ( ak+2a+b )( ak+a+b ) = k0 ( 1 ) k+1 ( ak+b ) 2 ( ak+2a+b )

由部分分式分解

( 1 ) k+1 ( ak+b ) 2 ( ak+2a+b ) = ( 1 ) k+1 2a ( ak+b ) 2 ( 1 ) k+1 4 a 2 ( ak+b ) + ( 1 ) k+1 4 a 2 ( ak+2a+b )

可以得到

k0 ( 1 ) k+1 ( ak+b ) 2 ( ak+2a+b ) = k0 { ( 1 ) k+1 2a ( ak+b ) 2 ( 1 ) k+1 4 a 2 ( ak+b ) + ( 1 ) k+1 4 a 2 ( ak+2a+b ) } = 1 2a k0 ( 1 ) k+1 ( ak+b ) 2 1 4 a 2 k2 ( 1 ) k ak+b + 1 4 a 2 k2 ( 1 ) k ak+b + 1 4ab( a+b ) = 1 2a k0 ( 1 ) k+1 ( ak+b ) 2 + 1 4ab( a+b )

则推出下面求和公式。

定理5 (含有 h k 2 ( a,b ) 的交错级数恒等式)

k1 ( 1 ) k ( 2ak+3a+2b ) h k 2 ( a,b ) ( ak+a+b )( ak+2a+b ) = 1 2a k0 ( 1 ) k+1 ( ak+b ) 2 + 1 4ab( a+b )

令参数 a=b=1 a=2,b=1 ,可以推出含有二阶调和数及其变换形式的表达式。

推论6 ( π 2 G的级数表达式)

k1 ( 1 ) k ( 2k+5 ) H k 2 ( k+2 )( k+3 ) = 1 8 π 2 24

k1 ( 1 ) k ( k+2 ) O k 2 ( 2k+3 )( 2k+5 ) = 1 96 G 16

定义序列 A k = ( 1 ) k ( ak+b )( ak+a+b ) ,计算其向后差分为

A k = 2 ( 1 ) k ( aka+b )( ak+a+b )

根据公式(2),得到如下恒等式

k1 2 ( 1 ) k h k 2 ( a,b ) ( aka+b )( ak+a+b ) = k0 ( 1 ) k+1 ( ak+b ) 3 ( ak+a+b )

由部分分式分解

( 1 ) k+1 ( ak+b ) 3 ( ak+a+b ) = ( 1 ) k+1 a ( ak+b ) 3 ( 1 ) k+1 a 2 ( ak+b ) 2 + ( 1 ) k+1 a 3 ( ak+b ) ( 1 ) k+1 a 3 ( ak+a+b )

可以得到

k0 ( 1 ) k+1 ( ak+b ) 3 ( ak+a+b ) = k0 { ( 1 ) k+1 a ( ak+b ) 3 ( 1 ) k+1 a 2 ( ak+b ) 2 + ( 1 ) k+1 a 3 ( ak+b ) ( 1 ) k+1 a 3 ( ak+a+b ) } = 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 1 a 3 k0 ( 1 ) k+1 ak+b 1 a 3 k1 ( 1 ) k ak+b = 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 2 a 3 k1 ( 1 ) k+1 ak+b 1 a 3 b

则推出下面求和公式。

定理7 (含有 h k 2 ( a,b ) 的交错级数恒等式)

k1 2 ( 1 ) k h k 2 ( a,b ) ( aka+b )( ak+a+b ) = 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 2 a 3 k1 ( 1 ) k+1 ak+b 1 a 3 b

令参数 a=b=1 a=2,b=1 ,可以推出含有二阶调和数及其变换形式的表达式。

推论8 ( π π 2 π 3 ζ( 3 ) ln2 G级数表达式)

k1 ( 1 ) k H k 2 k( k+2 ) = 1 2 + π 2 24 ln2 3 8 ζ( 3 )

k1 ( 1 ) k O k 2 ( 2k1 )( 2k+3 ) = 1 16 + G 8 π 3 128 π 32

4. 含有 h k 3 ( a,b ) 的无穷级数恒等式

在本节中,固定序列 { B k } B k = h k 3 ( a,b ) ,利用Abel分部求和法获得两个含有三阶广义调和数的无

穷级数恒等式和三个含有三阶广义调和数的交错级数恒等式。

定义序列 A k = 1 ( 2ak+a+2b )( 2aka+2b ) ,计算其向后差分为

A k = 4a ( 2ak+a+2b )( 2ak3a+2b )( 2aka+2b )

根据公式(2),得到如下恒等式

k1 4a h k 3 ( a,b ) ( 2ak+a+2b )( 2ak3a+2b )( 2aka+2b ) = k0 1 ( ak+b ) 3 ( 2ak+a+2b )( 2aka+2b )

由部分分式分解

1 ( ak+b ) 3 ( 2ak+a+2b )( 2aka+2b ) = 1 a 2 ( ak+b ) 3 4 a 4 ( ak+b ) + 4 a 4 ( 2aka+2b ) + 4 a 4 ( 2ak+a+2b )

可以得到

k0 1 ( ak+b ) 3 ( 2ak+a+2b )( 2aka+2b ) = k0 { 1 a 2 ( ak+b ) 3 4 a 4 ( ak+b ) + 4 a 4 ( 2aka+2b ) + 4 a 4 ( 2ak+a+2b ) } = 1 a 2 k0 1 ( ak+b ) 3 4 a 4 k0 1 ak+b + 4 a 4 k0 1 2aka+2b + 4 a 4 k1 1 2aka+2b = 1 a 2 k0 1 ( ak+b ) 3 8 a 4 k0 ( 1 ) k ak+2b + 4 a 4 ( 2ba )

则推出下面求和公式。

定理9 (含有 h k 3 ( a,b ) 的无穷级数恒等式)

k1 4a h k 3 ( a,b ) ( 2ak+a+2b )( 2ak3a+2b )( 2aka+2b ) = 1 a 2 k0 1 ( ak+b ) 3 8 a 4 k0 ( 1 ) k ak+2b + 4 a 4 ( 2ba )

令参数 a=b=1 ,可以推出下面含有三阶调和数的表达式。

推论10 ( ζ( 3 ) ln2 级数表达式)

k1 H k 3 ( 2k1 )( 2k+1 )( 2k+3 ) =1+2ln2 1 4 ζ( 3 )

在定理9中,令参数 a=2,b=1 ,利用 ln( 1+x ) 的泰勒展开式和 ζ( s ) 定义式计算如下极限

lim a2 b1 { 1 a 2 k0 1 ( ak+b ) 3 8 a 4 k0 ( 1 ) k ak+2b + 4 a 4 ( 2ba ) + 4a b 3 ( a2b )( a+2b )( 3a+2b ) } = lim a2 b1 { 4 a 4 +8 a 3 b+16 a 2 b 2 +20a b 3 +8 b 4 a 4 b 3 ( a+2b )( 3a+2b ) 1 a 2 k0 1 ( ak+b ) 3 8 a 4 k0 ( 1 ) k ak+2b } = 15 32 1 4 ln2 7 32 ζ( 3 )

能够推导出下面三阶调和数变换形式的计算公式。

推论11 ( ζ( 3 ) ln2 级数表达式)

k2 O k 3 k( k1 )( k+1 ) = 15 4 2ln2 7 4 ζ( 3 )

定义序列 A k = 1 2ak+a+2b ,计算其向后差分为

A k = 2a ( 2aka+2b )( 2ak+a+2b )

根据公式(2),得到如下恒等式

k1 2a h k 3 ( a,b ) ( 2aka+2b )( 2ak+a+2b ) = k0 1 ( ak+b ) 3 ( 2ak+a+2b )

由部分分式分解

1 ( ak+b ) 3 ( 2ak+a+2b ) = 1 a ( ak+b ) 3 2 a 2 ( ak+b ) 2 + 4 a 3 ( ak+b ) 8 a 3 ( 2ak+a+2b )

可以得到

k0 1 ( ak+b ) 3 ( 2ak+a+2b ) = k0 { 1 a ( ak+b ) 3 2 a 2 ( ak+b ) 2 + 4 a 3 ( ak+b ) 8 a 3 ( 2ak+a+2b ) } = 1 a k0 1 ( ak+b ) 3 2 a 2 k0 1 ( ak+b ) 2 + 8 a 3 k0 1 2ak+2b 8 a 3 k0 1 2ak+a+2b = 1 a k0 1 ( ak+b ) 3 2 a 2 k0 1 ( ak+b ) 2 + 8 a 3 k0 ( 1 ) k ak+2b

则推出下面求和公式。

定理12 (含有 h k 3 ( a,b ) 的无穷级数恒等式)

k1 2a h k 3 ( a,b ) ( 2aka+2b )( 2ak+a+2b ) = 1 a k0 1 ( ak+b ) 3 2 a 2 k0 1 ( ak+b ) 2 + 8 a 3 k0 ( 1 ) k ak+2b

令参数 a=b=1 a=2,b=1 ,可以推出含有三阶调和数及其变换形式的表达式。

推论13 ( ζ( 3 ) π 2 ln2 级数表达式)

k1 H k 3 ( 2k+1 )( 2k+3 ) =4 π 2 6 4ln2+ 1 2 ζ( 3 )

k1 O k 3 k( k+1 ) = π 2 4 +2ln2+ 7 4 ζ( 3 )

接下来,将给出含有三阶调和数及其变换形式的交错级数求和公式。

定义序列 A k = ( 1 ) k ( ak+a+b )( aka+b ) ,计算其向后差分为

A k = ( 1 ) k ( 2 a 2 k 2 2 a 2 k+4abk a 2 2ab+2 b 2 ) ( ak2a+b )( aka+b )( ak+b )( ak+a+b )

根据公式(2),得到如下恒等式

k1 ( 1 ) k ( 2 a 2 k 2 2 a 2 k+4abk a 2 2ab+2 b 2 ) h k 3 ( a,b ) ( ak2a+b )( aka+b )( ak+b )( ak+a+b ) = k0 ( 1 ) k+1 ( ak+b ) 3 ( ak+a+b )( aka+b )

由部分分式分解

( 1 ) k+1 ( ak+b ) 3 ( ak+a+b )( aka+b ) = ( 1 ) k a 2 ( ak+b ) 3 + ( 1 ) k a 4 ( ak+b ) + ( 1 ) k+1 2 a 4 ( ak+a+b ) + ( 1 ) k+1 2 a 4 ( aka+b )

可以得到

k0 ( 1 ) k+1 ( ak+b ) 3 ( ak+a+b )( aka+b ) = k0 { ( 1 ) k a 2 ( ak+b ) 3 + ( 1 ) k a 4 ( ak+b ) + ( 1 ) k+1 2 a 4 ( ak+a+b ) + ( 1 ) k+1 2 a 4 ( aka+b ) } = 1 a 2 k0 ( 1 ) k ( ak+b ) 3 + 1 a 4 k0 ( 1 ) k ak+b + 1 2 a 4 k0 ( 1 ) k+1 aka+b + 1 2 a 4 k2 ( 1 ) k+1 aka+b = 1 a 2 k0 ( 1 ) k ( ak+b ) 3 + 2 a 4 k2 ( 1 ) k+1 aka+b + 3a2b 2 a 4 b( ab )

则推出下面求和公式。

定理14 (含有 h k 3 ( a,b ) 的交错级数恒等式)

k1 ( 1 ) k ( 2 a 2 k 2 2 a 2 k+4abk a 2 2ab+2 b 2 ) h k 3 ( a,b ) ( ak2a+b )( aka+b )( ak+b )( ak+a+b ) = 1 a 2 k0 ( 1 ) k ( ak+b ) 3 + 2 a 4 k2 ( 1 ) k+1 aka+b + 3a2b 2 a 4 b( ab )

令参数 a=2,b=1 ,可以推出含有三阶调和数及其变换形式的表达式。

推论15 ( π π 3 的级数表达式)

k1 ( 1 ) k ( 8 k 2 6 ) O k 3 ( 2k3 )( 2k1 )( 2k+1 )( 2k+3 ) = π 32 + π 3 128

在定理14中,令参数 a=b=1 ,利用 ln( 1+x ) 的泰勒展开式和 ζ( s ) 定义式计算如下极限

lim a1 b1 { 1 a 2 k0 ( 1 ) k ( ak+b ) 3 + 2 a 4 k2 ( 1 ) k+1 aka+b + 3a2b 2 a 4 b( ab ) a 2 +2ab+2 b 2 ( ab )( a+b )( 2a+b ) b 4 } = lim a1 b1 { 2 a 5 2 a 4 b6 a 3 b 2 +5a b 4 +2 b 5 2 a 4 b 4 ( a+b )( 2a+b ) 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 3 + 2 a 4 k2 ( 1 ) k+1 aka+b } = 23 12 +2ln2+ 3 4 ζ( 3 )

可以推导出下面含有三阶调和数变换形式的计算公式。

推论16 ( ζ( 3 ) ln2 的级数表达式)

k2 ( 1 ) k ( 2 k 2 +2k1 ) H k 3 k( k1 )( k+1 )( k+2 ) = 23 12 +2ln2+ 3 4 ζ( 3 )

定义序列 A k = ( 1 ) k ak+a+b ,计算其向后差分为

A k = ( 1 ) k ( 2ak+a+2b ) ( ak+b )( ak+a+b )

根据公式(2),得到如下恒等式

k1 ( 1 ) k ( 2ak+a+2b ) h k 3 ( a,b ) ( ak+b )( ak+a+b ) = k0 ( 1 ) k+1 ( ak+b ) 3 ( ak+a+b )

由部分分式分解

( 1 ) k+1 ( ak+b ) 3 ( ak+a+b ) = ( 1 ) k+1 a ( ak+b ) 3 ( 1 ) k+1 a 2 ( ak+b ) 2 + ( 1 ) k+1 a 3 ( ak+b ) ( 1 ) k+1 a 3 ( ak+a+b )

可以得到

k0 ( 1 ) k+1 ( ak+b ) 3 ( ak+a+b ) = k0 { ( 1 ) k+1 a ( ak+b ) 3 ( 1 ) k+1 a 2 ( ak+b ) 2 + ( 1 ) k+1 a 3 ( ak+b ) ( 1 ) k+1 a 3 ( ak+a+b ) } = 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 1 a 3 k0 ( 1 ) k+1 ak+b + 1 a 3 k1 ( 1 ) k+1 ak+b = 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 2 a 3 k1 ( 1 ) k+1 ak+b 1 a 3 b

则推出下面求和公式。

定理17 (含有 h k 3 ( a,b ) 的交错级数恒等式)

k1 ( 1 ) k ( 2ak+a+2b ) h k 3 ( a,b ) ( ak+b )( ak+a+b ) = 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 2 a 3 k1 ( 1 ) k+1 ak+b 1 a 3 b

令参数 a=b=1 a=2,b=1 ,可以推出含有三阶调和数及其变换形式的表达式。

推论18 ( π π 2 π 3 ζ( 3 ) ln2 G的级数表达式)

k1 ( 1 ) k ( 2k+3 ) H k 3 ( k+1 )( k+2 ) =1+ π 2 12 2ln2 3 4 ζ( 3 )

k1 ( 1 ) k ( k+1 ) O k 3 ( 2k+1 )( 2k+3 ) = 1 32 + G 16 π 64 π 3 256

定义序列 A k = ( 1 ) k aka+b ,计算其向后差分为

A k = ( 1 ) k ( 2ak3a+2b ) ( ak2a+b )( aka+b )

根据公式(2),得到如下恒等式

k1 ( 1 ) k ( 2ak3a+2b ) h k 3 ( a,b ) ( ak2a+b )( aka+b ) = k0 ( 1 ) k+1 ( ak+b ) 3 ( aka+b )

由部分分式分解

( 1 ) k+1 ( ak+b ) 3 ( aka+b ) = ( 1 ) k a ( ak+b ) 3 + ( 1 ) k a 2 ( ak+b ) 2 + ( 1 ) k a 3 ( ak+b ) ( 1 ) k a 3 ( aka+b )

可以得到

k0 ( 1 ) k+1 ( ak+b ) 3 ( aka+b ) = k0 { ( 1 ) k a ( ak+b ) 3 + ( 1 ) k a 2 ( ak+b ) 2 + ( 1 ) k a 3 ( ak+b ) ( 1 ) k a 3 ( aka+b ) } = 1 a k0 ( 1 ) k ( ak+b ) 3 + 1 a 2 k0 ( 1 ) k ( ak+b ) 2 1 a 3 k1 ( 1 ) k aka+b 1 a 3 k0 ( 1 ) k aka+b = 1 a k0 ( 1 ) k ( ak+b ) 3 + 1 a 2 k0 ( 1 ) k ( ak+b ) 2 2 a 3 k1 ( 1 ) k aka+b + 1 a 3 ( ab )

则推出下面求和公式。

定理19 (含有 h k 3 ( a,b ) 的交错级数恒等式)

k1 ( 1 ) k ( 2ak3a+2b ) h k 3 ( a,b ) ( ak2a+b )( aka+b ) = 1 a k0 ( 1 ) k ( ak+b ) 3 + 1 a 2 k0 ( 1 ) k ( ak+b ) 2 2 a 3 k1 ( 1 ) k aka+b + 1 a 3 ( ab )

令参数 a=2,b=1 ,可以推出含有三阶调和数变换形式的表达式。

推论20 ( π π 3 G的级数表达式)

k1 ( 1 ) k ( k1 ) O k 3 ( 2k3 )( 2k1 ) = 1 32 + G 16 + π 64 + π 3 256

在定理19中,令参数 a=b=1 ,利用Euler和公式, ln( 1+x ) 的泰勒展开式和 ζ( s ) 定义式计算如下极限

lim a1 b1 { k0 ( 1 ) k+1 ( ak+b ) 3 ( aka+b ) a2b b 4 ( ba ) } = lim a1 b1 { a 3 a 2 ba b 2 b 3 a 3 b 4 1 a k0 ( 1 ) k+1 ( ak+b ) 3 1 a 2 k0 ( 1 ) k+1 ( ak+b ) 2 + 2 a 3 k1 ( 1 ) k+1 aka+b } =2+ π 2 12 +2ln2+ 3 4 ζ( 3 )

可以推导出下面含有三阶调和数的求和公式。

推论21 ( ζ( 3 ) π 2 ln2 的级数表达式)

k2 ( 1 ) k ( 2k1 ) H k 3 k( k1 ) =2+ π 2 12 +2ln2+ 3 4 ζ( 3 )

5. 结论

本文利用组合分析中的Abel分部求和法,固定差分对中的一个序列 { B k } B k = h k m ( a,b ) ,再选取恰当序列 { A k } ,利用主要公式(2)获得一些新的含有二阶广义调和数和三阶广义调和数的无穷级数与交错级数求和公式,同时建立一些新的关于 π π 2 π 3 ζ( 3 ) ,Catalan常数和 ln2 的无穷级数表达式。

该研究将Abel分部求和法应用到高阶广义调和数恒等式的证明中,对调和数的研究领域起到一定的推动作用。文中所得结果不仅推动圆周率近似值研究的进步,获得的Riemann Zeta函数 ζ( 3 ) 也在素数分布理论中扮演着重要角色,而且Catalan常数和对数的表达式在组合计数和组合分析中都有重要应用价值。此外在物理学中,量子场论里大规模高阶扰动运算中的Feynman积分计算需要利用很多调和数和的乘积展开式,而且一些含有调和数的求和公式可以对粒子物理学中的散射过程进行高阶校正的计算。因此,本文的研究成果将在数学和物理多个领域都有重要应用。最后,对本文研究内容感兴趣的读者也可以进一步尝试证明更多含有高阶广义调和数的组合恒等式。

基金项目

辽宁省教育厅科学研究项目(项目编号JDL2019028)。

参考文献

[1] Berndt, B.C. (1985) Ramanujan’s Notebooks, Part I. Springer.
[2] De Doelder, P.J. (1991) On Some Series Containing and for Certain Values of X and Y. Journal of Computational and Applied Mathematics, 37, 125-141.
https://doi.org/10.1016/0377-0427(91)90112-w
[3] Borwein, D. and Borwein, J.M. (1995) On an Intriguing Integral and Some Series Related to . Proceedings of the American Mathematical Society, 123, 1191.
https://doi.org/10.2307/2160718
[4] Flajolet, P. and Salvy, B. (1998) Euler Sums and Contour Integral Representations. Experimental Mathematics, 7, 15-35.
https://doi.org/10.1080/10586458.1998.10504356
[5] Chen, W.Y.C., Hou, Q. and Jin, H. (2011) The Abel-Zeilberger Algorithm. The Electronic Journal of Combinatorics, 18, Article No. P17.
https://doi.org/10.37236/2013
[6] Wei, C. and Gong, D. (2014) The Derivative Operator and Harmonic Number Identities. The Ramanujan Journal, 34, 361-371.
https://doi.org/10.1007/s11139-013-9510-2
[7] Wei, C., Yan, Q. and Gong, D. (2015) A Family of Summation Formulae Involving Harmonic Numbers. Integral Transforms and Special Functions, 26, 667-677.
https://doi.org/10.1080/10652469.2015.1034124
[8] Chu, W. and De Donno, L. (2005) Hypergeometric Series and Harmonic Number Identities. Advances in Applied Mathematics, 34, 123-137.
https://doi.org/10.1016/j.aam.2004.05.003
[9] Genčev, M. (2011) Binomial Sums Involving Harmonic Numbers. Mathematica Slovaca, 61, 215-226.
https://doi.org/10.2478/s12175-011-0006-5
[10] 刘红梅. 含有中心二项式系数以及广义调和数的无穷级数恒等式[J]. 数学物理学报, 2020, 40A(3): 631-640.
[11] 闫庆伦, 王照芬, 米娟. 高阶shifted调和数的有限求和[J]. 华东师范大学学报(自然科学版), 2021(6): 24-32.
[12] Wenchang, C. (2006) Abel’s Method on Summation by Parts and Hypergeometric Series. Journal of Difference Equations and Applications, 12, 783-798.
https://doi.org/10.1080/10236190600704096
[13] Chu, W. (2007) Abel’s Lemma on Summation by Parts and Basic Hypergeometric Series. Advances in Applied Mathematics, 39, 490-514.
https://doi.org/10.1016/j.aam.2007.02.001
[14] Chu, W. (2011) Infinite Series Identities on Harmonic Numbers. Results in Mathematics, 61, 209-221.
https://doi.org/10.1007/s00025-010-0089-2
[15] Wang, X. and Chu, W. (2018) Infinite Series Identities Involving Quadratic and Cubic Harmonic Numbers. Publicacions Matemàtiques, 62, 285-300.
https://doi.org/10.5565/publmat6211813
[16] Chen, K. and Chen, Y. (2020) Infinite Series Containing Generalized Harmonic Functions. Notes on Number Theory and Discrete Mathematics, 26, 85-104.
https://doi.org/10.7546/nntdm.2020.26.2.85-104
[17] Sondow, J. and Weisstein, W.E. Harmonic Number. Wolfram MathWorld.
http://mathworld.wolfram.com/HarmonicNumber.html
[18] Sofo, A. (2011) Harmonic Number Sums in Higher Powers. Journal of Mathematical Analysis, 2, 15-22.
[19] Chu, W. (2006) Bailey’s Very Well-Poised-Series Identity. Journal of Combinatorial Theory, Series A, 113, 966-979.
https://doi.org/10.1016/j.jcta.2005.08.009