拉梅函数

来自testwiki
跳转到导航 跳转到搜索
Lame function Maple animation plot

拉梅函数(Lame functions)是下列拉梅方程的解:[1][2][3]

雅可比形式

d2wdz2+(A+v(v+1)k2sn2(z,k))w=0+ 此拉梅方程的正则奇点在复数平面的2pK+(2q+1)iK 其中 p,q ∈Z,K代表模数为k的完全椭圆积分,K'代表模数为k=1k2的完全椭圆积分。

其中 k,v 都是实数,并且 0<k<1,

代数形式

作雅可比橢圓函數变数替换s=sn2(z,k)得拉梅方程的代数形式:

d2Λds2+12(12+1s1+1sh)dΛdsn(n+1)s+H4s(s1)(sh)Λ=0

h=k2=a2c2a2b2,

H=hA

h>1 此傅克型方程有四个正则奇点0,1,h,

魏尔斯特拉斯形式[3]

d2Λdz2+[Hn(n+1)(z)]Λ=0

其中魏尔斯特拉斯函数

三角函数形式

在雅可比形式的拉梅方程中做代换[4] snz=cosζ

ζ=12πamz

可得

[1(kcosζ)2]d2ΛdΛ2+k2cosζsinζdΛdζ+[hn(n+1)(kcosζ)2]Λ=0

在上列方程组 h,k,n等是实数或复数常数,而各变量为复数。

拉梅方程的本征值

对于给定的参数v,k,存在四套实数本征值h,令拉梅方程的奇数解或偶数解有2K或4K周期[5]

本征值 h 奇偶 周期
av2m(k2) 2K
av2m+1(k2) 4K
bv2m(k2) 4K
bv2m+1(k2) 2K

拉梅函数

与每一个本征值对应的本征函数,称为v阶拉梅函数,其记法及周期性列表于下:[6]

本征值 h 奇偶 周期 本征函数(拉梅函数)
av2m(k2) 2K Ecv2m(z,k2)
av2m+1(k2) 4K Ecv2m+1(z,k2)
bv2m(k2) 4K Esv2m+1(z,k2)
bv2m+1(k2) 2K Esv2m+2(z,k2)

其中2m,2m+1,2m+2代表在(0,2K)区间内的零点数。

拉梅函数是Heun函数的特例

Heun方程 gh:=d2(y(z)dz2+(γz+δz1+ϵza)d(y(z)dz+(αβzq)y(z)/(z(z1)(za))=0


令=γ=1/2,δ=1/2,ϵ=1/2,q=(1/4)ah,α=1/4,β=v(v+1)

则化为拉梅方程

d2(y(z)dz2+(1/2(1/z+1/(z1)+1/(za)))d(y(z)dz+(1/4)(ahν(ν+1)z)y(z)/(z(z1)(za))=0

拉梅方程的Heun函数解

由于拉梅方程式是Heun方程的特例,因此拉梅方程可以用HeunG函数表示[7] y(z)=C1HeunG(a,(1/4)ah,(1/2)ν,(1/2)ν+1/2,1/2,1/2,z)

+C2(z)HeunG(a,1/4+(1/4(h+1))a,1+(1/2)ν,1/2(1/2)ν,3/2,1/2,z) 其中二个HeunG函数是线性无关的。

拉梅函数的幂级数展开

拉梅函数可以展开成幂级数形式[8]


y(z)=v=0avzρ+v

其中ρ只能取0,1/2

例子

y(z)=1.2+2.3(z).600hz(.383(ah1.a1.))z(3/2)/a+(0.500e1(4.ah4.h+ah2+2.ν2+2.ν))z2/a+(0.192e1(10.a2h+9.a2+6.a10.ah+9.+a2h2+6.νa+6.ν2a))z(5/2)/a2+O(z3)

参考文献

  1. 王竹溪 第572页
  2. Whittaker p554
  3. 3.0 3.1 Erdelyi p55
  4. Erdelyi p 56
  5. Frank Oliver p685
  6. Frank, p684
  7. Frank Oliver,p713
  8. 王竹溪 第573页
  • 王竹溪 郭敦仁 《特殊函数概论》 北京大学出版 2000
  • Whittaker and Watson, A Course of Modern Analysis 1920, Cambridge University Press
  • Erdelyi, Higher Transcendental Functions Vol III