椭圆积分

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

积分学中,椭圆积分最初出现于椭圆弧长有关的问题中。Template:Le欧拉是最早的研究者。现代数学将椭圆积分定义为可以表达为如下形式的任何函数 f的积分

f(x)=cxR[t,P(t)] dt

其中R是其两个参数的有理函数P是一个无重根的34多项式,而c是一个常数。

通常,椭圆积分不能用基本函数表达。这个一般规则的例外出现在P有重根的时候,或者是R,(x,y)没有y的奇数幂时。但是,通过适当的简化公式,每个椭圆积分可以变为只涉及有理函数和三个经典形式的积分。(也即,第一,第二,和第三类的椭圆积分)。

除下面给出的形式之外,椭圆积分也可以表达为勒让德形式Carlson对称形式。通过对施瓦茨-克里斯托费尔映射的研究可以加深对椭圆积分理论的理解。历史上,椭圆函数是作为椭圆积分的逆函数被发现的,特别是这一个:F[sn(z;k);k]=z其中sn雅可比正弦椭圆函数

记法

椭圆积分通常表述为不同变量的函数。这些变量完全等价(它们给出同样的椭圆积分),但是它们看起来很不相同。很多文献使用单一一种标准命名规则。在定义积分之前,先来检视一下这些变量的命名常规:

上述三种常规完全互相确定。规定其中一个和规定另外一个一样。椭圆积分也依赖于另一个变量,可以有如下几种不同的设定方法:

规定其中一个决定另外两个。这样,它们可以互换地使用。注意u也依赖于m。其它包含u的关系有

cosϕ=cnu

1msin2ϕ=dnu.

后者有时称为δ幅度并写作Δ(ϕ)=dnu。有时文献也称之为补参数,补模或者补模角。这些在四分周期中有进一步的定义。

第一类不完全椭圆积分

第一类不完全椭圆积分 F定义为

F(ϕα)=F(ϕ|m)=0ϕdθ1(sinθsinα)2.

与此等价,用雅可比的形式,可以设 x=sinϕ,t=sinθ;则

F(ϕα)=F(x;k)=0xdt(1t2)(1k2t2)

其中,假定任何有竖直条出现的地方,紧跟竖直条的变量是(如上定义的)参数;而且,当反斜杠出现的时候,跟着出现的是模角。 在这个意义下,F(sinϕ;sinα)=F(ϕ|sin2α)=F(ϕα),这里的记法来自标准参考书Abramowitz and Stegun

但是,还有许多不同的用于椭圆积分的记法。取值为椭圆积分的函数没有(像平方根正弦误差函数那样的)标准和唯一的名字。甚至关于该领域的文献也常常采用不同的记法。Gradstein, Ryzhik[1] Template:Wayback, Eq.(8.111)]采用F(ϕ,k)。该记法和这里的F(ϕ|k2);以及下面的E(ϕ,k)=E(ϕ|k2)等价。

和上面的不同对应的是,如果从Mathematica语言翻译代码到Maple语言,必须将EllipticK函数的参数用它的平方根代替。反过来,如果从Maple翻到Mathematica,则参数应该用它的平方代替。Maple中的EllipticK(x)几乎和Mathematica中的EllipticK[x2]相等;至少当0<x<1时是相等的。

注意

F(x;k)=u

其中u如上文所定义:由此可见,雅可比椭圆函数是椭圆积分的逆。

加法公式

φ1,φ2]π2;π2[,
F(φ1,k)+F(φ2,k)=F(arctan(tanφ11k2sin2φ2)+arctan(tanφ21k2sin2φ1),k)
arctan(tanφ11k2sin2φ2)+arctan(tanφ21k2sin2φ1)[π/2;π/2]
F(φ1,k)+F(φ2,k)=F(arcsincosφ11k2sin2φ1sinφ2+cosφ21k2sin2φ2sinφ11k2sin2φ1sin2φ2,k)
arctan(tanφ11k2sin2φ2)+arctan(tanφ21k2sin2φ1)[0;π]
F(φ1,k)+F(φ2,k)=F(arccoscosφ1cosφ2sinφ1sinφ21k2sin2φ11k2sin2φ11k2sin2φ21k2sin2φ1sin2φ2,k)

性质

F(x+nπ;k)=F(x;k)+2nK(k)
F(x+nπ2;k)=nK(k)
n
F(x;k)=F(x;k)
F(x;0)=x
F(0;k)=F(x;k)
F(x;1)=arctanhsinx
π2<(x)<π2

第一类不完全椭圆积分的导数

ddxF(x;k)=11k2sin2x
ddkF(x;k)=E(x;k)2k(1k)F(x;k)2ksin2x4(1k)1ksin2x

第二类不完全椭圆积分

第二类不完全椭圆积分 E

E(ϕα)=E(ϕ|m)=0ϕE(θ) dθ=0ϕ1(sinθsinα)2 dθ.

与此等价,采用另外一个记法(作变量替换t=sinθ),

E(x;k)=0x1k2t21t2 dt.

其它关系包括

E(ϕ|m)=0udn2wdw=um0usn2wdw=(1m)u+m0ucn2wdw.
E(ϕ|k2)=(1k2)0ϕdθ(1k2sin2θ)1k2sin2θ+k2sinθcosθ1k2sin2θ

加法公式

φ1,φ2]π2;π2[,
E(φ1,k)+E(φ2,k)=[E(arctan(tanφ11k2sin2φ2)+arctan(tanφ21k2sin2φ1),k)+k2sinφ1sinφ2(cosφ11k2sin2φ1sinφ2+cosφ21k2sin2φ2sinφ1)1k2sinφ1sinφ2]

性质

E(ϕ+nπ;k)=E(ϕ;k)+2nE(k)
E(ϕ;k)=E(ϕ;k)

第二类不完全椭圆积分的导数

ddϕE(ϕ;k)=1k2sin2ϕ
ddkE(ϕ;k)=E(ϕ;k)F(ϕ;k)2k
dndknE(ϕ;k)=π2kn2F1(12,12;1n;k)πcosϕ2k2nF2×1×01×3×2[12;12,12,1;12,1;1,32;1n;;k2cosϕ,cos2ϕ]+πm1ncosϕ8F3×1×12×1×1[12,32,2;12,1;2,2n;1n;32;32;k2cos2ϕ,k2]

第三类不完全椭圆积分

第三类不完全椭圆积分Π

Π(n;ϕ|m)=0ϕdθ(1nsin2θ)1(sinθsinoε)2,

或者

Π(n;ϕ|m)=0sinϕdt(1nt2)(1k2t2)(1t2),

或者

Π(n;ϕ|m)=0F(ϕ|m)dw1nsn2(w|m).

数字n称为特征数,可以取任意值,和其它参数独立。但是要注意Π(1;π2|m)对于任意m是无穷的。

加法公式

Π(n;ϕ1,k)+Π(n;ϕ2,k)=Π[n;arccoscosϕ1cosϕ2sinϕ1sinϕ2(1k2sin2ϕ1)(1k2sin2ϕ2)1k2sin2ϕ1sin2ϕ2,k]n(1n)(nk2)arctan(1n)n(nk2)sinarccoscosϕ1cosϕ2sinϕ1sinϕ2(1k2sin2ϕ1)(1k2sin2ϕ2)1k2sin2ϕ1sin2ϕ2sinϕ1sinϕ2ncosϕ1cosϕ2nsinϕ1sinϕ2(1k2sin2ϕ1)(1k2sin2ϕ2)1k2sin2ϕ1sin2ϕ21k2sin2arccoscosϕ1cosϕ2sinϕ1sinϕ2(1k2sin2ϕ1)(1k2sin2ϕ2)1k2sin2ϕ1sin2ϕ2sinϕ1sinϕ2+1nsin2arccoscosϕ1cosϕ2sinϕ1sinϕ2(1k2sin2ϕ1)(1k2sin2ϕ2)1k2sin2ϕ1sin2ϕ2

第三类不完全椭圆积分的导数

nΠ(n;ϕ,k)=12(k2n)(n1)[E(ϕ;k)+(k2n)F(ϕ;k)n+(n2k2)Π(n;ϕ,k)nn1k2sinϕsin2ϕ2(1nsin2ϕ)]
mnmΠ(n;ϕ,k)=sinϕnmq=0q!(nsin2ϕ)q(2q+1)Γ(qm+1)F1(q+12,12,12;q+32;sin2ϕ,k2sin2ϕ)
ϕΠ(n;ϕ,k)=1(1k2sin2ϕ)
kΠ(n;ϕ,k)=knk2[E(ϕ;k)k21+Π(n;ϕ,k)k2sin2ϕ2(k21)1k2sin2ϕ]

特殊值

Π(n;ϕ,1)=12n2[nln1+nsinϕ1nsinϕ2ln(secϕ+tanϕ)]
π2(ϕ)π2
Π(0;ϕ,k)=F(ϕ,k)
Π(n;ϕ,0)=arctanh(n1tanϕ)n1
π2(ϕ)π2
Π(n;ϕ,n)=11n[E(ϕ,n)nsin2ϕ21nsin2ϕ]
Π(n;1k,k)=1kΠ(nk2,1k)
Π(1;ϕ,k)=1k2sin2ϕtanϕE(ϕ,k)1k2+F(ϕ,k)

第一类完全椭圆积分

第一类完全椭圆积分K(k)

如果幅度为π2或者x=1,则称椭圆积分为完全的。 第一类完全椭圆积分K可以定義为

K(k)=0π2dθ1k2sin2θ

或者

K(k)=01dt(1t2)(1k2t2).

它是第一类不完全椭圆积分的特例:

K(k)=F(1;k)=F(π2|k2)

这个特例可以表达为幂级数

K(k)=π2n=0[(2n)!22nn!2]2k2n

它等价于

K(k)=π2{1+(12)2k2+(1324)2k4++[(2n1)!!(2n)!!]2k2n+}.

其中n!!表示双阶乘。利用高斯的超几何函数,第一类完全椭圆积分可以表达为

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

第一类完全椭圆积分有时称为四分周期。它可以利用算术几何平均值來快速计算。

K(k)=π2agm(1,1k2).

复数值

[K(x+yi)]=π2F2×1×14×0×0[34,34,54,54,;;;1,32;12;32;y2,x2]+π8xF2×1×14×0×0[14,14,34,34,;;;1,12;12;12;y2,x2]
[K(x+yi)]=π8yF2×1×14×0×0[34,54,34,54,;;;1,32;32;12;y2,x2]+964πxyF2×1×14×0×0[54,74,74,54,;;;2,32;32;32;y2,x2]


特殊值

K(±)=0
K(±i)=0
K(0)=π2
K(1)=
K(22)=8πΓ2(14)π
K(17122)=(4+22)πΓ2(14)π
K(624)=43348πΓ3(13)
K(6+24)=432748πΓ3(13)
K(i)=2π8πΓ2(14)
K(2)=42ππΓ2(14)+42ππΓ2(14)i
K(ik)=1k2+1K(k2k2+1)

其中

Γ(14)3.62561
Γ(13)2.67893

第一类完全椭圆积分满足

E(k)K(k)+E(k)K(k)K(k)K(k)=π2

导数

ddkKn(k)=nKn1(k)E(k)2k(1k)nKn(k)2k

漸近表示

K(k2)π2+π8k21k2π16k41k2

這個近似在k<1/2時相對誤差小於Template:Math,若只保留前兩項則誤差在k<1/2時小於0.01

微分方程

此函數滿足以下微分方程

ddk[k(1k2)dK(k)dk]=kK(k)

此微分方程之另一解為K(1k2),此解滿足以下關係。

ddkK(1k2)=E(k)k(1k2)K(k)k.

第二类完全椭圆积分

第二类完全椭圆积分E(k)

第二类完全椭圆积分 E可以定义为

E(k)=0π21k2sin2θ dθ

或者

E(k)=011k2t21t2 dt.

它是第二类不完全椭圆积分的特殊情况:

E(k)=E(1;k)=E(π2|k2)

它可以用幂级数表达

E(k)=π2n=0[(2n)!22nn!2]2k2n12n

也就是

E(k)=π2{1(12)2k21(1324)2k43[(2n1)!!(2n)!!]2k2n2n1}.

高斯超几何函数表示的话,第二类完全椭圆积分可以写作

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

有如下性质

E(nπ2;k)=nE(k)
n


复数值

E(x+yi)={π2F2×1×14×0×0[34,54,14,34,;;;1,32;12;32;y2,x2]π8xF2×1×14×0×0[14,34,14,14,;;;1,12;12;12;y2,x2]}+i{π8yF2×1×14×0×0[34,54,14,34,;;;1,32;12;32;y2,x2]364πxyF2×1×14×0×0[54,74,34,54,;;;2,32;32;32;y2,x2]}

特殊值

E(0)=π2
E(1)=1
E()=i
E()=
E(i)=(2222i)
E(i)=2π2πΓ2(34)+2ππ24πΓ2(34)=π2πΓ2(14)+2π8πΓ2(14)
E(i)=(22+22i)
E(22)=π32Γ(14)2+18πΓ(14)2
E(624)=23 343Γ3(13)π2+43(334+274)48πΓ3(13)
E(6+24)=23 2743Γ3(13)π2+43(27434)16πΓ3(13)
E(21)=π8[Γ(18)Γ(58)+Γ(58)Γ(98)]
E(2)=12πΓ2(34)+12πΓ2(34)i


其中

Γ(18)7.53394
Γ(58)1.43452
Γ(98)0.94174
Γ(34)1.22541

导数、積分及微分方程

ddkE(k)=E(k)K(k)k
E(k)dk=23[kK(k)K(k)+kE(k)+E(k)]
(k21)ddk[kdE(k)dk]=kE(k)

此微分方程之另解為E(1k2)K(1k2)

第三类完全椭圆积分

不同n值的第三类完全椭圆积分Π(n,k)

第三类完全椭圆积分Π可以定义为

Π(n,k)=0π2 dθ(1nsin2θ)1k2sin2θ

注意有时第三类椭圆积分被定义为带相反符号的n,也即

Π(n,k)=0π2 dθ(1+nsin2θ)1k'2sin2θ.

阿佩尔函数可表示为

Π(m,n)=π2F1(12;1,12;1;m,n)

第三类完全椭圆积分和第一类椭圆积分之间的关系

Π[(1+x)(13x)(1x)(1+3x),(1+x)3(13x)(1x)3(1+3x)]1+3x6xK[(1+x)3(13x)(1x)3(1+3x)]=

{0for 0<x<1π(x1)(x1)(1+3x)12xfor x<0,x>1

K(22)=π4πΓ2(14)=36392Π(12332,12)

=3+6392Π(1+2332,12)π2+3+7+3893

偏导数

nΠ(n,k)=12(k2n)(n1)[E(k)+(k2n)K(k)n+(n2k2)Π(n,k)n]
kΠ(n,k)=knk2[E(k)k21+Π(n,k)]

特殊值

Π(0,0)=π2
Π(n,0)=π21n
Π(n,1)=sgnn1
Π(n,n)=E(n)1n
Π(0,n)=K(n)
Π(±,n)=0
Π(n,±)=0

函數關係

勒讓得關係指出了第一类和第二类完全椭圆积分之间的联系:

K(k)E(1k2)+E(k)K(1k2)K(k)K(1k2)=π2.

参看

参考