根据前面的知识,对拉普拉斯方程和亥姆霍兹方程进行分离变量,得到球函数方程:
sinθ1∂θ∂(sinθ∂θ∂Y)+sin2θ1∂φ2∂2Y+l(l+1)Y=0
球函数方程的解Y(θ,φ)称为球函数,即定义在半径r的球面上的函数,然后再分离变量,得到:
Y(θ,φ)=(Acosmφ+Bsinmφ)Θ(θ)
其中Θ(θ)需要从l阶连带的勒让德方程:
(1−x2)dx2d2Θ−2xdxdΘ+[l(l+1)−1−x2m2]Θ=0
中解出,式中采用变量代换:
x=cosθ
一、轴对称球函数#
1、勒让德多项式#
当m=0时,勒让德方程为:
(1−x2)dx2d2Θ−2xdxdΘ+l(l+1)Θ=0使用之前的方法,求出对应级数退化的多项式,将他们分别乘以适当的系数,使得最高次幂项xl的系数为:
al=2l(l!)2(2l)!反用系数递推公式,改写为:
ak=(k−l)(k+l+1)(k+2)(k+1)ak+2就可以把系数一一计算出来,得:
al−2n=(−1)nn!2l(l−n)!(l−2n)!(2l−2n)!将n记为k,求得l阶勒让德多项式的具体表达式:
Pl(x)=k=0∑[l/2](−1)k2lk!(l−k)!(l−2k)!(2l−2k)!xl−2k其中[l/2]表示不大于l/2的最大整数。由此写出前几个勒让德多项式:
P0(x)=1P1(x)=cosθP2(x)=21(3x2−1)=41(3cos2θ+1)现在计算Pl(0),这应当等于多项式Pl(x)的常数项,如l=2n+1,则P2n+1(x)只含奇次幂项,不含常数项,所以:
P2n+1(0)=0如果l=2n,则P2n(x)含有常数项,即:
P2n(0)=(−1)n2nn!2nn!(2n)!=(−1)n[(2n)!!]2(2n)!=(−1)n(2n)!!(2n−1)!!2、勒让德多项式的微分表示#
用二项式定理把(x2−1)l展开:
2ll!1(x2−1)l=2ll!1k=0∑l(l−k)!k!l!(x2)l−k(−1)k=k=0∑l(−1)k2lk!(l−k)!1x2l−2k把上式求导l次,凡是幂次2l−2k低于l的项在求导过程中都变成0,所以只需要保留幂次2l−2k≥的项,即k≤2l,这样:
2ll!1dxldl(x2−1)l=k=0∑[l/2](−1)k2lk!(l−k)!(2l−2k)(2l−2k−1)⋯(l−2k+1)xl−2k=k=0∑[l/2](−1)k2lk!(l−k)!(l−2k)!(2l−2k)!xl−2k=Pl(x)于是我们得到了勒让德多项式的微分形式:
Pl(x)=2ll!1dxldl(x2−1)l该公式叫做罗德里格斯公式。
3、勒让德多项式的积分形式#
利用高阶导数的柯西公式,微分表示可以改写为复变函数上的环路积分:
Pl(x)=2πi12l1∮C(z−x)l+1(z2−1)ldzC为z平面上围绕z=x点的任意闭合回路。该公式叫做施列夫利积分。该积分还可以进一步化为定积分形式,取C为圆周,圆心在z=x,半径为∣x2−1∣,在C上,z−x=x2−1eiψ,dz=ix2−1eiψdψ,代入到上式,得:
Pl(x)=2πi12l1∫−ππ(x2−1)l+1(eiψ)l+1[(x+x2−1eiψ)2−1]l⋅[ix2−1eiψ]dψ化简得到:
Pl(x)=π1∫0π[x+i1−x2cosψ]ldψ这叫做拉普拉斯积分,由此可以看出:
∣Pl(x)∣⩽π1∫0π∣cosθ+isinθcosψ∣ldψ=π1∫0π[cos2θ+sin2θcos2ψ]1/2dψ⩽π1∫0π[cos2θ+sin2θ]l/2dψ=π1∫0πdψ=1即:
∣Pl(x)∣≤1(−1≤x≤1)4、勒让德多项式的正交关系#
作为施图姆-刘维尔本征值问题的正交关系特例,不同阶的勒让德多项式在(−1,+1)上正交,即:
∫−1+1Pk(x)Pl(x)dx=0(l=k)如果从x回到原来的变数θ,则:
∫0πPk(cosθ)Pl(cosθ)sinθdθ(l=k)5、勒让德多项式的模#
如果前一部分中k=l,则可以计算其模Nl:
Nl2=∫−1+1[Pl(x)]2dx将上式的Pl(x)使用微分表示,并使用分部积分:
Nl2=22l(l!)21∫−1+1dxldl(x2−1)ldxd[dxl−1dl−1(x2−1)l]dx=22l(l!)21[dxldl(x2−1)ldxl−1dl−1(x2−1)l]−1+1−22l(l!)21∫−1+1dxl−1dl−1(x2−1)ldxd[dxldl(x2−1)l]dx这里的(x2−1)l是2l次多项式,它的2l阶导数也就是最高幂项x2l的二阶导数,即(2l)!,于是:
Nl2=(−1)l22l(l!)2(2l)!∫−1+1(x−1)l(x+1)ldx分部积分,得到:
Nl2=(−1)l22l(l!)2(2l)!⋅(−1)l+1l∫−11(x−1)l−1(x+1)l+1dx观察分部积分过程,每次分部积分都会使(x−1)的幂次降低一次,(x+1)的幂次升高一次,且积分乘上一个相应的常数因子,共操作l次,得到勒让德多项式的模:
Nl=2l+12,(l∈N)6、广义傅里叶级数#
根据上面的证明,可知勒让德多项式Pl(x)是完备的,可以作为广义傅里叶级数展开的基,定义在x的区间[−1,+1]上的函数f(x)或者定义在θ区间[0,π]上的函数f(θ)都可以展开为广义傅里叶级数。
⎩⎨⎧f(x)=l=0∑∞flPl(x)fl=22l+1∫−1+1f(x)Pl(x)dx或:
⎩⎨⎧f(θ)=l=0∑∞flPl(cosθ),fl=22l+1∫0πf(θ)Pl(cosθ)sinθdθ7、拉普拉斯方程的轴对称定解问题#
拉普拉斯方程的定解问题,如果有对称轴,自然就取这对称轴为球坐标的极轴,取这对称轴为球坐标的极轴,由此问题与φ无关,只需要用m=0的轴对称球函数。如果定解问题只在半球区域内有意义,则应当根据边界条件类型,选择使用奇延拓或者偶延拓,然后使用勒让德方法。
8、母函数#
在单位球北极放置一个4πε0单位的正电荷,则在球内任意一点M的静电势为:
d1=1−2rcosθ+r21静电势d1遵从拉普拉斯方程,且以球坐标系的极轴作为对称轴,因此应当满足拉普拉斯方程一般解的形式,即:
d1=l=0∑∞(Alrl+Blrl+11)Pl(cosθ)研究球心处的电势,应该有限,所以Bl=0,故:
1−2rcos2θ+r21=l=0∑∞AlrlPl(cosθ)可以利用广义傅里叶方法展开,更简单地,可以令θ=0:
1−r1=i=0∑∞Alrl上式左边使用泰勒展开,配凑系数相等得到Al=1,于是:
1−2rcos2θ+r21=l=0∑∞rlPl(cosθ),r<1同样的方法,计算出球外的静电势,结果为:
1−2rcosθ+r21=l=0∑∞rl+11Pl(cosθ),r>1于是,1−2rcosθ+r21被叫做勒让德多项式的母函数(或者是生成函数),用半径R的球代替单位球,则:
R2−2rRcosθ+r21=⎩⎨⎧l=0∑∞Rl+11rlPl(cosθ)l=0∑∞Rlrl+11Pl(cosθ)(r<R)(r>R)9、递推公式#
母函数可以推导出勒让德多项式的递推公式,将上式改写为:
1−2rx+r21=l=0∑∞rlPl(x)两边对r求导,并化简得:
1−2rx+r2x−r=(1−2rx+r2)l=0∑∞lrl−1Pl(x)将该写的未求导的公式代入到求导后的公式,得:
(x−r)l=0∑∞rlPl(x)=(1−2rx+r2)l=0∑∞lrl−1Pl(x)观察两边rk的系数,得到:
(k+1)Pk+1(x)−(2k+1)xPk(x)+kPk−1(x)=0(k≥1)利用这个式子可以从Pk−1(x),Pk(x)推导出Pk+1,这个称为勒让德多项式的递推公式。如果将改写的式子对x求导,采取类似的方法化简,比较系数后,可以得到另外一个递推公式:
Pk(x)=Pk+1′(x)−2xPk′(x)+Pk−1′(x)(k≥1)类似的递推公式还有:
Pk+1′(x)=(k+1)Pk(x)+xPk′(x)kPk(x)=xPk′(x)−Pk−1′(x)(k≥1)(x2−1)Pk′(x)=kxPk(x)−kPk−1(x)(k≥1)(2k+1)Pk(x)=Pk+1′(x)−Pk−1′(x)(k≥1)
二、连带勒让德函数#
1、连带勒让德函数#
l阶连带勒让德函数的表达式为:
(1−x2)dx2d2Θ−2xdxdΘ+[l(l+1)−1−x2m2]Θ=0(x=cosθ)x0=0是连带勒让德方程的常点,可以通过一般级数解法求得级数解,但是递推公式极为复杂,难以求出系数的一般形式。所以考虑其他方法,作变换:
Θ=(1−x2)2my(x)在这个这个变换下:
dxdΘ=(1−x2)2my′−m(1−x2)2m−xydx2d2Θ=(1−x2)2my′′−2m(1−x2)2m−1xy′−m(1−x2)2m−1y+m(m−2)(1−x2)2m−2x2y将上面的式子代入到连带勒让德方程中,化为y(x)的微分方程:
(1−x2)y′′−2(m+1)xy′+[l(l+1)−m(m+1)]y=0上式可以使用级数解法直接求解,但是还有更简单的方法。观察上面的微分方程,就是勒让德方程逐项求导m次后得到的方程。
证明:
上面的微分方程是勒让德方程逐项求导m次后得到的方程。根据莱布尼茨求导法则,有:
(uv)(m)=uv(m)+1!mu′v(m−1)+2!m(m−1)u′′v(m−2)+⋯++k!m(m−1)(m−2)⋯(m−k+1)u(k)v(m−k)+⋯+u(m)v,
将勒让德方程:
(1−x2)P′′−2xP′+l(l+1)P=0
对x求导m次,结果为:
[(1−x2)P(m+2)−m2xP(m+1)−2m(m−1)2P(m)]−2[xP(m+1)+mP(m)]+l(l+1)P(m)=0
化简得:
(1−x2)P(m)′′−2(m+1)xP(m)′+[l(l+1)−m(m+1)]P(m)=0
即为l阶连带勒让德方程。其中:
y(x)=P(m)(x)
前面已经学过勒让德方程和自然边界条件(在x=±1有限)构成本征值问题,本征值为l(l+1),而l为整数,本征函数则是勒让德多项式Pl(x),那么l阶连带勒让德方程与自然边界条件也构成本征值问题,本征值为l(l+1),本征函数为Pl(x)的m阶导数,即:
y(x)=Pl(m)(x)带回原连带勒让德方程,得到:
Plm(x)=Θ=(1−x2)2mPl(m)(x)其中Plm(x)为连带勒让德函数。由于Pl(x)为l次多项式,最多只能求导l次,所以l≥m,当m=0时,方程和结果均退化为普通的l阶勒让德方程和勒让德多项式。
2、连带勒让德函数的微分表示#
由勒让德多项式的微分表示可以得到连带勒让德函数的微分表示:
Plm(x)=2ll!(1−x2)2mdxl+mdl+m(x2−1)l这也叫做罗德里格斯公式,观察公式可知,当l−m为偶数时,Plm(x)为偶函数,否则为奇函数。
在之前的推导中,一直是把m看作是正整数,但是在连带勒让德方程中只出现了m2,所以可以将−m代入,连带勒让德方程仍然不变,此时得到连带勒让德函数为:
Pl−m(x)=2ll!(1−x2)−2mdxl−mdl−m(x2−1)l这也是连带勒让德方程的解且满足自然边界条件,但是他实际上与Plm是相关的,他们只相差一个常数,将两者作比,得到:
Plm(x)=(−1)m(l−m)!(l+m)!Pl−m(x)3、连带勒让德函数的积分表示#
按照柯西公式,微分表达式可以表示成环路积分:
Plm(x)=2l(1−x2)2m2πi1l!(l+m)!∮C(z−x)l+m+1(z2−1)ldzC为z平面上围绕z=x的任意闭合回路,这也叫作施列夫利积分。同样地,这也能进一步表示为定积分,取C为圆周,圆心在z=x,半径为∣x2−1∣,在C上,z−x=x2−1eiψ,dz=ix2−1eiψdψ,代入到上式,得:
Plm(x)=2πiml!(l+m)!∫−ππe−imψ[x+x2−121(e−iψ+eiψ)]ldψ这也叫作拉普拉斯积分。
4、连带勒让德函数的正交关系#
作为施图姆-刘维尔本征值问题的正交关系特例,不同阶的勒让德多项式在(−1,+1)上正交,即:
∫−1+1Pkm(x)Plm(x)dx=0(l=k)如果从x回到原来的变数θ,则:
∫0πPkm(cosθ)Plm(cosθ)sinθdθ(l=k)5、连带勒让德函数的模#
计算连带勒让德函数的模时,可以利用Plm(x)与Pl−m(x)的关系,有:
(Nlm)2=∫−1+1[Plm(x)]2dx=(−1)m(l−m)!(l+m)!∫−1+1Pl−m(x)Plm(x)dx=(−1)m(l−m)!(l+m)!22l(l!)21∫−1+1dxl−mdl−m(x2−1)l×dxl+mdl+m(x2−1)ldx仿照勒让德多项式计算模的方法,进行连续地分部积分,得到:
Nlm=(l−m)!(2l+1)(l+m)!26、广义傅里叶级数#
根据上面的证明,可知勒让德函数Plm(x)是完备的,可以作为广义傅里叶级数展开的基,定义在x的区间[−1,+1]上的函数f(x)或者定义在θ区间[0,π]上的函数f(θ)都可以展开为广义傅里叶级数。
⎩⎨⎧f(x)=l=0∑∞flPlm(x)fl=22l+1(l+m)!(l−m)!∫−1+1f(x)Plm(x)dx或:
⎩⎨⎧f(θ)=l=0∑∞flPlm(cosθ),fl=22l+1(l+m)!(l−m)!∫0πf(θ)Plm(cosθ)sinθdθ7、连带勒让德函数的递推公式#
连带勒让德函数的基本递推公式有四个:
(2k+1)xPkm(x)=(k+m)Pk−1m(x)+(k−m+1)Pk+1m(x)(k≥1)(2k+1)(1−x2)1/2Pkm(x)=Pk+1m+1(x)−Pk−1m+1(x)(k≥1)(2k+1)(1−x2)1/2Pkm(x)=(k+m)(k+m−1)Pk−1m−1(x)−(k−m+2)(k−m+1)Pk+1m−1(x)(k≥1)(2k+1)(1−x2)dxdPkm(x)=(k+1)(k+m)Pk−1m(x)−k(k−m+1)Pk+1m(x)(k≥1)
三、一般的球函数#
1、球函数#
球函数方程的解叫做球函数:
Ylm(θ,φ)=Plm(cosθ)(Asinmφ+Bcosmφ)(m=0,1,2,⋯,l;l=0,1,2,3,⋯)线性独立的l阶球函数共有2l+1个,这是因为对于m=0,有一个球函数Pl(x),当m=0,时各有两个球函数Plm(x)cosmφ,Plm(x)sinmφ。根据欧拉公式,可以组合为:
Ylm(θ,φ)=Pl∣m∣(cosθ)eimφ(m=−l,⋯,0,1,⋯,l;l=0,1,2,⋯)2、球函数的正交关系#
球函数中的任意两个在球面S上(即0≤θ≤π,0≤φ≤2π)正交:
∬SYlm(θ,φ)Ykn(θ,φ)sinθdθdφ=0(m=n或l=k)3、球函数的模#
使用三角函数表示的球函数进行积分:
(Nlm)2=∬S[Ylm(θ,φ)]2sinθdθdφ=∫0π[Plm(cosθ)]2sinθdθ∫02π(Asin2mφ+Bcos2mφ)dφ=∫−1+1[Plm(x)]2dx∫02π(Asin2mφ+Bcos2mφ)dφ前一个积分值已经得出,后一个积分不会积建议重修微积分。代入得到球函数的模:
Nlm=2l+12πδm(l−m)!(l+m)!其中δm={2,m=01,m=1,2,3,⋯。
因为:
∫02πeimφ[eimφ]∗dφ=∫02π1dφ=2π所以使用复数形式的球函数计算模的平方:
(Nlm)2=∬SYlm(θ,φ)[Ylm(θ,φ)]∗sinθdθdφ=∫0π[Pl∣m∣(cosθ)]2sinθdθ⋅∫02πeimφ[eimφ]∗dφ=2l+12⋅(l−∣m∣)!(l+∣m∣)!⋅2π得到:
Nlm=2l+14π⋅(l−∣m∣)!(l+∣m∣)!4、球面上的函数的广义傅里叶级数#
定义在球面S(0≤θ≤π,0≤φ≤2π)上的函数f(θ,φ)可以用球函数展开成二重广义傅里叶级数。现在以三角函数形式的球函数为基,把函数f进行展开:
f(θ,φ)=m=0∑∞[Am(θ)cosmφ+Bm(θ)sinmφ]这里θ作为参数出现在傅里叶系数中:
⎩⎨⎧Am(θ)Bm(θ)=πδm1∫02πf(θ,φ)cosmφdφ=π1∫02πf(θ,φ)sinmφdφ然后以Plm(cosθ)为基,继续把Am(θ),Bm(θ)展开,有:
⎩⎨⎧Am(θ)=l=m∑∞AlmPlm(cosθ)Bm(θ)=l=m∑∞BlmPlm(cosθ)然后利用广义傅里叶展开,得到:
⎩⎨⎧AlmBlm=22l+1⋅(l+m)!(l−m)!∫0πAm(θ)Plm(cosθ)sinθdθ=2πδm2l+1⋅(l+m)!(l−m)!∫0π∫02πf(θ,φ)Plm(cosθ)cosmφsinθdθdφ=22l+1⋅(l+m)!(l−m)!∫0πBm(θ)Plm(cosθ)sinθdθ=2π2l+1⋅(l+m)!(l−m)!∫0π∫02πf(θ,φ)Plm(cosθ)sinmφsinθdθdφ带回到f(θ,φ)在球面S上的展开式:
f(θ,φ)=m=0∑∞l=m∑∞[Almcosmφ+Blmsinmφ]Plm(cosθ)这里两个求和符号可以交换。如果想要展开复数形式的广义傅里叶级数,则有:
f(θ,φ)=l=0∑∞m=−l∑lClmPl∣m∣(cosθ)eimφ其中系数Clm的计算公式为:
Clm=4π2l+1(l+∣m∣)!(l−∣m∣)!∫0π∫02πf(θ,φ)Pl∣m∣(cosθ)[eimφ]∗sinθdθdφ5、正交归一化的球函数#
物理学常常使用正交归一化的球函数,定义为:
Ylm(θ,φ)=Nlm1Ylm(θ,φ)=4π2l+1⋅(l+∣m∣)!(l−∣m∣)!Pl∣m∣(cosθ)eimφ(l=0,1,2,⋯;m=−l,−l+1,⋯,0,⋯,l−1,l)在这个正交归一化的球函数下,其正交归一,球面上的函数f(θ,φ)可以使用正交归一球函数展开,即:
f(θ,φ)=l=0∑∞m=−l∑lClmYlm(θ,φ)广义傅里叶系数Clm为:
Clm=∫02π∫0πf(θ,φ)Ylm∗(θ,φ)sin2θdθdφ