马丢函数

来自testwiki
imported>IHEPer2024年5月23日 (四) 10:18的版本 傅立叶展开式
(差异) ←上一版本 | 最后版本 (差异) | 下一版本→ (差异)
跳转到导航 跳转到搜索
MathieuCE 3D
MathieuSE 3D

马丟函数Template:Lang-fr)是1868年法國數學家Template:Tsl因研究数学物理所推得的特殊函數,下列马丟方程的解析解:

d2ydx2+[a2qcos(2x)]y=0.

马丟方程有两个线性无关的解:

奇数解

MathieuCE(n, q, x),或记为wI(n,q,x),

偶数解

MathieuSE(n, q, x).或记为wII(n,q,x) 称为基本解[1]

周期性

马丟函数 MathieuC(a,q,z) 或 MathieuS(a,q,z) 只有一个是周期为 π2π的周期解,另一个不是。

马丟函数 MathieuC(a,q,z) 和 MathieuS(a,q,z) 两者都有是周期为2nπ(n≥2)的周期函数。 [1]

正交性

  • 02*πcem(x,q)*cen(x,q)dx=0
  • 02*πcem(x,q)*sen(x,q)dx=0
  • 02*πsem(x,q)*sen(x,q)dx=0

特征方程

Mathieu Eigen value a(n,q)
Mathieu eigenvalue b(n,q)

马丟方程的特征方程是[1]

cos(π*v)=wI(a,q,π)

cos(π*v)=wII(b,q,π)

对于给定的v,q, 上列特征方程给出无穷多个a、b解称为特征值。

特征值的展开

马丟函数体特征值可展开成级数:[2]

a0(q)=(1/2)*z2+(7/128)*z4(29/2304)*z6+(68687/18874368)*z8+𝒪(z10) a1(q)=1+z(1/8)*z2(1/64)*z3(1/1536)*z4+(11/36864)*z5+(49/589824)*z6+(55/9437184)*z7(83/35389440)*z8(12121/15099494400)*z9+𝒪(z10) a2(q)=4+(5/12)*z2(763/13824)*z4+(1002401/79626240)*z6(1669068401/458647142400)*z8+𝒪(z10) a3(q)=9+(1/16)*z2+(1/64)*z3+(13/20480)*z4(5/16384)*z5(1961/23592960)*z6(609/104857600)*z7+(4957199/2113929216000)*z8+(872713/1087163596800)*z9+𝒪(z10)

b1(q)=1z(1/8)*z2+(1/64)*z3(1/1536)*z4(11/36864)*z5+(49/589824)*z6(55/9437184)*z7(83/35389440)*z8+(12121/15099494400)*z9+𝒪(z10) b2(q)=4(1/12)*z2+(5/13824)*z4(289/79626240)*z6+(21391/458647142400)*z8+𝒪(z10) b3(q)=9+(1/16)*z2(1/64)*z3+(13/20480)*z4+(5/16384)*z5(1961/23592960)*z6+(609/104857600)*z7+(4957199/2113929216000)*z8(872713/1087163596800)*z9+𝒪(z10) b4(q)=16+(1/30)*z2(317/864000)*z4+(10049/2721600000)*z6(93824197/2006581248000000)*z8+𝒪(z10) b5(q)=25+(1/48)*z2+(11/774144)*z4(1/147456)*z5+(37/891813888)*z6+(7/339738624)*z7+(63439/201364441399296)*z8+(1/2130840649728)*z9+𝒪(z10)

级数展开

马丟函数ce,se的级数展开[3]

ce0(z,q)=1(1/2)*cos(2*z)*q+(1/16+(1/32)*cos(4*z))*q2+((11/128)*cos(2*z)(1/1152)*cos(6*z))*q3+O(q4) ce1(z,q)=cos(z)(1/8)*cos(3*z)*q+((1/128)*cos(z)(1/64)*cos(3*z)+(1/192)*cos(5*z))*q2+((1/512)*cos(z)+(1/3072)*cos(3*z)+(1/1152)*cos(5*z)(1/9216)*cos(7*z))*q3+O(q4) ce2(z,q)=cos(2*z)+(1/4(1/12)*cos(4*z))*q+((19/288)*cos(2*z)+(1/384)*cos(6*z))*q2+(49/1152+(11/4608)*cos(4*z)(1/23040)*cos(8*z))*q3+O(q4) ce3(z,q)=cos(3*z)+((1/8)*cos(z)(1/16)*cos(5*z))*q+((5/512)*cos(3*z)+(1/64)*cos(z)+(1/640)*cos(7*z))*q2+((1/512)*cos(3*z)(1/4096)*cos(z)+(11/40960)*cos(5*z)(1/46080)*cos(9*z))*q3+O(q4) ce4(z,q)=cos(4*z)+((1/12)*cos(2*z)(1/20)*cos(6*z))*q+((17/3600)*cos(4*z)+1/192+(1/960)*cos(8*z))*q2+((7/28800)*cos(2*z)+(29/288000)*cos(6*z)(1/80640)*cos(10*z))*q3+O(q4)

se1(z,q)=sin(z)(1/8)*sin(3*z)*q+((1/128)*sin(z)+(1/64)*sin(3*z)+(1/192)*sin(5*z))*q2+((1/512)*sin(z)+(1/3072)*sin(3*z)(1/1152)*sin(5*z)(1/9216)*sin(7*z))*q3+O(q4) se2(z,q)=sin(2*z)(1/12)*sin(4*z)*q+((1/288)*sin(2*z)+(1/384)*sin(6*z))*q2+((1/1536)*sin(4*z)(1/23040)*sin(8*z))*q3+O(q4) se3(z,q)=sin(3*z)+((1/8)*sin(z)(1/16)*sin(5*z))*q+((5/512)*sin(3*z)(1/64)*sin(z)+(1/640)*sin(7*z))*q2+((1/512)*sin(3*z)(1/4096)*sin(z)+(11/40960)*sin(5*z)(1/46080)*sin(9*z))*q3+O(q4) se4(z,q)=sin(4*z)+((1/12)*sin(2*z)(1/20)*sin(6*z))*q+((17/3600)*sin(4*z)+(1/960)*sin(8*z))*q2+((1/1600)*sin(2*z)+(29/288000)*sin(6*z)(1/80640)*sin(10*z))*q3+O(q4) se5(z,q)=sin(5*z)+((1/16)*sin(3*z)(1/24)*sin(7*z))*q+((13/4608)*sin(5*z)+(1/384)*sin(z)+(1/1344)*sin(9*z))*q2+((7/73728)*sin(3*z)+(13/258048)*sin(7*z)(1/9216)*sin(z)(1/129024)*sin(11*z))*q3+O(q4)

傅立叶展开式

马丟函数的傅立叶展开:[3]

  • MathieuCE(2n,q,x)=m=0A2m2n(q)cos(2mx)
  • MathieuCE(2n+1,q,x)=m=0A2m+12n+1(q)cos[(2m+1)x]
  • MathieuSE(2n+1,q,x)=m=0B2m+12n+1(q)sin[(2m+1)x]
  • MathieuSE(2n+2,q,x)=m=0B2m+22n+2(q)sin[(2m+2)x]

其中系数A,B满足下列递归关系:[3]


aA0=qA2

(a4)A2=q(2A0+A4)

(a4m2)A2m=q(A2m2+A2m+2)

(a1+q)B1=qB3

(a(2m+1)2)B2m+1=q(B2m1+B2m+1)

关系式

马丟方程的基本解WIWII满足下列关系:[3]:

|wI(n,q,0)wII(n,q,0)wi'(n,q,0)wII'(n,q,0)|= |1001|

郎斯基行列式: W[wI,wII]=1

wI(a,q,z+π)=wI(a,q,π)*wI(a,q,z)+wI'(a,q,π)*wII(a,q,z) wI(a,q,zπ)=wI(a,q,π)*wI(a,q,z)wI'(a,q,π)*wII(a,q,z) wII(a,q,z+π)=wII(a,q,π)*wII(a,q,z)+wII'(a,q,π)*wII(a,q,z) wII(a,q,zπ)=wII(a,q,π)*wII(a,q,z)wI'(a,q,π)*wII(a,q,z)


wI(z)=wI(z) wII(z)=wII(z)

特例

  • CE(a,0,z)=cos(az)
  • SE(a,0,z)=sin(az)
  • MathieuA(1,0)=1
  • MathieuA(a,0)=a2
  • MathieuB(a,0)=a2
  • MathieuFloquet(a,0,z)=exp(I*sqrt(a)*z)

夫洛开解

Mathieu Floquet

马丟函数中,如果f(x) 是一个周期为ω的解,并满足下列条件

f(x+ω)=σ*f(x),其中σ与x 无关,则此解称为夫洛开解。

级数展开

MF(1,1,z)=.7992.5734*I+(.9134+.6553*I)*z+(.3996.2867*I)*z2+(.1523+.1092*I)*z3+(.2331+.1673*I)*z4+O(z5) MF(1,2,z)=.7643.4526*I+(1.167+.6910*I)*z+(1.146.6789*I)*z2+(.5835+.3455*I)*z3+(.2229+.1320*I)*z4+O(z5) MF(1,3,z)=.6841.3703*I+(1.318+.7135*I)*z+(1.710.9258*I)*z2+(1.098+.5946*I)*z3+(0.2851e10.1543e1*I)*z4+O(z5)

参考文献

  1. 1.0 1.1 1.2 王竹溪 郭敦仁 603 引用错误:<ref>标签无效;同一name(名称)“W”以不同内容定义了多次
  2. Frank p659
  3. 3.0 3.1 3.2 3.3 Frank p660 引用错误:<ref>标签无效;同一name(名称)“F”以不同内容定义了多次
  • 王竹溪 郭敦仁 《特殊函数概论》 第十二章 马丟函数 北京大学出版社 2000
  • Frank J Oliver NIST Handbook of Mathematical Functions,Cambridge University PRESS, 2010