関連ルジャンドル多項式の直交性
📂関数関連ルジャンドル多項式の直交性
定理
区間[−1,1]で固定されたmに対する関連ルジャンドル多項式は、直交集合を形成する。
∫−11Plm(x)Pkm(x)dx=2l+12(l−m)!(l+m)!δlk
x=cosθの場合、
∫0πPlm(cosθ)Pkm(cosθ)sinθdθ=2l+12(l−m)!(l+m)!δlk
関連ルジャンドル多項式
Plm(x)=(1−x2)2m2ll!1dxl+mdl+m(x2−1)l
証明
便宜上、簡単にPlm=Plm(x)と表記しよう。関連ルジャンドル微分方程式は以下の通り。
dxd[(1−x2)Plm′]+[l(l+1)−1−x2m2]Plm=0
場合 1: l=k
証明方法は、ルジャンドル多項式の直交性を示すのと同じだ。(1)をlとkに対して書くと、
dxd[(1−x2)Plm′]+[l(l+1)−1−x2m2]Plm=0dxd[(1−x2)Pkm′]+[k(k+1)−1−x2m2]Pkm=0
(lに対する式)にPkmを掛け、(kに対する式)にPlmを掛けて互いに引くと、以下のようになる。
Pkmdxd[(1−x2)Plm′]−Plmdxd[(1−x2)Pkm′]+[l(l+1)−k(k+1)]PlmPkm=0
ここで第一項、第二項を次のように整理できる。
Pkmdxd[(1−x2)Plm′]−Plmdxd[(1−x2)Pkm′]=Pkm(1−x2)′Plm′+Pkm(1−x2)Plm′′−Plm(1−x2)′Pkm′−Plm(1−x2)Pkm′′=Pkm(1−x2)′Plm′+Pkm(1−x2)Plm′′−Plm(1−x2)′Pkm′−Plm(1−x2)Pkm′′+Pkm′(1−x2)Plm′−Pkm′(1−x2)Plm′=dxd[(1−x2)Plm′Pkm−(1−x2)PlmPkm′]=dxd[(1−x2)(Plm′Pkm−PlmPkm′)]
これを(2)に代入すると、
dxd[(1−x2)(Plm′Pkm−PlmPkm′)]+[l(l+1)−k(k+1)]PlmPkm=0
両辺を区間[−1,1]で積分すると、以下を得る。
[(1−x2)(Plm′Pkm−PlmPkm′)]−11+[l(l+1)−k(k+1)]∫−11PlmPkmdx=0
最初の項は0なので、以下の通り。
[l(l+1)−k(k+1)]∫−11PlmPkmdx=0
この時、l=kなのでl(l+1)−k(k+1)=0となる。したがって、以下の通り。
∫−11Plm(x)Pkm(x)dx=0
■
場合 2: l=k
補助定理
次が成り立つ。
dxl+mdl+m(x2−1)l=(l−m)!(l+m)!(x2−1)−mdxl−mdl−m(x2−1)l(3)
上記の公式を関連ルジャンドル多項式に適用すると、
Plm(x)=(1−x2)2m2ll!1dxl+mdl+m(x2−1)l=(1−x2)2m2ll!1(l−m)!(l+m)!(x2−1)−mdxl−mdl−m(x2−1)l=2ll!(−1)m(l−m)!(l+m)!(1−x2)−2mdxl−mdl−m(x2−1)l
上記の式の両辺を二乗すると、以下を得る。
[Plm(x)]2=22l(l!)21[(l−m)!(l+m)!]2(1−x2)−mdxl−mdl−m(x2−1)ldxl−mdl−m(x2−1)l
上記の式に(3)を代入すると、
[Plm(x)]2=22l(l!)21[(l−m)!(l+m)!]2(1−x2)−m[(l+m)!(l−m)!(x2−1)mdxl+mdl+m(x2−1)l]dxl−mdl−m(x2−1)l=22l(l!)2(−1)m(l−m)!(l+m)!dxl+mdl+m(x2−1)ldxl−mdl−m(x2−1)l
今、両辺を区間[−1,1]で積分すると以下のようになる。
∫−11[Plm(x)]2dx=22l(l!)2(−1)m(l−m)!(l+m)!∫−11[dxl+mdl+m(x2−1)ldxl−mdl−m(x2−1)l]dx(4)
右辺の積分部分だけを見る。部分積分で解くと、次のようになる。
∫−11[dxl+mdl+m(x2−1)ldxl−mdl−m(x2−1)l]dx=∫−11[dxl+m−1dl+m−1(x2−1)l]′dxl−mdl−m(x2−1)ldx=[dxl+m−1dl+m−1(x2−1)ldxl−mdl−m(x2−1)l]−11−∫−11dxl+m−1dl+m−1(x2−1)ldxl−m+1dl−m+1(x2−1)ldx
ここで第一項は0である。(x2−1)lは2l次の多項式であり、∣m∣<lなのでl+m−1とl−mはlより小さく、少なくとも(x2−1)が微分されずに残るからだ。ここに±1を代入すると0になる。残りの項を再び部分積分で解いてみると、
−∫−11dxl+m−1dl+m−1(x2−1)ldxl−m+1dl−m+1(x2−1)ldx=[−dxl+m−2dl+m−2(x2−1)ldxl−m+1dl−m+1(x2−1)l]+∫−11dxl+m−2dl+m−2(x2−1)ldxl−m+2dl−m+2(x2−1)ldx
ここで最初の項は上記の理由で0である。このような部分積分をm回繰り返すと、以下の式を得る。
∫−11dxl+mdl+m(x2−1)ldxl−mdl−m(x2−1)ldx=(−1)m∫−11dxldl(x2−1)ldxldl(x2−1)ldx
したがって、(4)は以下の通り。
∫−11[Plm(x)]2dx=22l(l!)21(l−m)!(l+m)!∫−11[dxldl(x2−1)l]2dx
ロドリゲスの公式
Pl(x)=2ll!1dxldl(x2−1)l
それにより、ロドリゲスの公式により、以下が導かれる。
∫−11[Plm(x)]2dx=22l(l!)21(l−m)!(l+m)!22l(l!)2∫−11[Pl(x)]2dx=(l−m)!(l+m)!∫−11[Pl(x)]2dx
最後に、ルジャンドル多項式の直交性により、最終的に以下が得られる。
∫−11[Plm(x)]2dx=2l+12(l−m)!(l+m)!
■