超几何函数

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA Template:Otheruses 在数学中,高斯超几何函数普通超几何函数2F1(a,b;c;z)是一个用超几何级数定义的函数,很多特殊函数都是它的特例或极限。所有具有三个Template:Link-en的二阶线性常微分方程的解都可以用超几何函数表示。

超几何级数

c为正整数时,对于|z| < 1,超几何函数可用如下幂级数定义

2F1(a,b;c;z)=n=0a(n)b(n)c(n)znn!

其中  x(n)遞進階乘,定义为:

q(n)={1if n=0q(q+1)(q+n1)if n>0

ab0或负整数时级数只有有限项,另有避免这种情况出现的正则超几何函数。

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

2F1(a,b;c;z)=Γ(ba)Γ(c)(z)aΓ(b)Γ(ca)k=0(a)k(ac+1)kzkk!(ab+1)k+Γ(ab)Γ(c)(z)bΓ(a)Γ(cb)k=0(b)k(bc+1)kzkk!(a+b+1)k/;|z|1ab

特殊情形

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

ln(1+z)=z2F1(1,1;2;z).
(1z)a=2F1(a,1;1;z)
arcsinz=z2F1(12,12;32;z2)

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

M(a,c,z)=limb2F1(a,b;c;b1z)

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

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

2F1(a,1a;c;z)=Γ(c)z1c2(1z)c12Pa1c(12z)

很多多项式,例如贾可比多项式 PTemplate:Su及其特殊情形勒让德多项式车比雪夫多项式Gegenbauer多项式都能用超几何函数表示

2F1(n,α+1+β+n;α+1;x)=n!(α+1)nPn(α,β)(12x) 其它特殊情形还包括Krawtchouk多项式Meixner多项式Meixner–Pollaczek多项式

Template:Link-en有时能表示成参数a, b, c是1, 1/2, 1/3, ... 或 0 的超几何函数之比的反函数。例如,若

τ=i2F1(12,12;1;1z)2F1(12,12;1;z)

z=κ2(τ)=θ2(τ)4θ3(τ)4

τ的椭圆模函数.

不完全 Beta 函数 Bx(p,q) 表示成

Bx(p,q)=xpp2F1(p,1q;p+1;x)

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

K(k)=π22F1(12,12;1;k2)
E(k)=π22F1(12,12;1;k2)

超几何方程

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

z(zddz+a)(zddz+b)w=zddz(zddz+c1)w,w(z)=2F1(a,b;c;z).

展开后,得

z(1z)d2wdz2+[c(a+b+1)z]dwdzabw=0.

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

正则奇点 0 附近的解

超几何方程的Template:Link-en

ρ(ρ1)+cρ=0

它的两个指标 Template:Mvar 是 0 和 1-Template:Mvar

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

2F1(a,b;c;z) and z1c2F1(1+ac,1+bc;2c;z)

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

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

2F1(a,b;c;z) and G2,22,0(1a,1b;0,c1;z), if c+
z1c2F1(1+ac,1+bc;2c;z) and G2,22,0(1a,1b;0,1c;z), if c0

正则奇点 1 附近的解

只需作代换 Template:Mvar=1-Template:Mvar,方程变为:

t(1t)d2wdt2+[1+a+bc(a+b+1)t]dwdtabw=0.

Template:Mvar 不是整数时,两个线性无关的正则特解为:

2F1(a,b;1+a+bc;1z) and (1z)cab2F1(cb,ca;1ab+c;1z)

正则奇点 ∞ 附近的解

Template:Mvar 不是整数时,两个线性无关的正则特解为:

za2F1(a,1+ac;1+ab;z1) and zb2F1(b,1+bc;1+ba;z1).

李代数参数与连接关系

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

Fα,β,μ(z)=2F1(a,b;c;z)
α=c1,β=a+bc,μ=ba
a=1+α+βμ2,b=1+α+β+μ2,c=1+α

参数 Template:Mvar 称为李代数参数。

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

At 0:Fα,β,μ(z) and zαFα,β,μ(z)
At 1:Fβ,α,μ(1z) and (1z)βFβ,α,μ(1z)
At :(z)1αβ+μ2Fμ,β,α(z1) and (z)1αβμ2Fμ,β,α(z1)

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

引入记号:

G(m;n,p)=πsinmπΓ(n)Γ(p)=G(m;p,n)
𝐅α,β,μ(z)=1Γ(1+α)Fα,β,μ(z)

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

𝐅β,α,μ(1z)=G(α;aα,bα)𝐅α,β,μ(z)+G(α;a,b)zα𝐅α,β,μ(z),
(1z)β𝐅β,α,μ(1z)=G(α;1a,1b)𝐅α,β,μ(z)+G(α;bβ,aβ)zα𝐅α,β,μ(z)=G(α;1a,1b)𝐅α,β,μ(z)+G(α;1(aα),1(bα))zα𝐅α,β,μ(z);
(z)a𝐅μ,β,α(z1)=G(α;1b,aα)𝐅α,β,μ(z)+G(α;a,aβ)zα𝐅α,β,μ(z),
(z)b𝐅μ,β,α(z1)=G(α;1a,bα)𝐅α,β,μ(z)+G(α;b,bβ)zα𝐅α,β,μ(z)=G(α;1a,1(aβ))𝐅α,β,μ(z)+G(α;b,1(aα))zα𝐅α,β,μ(z).

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

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

积分表示

B(a,ca)2F1(a,b;c;z)=1tbc(t1)ca1(tz)bdt,(c)>(a)>0,|arg(1z)|<π

式中的 Β 是beta函数

证明

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

p(a,b,c;t,z)=tbc(t1)ca1(tz)b2,w(a,b,c;t,z)=(tz)2p(a,b,c;t,z);

wz=b(tz)p(a,b,c;t,z),2wz2=b(b+1)p(a,b,c;t,z)
z(1z)2wz2+[c(a+b+1)z]wzabw=bp(a,b,c;t,z){z(1z)(b+1)+[c(a+b+1)z](tz)a(tz)2}=bp(a,b,c;t,z){at2+[c(ba+1)z]t+(bc+1)z}=bp(a,b,c;t,z){(bc+1)(t1)(tz)+(ca)t(tz)+(b1)t(t1)}=bt[t(t1)(tz)p(a,b,c;t,z)],

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

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

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

变换公式

分式线性变换

Pfaff 变换

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

2F1(a,b;c;z)=(1z)b2F1(ca,b;c;zz1),|arg(1z)|<π

Template:Mvar 的对称性自然有:

2F1(a,b;c;z)=(1z)a2F1(a,cb;c;zz1),|arg(1z)|<π
证明

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

u=zz1=1+1z1

z=uu1,(1z)a=(1u)a,dudz=(1u)2,d2udz2=(1u)3
z(1z)d2dz2[(1z)bw]+[c(a+b+1)z]ddz[(1z)bw]ab(1z)bw=(1z)b1{z[b(b+1)+2b(1z)ddz+(1z)2d2dz2]+[c(a+b+1)z][b+(1z)ddz]ab(1z)}w=(1z)b1{z(1z)2d2dz2+(1z)[c(ab+1)z]ddz+b(ca)}w=(1u)b+1{u(1u)d2du2+2uddu(1u)[c+(ab+1)(1u)1u]ddu+b(ca)}w=(1u)b+1{u(1u)d2du2+[c(ca+b+1)u]ddub(ca)}w

w=2F1(ca,b;c;u)

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

Euler 变换

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

2F1(a,b;c;z)=(1z)b2F1(ca,b;c;zz1)=(1z)b(1zz1)ac2F1(ca,cb;c;zz1zz11)=(1z)cab2F1(ca,cb;c;z),|arg(1z)|<π

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

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

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

例如 Euler 变换可以表示为:

Fα,β,μPfaffFα,μ,βFα,μ,βPfaffFα,β,μ

二次变换

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

2F1(a,b;2a;z)=(1z)b22F1(ab2,b2;a+12;z24z4),|arg(1z)|<π

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

证明

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

z(1z)d2dz2[(1z)b2w]+[c(a+b+1)z]ddz[(1z)b2w]ab(1z)b2w=(1z)b21{z[b2(b2+1)+b(1z)ddz+(1z)2d2dz2]+[c(a+b+1)z][b2+(1z)ddz]ab(1z)}w=(1z)b21{z(1z)2d2dz2+(1z)[c(a+1)z]ddz+b4[2(c2a)+(2ab)z]}w

c=2a,u=z24z4=14(z+111z)

1u=(z2)24(1z),dudz=z(z2)4(1z)2,d2udz2=12(1z)3
z(1z)2d2dz2+(1z)[c(a+1)z]ddz+b4[2(c2a)+(2ab)z]=z(1z)2d2dz2+(1z)[2a(a+1)z]ddz+b2(ab2)z=z3(z2)216(1z)2d2du2z2(1z)ddu+z(z2)(2aazz)4(1z)ddu+b2(ab2)z=z{u(1u)d2du2+[a+12(a+1)u]ddub2(ab2)}

w=2F1(ab2,b2;a+12;u)

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

其它例子

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

Fα,β,μ(z)=f(z)Fα,β,μ(g(z)),P(z)

其中 Template:Mvar(Template:Mvar),Template:Mvar(Template:Mvar) 是 Template:Mvar 的函数, Template:Mvar(Template:Mvar) 表示 Template:Mvar 要满足的约束。

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

李代数参数(左) 李代数参数(右) f(z) g(z) P(z)
α,μ,μ α2,μ,12 (112z)b (z2z)2 |arg(1z)|<π
μ,β,μ μ,β2,12 (1+z)b 4z(1+z)2 |z|<1
α,α,μ α,μ2,12 (12z)b 4z(z1)(12z)2 z<12

另外还有:

2Γ(12)Γ(a+b+12)Γ(a+12)Γ(b+12)F12,β,μ2(z)=Fβ,β,μ(1212z)+Fβ,β,μ(12+12z),|argz|<π,|arg(1z)|<π

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

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

三次及高次变换

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

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

特殊值

2F1(a,b;c;0)=1
2F1(a,b;c;1)=B(a,cab)B(a,ca)=Γ(c)Γ(cab)Γ(ca)Γ(cb),(c)>(a+b)

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

2F1(a,b;1+ab;1)=Γ(1+ab)Γ(1+12a)Γ(1+a)Γ(1+12ab)

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

2F1(a,b;12(1+a+b);12)=Γ(12)Γ(12(1+a+b))Γ(12(1+a))Γ(12(1+b)).
2F1(a,1a;c;12)=Γ(12c)Γ(12(1+c))Γ(12(c+a))Γ(12(1+ca)).

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

参考文献