一般ルジャンドル方程式の正準解
数学 において 、 ルジャンドル多項式は、 一般ルジャンドル方程式 の標準解である。
(
1
−
x
2
)
d
2
d
x
2
ポ
ℓ
メートル
(
x
)
−
2
x
d
d
x
ポ
ℓ
メートル
(
x
)
+
[
ℓ
(
ℓ
+
1
)
−
メートル
2
1
−
x
2
]
ポ
ℓ
メートル
(
x
)
=
0
、
{\displaystyle \left(1-x^{2}\right){\frac {d^{2}}{dx^{2}}}P_{\ell }^{m}(x)-2x{\frac {d}{dx}}P_{\ell }^{m}(x)+\left[\ell (\ell +1)-{\frac {m^{2}}{1-x^{2}}}\right]P_{\ell }^{m}(x)=0,}
または同等
d
d
x
[
(
1
−
x
2
)
d
d
x
ポ
ℓ
メートル
(
x
)
]
+
[
ℓ
(
ℓ
+
1
)
−
メートル
2
1
−
x
2
]
ポ
ℓ
メートル
(
x
)
=
0
、
{\displaystyle {\frac {d}{dx}}\left[\left(1-x^{2}\right){\frac {d}{dx}}P_{\ell }^{m}(x)\right]+\left[\ell (\ell +1)-{\frac {m^{2}}{1-x^{2}}}\right]P_{\ell }^{m}(x)=0,}
ここで、添え字 ℓ と m (整数) は、それぞれルジャンドル多項式の次数と位数と呼ばれます。この方程式は、 ℓ と m が 0 ≤ m ≤ ℓ の整数であるか、自明に同値な負の値である 場合に限り、 [−1, 1]上で非特異な非ゼロ解を持ちます。さらに m が偶数の場合、関数は 多項式 です。m が 0 で ℓ が整数の場合、これらの関数は ルジャンドル多項式 と同一です 。一般に、ℓ と m が整数の場合 、 正規 解は m が奇数のときは 多項式 ではありませんが、「ルジャンドル多項式」と呼ばれることがあります 。ℓ と m の 任意 の 実数値または複素数値を持つ関数の完全に一般的なクラスは、 ルジャンドル関数 です 。その場合、パラメータは通常ギリシャ文字でラベル付けされます。
ルジャンドル 常微分方程式は、 物理学 やその他の技術分野で頻繁に登場します 。特に、 球面座標で ラプラス方程式 (および関連する 偏微分方程式 )を解くときに発生します。ルジャンドル多項式は、 球面調和関数 の定義において重要な役割を果たします 。
非負整数パラメータの定義 ℓ そして メートル
これらの関数は と表記され、上付き文字は P のべき乗ではなく次数を示す。最も簡単な定義は、通常の ルジャンドル多項式 ( m ≥ 0)
の導関数である。
ポ
ℓ
メートル
(
x
)
{\displaystyle P_{\ell}^{m}(x)}
ポ
ℓ
メートル
(
x
)
=
(
−
1
)
メートル
(
1
−
x
2
)
メートル
/
2
d
メートル
d
x
メートル
(
ポ
ℓ
(
x
)
)
、
{\displaystyle P_{\ell}^{m}(x)=(-1)^{m}(1-x^{2})^{m/2}{\frac {d^{m}}{dx^{m}}}\left(P_{\ell}(x)\right),}
この式の(−1) m 因子は コンドン・ショートリー位相 として知られている。著者によっては省略している。この式で記述される関数が、パラメータ ℓ と m の指定された値を持つ一般ルジャンドル微分方程式を満たすことは、 ルジャンドル 方程式を Pℓ についてm 回微分することによってわかる 。 [1]
(
1
−
x
2
)
d
2
d
x
2
ポ
ℓ
(
x
)
−
2
x
d
d
x
ポ
ℓ
(
x
)
+
ℓ
(
ℓ
+
1
)
ポ
ℓ
(
x
)
=
0.
{\displaystyle \left(1-x^{2}\right){\frac {d^{2}}{dx^{2}}}P_{\ell }(x)-2x{\frac {d}{dx}}P_{\ell }(x)+\ell (\ell +1)P_{\ell }(x)=0.}
さらに、ロドリゲスの公式 に よれば 、
P
ポ
ℓ
(
x
)
=
1
2
ℓ
ℓ
!
d
ℓ
d
x
ℓ
[
(
x
2
−
1
)
ℓ
]
、
{\displaystyle P_{\ell }(x)={\frac {1}{2^{\ell }\,\ell !}}\ {\frac {d^{\ell }}{dx^{\ell }}}\left[(x^{2}-1)^{\ell }\right],}
メートル リットル という形で表現できる。
ポ
ℓ
メートル
(
x
)
=
(
−
1
)
メートル
2
ℓ
ℓ
!
(
1
−
x
2
)
メートル
/
2
d
ℓ
+
メートル
d
x
ℓ
+
メートル
(
x
2
−
1
)
ℓ
。
{\displaystyle P_{\ell }^{m}(x)={\frac {(-1)^{m}}{2^{\ell }\ell !}}(1-x^{2}) ^{m/2}\ {\frac {d^{\ell +m}}{dx^{\ell +m}}}(x^{2}-1)^{\ell }.}
この式により、 m の範囲を 次のように拡張できます: − ℓ ≤ m ≤ ℓ 。この式に ± mを代入して得られる P ℓ ± m の定義は 比例的です。確かに、の左辺と右辺の等指数の係数を等しくする
と、比例定数は
次のように
なります。
d
ℓ
−
メートル
d
x
ℓ
−
メートル
(
x
2
−
1
)
ℓ
=
c
l
メートル
(
1
−
x
2
)
メートル
d
ℓ
+
メートル
d
x
ℓ
+
メートル
(
x
2
−
1
)
ℓ
、
{\displaystyle {\frac {d^{\ell -m}}{dx^{\ell -m}}}(x^{2}-1)^{\ell }=c_{lm}(1-x^{2})^{m}{\frac {d^{\ell +m}}{dx^{\ell +m}}}(x^{2}-1)^{\ell },}
c
l
メートル
=
(
−
1
)
メートル
(
ℓ
−
メートル
)
!
(
ℓ
+
メートル
)
!
、
{\displaystyle c_{lm}=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}},}
ポ
ℓ
−
メートル
(
x
)
=
(
−
1
)
メートル
(
ℓ
−
メートル
)
!
(
ℓ
+
メートル
)
!
ポ
ℓ
メートル
(
x
)
。
{\displaystyle P_{\ell }^{-m}(x)=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}}P_{\エル }^{m}(x).}
代替表記
文献では以下の代替表記も使用されている: [2]
ポ
ℓ
メートル
(
x
)
=
(
−
1
)
メートル
ポ
ℓ
メートル
(
x
)
{\displaystyle P_{\ell m}(x)=(-1)^{m}P_{\ell }^{m}(x)}
ルジャンドル多項式は、
単純な単項式と 二項係数の一般化された形式 を使用して次のように記述することもできます。
[ 引用が必要 ]
ポ
l
メートル
(
x
)
=
(
−
1
)
メートル
⋅
2
l
⋅
(
1
−
x
2
)
メートル
/
2
⋅
∑
け
=
メートル
l
け
!
(
け
−
メートル
)
!
⋅
x
け
−
メートル
⋅
(
l
け
)
(
l
+
け
−
1
2
l
)
{\displaystyle P_{l}^{m}(x)=(-1)^{m}\cdot 2^{l}\cdot (1-x^{2})^{m/2}\cdot \sum _{k=m}^{l}{\frac {k!}{(km)!}}\cdot x^{km}\cdot {\binom {l}{k}}{\binom {\frac {l+k-1}{2}}{l}}}
直交性
ルジャンドルの従属多項式は一般に互いに直交しません。たとえば、 は と直交しません 。ただし、一部のサブセットは直交します。 0 ≤ m ≤ ℓと仮定すると、それらは m を 固定した場合の直交条件を満たします 。
ポ
1
1
{\displaystyle P_{1}^{1}}
ポ
2
2
{\displaystyle P_{2}^{2}}
∫
−
1
1
ポ
け
メートル
ポ
ℓ
メートル
d
x
=
2
(
ℓ
+
メートル
)
!
(
2
ℓ
+
1
)
(
ℓ
−
メートル
)
!
δ
け
、
ℓ
{\displaystyle \int _{-1}^{1}P_{k}^{m}P_{\ell }^{m}dx={\frac {2(\ell +m)!}{(2\ ell +1)(\ell -m)!}}\ \delta _{k,\ell }}
ここで δ k , ℓ は クロネッカーのデルタ です 。
また、これらはℓを 固定した場合の直交条件を満たす 。
∫
−
1
1
ポ
ℓ
メートル
ポ
ℓ
ん
1
−
x
2
d
x
=
{
0
もし
メートル
≠
ん
(
ℓ
+
メートル
)
!
メートル
(
ℓ
−
メートル
)
!
もし
メートル
=
ん
≠
0
∞
もし
メートル
=
ん
=
0
{\displaystyle \int _{-1}^{1}{\frac {P_{\ell }^{m}P_{\ell }^{n}}{1-x^{2}}}dx={\begin{cases}0&{\text{if }}m\neq n\\{\frac {(\ell +m)!}{m(\ell -m)!}}&{\text{if }}m=n\neq 0\\\infty &{\text{if }}m=n=0\end{cases}}}
ネガティブ メートル および/または否定的 ℓ
微分方程式は、 m の符号の変化に対して明らかに不変です 。
負のm に対する関数は正の m に対する関数に比例することが上で示されました 。
ポ
ℓ
−
メートル
=
(
−
1
)
メートル
(
ℓ
−
メートル
)
!
(
ℓ
+
メートル
)
!
ポ
ℓ
メートル
{\displaystyle P_{\ell }^{-m}=(-1)^{m}{\frac {(\ell -m)!}{(\ell +m)!}}P_{\ell }^ {m}}
(これはロドリゲスの公式の定義に由来する。この定義により、さまざまな漸化式が正または負の m に対しても機能するようになる。)
もし
|
メートル
|
>
ℓ
それから
ポ
ℓ
メートル
=
0.
{\displaystyle {\text{If}}\quad |m|>\ell \,\quad {\text{then}}\quad P_{\ell }^{m}=0.\,}
この微分方程式はℓから −ℓ − 1 への 変化に対しても不変であり 、 ℓが 負の場合の関数は次のように定義される。
ポ
−
ℓ
メートル
=
ポ
ℓ
−
1
メートル
、
(
ℓ
=
1
、
2
、
…
)
。
{\displaystyle P_{-\ell }^{m}=P_{\ell -1}^{m},\ (\ell =1,\,2,\,\dots ).}
パリティ
定義から、ルジャンドル関数は偶数か奇数かが次のように表せることが分かる。
ポ
ℓ
メートル
(
−
x
)
=
(
−
1
)
ℓ
−
メートル
ポ
ℓ
メートル
(
x
)
{\displaystyle P_{\ell}^{m}(-x)=(-1)^{\ell -m}P_{\ell}^{m}(x)}
最初のいくつかのルジャンドル関数
m = 0 のルジャンドル関数
m = 1 のルジャンドル関数
m = 2 のルジャンドル関数
m が負の値の場合も含め、最初のいくつかのルジャンドル関数は次 のようになります。
ポ
0
0
(
x
)
=
1
{\displaystyle P_{0}^{0}(x)=1}
ポ
1
−
1
(
x
)
=
−
1
2
ポ
1
1
(
x
)
ポ
1
0
(
x
)
=
x
ポ
1
1
(
x
)
=
−
(
1
−
x
2
)
1
/
2
{\displaystyle {\begin{aligned}P_{1}^{-1}(x)&=-{\tfrac {1}{2}}P_{1}^{1}(x)\\P_{1}^{0}(x)&=x\\P_{1}^{1}(x)&=-(1-x^{2})^{1/2}\end{aligned}}}
P
2
−
2
(
x
)
=
1
24
P
2
2
(
x
)
P
2
−
1
(
x
)
=
−
1
6
P
2
1
(
x
)
P
2
0
(
x
)
=
1
2
(
3
x
2
−
1
)
P
2
1
(
x
)
=
−
3
x
(
1
−
x
2
)
1
/
2
P
2
2
(
x
)
=
3
(
1
−
x
2
)
{\displaystyle {\begin{aligned}P_{2}^{-2}(x)&={\tfrac {1}{24}}P_{2}^{2}(x)\\P_{2}^{-1}(x)&=-{\tfrac {1}{6}}P_{2}^{1}(x)\\P_{2}^{0}(x)&={\tfrac {1}{2}}(3x^{2}-1)\\P_{2}^{1}(x)&=-3x(1-x^{2})^{1/2}\\P_{2}^{2}(x)&=3(1-x^{2})\end{aligned}}}
P
3
−
3
(
x
)
=
−
1
720
P
3
3
(
x
)
P
3
−
2
(
x
)
=
1
120
P
3
2
(
x
)
P
3
−
1
(
x
)
=
−
1
12
P
3
1
(
x
)
P
3
0
(
x
)
=
1
2
(
5
x
3
−
3
x
)
P
3
1
(
x
)
=
3
2
(
1
−
5
x
2
)
(
1
−
x
2
)
1
/
2
P
3
2
(
x
)
=
15
x
(
1
−
x
2
)
P
3
3
(
x
)
=
−
15
(
1
−
x
2
)
3
/
2
{\displaystyle {\begin{aligned}P_{3}^{-3}(x)&=-{\tfrac {1}{720}}P_{3}^{3}(x)\\P_{3}^{-2}(x)&={\tfrac {1}{120}}P_{3}^{2}(x)\\P_{3}^{-1}(x)&=-{\tfrac {1}{12}}P_{3}^{1}(x)\\P_{3}^{0}(x)&={\tfrac {1}{2}}(5x^{3}-3x)\\P_{3}^{1}(x)&={\tfrac {3}{2}}(1-5x^{2})(1-x^{2})^{1/2}\\P_{3}^{2}(x)&=15x(1-x^{2})\\P_{3}^{3}(x)&=-15(1-x^{2})^{3/2}\end{aligned}}}
P
4
−
4
(
x
)
=
1
40320
P
4
4
(
x
)
P
4
−
3
(
x
)
=
−
1
5040
P
4
3
(
x
)
P
4
−
2
(
x
)
=
1
360
P
4
2
(
x
)
P
4
−
1
(
x
)
=
−
1
20
P
4
1
(
x
)
P
4
0
(
x
)
=
1
8
(
35
x
4
−
30
x
2
+
3
)
P
4
1
(
x
)
=
−
5
2
(
7
x
3
−
3
x
)
(
1
−
x
2
)
1
/
2
P
4
2
(
x
)
=
15
2
(
7
x
2
−
1
)
(
1
−
x
2
)
P
4
3
(
x
)
=
−
105
x
(
1
−
x
2
)
3
/
2
P
4
4
(
x
)
=
105
(
1
−
x
2
)
2
{\displaystyle {\begin{aligned}P_{4}^{-4}(x)&={\tfrac {1}{40320}}P_{4}^{4}(x)\\P_{4}^{-3}(x)&=-{\tfrac {1}{5040}}P_{4}^{3}(x)\\P_{4}^{-2}(x)&={\tfrac {1}{360}}P_{4}^{2}(x)\\P_{4}^{-1}(x)&=-{\tfrac {1}{20}}P_{4}^{1}(x)\\P_{4}^{0}(x)&={\tfrac {1}{8}}(35x^{4}-30x^{2}+3)\\P_{4}^{1}(x)&=-{\tfrac {5}{2}}(7x^{3}-3x)(1-x^{2})^{1/2}\\P_{4}^{2}(x)&={\tfrac {15}{2}}(7x^{2}-1)(1-x^{2})\\P_{4}^{3}(x)&=-105x(1-x^{2})^{3/2}\\P_{4}^{4}(x)&=105(1-x^{2})^{2}\end{aligned}}}
これらの関数には、いくつかの繰り返しプロパティがあります。
(
ℓ
−
m
−
1
)
(
ℓ
−
m
)
P
ℓ
m
(
x
)
=
−
P
ℓ
m
+
2
(
x
)
+
P
ℓ
−
2
m
+
2
(
x
)
+
(
ℓ
+
m
)
(
ℓ
+
m
−
1
)
P
ℓ
−
2
m
(
x
)
{\displaystyle (\ell -m-1)(\ell -m)P_{\ell }^{m}(x)=-P_{\ell }^{m+2}(x)+P_{\ell -2}^{m+2}(x)+(\ell +m)(\ell +m-1)P_{\ell -2}^{m}(x)}
(
ℓ
−
m
+
1
)
P
ℓ
+
1
m
(
x
)
=
(
2
ℓ
+
1
)
x
P
ℓ
m
(
x
)
−
(
ℓ
+
m
)
P
ℓ
−
1
m
(
x
)
{\displaystyle (\ell -m+1)P_{\ell +1}^{m}(x)=(2\ell +1)xP_{\ell }^{m}(x)-(\ell +m)P_{\ell -1}^{m}(x)}
2
m
x
P
ℓ
m
(
x
)
=
−
1
−
x
2
[
P
ℓ
m
+
1
(
x
)
+
(
ℓ
+
m
)
(
ℓ
−
m
+
1
)
P
ℓ
m
−
1
(
x
)
]
{\displaystyle 2mxP_{\ell }^{m}(x)=-{\sqrt {1-x^{2}}}\left[P_{\ell }^{m+1}(x)+(\ell +m)(\ell -m+1)P_{\ell }^{m-1}(x)\right]}
1
1
−
x
2
P
ℓ
m
(
x
)
=
−
1
2
m
[
P
ℓ
−
1
m
+
1
(
x
)
+
(
ℓ
+
m
−
1
)
(
ℓ
+
m
)
P
ℓ
−
1
m
−
1
(
x
)
]
{\displaystyle {\frac {1}{\sqrt {1-x^{2}}}}P_{\ell }^{m}(x)={\frac {-1}{2m}}\left[P_{\ell -1}^{m+1}(x)+(\ell +m-1)(\ell +m)P_{\ell -1}^{m-1}(x)\right]}
1
1
−
x
2
P
ℓ
m
(
x
)
=
−
1
2
m
[
P
ℓ
+
1
m
+
1
(
x
)
+
(
ℓ
−
m
+
1
)
(
ℓ
−
m
+
2
)
P
ℓ
+
1
m
−
1
(
x
)
]
{\displaystyle {\frac {1}{\sqrt {1-x^{2}}}}P_{\ell }^{m}(x)={\frac {-1}{2m}}\left[P_{\ell +1}^{m+1}(x)+(\ell -m+1)(\ell -m+2)P_{\ell +1}^{m-1}(x)\right]}
1
−
x
2
P
ℓ
m
(
x
)
=
1
2
ℓ
+
1
[
(
ℓ
−
m
+
1
)
(
ℓ
−
m
+
2
)
P
ℓ
+
1
m
−
1
(
x
)
−
(
ℓ
+
m
−
1
)
(
ℓ
+
m
)
P
ℓ
−
1
m
−
1
(
x
)
]
{\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m}(x)={\frac {1}{2\ell +1}}\left[(\ell -m+1)(\ell -m+2)P_{\ell +1}^{m-1}(x)-(\ell +m-1)(\ell +m)P_{\ell -1}^{m-1}(x)\right]}
1
−
x
2
P
ℓ
m
(
x
)
=
−
1
2
ℓ
+
1
[
P
ℓ
+
1
m
+
1
(
x
)
−
P
ℓ
−
1
m
+
1
(
x
)
]
{\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m}(x)={\frac {-1}{2\ell +1}}\left[P_{\ell +1}^{m+1}(x)-P_{\ell -1}^{m+1}(x)\right]}
1
−
x
2
P
ℓ
m
+
1
(
x
)
=
(
ℓ
−
m
)
x
P
ℓ
m
(
x
)
−
(
ℓ
+
m
)
P
ℓ
−
1
m
(
x
)
{\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m+1}(x)=(\ell -m)xP_{\ell }^{m}(x)-(\ell +m)P_{\ell -1}^{m}(x)}
1
−
x
2
P
ℓ
m
+
1
(
x
)
=
(
ℓ
−
m
+
1
)
P
ℓ
+
1
m
(
x
)
−
(
ℓ
+
m
+
1
)
x
P
ℓ
m
(
x
)
{\displaystyle {\sqrt {1-x^{2}}}P_{\ell }^{m+1}(x)=(\ell -m+1)P_{\ell +1}^{m}(x)-(\ell +m+1)xP_{\ell }^{m}(x)}
1
−
x
2
d
d
x
P
ℓ
m
(
x
)
=
1
2
[
(
ℓ
+
m
)
(
ℓ
−
m
+
1
)
P
ℓ
m
−
1
(
x
)
−
P
ℓ
m
+
1
(
x
)
]
{\displaystyle {\sqrt {1-x^{2}}}{\frac {d}{dx}}{P_{\ell }^{m}}(x)={\frac {1}{2}}\left[(\ell +m)(\ell -m+1)P_{\ell }^{m-1}(x)-P_{\ell }^{m+1}(x)\right]}
(
1
−
x
2
)
d
d
x
P
ℓ
m
(
x
)
=
1
2
ℓ
+
1
[
(
ℓ
+
1
)
(
ℓ
+
m
)
P
ℓ
−
1
m
(
x
)
−
ℓ
(
ℓ
−
m
+
1
)
P
ℓ
+
1
m
(
x
)
]
{\displaystyle (1-x^{2}){\frac {d}{dx}}{P_{\ell }^{m}}(x)={\frac {1}{2\ell +1}}\left[(\ell +1)(\ell +m)P_{\ell -1}^{m}(x)-\ell (\ell -m+1)P_{\ell +1}^{m}(x)\right]}
(
x
2
−
1
)
d
d
x
P
ℓ
m
(
x
)
=
ℓ
x
P
ℓ
m
(
x
)
−
(
ℓ
+
m
)
P
ℓ
−
1
m
(
x
)
{\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)={\ell }xP_{\ell }^{m}(x)-(\ell +m)P_{\ell -1}^{m}(x)}
(
x
2
−
1
)
d
d
x
P
ℓ
m
(
x
)
=
−
(
ℓ
+
1
)
x
P
ℓ
m
(
x
)
+
(
ℓ
−
m
+
1
)
P
ℓ
+
1
m
(
x
)
{\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)=-(\ell +1)xP_{\ell }^{m}(x)+(\ell -m+1)P_{\ell +1}^{m}(x)}
(
x
2
−
1
)
d
d
x
P
ℓ
m
(
x
)
=
1
−
x
2
P
ℓ
m
+
1
(
x
)
+
m
x
P
ℓ
m
(
x
)
{\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)={\sqrt {1-x^{2}}}P_{\ell }^{m+1}(x)+mxP_{\ell }^{m}(x)}
(
x
2
−
1
)
d
d
x
P
ℓ
m
(
x
)
=
−
(
ℓ
+
m
)
(
ℓ
−
m
+
1
)
1
−
x
2
P
ℓ
m
−
1
(
x
)
−
m
x
P
ℓ
m
(
x
)
{\displaystyle (x^{2}-1){\frac {d}{dx}}{P_{\ell }^{m}}(x)=-(\ell +m)(\ell -m+1){\sqrt {1-x^{2}}}P_{\ell }^{m-1}(x)-mxP_{\ell }^{m}(x)}
役に立つアイデンティティ(最初の再帰の初期値):
P
ℓ
+
1
ℓ
+
1
(
x
)
=
−
(
2
ℓ
+
1
)
1
−
x
2
P
ℓ
ℓ
(
x
)
{\displaystyle P_{\ell +1}^{\ell +1}(x)=-(2\ell +1){\sqrt {1-x^{2}}}P_{\ell }^{\ell }(x)}
P
ℓ
ℓ
(
x
)
=
(
−
1
)
ℓ
(
2
ℓ
−
1
)
!
!
(
1
−
x
2
)
(
ℓ
/
2
)
{\displaystyle P_{\ell }^{\ell }(x)=(-1)^{\ell }(2\ell -1)!!(1-x^{2})^{(\ell /2)}}
P
ℓ
+
1
ℓ
(
x
)
=
x
(
2
ℓ
+
1
)
P
ℓ
ℓ
(
x
)
{\displaystyle P_{\ell +1}^{\ell }(x)=x(2\ell +1)P_{\ell }^{\ell }(x)}
二 重階 乗
で !!
3 つのルジャンドル多項式 (次数が同じ) の積の積分は、ルジャンドル多項式の積をルジャンドル多項式の線型級数に展開するときに必要な要素です。たとえば、クーロン 演算子の行列要素が必要な ハートリー・フォック 多様体の原子計算を行うときにこれが必要になります。これには、ゴーントの公式 [3] があります。
この公式は、次の仮定の下で使用されます。
1
2
∫
−
1
1
P
l
u
(
x
)
P
m
v
(
x
)
P
n
w
(
x
)
d
x
=
(
−
1
)
s
−
m
−
w
(
m
+
v
)
!
(
n
+
w
)
!
(
2
s
−
2
n
)
!
s
!
(
m
−
v
)
!
(
s
−
l
)
!
(
s
−
m
)
!
(
s
−
n
)
!
(
2
s
+
1
)
!
×
∑
t
=
p
q
(
−
1
)
t
(
l
+
u
+
t
)
!
(
m
+
n
−
u
−
t
)
!
t
!
(
l
−
u
−
t
)
!
(
m
−
n
+
u
+
t
)
!
(
n
−
w
−
t
)
!
{\displaystyle {\begin{aligned}{\frac {1}{2}}\int _{-1}^{1}P_{l}^{u}(x)P_{m}^{v}(x)P_{n}^{w}(x)dx={}&{}(-1)^{s-m-w}{\frac {(m+v)!(n+w)!(2s-2n)!s!}{(m-v)!(s-l)!(s-m)!(s-n)!(2s+1)!}}\\&{}\times \ \sum _{t=p}^{q}(-1)^{t}{\frac {(l+u+t)!(m+n-u-t)!}{t!(l-u-t)!(m-n+u+t)!(n-w-t)!}}\end{aligned}}}
度数は負でない整数である
l
,
m
,
n
≥
0
{\displaystyle l,m,n\geq 0}
3つの順序はすべて非負の整数である
u
,
v
,
w
≥
0
{\displaystyle u,v,w\geq 0}
u
{\displaystyle u}
3つの秩序の中で最大の秩序である
注文は合計すると
u
=
v
+
w
{\displaystyle u=v+w}
度合いは従う
m
≥
n
{\displaystyle m\geq n}
式に現れる他の量は次のように定義される。
2
s
=
l
+
m
+
n
{\displaystyle 2s=l+m+n}
p
=
max
(
0
,
n
−
m
−
u
)
{\displaystyle p=\max(0,\,n-m-u)}
q
=
min
(
m
+
n
−
u
,
l
−
u
,
n
−
w
)
{\displaystyle q=\min(m+n-u,\,l-u,\,n-w)}
積分はゼロとなる。
度数の合計は偶数なので 整数となる
s
{\displaystyle s}
三角形の条件は満たされる
m
+
n
≥
l
≥
m
−
n
{\displaystyle m+n\geq l\geq m-n}
ドンとレムス(2002) [4]は、 この式の導出を、任意の数のルジャンドル多項式の積の積分に一般化した。
超幾何関数による一般化
これらの関数は、実際には一般的な複雑なパラメータと引数に対して定義できます。
P
λ
μ
(
z
)
=
1
Γ
(
1
−
μ
)
[
1
+
z
1
−
z
]
μ
/
2
2
F
1
(
−
λ
,
λ
+
1
;
1
−
μ
;
1
−
z
2
)
{\displaystyle P_{\lambda }^{\mu }(z)={\frac {1}{\Gamma (1-\mu )}}\left[{\frac {1+z}{1-z}}\right]^{\mu /2}\,_{2}F_{1}(-\lambda ,\lambda +1;1-\mu ;{\frac {1-z}{2}})}
ここで は ガンマ関数 であり 、 は 超幾何関数 である。
Γ
{\displaystyle \Gamma }
2
F
1
{\displaystyle _{2}F_{1}}
2
F
1
(
α
,
β
;
γ
;
z
)
=
Γ
(
γ
)
Γ
(
α
)
Γ
(
β
)
∑
n
=
0
∞
Γ
(
n
+
α
)
Γ
(
n
+
β
)
Γ
(
n
+
γ
)
n
!
z
n
,
{\displaystyle \,_{2}F_{1}(\alpha ,\beta ;\gamma ;z)={\frac {\Gamma (\gamma )}{\Gamma (\alpha )\Gamma (\beta )}}\sum _{n=0}^{\infty }{\frac {\Gamma (n+\alpha )\Gamma (n+\beta )}{\Gamma (n+\gamma )\ n!}}z^{n},}
より一般的な方法で定義される関数は、ルジャンドル関数 と呼ばれます 。ルジャンドル関数は、前と同じ微分方程式を満たします。
(
1
−
z
2
)
y
″
−
2
z
y
′
+
(
λ
[
λ
+
1
]
−
μ
2
1
−
z
2
)
y
=
0.
{\displaystyle (1-z^{2})\,y''-2zy'+\left(\lambda [\lambda +1]-{\frac {\mu ^{2}}{1-z^{2}}}\right)\,y=0.\,}
これは2階微分方程式なので、2番目の解は次 のように定義されます。
Q
λ
μ
(
z
)
{\displaystyle Q_{\lambda }^{\mu }(z)}
Q
λ
μ
(
z
)
=
π
Γ
(
λ
+
μ
+
1
)
2
λ
+
1
Γ
(
λ
+
3
/
2
)
1
z
λ
+
μ
+
1
(
1
−
z
2
)
μ
/
2
2
F
1
(
λ
+
μ
+
1
2
,
λ
+
μ
+
2
2
;
λ
+
3
2
;
1
z
2
)
{\displaystyle Q_{\lambda }^{\mu }(z)={\frac {{\sqrt {\pi }}\ \Gamma (\lambda +\mu +1)}{2^{\lambda +1}\Gamma (\lambda +3/2)}}{\frac {1}{z^{\lambda +\mu +1}}}(1-z^{2})^{\mu /2}\,_{2}F_{1}\left({\frac {\lambda +\mu +1}{2}},{\frac {\lambda +\mu +2}{2}};\lambda +{\frac {3}{2}};{\frac {1}{z^{2}}}\right)}
P
λ
μ
(
z
)
{\displaystyle P_{\lambda }^{\mu }(z)}
そして、 どちらも前述のさまざまな漸化式に従います。
Q
λ
μ
(
z
)
{\displaystyle Q_{\lambda }^{\mu }(z)}
角度による再パラメータ化
これらの関数は、引数が角度に関して再パラメータ化され、次のようになる場合に最も役立ちます 。
x
=
cos
θ
{\displaystyle x=\cos \theta }
P
ℓ
m
(
cos
θ
)
=
(
−
1
)
m
(
sin
θ
)
m
d
m
d
(
cos
θ
)
m
(
P
ℓ
(
cos
θ
)
)
{\displaystyle P_{\ell }^{m}(\cos \theta )=(-1)^{m}(\sin \theta )^{m}\ {\frac {d^{m}}{d(\cos \theta )^{m}}}\left(P_{\ell }(\cos \theta )\right)}
関係式を使用すると 、上記のリストは次のようにパラメータ化された最初のいくつかの多項式を生成します。
(
1
−
x
2
)
1
/
2
=
sin
θ
{\displaystyle (1-x^{2})^{1/2}=\sin \theta }
P
0
0
(
cos
θ
)
=
1
P
1
0
(
cos
θ
)
=
cos
θ
P
1
1
(
cos
θ
)
=
−
sin
θ
P
2
0
(
cos
θ
)
=
1
2
(
3
cos
2
θ
−
1
)
P
2
1
(
cos
θ
)
=
−
3
cos
θ
sin
θ
P
2
2
(
cos
θ
)
=
3
sin
2
θ
P
3
0
(
cos
θ
)
=
1
2
(
5
cos
3
θ
−
3
cos
θ
)
P
3
1
(
cos
θ
)
=
−
3
2
(
5
cos
2
θ
−
1
)
sin
θ
P
3
2
(
cos
θ
)
=
15
cos
θ
sin
2
θ
P
3
3
(
cos
θ
)
=
−
15
sin
3
θ
P
4
0
(
cos
θ
)
=
1
8
(
35
cos
4
θ
−
30
cos
2
θ
+
3
)
P
4
1
(
cos
θ
)
=
−
5
2
(
7
cos
3
θ
−
3
cos
θ
)
sin
θ
P
4
2
(
cos
θ
)
=
15
2
(
7
cos
2
θ
−
1
)
sin
2
θ
P
4
3
(
cos
θ
)
=
−
105
cos
θ
sin
3
θ
P
4
4
(
cos
θ
)
=
105
sin
4
θ
{\displaystyle {\begin{aligned}P_{0}^{0}(\cos \theta )&=1\\[8pt]P_{1}^{0}(\cos \theta )&=\cos \theta \\[8pt]P_{1}^{1}(\cos \theta )&=-\sin \theta \\[8pt]P_{2}^{0}(\cos \theta )&={\tfrac {1}{2}}(3\cos ^{2}\theta -1)\\[8pt]P_{2}^{1}(\cos \theta )&=-3\cos \theta \sin \theta \\[8pt]P_{2}^{2}(\cos \theta )&=3\sin ^{2}\theta \\[8pt]P_{3}^{0}(\cos \theta )&={\tfrac {1}{2}}(5\cos ^{3}\theta -3\cos \theta )\\[8pt]P_{3}^{1}(\cos \theta )&=-{\tfrac {3}{2}}(5\cos ^{2}\theta -1)\sin \theta \\[8pt]P_{3}^{2}(\cos \theta )&=15\cos \theta \sin ^{2}\theta \\[8pt]P_{3}^{3}(\cos \theta )&=-15\sin ^{3}\theta \\[8pt]P_{4}^{0}(\cos \theta )&={\tfrac {1}{8}}(35\cos ^{4}\theta -30\cos ^{2}\theta +3)\\[8pt]P_{4}^{1}(\cos \theta )&=-{\tfrac {5}{2}}(7\cos ^{3}\theta -3\cos \theta )\sin \theta \\[8pt]P_{4}^{2}(\cos \theta )&={\tfrac {15}{2}}(7\cos ^{2}\theta -1)\sin ^{2}\theta \\[8pt]P_{4}^{3}(\cos \theta )&=-105\cos \theta \sin ^{3}\theta \\[8pt]P_{4}^{4}(\cos \theta )&=105\sin ^{4}\theta \end{aligned}}}
上記の直交関係は、次のように定式化されます。m を 固定すると、 は直交し、 に対して θ でパラメータ化され 、重みは になります 。
P
ℓ
m
(
cos
θ
)
{\displaystyle P_{\ell }^{m}(\cos \theta )}
[
0
,
π
]
{\displaystyle [0,\pi ]}
sin
θ
{\displaystyle \sin \theta }
∫
0
π
P
k
m
(
cos
θ
)
P
ℓ
m
(
cos
θ
)
sin
θ
d
θ
=
2
(
ℓ
+
m
)
!
(
2
ℓ
+
1
)
(
ℓ
−
m
)
!
δ
k
,
ℓ
{\displaystyle \int _{0}^{\pi }P_{k}^{m}(\cos \theta )P_{\ell }^{m}(\cos \theta )\,\sin \theta \,d\theta ={\frac {2(\ell +m)!}{(2\ell +1)(\ell -m)!}}\ \delta _{k,\ell }}
また、 ℓを 固定した場合:
∫
0
π
P
ℓ
m
(
cos
θ
)
P
ℓ
n
(
cos
θ
)
csc
θ
d
θ
=
{
0
if
m
≠
n
(
ℓ
+
m
)
!
m
(
ℓ
−
m
)
!
if
m
=
n
≠
0
∞
if
m
=
n
=
0
{\displaystyle \int _{0}^{\pi }P_{\ell }^{m}(\cos \theta )P_{\ell }^{n}(\cos \theta )\csc \theta \,d\theta ={\begin{cases}0&{\text{if }}m\neq n\\{\frac {(\ell +m)!}{m(\ell -m)!}}&{\text{if }}m=n\neq 0\\\infty &{\text{if }}m=n=0\end{cases}}}
θに関して、 解は
P
ℓ
m
(
cos
θ
)
{\displaystyle P_{\ell }^{m}(\cos \theta )}
d
2
y
d
θ
2
+
cot
θ
d
y
d
θ
+
[
λ
−
m
2
sin
2
θ
]
y
=
0
{\displaystyle {\frac {d^{2}y}{d\theta ^{2}}}+\cot \theta {\frac {dy}{d\theta }}+\left[\lambda -{\frac {m^{2}}{\sin ^{2}\theta }}\right]\,y=0\,}
より正確には、整数 m 0 が与えられた場合、上記の式は、 ℓ に対して
整数 ≥ m の場合にのみ特異でない解を持ち、それらの解は に比例します
。
≥
{\displaystyle \geq }
λ
=
ℓ
(
ℓ
+
1
)
{\displaystyle \lambda =\ell (\ell +1)\,}
P
ℓ
m
(
cos
θ
)
{\displaystyle P_{\ell }^{m}(\cos \theta )}
物理学への応用: 球面調和関数
物理学 では、 球面 対称性 が関係する 角度に関するルジャンドル多項式が頻繁に登場する。 球面座標 の余緯度角は、上で使用した 角度である。経度角 は、乗数として現れる。これらは、 球面調和関数 と呼ばれる関数の集合を構成する。これらの関数は、 リー群 SO(3)の作用下での 2 球面 の対称性を表現する 。 [ 要出典 ]
θ
{\displaystyle \theta }
ϕ
{\displaystyle \phi }
これらの関数が便利なのは、球面上の 方程式を解く上で中心的な役割を果たすからです
。球面座標θ(余緯度)とφ(経度)では、 ラプラシアン は
∇
2
ψ
+
λ
ψ
=
0
{\displaystyle \nabla ^{2}\psi +\lambda \psi =0}
∇
2
ψ
=
∂
2
ψ
∂
θ
2
+
cot
θ
∂
ψ
∂
θ
+
csc
2
θ
∂
2
ψ
∂
ϕ
2
.
{\displaystyle \nabla ^{2}\psi ={\frac {\partial ^{2}\psi }{\partial \theta ^{2}}}+\cot \theta {\frac {\partial \psi }{\partial \theta }}+\csc ^{2}\theta {\frac {\partial ^{2}\psi }{\partial \phi ^{2}}}.}
偏 微分方程式
∂
2
ψ
∂
θ
2
+
cot
θ
∂
ψ
∂
θ
+
csc
2
θ
∂
2
ψ
∂
ϕ
2
+
λ
ψ
=
0
{\displaystyle {\frac {\partial ^{2}\psi }{\partial \theta ^{2}}}+\cot \theta {\frac {\partial \psi }{\partial \theta }}+\csc ^{2}\theta {\frac {\partial ^{2}\psi }{\partial \phi ^{2}}}+\lambda \psi =0}
を変数分離 法で解くと 、φ依存部分 または 整数m≥0に対しては次式が得られる。θ依存部分については次式が
得られる。
sin
(
m
ϕ
)
{\displaystyle \sin(m\phi )}
cos
(
m
ϕ
)
{\displaystyle \cos(m\phi )}
d
2
y
d
θ
2
+
cot
θ
d
y
d
θ
+
[
λ
−
m
2
sin
2
θ
]
y
=
0
{\displaystyle {\frac {d^{2}y}{d\theta ^{2}}}+\cot \theta {\frac {dy}{d\theta }}+\left[\lambda -{\frac {m^{2}}{\sin ^{2}\theta }}\right]\,y=0\,}
解は および
です 。
P
ℓ
m
(
cos
θ
)
{\displaystyle P_{\ell }^{m}(\cos \theta )}
ℓ
≥
m
{\displaystyle \ell {\geq }m}
λ
=
ℓ
(
ℓ
+
1
)
{\displaystyle \lambda =\ell (\ell +1)}
したがって、方程式
∇
2
ψ
+
λ
ψ
=
0
{\displaystyle \nabla ^{2}\psi +\lambda \psi =0}
の場合にのみ非特異分離解を持ち 、それらの解は
λ
=
ℓ
(
ℓ
+
1
)
{\displaystyle \lambda =\ell (\ell +1)}
P
ℓ
m
(
cos
θ
)
cos
(
m
ϕ
)
0
≤
m
≤
ℓ
{\displaystyle P_{\ell }^{m}(\cos \theta )\ \cos(m\phi )\ \ \ \ 0\leq m\leq \ell }
そして
P
ℓ
m
(
cos
θ
)
sin
(
m
ϕ
)
0
<
m
≤
ℓ
.
{\displaystyle P_{\ell }^{m}(\cos \theta )\ \sin(m\phi )\ \ \ \ 0<m\leq \ell .}
ℓ の各選択に対して、 m のさまざまな値と 正弦および余弦の選択に対して 2ℓ + 1 個の関数があります。 球面上で積分すると、
それらはすべて ℓ と m の両方で直交します。
解は通常、複素指数 で表されます 。
Y
ℓ
,
m
(
θ
,
ϕ
)
=
(
2
ℓ
+
1
)
(
ℓ
−
m
)
!
4
π
(
ℓ
+
m
)
!
P
ℓ
m
(
cos
θ
)
e
i
m
ϕ
−
ℓ
≤
m
≤
ℓ
.
{\displaystyle Y_{\ell ,m}(\theta ,\phi )={\sqrt {\frac {(2\ell +1)(\ell -m)!}{4\pi (\ell +m)!}}}\ P_{\ell }^{m}(\cos \theta )\ e^{im\phi }\qquad -\ell \leq m\leq \ell .}
関数は 球面調和関数 であり 、平方根の量は正規化因子である。正と負の m のルジャンドル関数の関係を思い出すと、球面調和関数が恒等式を満たすことが簡単に示される [5]
Y
ℓ
,
m
(
θ
,
ϕ
)
{\displaystyle Y_{\ell ,m}(\theta ,\phi )}
Y
ℓ
,
m
∗
(
θ
,
ϕ
)
=
(
−
1
)
m
Y
ℓ
,
−
m
(
θ
,
ϕ
)
.
{\displaystyle Y_{\ell ,m}^{*}(\theta ,\phi )=(-1)^{m}Y_{\ell ,-m}(\theta ,\phi ).}
球面調和関数は、フーリエ級数 の意味で完全な正規直交関数の集合を形成します 。測地学、地磁気学、スペクトル解析の分野の研究者は、ここで示されているものとは異なる位相係数と正規化係数を使用します ( 球面調和関数 を参照)。
3次元球対称偏微分方程式を球座標の変数分離法で解く場合、放射状部分を除いた残りの部分は典型的には次のような形になる。
∇
2
ψ
(
θ
,
ϕ
)
+
λ
ψ
(
θ
,
ϕ
)
=
0
,
{\displaystyle \nabla ^{2}\psi (\theta ,\phi )+\lambda \psi (\theta ,\phi )=0,}
したがって、解は球面調和関数になります。
一般化
ルジャンドル多項式は 超幾何級数 と密接な関係がある。球面調和関数の形で、ルジャンドル多項式は リー群 SO(3)の作用下での 二球面 の対称性を表現する。SO(3) 以外にも多くのリー群があり、ルジャンドル多項式の類似の一般化は半単純リー群と リーマン対称空間 の対称性を表現するために存在する。大まかに言えば、対称空間上の ラプラシアン を定義することができ、ラプラシアンの固有関数は球面調和関数の他の設定への一般化として考えることができる。
参照
注釈と参考文献
^ クーラント&ヒルベルト 1953、V、§10。
^ アブラモウィッツ、ミルトン 、 ステグン、アイリーン ・アン編 (1983) [1964 年 6 月]。「第 8 章」。 数式、グラフ、および数学表付き数学関数ハンドブック 。応用数学シリーズ。第 55 巻 (訂正を含む第 10 回目の初版 (1972 年 12 月) の追加訂正を含む第 9 回目の再版、初版)。ワシントン DC、ニューヨーク: 米国商務省、国立標準局、ドーバー出版。p. 332。ISBN 978-0-486-61272-0 LCCN 64-60036。MR 0167642。LCCN 65-12253 。
^ ジョン・C・スレーター著 『原子構造の量子理論 』マグロウヒル社(ニューヨーク、1960年)第1巻309ページより。原著はJA・ゴーント著『 ロンドン王立協会哲学論文集 』A228:151(1929年)を引用している。
^ Dong SH, Lemus R., (2002)、「3つのルジャンドル多項式の重なり積分」、Appl. Math. Lett. 15、541-546。
^ この恒等式は、球面調和関数を ウィグナー D 行列 に関連付け、後者の時間反転特性を使用することによっても示されます。 ± m のルジャンドル関数間の関係は、球面調和関数の複素共役恒等式から証明できます。
Arfken, GB; Weber, HJ (2001)、 物理学者のための数学的手法 、Academic Press、 ISBN 978-0-12-059825-0 ; セクション 12.5 (異なる符号規則を使用します。)
Belousov, SL (1962)、 「正規化準ルジャンドル多項式の表」 、数学表、第 18 巻、Pergamon Press 。
コンドン、EU; ショートリー、GH (1970)、 原子スペクトルの理論 、ケンブリッジ、イギリス: ケンブリッジ大学出版局、 OCLC 5388084 ; 第3章。
クーラント、リチャード 、 ヒルベルト、デイヴィッド (1953)、 数理物理学の方法、第1巻 、ニューヨーク:インターサイエンス出版社 。
Dunster, TM (2010)、「Legendre と関連関数」、 Olver, Frank WJ 、Lozier, Daniel M.、Boisvert, Ronald F.、Clark, Charles W. (編)、 NIST Handbook of Mathematical Functions 、Cambridge University Press、 ISBN 978-0-521-19225-5 、 MR 2723248 。
エドモンズ、AR(1957)、 量子力学における角運動量 、プリンストン大学出版、 ISBN 978-0-691-07912-7 ; 第2章。
ヒルデブランド、FB (1976)、 応用のための高度な微積分 、プレンティスホール、 ISBN 978-0-13-011189-0 。
Koornwinder, Tom H.; Wong, Roderick SC; Koekoek, Roelof; Swarttouw, René F. (2010)、「直交多項式」、 Olver, Frank WJ ; Lozier, Daniel M.; Boisvert, Ronald F.; Clark, Charles W. (eds.)、 NIST Handbook of Mathematical Functions 、Cambridge University Press、 ISBN 978-0-521-19225-5 、 MR 2723248 。
Schach, SR (1973) 積分位数と次数のルジャンドル関連関数の新しい恒等式 、応用数学協会数学解析誌、1976 年、第 7 巻、第 1 号、pp. 59–69
外部リンク
MathWorld の準ルジャンドル多項式
MathWorld のルジャンドル多項式
DLMF におけるルジャンドル関数と関連関数