超几何函数



在数学中,高斯超几何函数普通超几何函数2F1(a,b;c;z)是一个用超几何级数定义的函数,很多特殊函数都是它的特例或极限。所有具有三个正则奇点英语Regular singular point的二阶线性常微分方程的解都可以用超几何函数表示。




目录





  • 1 超几何级数


  • 2 特殊情形


  • 3 超几何方程

    • 3.1 正则奇点 0 附近的解


    • 3.2 正则奇点 1 附近的解


    • 3.3 正则奇点 ∞ 附近的解


    • 3.4 李代数参数与连接关系



  • 4 积分表示

    • 4.1 证明



  • 5 变换公式

    • 5.1 分式线性变换

      • 5.1.1 Pfaff 变换

        • 5.1.1.1 证明



      • 5.1.2 Euler 变换



    • 5.2 二次变换

      • 5.2.1 证明


      • 5.2.2 其它例子



    • 5.3 三次及高次变换



  • 6 特殊值

    • 6.1 z=0


    • 6.2 z=1


    • 6.3 z=-1


    • 6.4 z=1/2



  • 7 参考文献




超几何级数


cdisplaystyle cc不是0,-1,-2...时,对于|z| < 1,超几何函数可用如下幂级数定义


2F1(a,b;c;z)=∑n=0∞a(n)b(n)c(n)znn!displaystyle ,_2F_1(a,b;c;z)=sum _n=0^infty a^(n)b^(n) over c^(n),z^n over n!,_2F_1(a,b;c;z)=sum _n=0^infty a^(n)b^(n) over c^(n),z^n over n!


其中  x(n)displaystyle x^(n) x^(n) 是Pochhammer符号,定义为:


q(n)={1if n=0q(q+1)⋯(q+n−1)if n>0{displaystyle q^(n)=leftbeginarrayll1&mboxif n=0\q(q+1)cdots (q+n-1)&mboxif n>0endarrayright.{displaystyle q^(n)=leftbeginarrayll1&mboxif n=0\q(q+1)cdots (q+n-1)&mboxif n>0endarrayright.

ab0或负整数时级数只有有限项。


对于满足|z| ≥ 1 的复数z,超几何函数可以通过将上述在单位圆内定义的函数沿着避开支点01的任意路径做解析延拓来得到。



特殊情形


很多普通的数学函数可以用超几何函数或它的极限表示出来,一些典型的例子如下:



ln⁡(1+z)=z2F1(1,1;2;−z)displaystyle ln(1+z)=z,_2F_1(1,1;2;-z)ln(1+z)=z,_2F_1(1,1;2;-z).

(1−z)−a=2F1(a,1;1;z)displaystyle (1-z)^-a=,_2F_1(a,1;1;z)(1-z)^-a=,_2F_1(a,1;1;z)

arcsin⁡z=z2F1(12,12;32;z2)displaystyle arcsin z=z,_2F_1left(tfrac 12,tfrac 12;tfrac 32;z^2right)arcsin z=z,_2F_1left(tfrac 12,tfrac 12;tfrac 32;z^2right)

合流超几何函数(Kummer函数)可以用超几何函数的极限表示如下


M(a,c,z)=limb→∞2F1(a,b;c;b−1z)displaystyle M(a,c,z)=lim _brightarrow infty _2F_1(a,b;c;b^-1z)M(a,c,z)=lim _brightarrow infty _2F_1(a,b;c;b^-1z)

因此,所有合流超几何函数的特例,例如贝塞尔函数都可以表示成超几何函数的极限。


勒让德函数是有3个正则奇点的二阶线性常微分方程的解,可以用以不同的形式用超几何函数表示,例如


2F1(a,1−a;c;z)=Γ(c)z1−c2(1−z)c−12P−a1−c(1−2z)displaystyle _2F_1(a,1-a;c;z)=Gamma (c)z^tfrac 1-c2(1-z)^tfrac c-12P_-a^1-c(1-2z)_2F_1(a,1-a;c;z)=Gamma (c)z^tfrac 1-c2(1-z)^tfrac c-12P_-a^1-c(1-2z)


很多多项式,例如贾可比多项式 P(α,β)
n
及其特殊情形勒让德多项式, 车比雪夫多项式, Gegenbauer多项式都能用超几何函数表示


2F1(−n,α+1+β+n;α+1;x)=n!(α+1)nPn(α,β)(1−2x)displaystyle _2F_1(-n,alpha +1+beta +n;alpha +1;x)=frac n!(alpha +1)_nP_n^(alpha ,beta )(1-2x)_2F_1(-n,alpha +1+beta +n;alpha +1;x)=frac n!(alpha +1)_nP_n^(alpha ,beta )(1-2x)
其它特殊情形还包括Krawtchouk多项式, Meixner多项式, Meixner–Pollaczek多项式。


椭圆模函数英语Elliptic modular function有时能表示成参数a, b, c是1, 1/2, 1/3, ... 或 0 的超几何函数之比的反函数。例如,若


τ=i2F1(12,12;1;1−z)2F1(12,12;1;z)displaystyle tau =rm ifrac _2F_1(frac 12,frac 12;1;1-z)_2F_1(frac 12,frac 12;1;z)tau =rm ifrac _2F_1(frac 12,frac 12;1;1-z)_2F_1(frac 12,frac 12;1;z)



z=κ2(τ)=θ2(τ)4θ3(τ)4displaystyle z=kappa ^2(tau )=frac theta _2(tau )^4theta _3(tau )^4z=kappa ^2(tau )=frac theta _2(tau )^4theta _3(tau )^4

τ的椭圆模函数.


不完整的beta函数 Bx(p,q) 表示成


Bx(p,q)=xpp2F1(p,1−q;p+1;x)displaystyle B_x(p,q)=frac x^pp_2F_1(p,1-q;p+1;x)B_x(p,q)=frac x^pp_2F_1(p,1-q;p+1;x)

完整的椭圆积分 KE 如下给出


K(k)=π22F1(12,12;1;k2)displaystyle K(k)=tfrac pi 2,_2F_1left(tfrac 12,tfrac 12;1;k^2right)K(k)=tfrac pi 2,_2F_1left(tfrac 12,tfrac 12;1;k^2right)
E(k)=π22F1(−12,12;1;k2)displaystyle E(k)=tfrac pi 2,_2F_1left(-tfrac 12,tfrac 12;1;k^2right)E(k)=tfrac pi 2,_2F_1left(-tfrac 12,tfrac 12;1;k^2right)


超几何方程


超几何函数满足的微分方程称为超几何方程,其形式为(参见广义超几何函数)



z(zddz+a)(zddz+b)w=zddz(zddz+c−1)w,w(z)=2F1(a,b;c;z)displaystyle zleft(zfrac rm drm dz+aright)left(zfrac rm drm dz+bright)w=zfrac rm drm dzleft(zfrac rm drm dz+c-1right)w,quad w(z)=_2F_1(a,b;c;z)zleft(zfrac rm drm dz+aright)left(zfrac rm drm dz+bright)w=zfrac rm drm dzleft(zfrac rm drm dz+c-1right)w,quad w(z)=_2F_1(a,b;c;z).

展开后,得


z(1−z)d2wdz2+[c−(a+b+1)z]dwdz−abw=0.displaystyle z(1-z)frac mathrm d ^2wmathrm d z^2+left[c-(a+b+1)zright]frac mathrm d wmathrm d z-abw=0.z(1-z)frac mathrm d^2wmathrm dz^2+left[c-(a+b+1)zright]frac mathrm dwmathrm dz-abw=0.

它有三个正则奇点:0, 1, ∞.



正则奇点 0 附近的解


超几何方程的指标方程英语Frobenius method


ρ(ρ−1)+cρ=0displaystyle rho (rho -1)+crho =0rho (rho -1)+crho =0

它的两个指标 ρ 是 0 和 1-c


c不是整数时,超几何方程在 0 附近的两个线性无关的正则特解为:


2F1(a,b;c;z) and z1−c2F1(1+a−c,1+b−c;2−c;z)displaystyle ,_2F_1(a,b;c;z)text and z^1-c,_2F_1(1+a-c,1+b-c;2-c;z),_2F_1(a,b;c;z)text and z^1-c,_2F_1(1+a-c,1+b-c;2-c;z)

c 为 1 时,方程只有一个正则解。当 c 为其余整数时,另一个线性无关的正则特解涉及对数项。


事实上,当 c 为整数时,另一个线性无关的特解总可以选取为 Meijer G-函数:


2F1(a,b;c;z) and G2,22,0(1−a,1−b;0,c−1;z), if c∈Z+displaystyle ,_2F_1(a,b;c;z)text and ,G_2,2^2,0(1-a,1-b;0,c-1;z),text if cin mathbb Z ^+,_2F_1(a,b;c;z)text and ,G_2,2^2,0(1-a,1-b;0,c-1;z),text if cin mathbb Z^+

z1−c2F1(1+a−c,1+b−c;2−c;z) and G2,22,0(1−a,1−b;0,1−c;z), if c∈Z0−displaystyle ,z^1-c,_2F_1(1+a-c,1+b-c;2-c;z)text and ,G_2,2^2,0(1-a,1-b;0,1-c;z),text if cin mathbb Z _0^-,z^1-c,_2F_1(1+a-c,1+b-c;2-c;z)text and ,G_2,2^2,0(1-a,1-b;0,1-c;z),text if cin mathbb Z_0^-


正则奇点 1 附近的解


只需作代换 t=1-z,方程变为:


t(1−t)d2wdt2+[1+a+b−c−(a+b+1)t]dwdt−abw=0.displaystyle t(1-t)frac d^2wdt^2+left[1+a+b-c-(a+b+1)tright]frac dwdt-abw=0.t(1-t)frac d^2wdt^2+left[1+a+b-c-(a+b+1)tright]frac dwdt-abw=0.

a+b-c 不是整数时,两个线性无关的正则特解为:


2F1(a,b;1+a+b−c;1−z) and (1−z)c−a−b2F1(c−b,c−a;1−a−b+c;1−z)displaystyle ,_2F_1(a,b;1+a+b-c;1-z)text and (1-z)^c-a-b,_2F_1(c-b,c-a;1-a-b+c;1-z),_2F_1(a,b;1+a+b-c;1-z)text and (1-z)^c-a-b,_2F_1(c-b,c-a;1-a-b+c;1-z)


正则奇点 ∞ 附近的解


a-b 不是整数时,两个线性无关的正则特解为:


z−a2F1(a,1+a−c;1+a−b;z−1) and z−b2F1(b,1+b−c;1+b−a;z−1).displaystyle z^-a,_2F_1left(a,1+a-c;1+a-b;z^-1right)text and z^-b,_2F_1left(b,1+b-c;1+b-a;z^-1right).z^-a,_2F_1left(a,1+a-c;1+a-b;z^-1right)text and z^-b,_2F_1left(b,1+b-c;1+b-a;z^-1right).


李代数参数与连接关系


在讨论超几何方程的解的连接关系的时候,采用另外一套参数[1]会更加方便。这组参数是根据方程在三个正则奇点处的指标之差来定义的。


Fα,β,μ(z)=2F1(a,b;c;z)displaystyle F_alpha ,beta ,mu (z)=_2F_1(a,b;c;z)F_alpha ,beta ,mu (z)=_2F_1(a,b;c;z)

α=c−1,β=a+b−c,μ=b−adisplaystyle alpha =c-1,beta =a+b-c,mu =b-aalpha =c-1,beta =a+b-c,mu =b-a

a=1+α+β−μ2,b=1+α+β+μ2,c=1+αdisplaystyle a=frac 1+alpha +beta -mu 2,b=frac 1+alpha +beta +mu 2,c=1+alpha a=frac 1+alpha +beta -mu 2,b=frac 1+alpha +beta +mu 2,c=1+alpha

参数 α,β,γ 称为李代数参数。


运用李代数参数,超几何方程在三个正则奇点处的正则解可以分别表示为:


At 0:Fα,β,μ(z) and z−αF−α,β,−μ(z)displaystyle textAt 0:F_alpha ,beta ,mu (z)text and z^-alpha F_-alpha ,beta ,-mu (z)textAt 0:F_alpha ,beta ,mu (z)text and z^-alpha F_-alpha ,beta ,-mu (z)

At 1:Fβ,α,μ(1−z) and (1−z)−βF−β,α,−μ(1−z)displaystyle textAt 1:F_beta ,alpha ,mu (1-z)text and (1-z)^-beta F_-beta ,alpha ,-mu (1-z)textAt 1:F_beta ,alpha ,mu (1-z)text and (1-z)^-beta F_-beta ,alpha ,-mu (1-z)

At ∞:(−z)−1−α−β+μ2F−μ,β,−α(z−1) and (−z)−1−α−β−μ2Fμ,β,α(z−1)displaystyle textAt infty :(-z)^frac -1-alpha -beta +mu 2F_-mu ,beta ,-alpha (z^-1)text and (-z)^frac -1-alpha -beta -mu 2F_mu ,beta ,alpha (z^-1)textAt infty :(-z)^frac -1-alpha -beta +mu 2F_-mu ,beta ,-alpha (z^-1)text and (-z)^frac -1-alpha -beta -mu 2F_mu ,beta ,alpha (z^-1)

从上面的表达式可见,李代数参数比起通常用的参数 a,b,c 的优势在于能够体现不同区域的解之间的对称性。


引入记号:


G(m;n,p)=πsin⁡mπΓ(n)Γ(p)=G(m;p,n)displaystyle G(m;n,p)=frac pi sin mpi Gamma (n)Gamma (p)=G(m;p,n)G(m;n,p)=frac pi sin mpi Gamma (n)Gamma (p)=G(m;p,n)

Fα,β,μ(z)=1Γ(1+α)Fα,β,μ(z)displaystyle mathbf F _alpha ,beta ,mu (z)=frac 1Gamma (1+alpha )F_alpha ,beta ,mu (z)mathbf F_alpha ,beta ,mu (z)=frac 1Gamma (1+alpha )F_alpha ,beta ,mu (z)

则超几何方程在不同区域的解的连接关系可以表示为:


Fβ,α,μ(1−z)=G(−α;a−α,b−α)Fα,β,μ(z)+G(α;a,b)z−αF−α,β,−μ(z),displaystyle mathbf F _beta ,alpha ,mu (1-z)=G(-alpha ;a-alpha ,b-alpha )mathbf F _alpha ,beta ,mu (z)+G(alpha ;a,b)z^-alpha mathbf F _-alpha ,beta ,-mu (z),mathbf F_beta ,alpha ,mu (1-z)=G(-alpha ;a-alpha ,b-alpha )mathbf F_alpha ,beta ,mu (z)+G(alpha ;a,b)z^-alpha mathbf F_-alpha ,beta ,-mu (z),

(1−z)−βF−β,α,−μ(1−z)=G(−α;1−a,1−b)Fα,β,μ(z)+G(α;b−β,a−β)z−αF−α,β,−μ(z)=G(−α;1−a,1−b)Fα,β,μ(z)+G(α;1−(a−α),1−(b−α))z−αF−α,β,−μ(z);displaystyle beginarrayrcl(1-z)^-beta mathbf F _-beta ,alpha ,-mu (1-z)&=&G(-alpha ;1-a,1-b)mathbf F _alpha ,beta ,mu (z)+G(alpha ;b-beta ,a-beta )z^-alpha mathbf F _-alpha ,beta ,-mu (z)\&=&G(-alpha ;1-a,1-b)mathbf F _alpha ,beta ,mu (z)+G(alpha ;1-(a-alpha ),1-(b-alpha ))z^-alpha mathbf F _-alpha ,beta ,-mu (z);endarraybeginarrayrcl(1-z)^-beta mathbf F_-beta ,alpha ,-mu (1-z)&=&G(-alpha ;1-a,1-b)mathbf F_alpha ,beta ,mu (z)+G(alpha ;b-beta ,a-beta )z^-alpha mathbf F_-alpha ,beta ,-mu (z)\&=&G(-alpha ;1-a,1-b)mathbf F_alpha ,beta ,mu (z)+G(alpha ;1-(a-alpha ),1-(b-alpha ))z^-alpha mathbf F_-alpha ,beta ,-mu (z);endarray

(−z)−aF−μ,β,−α(z−1)=G(−α;1−b,a−α)Fα,β,μ(z)+G(α;a,a−β)z−αF−α,β,−μ(z),displaystyle (-z)^-amathbf F _-mu ,beta ,-alpha (z^-1)=G(-alpha ;1-b,a-alpha )mathbf F _alpha ,beta ,mu (z)+G(alpha ;a,a-beta )z^-alpha mathbf F _-alpha ,beta ,-mu (z),(-z)^-amathbf F_-mu ,beta ,-alpha (z^-1)=G(-alpha ;1-b,a-alpha )mathbf F_alpha ,beta ,mu (z)+G(alpha ;a,a-beta )z^-alpha mathbf F_-alpha ,beta ,-mu (z),

(−z)−bFμ,β,α(z−1)=G(−α;1−a,b−α)Fα,β,μ(z)+G(α;b,b−β)z−αF−α,β,−μ(z)=G(−α;1−a,1−(a−β))Fα,β,μ(z)+G(α;b,1−(a−α))z−αF−α,β,−μ(z).displaystyle beginarrayrcl(-z)^-bmathbf F _mu ,beta ,alpha (z^-1)&=&G(-alpha ;1-a,b-alpha )mathbf F _alpha ,beta ,mu (z)+G(alpha ;b,b-beta )z^-alpha mathbf F _-alpha ,beta ,-mu (z)\&=&G(-alpha ;1-a,1-(a-beta ))mathbf F _alpha ,beta ,mu (z)+G(alpha ;b,1-(a-alpha ))z^-alpha mathbf F _-alpha ,beta ,-mu (z).endarraybeginarrayrcl(-z)^-bmathbf F_mu ,beta ,alpha (z^-1)&=&G(-alpha ;1-a,b-alpha )mathbf F_alpha ,beta ,mu (z)+G(alpha ;b,b-beta )z^-alpha mathbf F_-alpha ,beta ,-mu (z)\&=&G(-alpha ;1-a,1-(a-beta ))mathbf F_alpha ,beta ,mu (z)+G(alpha ;b,1-(a-alpha ))z^-alpha mathbf F_-alpha ,beta ,-mu (z).endarray

分别对比两组式子最后一个等号之后的部分,可以看出每组的两个式子之间的对称性。


完整的连接关系表称为 Kummer 表,上面四式是 Kummer 表的一部分。



积分表示


B(a,c−a)2F1(a,b;c;z)=∫1∞tb−c(t−1)c−a−1(t−z)−bdt,ℜ(c)>ℜ(a)>0,|arg⁡(1−z)|<π<pi mathrmB (a,c-a)_2F_1(a,b;c;z)=int _1^infty t^b-c(t-1)^c-a-1(t-z)^-bmathrm dt,Re (c)>Re (a)>0,|arg(1-z)|<pi

式中的 Β 是beta函数。



证明


可以证明等号右边的表达式是超几何方程的解。再考虑这个解在 z=0 附近的性质,可以确定它的具体形式。




p(a,b,c;t,z)=tb−c(t−1)c−a−1(t−z)−b−2,w(a,b,c;t,z)=(t−z)2p(a,b,c;t,z);displaystyle p(a,b,c;t,z)=t^b-c(t-1)^c-a-1(t-z)^-b-2,quad w(a,b,c;t,z)=(t-z)^2p(a,b,c;t,z);p(a,b,c;t,z)=t^b-c(t-1)^c-a-1(t-z)^-b-2,quad w(a,b,c;t,z)=(t-z)^2p(a,b,c;t,z);



∂w∂z=b(t−z)p(a,b,c;t,z),∂2w∂z2=b(b+1)p(a,b,c;t,z)displaystyle frac partial wpartial z=b(t-z)p(a,b,c;t,z),quad frac partial ^2wpartial z^2=b(b+1)p(a,b,c;t,z)frac partial wpartial z=b(t-z)p(a,b,c;t,z),quad frac partial ^2wpartial z^2=b(b+1)p(a,b,c;t,z)

z(1−z)∂2w∂z2+[c−(a+b+1)z]∂w∂z−abw=bp(a,b,c;t,z)z(1−z)(b+1)+[c−(a+b+1)z](t−z)−a(t−z)2=bp(a,b,c;t,z)−at2+[c−(b−a+1)z]t+(b−c+1)z=bp(a,b,c;t,z)(b−c+1)(t−1)(t−z)+(c−a)t(t−z)+(−b−1)t(t−1)=b∂∂t[t(t−1)(t−z)p(a,b,c;t,z)],displaystyle beginarraycl&z(1-z)frac partial ^2wpartial z^2+left[c-(a+b+1)zright]frac partial wpartial z-abw\=&bp(a,b,c;t,z)leftz(1-z)(b+1)+[c-(a+b+1)z](t-z)-a(t-z)^2right\=&bp(a,b,c;t,z)left-at^2+[c-(b-a+1)z]t+(b-c+1)zright\=&bp(a,b,c;t,z)left(b-c+1)(t-1)(t-z)+(c-a)t(t-z)+(-b-1)t(t-1)right\=&bfrac partial partial t[t(t-1)(t-z)p(a,b,c;t,z)],endarraybeginarraycl&z(1-z)frac partial ^2wpartial z^2+left[c-(a+b+1)zright]frac partial wpartial z-abw\=&bp(a,b,c;t,z)leftz(1-z)(b+1)+[c-(a+b+1)z](t-z)-a(t-z)^2right\=&bp(a,b,c;t,z)left-at^2+[c-(b-a+1)z]t+(b-c+1)zright\=&bp(a,b,c;t,z)left(b-c+1)(t-1)(t-z)+(c-a)t(t-z)+(-b-1)t(t-1)right\=&bfrac partial partial t[t(t-1)(t-z)p(a,b,c;t,z)],endarray

上式中的第二、三个等号可以通过直接展开大括号内的多项式乘积得到。上式两边分别对 t 从 1 到无穷大进行积分,等号右边为 0,于是我们证明了上面的积分表达式的确是超几何方程的解。


另一方面,利用二项式定理,积分表达式等号右边的部分可以按 z 展开成幂级数,故可知等号右边应取 C 2F1(a,b,c;z) 的形式(因为另一个线性无关的特解无法展开成幂级数),其中 C 为待定的常数。


对比积分表达式在 z=0 处的值与 Β 函数的定义,即可确定常数 C



变换公式



分式线性变换



Pfaff 变换


Pfaff 变换将正则奇点 1 和 ∞ 交换(也就是将李代数参数中的 βμ 对换):


2F1(a,b;c;z)=(1−z)−b2F1(c−a,b;c;zz−1),|arg⁡(1−z)|<π<pi _2F_1(a,b;c;z)=(1-z)^-b,_2F_1(c-a,b;c;tfrac zz-1),quad |arg(1-z)|<pi

a,b 的对称性自然有:


2F1(a,b;c;z)=(1−z)−a2F1(a,c−b;c;zz−1),|arg⁡(1−z)|<πarg(1-z)_2F_1(a,b;c;z)=(1-z)^-a,_2F_1(a,c-b;c;tfrac zz-1),quad |arg(1-z)|<pi

证明

Pfaff 变换可以根据超几何方程得到。事实上,令


u=zz−1=1+1z−1displaystyle u=tfrac zz-1=1+tfrac 1z-1u=tfrac zz-1=1+tfrac 1z-1



z=uu−1,(1−z)a=(1−u)−a,dudz=−(1−u)2,d2udz2=−(1−u)3displaystyle z=tfrac uu-1,quad (1-z)^a=(1-u)^-a,quad tfrac mathrm d umathrm d z=-(1-u)^2,quad tfrac mathrm d ^2umathrm d z^2=-(1-u)^3z=tfrac uu-1,quad (1-z)^a=(1-u)^-a,quad tfrac mathrm dumathrm dz=-(1-u)^2,quad tfrac mathrm d^2umathrm dz^2=-(1-u)^3
z(1−z)d2dz2[(1−z)−bw]+[c−(a+b+1)z]ddz[(1−z)−bw]−ab(1−z)−bw=(1−z)−b−1z[b(b+1)+2b(1−z)ddz+(1−z)2d2dz2]+[c−(a+b+1)z][b+(1−z)ddz]−ab(1−z)w=(1−z)−b−1z(1−z)2d2dz2+(1−z)[c−(a−b+1)z]ddz+b(c−a)w=(1−u)b+1−u(1−u)d2du2+2uddu−(1−u)[c+(a−b+1)(1−u)−1u]ddu+b(c−a)w=−(1−u)b+1u(1−u)d2du2+[c−(c−a+b+1)u]ddu−b(c−a)wdisplaystyle beginarraycl&z(1-z)tfrac mathrm d ^2mathrm d z^2[(1-z)^-bw]+left[c-(a+b+1)zright]tfrac mathrm d mathrm d z[(1-z)^-bw]-ab(1-z)^-bw\=&(1-z)^-b-1leftz[b(b+1)+2b(1-z)tfrac mathrm d mathrm d z+(1-z)^2tfrac mathrm d ^2mathrm d z^2]+[c-(a+b+1)z][b+(1-z)tfrac mathrm d mathrm d z]-ab(1-z)rightw\=&(1-z)^-b-1leftz(1-z)^2tfrac mathrm d ^2mathrm d z^2+(1-z)[c-(a-b+1)z]tfrac mathrm d mathrm d z+b(c-a)rightw\=&(1-u)^b+1left-u(1-u)tfrac mathrm d ^2mathrm d u^2+2utfrac mathrm d mathrm d u-(1-u)[c+(a-b+1)(1-u)^-1u]tfrac mathrm d mathrm d u+b(c-a)rightw\=&-(1-u)^b+1leftu(1-u)tfrac mathrm d ^2mathrm d u^2+[c-(c-a+b+1)u]tfrac mathrm d mathrm d u-b(c-a)rightwendarraybeginarraycl&z(1-z)tfrac mathrm d^2mathrm dz^2[(1-z)^-bw]+left[c-(a+b+1)zright]tfrac mathrm dmathrm dz[(1-z)^-bw]-ab(1-z)^-bw\=&(1-z)^-b-1leftz[b(b+1)+2b(1-z)tfrac mathrm dmathrm dz+(1-z)^2tfrac mathrm d^2mathrm dz^2]+[c-(a+b+1)z][b+(1-z)tfrac mathrm dmathrm dz]-ab(1-z)rightw\=&(1-z)^-b-1leftz(1-z)^2tfrac mathrm d^2mathrm dz^2+(1-z)[c-(a-b+1)z]tfrac mathrm dmathrm dz+b(c-a)rightw\=&(1-u)^b+1left-u(1-u)tfrac mathrm d^2mathrm du^2+2utfrac mathrm dmathrm du-(1-u)[c+(a-b+1)(1-u)^-1u]tfrac mathrm dmathrm du+b(c-a)rightw\=&-(1-u)^b+1leftu(1-u)tfrac mathrm d^2mathrm du^2+[c-(c-a+b+1)u]tfrac mathrm dmathrm du-b(c-a)rightwendarray



w=2F1(c−a,b;c;u)displaystyle w=_2F_1(c-a,b;c;u)w=_2F_1(c-a,b;c;u)

w(u) 满足的超几何方程知等号右边为 0,再考虑函数 (1-z)-bw(z)z=0 附近的性质即可得到 Pfaff 变换的公式。



Euler 变换


Pfaff 变换可以导出 Euler 变换,它将李代数参数 β 变成 -β


2F1(a,b;c;z)=(1−z)−b2F1(c−a,b;c;zz−1)=(1−z)−b(1−zz−1)a−c2F1(c−a,c−b;c;zz−1zz−1−1)=(1−z)c−a−b2F1(c−a,c−b;c;z),|arg⁡(1−z)|<πdisplaystyle <pi endarray<pi endarray

Pfaff 变换和 Euler 变换都是分式线性变换的例子,这得名于等式两边的超几何函数的宗量的联系,参见莫比乌斯变换。


将上面提到的四个连接关系与 Pfaff 变换及 Euler 变换组合起来,就得到完整的 Kummer 表。


给定一组李代数参数(α,β,μ),(±α,±β,±μ) 及其轮换对应着 24 个不同但彼此关联的超几何函数(Fα,β,μ 恒等于 Fα,β,),利用前面提到的四个连接关系和 Pfaff 变换,它们中的任意一个可以通过任意另外两个表出。


例如 Euler 变换可以表示为:


Fα,β,μ→PfaffFα,μ,β≡Fα,μ,−β→PfaffFα,−β,μdisplaystyle F_alpha ,beta ,mu xrightarrow textPfaffF_alpha ,mu ,beta equiv F_alpha ,mu ,-beta xrightarrow textPfaffF_alpha ,-beta ,mu F_alpha ,beta ,mu xrightarrow textPfaffF_alpha ,mu ,beta equiv F_alpha ,mu ,-beta xrightarrow textPfaffF_alpha ,-beta ,mu


二次变换


下面是一个二次变换的例子:


2F1(a,b;2a;z)=(1−z)−b22F1(a−b2,b2;a+12;z24z−4),|arg⁡(1−z)|<πarg(1-z)_2F_1(a,b;2a;z)=(1-z)^-tfrac b2,_2F_1(a-tfrac b2,tfrac b2;a+tfrac 12;tfrac z^24z-4),quad |arg(1-z)|<pi

二次变换得名于等号两边超几何函数宗量的联系(一个二次函数和一个莫比乌斯变换的组合)。



证明


仿照上面 Pfaff 变换的证明,有:


z(1−z)d2dz2[(1−z)−b2w]+[c−(a+b+1)z]ddz[(1−z)−b2w]−ab(1−z)−b2w=(1−z)−b2−1z[b2(b2+1)+b(1−z)ddz+(1−z)2d2dz2]+[c−(a+b+1)z][b2+(1−z)ddz]−ab(1−z)w=(1−z)−b2−1z(1−z)2d2dz2+(1−z)[c−(a+1)z]ddz+b4[2(c−2a)+(2a−b)z]wdisplaystyle beginarraycl&z(1-z)tfrac mathrm d ^2mathrm d z^2[(1-z)^-tfrac b2w]+left[c-(a+b+1)zright]tfrac mathrm d mathrm d z[(1-z)^-tfrac b2w]-ab(1-z)^-tfrac b2w\=&(1-z)^-tfrac b2-1leftz[tfrac b2(tfrac b2+1)+b(1-z)tfrac mathrm d mathrm d z+(1-z)^2tfrac mathrm d ^2mathrm d z^2]+[c-(a+b+1)z][tfrac b2+(1-z)tfrac mathrm d mathrm d z]-ab(1-z)rightw\=&(1-z)^-tfrac b2-1leftz(1-z)^2tfrac mathrm d ^2mathrm d z^2+(1-z)[c-(a+1)z]tfrac mathrm d mathrm d z+tfrac b4[2(c-2a)+(2a-b)z]rightw\endarraybeginarraycl&z(1-z)tfrac mathrm d^2mathrm dz^2[(1-z)^-tfrac b2w]+left[c-(a+b+1)zright]tfrac mathrm dmathrm dz[(1-z)^-tfrac b2w]-ab(1-z)^-tfrac b2w\=&(1-z)^-tfrac b2-1leftz[tfrac b2(tfrac b2+1)+b(1-z)tfrac mathrm dmathrm dz+(1-z)^2tfrac mathrm d^2mathrm dz^2]+[c-(a+b+1)z][tfrac b2+(1-z)tfrac mathrm dmathrm dz]-ab(1-z)rightw\=&(1-z)^-tfrac b2-1leftz(1-z)^2tfrac mathrm d^2mathrm dz^2+(1-z)[c-(a+1)z]tfrac mathrm dmathrm dz+tfrac b4[2(c-2a)+(2a-b)z]rightw\endarray



c=2a,u=z24z−4=14(z+1−11−z)displaystyle c=2a,quad u=tfrac z^24z-4=tfrac 14(z+1-tfrac 11-z)c=2a,quad u=tfrac z^24z-4=tfrac 14(z+1-tfrac 11-z)



1−u=(z−2)24(1−z),dudz=z(z−2)4(1−z)2,d2udz2=−12(1−z)3displaystyle 1-u=tfrac (z-2)^24(1-z),quad tfrac mathrm d umathrm d z=tfrac z(z-2)4(1-z)^2,quad tfrac mathrm d ^2umathrm d z^2=-tfrac 12(1-z)^31-u=tfrac (z-2)^24(1-z),quad tfrac mathrm dumathrm dz=tfrac z(z-2)4(1-z)^2,quad tfrac mathrm d^2umathrm dz^2=-tfrac 12(1-z)^3
z(1−z)2d2dz2+(1−z)[c−(a+1)z]ddz+b4[2(c−2a)+(2a−b)z]=z(1−z)2d2dz2+(1−z)[2a−(a+1)z]ddz+b2(a−b2)z=z3(z−2)216(1−z)2d2du2−z2(1−z)ddu+z(z−2)(2a−az−z)4(1−z)ddu+b2(a−b2)z=−zu(1−u)d2du2+[a+12−(a+1)u]ddu−b2(a−b2)displaystyle beginarraycl&z(1-z)^2tfrac mathrm d ^2mathrm d z^2+(1-z)[c-(a+1)z]tfrac mathrm d mathrm d z+tfrac b4[2(c-2a)+(2a-b)z]\=&z(1-z)^2tfrac mathrm d ^2mathrm d z^2+(1-z)[2a-(a+1)z]tfrac mathrm d mathrm d z+tfrac b2(a-tfrac b2)z\=&tfrac z^3(z-2)^216(1-z)^2tfrac mathrm d ^2mathrm d u^2-tfrac z2(1-z)tfrac mathrm d mathrm d u+tfrac z(z-2)(2a-az-z)4(1-z)tfrac mathrm d mathrm d u+tfrac b2(a-tfrac b2)z\=&-zleftu(1-u)tfrac mathrm d ^2mathrm d u^2+[a+tfrac 12-(a+1)u]tfrac mathrm d mathrm d u-tfrac b2(a-tfrac b2)rightendarraybeginarraycl&z(1-z)^2tfrac mathrm d^2mathrm dz^2+(1-z)[c-(a+1)z]tfrac mathrm dmathrm dz+tfrac b4[2(c-2a)+(2a-b)z]\=&z(1-z)^2tfrac mathrm d^2mathrm dz^2+(1-z)[2a-(a+1)z]tfrac mathrm dmathrm dz+tfrac b2(a-tfrac b2)z\=&tfrac z^3(z-2)^216(1-z)^2tfrac mathrm d^2mathrm du^2-tfrac z2(1-z)tfrac mathrm dmathrm du+tfrac z(z-2)(2a-az-z)4(1-z)tfrac mathrm dmathrm du+tfrac b2(a-tfrac b2)z\=&-zleftu(1-u)tfrac mathrm d^2mathrm du^2+[a+tfrac 12-(a+1)u]tfrac mathrm dmathrm du-tfrac b2(a-tfrac b2)rightendarray



w=2F1(a−b2,b2;a+12;u)displaystyle w=,_2F_1(a-tfrac b2,tfrac b2;a+tfrac 12;u)w=,_2F_1(a-tfrac b2,tfrac b2;a+tfrac 12;u)

仿照上面关于 Pfaff 变换的讨论,可得二次变换的公式。



其它例子


运用李代数参数,一般的二次变换可以表示为


Fα,β,μ(z)=f(z)Fα′,β′,μ′(g(z)),P(z)displaystyle F_alpha ,beta ,mu (z)=f(z)F_alpha ',beta ',mu '(g(z)),quad P(z)F_alpha ,beta ,mu (z)=f(z)F_alpha ',beta ',mu '(g(z)),quad P(z)

其中 f(z),g(z) 是 z 的函数, P(z) 表示 z 要满足的约束。


下表给出了一些二次变换。






















李代数参数(左)李代数参数(右)f(z)displaystyle f(z)f(z)g(z)displaystyle g(z)g(z)
P(z)displaystyle P(z)P(z)
α,μ,μdisplaystyle alpha ,mu ,mu alpha ,mu ,mu α2,μ,12displaystyle tfrac alpha 2,mu ,tfrac 12tfrac alpha 2,mu ,tfrac 12(1−12z)−bdisplaystyle (1-tfrac 12z)^-b(1-tfrac 12z)^-b(z2−z)2displaystyle left(tfrac z2-zright)^2left(tfrac z2-zright)^2
|arg⁡(1−z)|<πarg(1-z)|arg(1-z)|<pi
μ,β,μdisplaystyle mu ,beta ,mu mu ,beta ,mu μ,β2,12displaystyle mu ,tfrac beta 2,tfrac 12mu ,tfrac beta 2,tfrac 12(1+z)−bdisplaystyle (1+z)^-b(1+z)^-b4z(1+z)2displaystyle tfrac 4z(1+z)^2tfrac 4z(1+z)^2
|z|<1displaystyle |z|<1
α,α,μdisplaystyle alpha ,alpha ,mu alpha ,alpha ,mu α,μ2,12displaystyle alpha ,tfrac mu 2,tfrac 12alpha ,tfrac mu 2,tfrac 12(1−2z)−bdisplaystyle (1-2z)^-b(1-2z)^-b4z(z−1)(1−2z)2displaystyle tfrac 4z(z-1)(1-2z)^2tfrac 4z(z-1)(1-2z)^2
ℜz<12displaystyle Re z<tfrac 12Re z<tfrac 12

另外还有:


2Γ(12)Γ(a+b+12)Γ(a+12)Γ(b+12)F−12,β,μ2(z)=Fβ,β,μ(12−12z)+Fβ,β,μ(12+12z),|arg⁡z|<π,|arg⁡(1−z)|<π<pi tfrac 2Gamma (tfrac 12)Gamma (tfrac a+b+12)Gamma (tfrac a+12)Gamma (tfrac b+12)F_-tfrac 12,beta ,tfrac mu 2(z)=F_beta ,beta ,mu left(tfrac 12-tfrac 12sqrt zright)+F_beta ,beta ,mu left(tfrac 12+tfrac 12sqrt zright),quad |arg z|<pi ,|arg(1-z)|<pi

将它们与 Kummer 表组合起来,就得到所有的含有两个独立参变量的二次变换关系式。例如上面的例子可以通过组合第一行中的变换与 Pfaff 变换得到。


另外还有一些只含有一个独立参变量的二次变换关系式。



三次及高次变换


若一组李代数参数满足下列条件:有两个是 ±1/3,或者三个参数的绝对值相等,则有一个三次变换的公式将它与另一个超几何函数联系起来。


另外有一些 4 次和 6 次变换的公式。其它次数的变换公式只有当参数取特定有理数值时存在。参见Goursat (1881)。



特殊值



z=0


2F1(a,b;c;0)=1displaystyle _2F_1(a,b;c;0)=1_2F_1(a,b;c;0)=1


z=1


2F1(a,b;c;1)=B(a,c−a−b)B(a,c−a)=Γ(c)Γ(c−a−b)Γ(c−a)Γ(c−b),ℜ(c)>ℜ(a+b)displaystyle _2F_1(a,b;c;1)=tfrac mathrm B (a,c-a-b)mathrm B (a,c-a)=tfrac Gamma (c)Gamma (c-a-b)Gamma (c-a)Gamma (c-b),quad Re (c)>Re (a+b)_2F_1(a,b;c;1)=tfrac mathrmB (a,c-a-b)mathrmB (a,c-a)=tfrac Gamma (c)Gamma (c-a-b)Gamma (c-a)Gamma (c-b),quad Re (c)>Re (a+b)

这称为高斯原理(Gauss's theorem),可以由超几何函数的积分表示得到。范德蒙恒等式是它的特殊情形。



z=-1


2F1(a,b;1+a−b;−1)=Γ(1+a−b)Γ(1+12a)Γ(1+a)Γ(1+12a−b)displaystyle _2F_1(a,b;1+a-b;-1)=frac Gamma (1+a-b)Gamma (1+tfrac 12a)Gamma (1+a)Gamma (1+tfrac 12a-b)_2F_1(a,b;1+a-b;-1)=frac Gamma (1+a-b)Gamma (1+tfrac 12a)Gamma (1+a)Gamma (1+tfrac 12a-b)

这可以通过组合上表中的第二个二次变换和 Pfaff 变换,并利用 z=1 时的特殊值得到。



z=1/2


2F1(a,b;12(1+a+b);12)=Γ(12)Γ(12(1+a+b))Γ(12(1+a))Γ(12(1+b)).displaystyle _2F_1left(a,b;tfrac 12left(1+a+bright);tfrac 12right)=frac Gamma (tfrac 12)Gamma (tfrac 12left(1+a+bright))Gamma (tfrac 12left(1+a)right)Gamma (tfrac 12left(1+bright))._2F_1left(a,b;tfrac 12left(1+a+bright);tfrac 12right)=frac Gamma (tfrac 12)Gamma (tfrac 12left(1+a+bright))Gamma (tfrac 12left(1+a)right)Gamma (tfrac 12left(1+bright)).
2F1(a,1−a;c;12)=Γ(12c)Γ(12(1+c))Γ(12(c+a))Γ(12(1+c−a)).displaystyle _2F_1left(a,1-a;c;tfrac 12right)=frac Gamma (tfrac 12c)Gamma (tfrac 12left(1+cright))Gamma (tfrac 12left(c+aright))Gamma (tfrac 12left(1+c-aright))._2F_1left(a,1-a;c;tfrac 12right)=frac Gamma (tfrac 12c)Gamma (tfrac 12left(1+cright))Gamma (tfrac 12left(c+aright))Gamma (tfrac 12left(1+c-aright)).

上面两式分别被称为高斯第二求和原理与 Balley 原理。它们都可以通过组合第三个二次变换和 Pfaff 变换,并利用 z=1 时的特殊值得到。



参考文献



  • Hazewinkel, Michiel (编), Hypergeometric function, 数学百科全书, Springer, 2001, ISBN 978-1-55608-010-4 

  • John Pearson, Computation of Hypergeometric Functions (University of Oxford, MSc Thesis)

  • Marko Petkovsek, Herbert Wilf and Doron Zeilberger, The book "A = B" (freely downloadable)

  • 埃里克·韦斯坦因. Hypergeometric Function. MathWorld. 


  • Olde Daalhuis, A. B., Hypergeometric Function, (编) Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W., NIST Handbook of Mathematical Functions, Cambridge University Press, 2010, ISBN 978-0521192255, MR2723248 


  • Goursat, Édouard. Sur l'équation différentielle linéaire, qui admet pour intégrale la série hypergéométrique. Annales Scientifiques de l'École Normale Supérieure. 1881, 10: 3–142 [2008-10-16] (法语). 


  • ^ Dereziński, Jan. Hypergeometric type functions and their symmetries. arXiv:1305.3113. 








  • Popular posts from this blog

    【情報】本週珍珠商品重點:煉金時裝 + 艾港勞工宿舍!!

    京昆高速公路

    【攻略】陳戈-謝勒汗智慧的古書 (完成)