超几何函数
F(a1,a2,⋯,amb1,b2,⋯,bn∣∣∣∣∣∣ z)=k⩾0∑b1k⋯bnka1k⋯amk⋅k!zk
其中 k!ak=(ka+k−1)
超几何函数的特殊情形
F(∣ z)=k⩾0∑k!zk=ez,为了使记号整齐,常常写成 F(11∣∣∣∣∣∣ z)=k⩾0∑k!zk=ez
又因为 1k=k!,所以我们有
F(1,11∣∣∣∣∣∣ z)=k⩾0∑zk=1−z1
另外,根据 k!(ka+k−1)=ak 可以得到封闭形式
F(a,11∣∣∣∣∣∣ z)=k⩾0∑akk!zk=k∑(ka+k−1)zk=(1−z)a1
最后一步是根据上指标反转,如果令 a←−a, z←−z,可以得到
F(−a,11∣∣∣∣∣∣ −z)=(1+z)a
修正贝塞尔函数
F(1b,1∣∣∣∣∣∣ z)=k⩾0∑(b−1+k)!(b−1)!k!zk=Ib−1(2z)z(b−1)/2(b−1)!
其中 Ib−1 为阶是 b−1 的修正贝塞尔函数,当 b=1 时,有
F(1,11∣∣∣∣∣∣ z)=I0(2z)=k⩾0∑k!zk
合流超几何函数
m=n=1 时,超几何函数有特殊情形
F(ab∣∣∣∣∣∣ z)=k⩾0∑bkakk!zk=M(a,b,z)
高斯超几何级数
F(a,bc∣∣∣∣∣∣ z)=k⩾0∑ckakbkk!zk
其中一种特殊情况是,ln(1+x)=k=1∑∞(−1)k−1k!xk=k=0∑∞(−1)kk+1xk+1=xk=0∑∞k+1(−x)k
从而有
zF(1,12∣∣∣∣∣∣ −z)=zk⩾0∑(k+1)!k!k!k!(−z)k
=z−2z2+3z3−4z4+⋯=ln(1+z)
高斯超几何函数具有循环分解形式
F(a,bc∣∣∣∣∣∣ z)=1+1acbz(1+2a+1c+1b+1z(1+3a+2c+2b+2z(1+⋯)))
超几何函数的常用结论
(1)
(−a)k=0,当 k>a⩾0 且 a 为整数
(2)
Γ 函数
z!=∫0∞tze−tdt,Rz>−1
证明如下
令 s>0,则有 dxd(e−xxs)=−e−xxs+se−xxs−1
两边同时积分
[e−xxs]ϵl=−∫ϵle−xxsdx+s∫ϵle−xxs−1dx
如果记 Γ(s)=∫0∞e−xxs−1dx,s>0
当 l→∞,ϵ→0,[e−xxs]ϵl=0,所以有
Γ(s+1)=sΓ(s),又因为 Γ(1)=∫0∞e−xdx=−e−x∣∣∣0∞=1
由此推出 Γ 函数的两个性质
Γ(z+1)=z!
(−z)!Γ(z)=sinπzπ
(3)
广义阶乘幂
zw=(z−w)!z!
zw=Γ(z)Γ(z+w)
(4)
(x)n=(−1)n⋅(−x)n
超几何表示
对于超几何函数
F(a1,a2,⋯,amb1,b2,⋯,bn∣∣∣∣∣∣ z)=k⩾0∑tk,tk=b1k⋯bnka1k⋯amk⋅k!zk
其中第一项是 t0=1,其他项可以从相邻比值中得到
tktk+1=a1ka2k⋯amka1k+1a2k+1⋯amk+1⋅b1k+1b2k+1⋯bnk+1b1kb2k⋯bnk
=(k+b1)(k+b2)⋯(k+bn)(k+1)(k+a1)(k+a2)⋯(k+am)⋅z
这是一个有理函数,给定的某个级数,写成超几何级数
k⩾0∑tk=t0⋅F(a1,a2,⋯,amb1,b2,⋯,bn∣∣∣∣∣∣ z)
这里有几点需要注意
如果 tk 中包含有 (a+k)! 这样的因子,将 (a+k+1) 写到 tktk+1 的相应位置
就是说如果 (a+k)! 出现在分母上,就将 (a+k+1) 写到 tktk+1 的分母上,否则就写到分母上
如果 tk 中包含有 (a−k)! 这样的因子,将 (a−k) 写到 tktk+1 的相反位置
就是说如果 (a−k)! 出现在分子上,就将 (a−k) 相反地写到 tktk+1 的分母上
另外注意 tktk+1 对应的 (a−k),实际上写成超几何函数的参数是 (−1)⋅(k−a),参数为 −a
分母没有 (k+1) 的要补上
几个恒等式的超几何表示
(1)
k⩽n∑(kn−k)zk=1+4z1((21+1+4z)n+1−(21−1+4z)n+1),n⩾0
来观察 k⩽n∑(kn−k)zk=k∑(2kn+k),这里的化归,大致如下
(kn−k)=k!(n−2k)!(n−k)!⟹2k′←(n−2k),k=2n−2k′,进一步化简得到
(2k′)!(n/2−k′)!(n/2+k′)!,令 n←n/2,k←k′
所以对应成 tk=(2kn+k)
t0=1
tktk+1=(k+1)(k+1/2)(n+k+1)(k−n)⋅4−1,超几何表示为
F(1+n,−n1/2∣∣∣∣∣∣ −z/4)=F(1+2⋅⌈n/2⌉,−n1/2∣∣∣∣∣∣ −z/4)
(2)
k=0∑2n(kr)(2n−kr)(−1)k=(−1)n(nr)
其中 t0=(2nr)=(2n)!(r−2n)!r!,这样 F=(−1)n(nr)r!(2n)!(r−2n)!=(−1)nn!(2n)!(r−n)!(r−2n)!
展开超几何函数
tktk+1=(k+1)(r−2n+1+k)(r−k)(2n−k)=F(−r,−2nr−2n+1∣∣∣∣∣∣ −1)
令 r−2n+1=c,那么有恒等式
F(1−c−2n,−2nc∣∣∣∣∣∣ −1)=(−1)nn!(2n)!(c+n−1)!(c−1)!,n⩾0
下面证明一个恒等式
(−1)nn!(2n)!=b→−2nlimb!(b/2)!=x→−nlim(2x)!x!
根据 (−z)!Γ(z)=sinπzπ
(−2n−2ϵ)!(−n−ϵ)!Γ(2n+2ϵ)Γ(n+ϵ)=sin(n+ϵ)πsin(2n+2ϵ)π
用 sin(x+y) 的展开,ϵ→0limcosnπsinϵπcos2nπsin2ϵπ=(−1)n(2+O(ϵ))
最后一步用了洛必达法则,分子分母同时对 dϵ 求导,得到 π2πcosϵπcos2ϵπ=2
从而我们有
(−2n−2ϵ)!(−n−ϵ)!=2(−1)n⋅Γ(n)Γ(2n)=2(−1)n(n−1)!(2n−1)!=(−1)nn!(2n)!
由此,如果令 b=−2n,代入 (2) 的超几何恒等式中
右边是 (−1)b/2(−b/2)!(−b)!(c−1−b/2)!(c−1)!=b!(b/2)!⋅(c−1)b/2
由此可以得到库默尔公式
F(a,b1+b−a∣∣∣∣∣∣ −1)=b!(b/2)!(b−a)b/2
(3)
范德蒙卷积的超几何表示
k∑(kr)(n−ks)=(nr+s)
其中 t0=(ns),可以推出
F(−r,−ns−n+1∣∣∣∣∣∣ 1)=(nr+s)/(ns)=(s)!(r+s)!(s−n+r)!(s−n)!
=(s−n+1)r(s+1)r=Γ(s+1)Γ(s+1+r)⋅Γ(s−n+1+r)Γ(s−n+1)
令 a=−r, b=−n,那么上式可以写成
F(a,bc∣∣∣∣∣∣ 1)=Γ(c−a)Γ(c−b)Γ(c−a−b)Γ(c),b⩽0,Rc⩾Ra+Rb
当 b=−n 时,可以用阶乘幂表示
F(a,−nc∣∣∣∣∣∣ 1)=cn(c−a)n=(−c)n(a−c)n,n⩾0
(4)
k∑(a+ka+b)(b+kb+c)(c+kc+a)(−1)k=a!b!c!(a+b+c)!,a,b,c⩾0
针对这个式子进行超几何表示,会比较麻烦
F=a!b!c!(a+b+c)!/t0=(a+b)!(b+c)!(a+c)!(a!b!c!)⋅(a+b+c)!
tktk+1=(k+a+1)(k+b+1)(k+c+1)(k−a)(k−b)(k−c)
令 k′=c+k,即式子中出现的 k 用 k−c 代换
F(−c−a,−c−b,−2ca−c+1,b−c+1∣∣∣∣∣∣ 1)
再做一替换,令 c=n,a−n+1=a′,b−n+1=b′,超几何表示为
F(1−2n−a,1−2n−b,−2na,b∣∣∣∣∣∣ 1),这里的 a,b 实际上是 a′,b′
那超几何恒等式右边等于多少呢?
(a+b+1)c⋅a!(c+a)!b!(c+b)!c!=(a′+b′+2n−1)n(a+1)n⋅(b+1)nn!
=(a′+b′+2n−1)n(a′+n)n⋅(b′+n)nn!
还可以进一步化简
(a′+n)n⋅(b′+n)nn!=(1−2n−a′)⋯(−a′−n)⋅(1−2n−b′)⋯(−b′−n)n!
注意到 (1−2n−a′)⋯(−a′−n)=(1−2n−a′)n,(1−2n−a′) 又恰好是超几何函数的上参数
上式 =(1−2n−a′)n(1−2n−b′)n(−2n)nn!(−2n)n
注意到 n!(−2n)n=(−1)n(2n)nn!=(−1)n(2n)!
上式的分母可以写成 tn⋅n!⋅anbn
由此 F=(−1)nn!(2n)!tn⋅anbn(a+b+2n−1)n
上式还不够完美,因为分子中有一个 tn,如果 tn=1 是不是就完美了呢?
实际上,我们可以从 k=n 开始求级数,即 ∑[k=⟨n,n+1,⋯⟩]
也就是说,对超几何级数 k′=k+n,k′⩾n
此时有 tn=1,并且
F=a!b!c!(a+b+c)!/tn=(−1)nn!(2n)!(a+b+1)n(a−n+1)n⋅(b−n+1)n1
作代换 a−n+1=a′,b−n+1=b′,即可得超几何恒等式
F(1−2n−a,1−2n−b,−2na,b∣∣∣∣∣∣ 1)=(−1)nn!(2n)!anbn(a+b+2n−1)n,n⩾0
推广到复数,就是迪克逊公式
F(a,b,c1+c−a,1+c−b∣∣∣∣∣∣ 1)=c!(c/2)!(c−a−b)c/2(c−a)c/2(c−b)c/2,Ra+Rb<1+Rc/2
(5)
k∑(km−r+s)(n−kn+r−s)(m+nr+k)=(mr)(ns)
写成超几何形式,不难发现
t0=(nn+r−s)(m+nr),超几何恒等式的右边为
(m+nr)(mr)(nn+r−s)(ns)=(r−m)n(m+1)n⋅cn(b−c)n
超几何表示为
F(r−s−m,r+1,−nr−s+1,r−m−n+1∣∣∣∣∣∣ 1)
令 a=r−s−m,b=r+1,c=r−s+1,可以推出 Saalschütz 恒等式
F(a,b,−nc,a+b−c−n+1∣∣∣∣∣∣ 1)=(c)n(c−a)n⋅(c−a−b)n(c−b)n
=(−c)n(a+b−c)n(a−c)n(b−c)n,n⩾0
假设上参数为 ak,下参数为 bk,恒等式给出,超几何函数满足 b1+b2=a1+a2+a3+1 时的封闭形式
推广
(a1−b1)F(a1,a2,⋯,amb1+1,b2,⋯bn∣∣∣∣∣∣ z)=a1F(a1+1,a2,⋯,amb1+1,b2,⋯,bn∣∣∣∣∣∣ z)−b1F(a1,a2,⋯,amb1,b2,⋯,bn∣∣∣∣∣∣ z)
实际上,等价于证明,对于第 k 项而言,我们有
(a−b)(b+1)kak=a(b+1)k(a+1)k−bbkak
这个不等式很好证明,右边通分即可,分母为 bk(b+1)k,分子为 bk⋅(ak(a+k)−ak(b+k))
化简完就等于左边