ルジャンドル多項式

ルジャンドル多項式（ルジャンドルたこうしき、英: Legendre polynomial）とは、ルジャンドルの微分方程式を満たすルジャンドル関数のうち次数が非負整数のものを言う。直交多項式の一種である。

定義

{\mathrm {d}  \over \mathrm {d} x}\left[(1-x^{2}){\mathrm {d}  \over \mathrm {d} x}f(x)\right]+\lambda (\lambda +1)f(x)=0

（ $λ$ は任意の複素数とする）は標準的な冪級数法を用いて解けることが知られており、その解は一般にルジャンドル関数と呼ばれる（何れもアドリアン＝マリ・ルジャンドルに名を因む）。この方程式は $x = \pm1$ に確定特異点（英語版）を持つから、一般には原点の周りでの級数解の収束半径は 1 である。

n = 5 までのルジャンドル多項式のグラフ

$λ$ が非負整数 $n = 0, 1, 2, \dots$ のときの解は $x = \pm1$ の両点においても正則であり、かつ級数は途中で止まって多項式となる。さらに、 $x = 1$ において値 1 を取るという初期条件を課すと、解は一意に定まる。これを $n$ 次のルジャンドル多項式と呼び、普通は $P n (x)$ と記す^[1]。また、全ての非負整数についての $n$ 次のルジャンドル多項式全体が成す関数族を総称的にルジャンドル多項式と呼ぶ。ルジャンドル多項式は後述する関数空間の内積に関して直交系を成す。ただし、この内積についての各 $P n (x)$ の大きさは 1 ではないため (これは $P n (1) = 1$ という初期条件を課したためである)、正規直交系にはなっていない点は注意を要する。各ルジャンドル多項式 $P n (x)$ は $n$ 次多項式で、ロドリゲスの公式

P_{n}(x)={1 \over 2^{n}n!}{\mathrm {d} ^{n} \over \mathrm {d} x^{n}}\left[(x^{2}-1)^{n}\right]

で表すことができる。

ルジャンドル多項式がルジャンドルの微分方程式を満たすことは、恒等式

(x^{2}-1){\frac {\mathrm {d} }{\mathrm {d} x}}(x^{2}-1)^{n}=2nx(x^{2}-1)^{n}

の両辺を $n + 1$ 回微分して、高階微分に関する一般ライプニッツ則を適用すればわかる^[2]。各ルジャンドル多項式 $P n$ は以下のテイラー級数

{\frac {1}{\sqrt {1-2xt+t^{2}}}}=\sum _{n=0}^{\infty }P_{n}(x)t^{n}

(1)

の係数として定義することもできる^[3]。この母函数は物理学において多重極展開（英語版）に利用される。

帰納的定義

上記の式 (1) で与えられたテイラー展開の最初の 2 項から、最初の 2 つのルジャンドル多項式が

P_{0}(x)=1,P_{1}(x)=x

となることがわかる。残りの多項式を得るのには、上記のテイラー展開を直截に計算するよりも、ボネの漸化式

(n+1)P_{n+1}(x)=(2n+1)xP_{n}(x)-nP_{n-1}(x)

を用いるのが適当である。この漸化式は、式 (1) の両辺を $t$ に関して微分したものを整理して得られる等式

{\frac {x-t}{\sqrt {1-2xt+t^{2}}}}=(1-2xt+t^{2})\sum _{n=1}^{\infty }nP_{n}(x)t^{n-1}

の分母に現れる平方根を式 (1) で置き換えて、 $t$ の冪に対する係数比較（英語版）を行えば得られる。漸化式に初期条件としてすでに得られている $P 0, P 1$ を当てはめれば、全てのルジャンドル多項式が帰納的に生成される。

漸化式を解いて陽に表せば

{\begin{aligned}P_{n}(x)&={\frac {1}{2^{n}}}\sum _{k=0}^{n}{n \choose k}^{2}(x-1)^{n-k}(x+1)^{k}\\&=\sum _{k=0}^{n}{n \choose k}{-n-1 \choose k}\left({\frac {1-x}{2}}\right)^{k}\\&=2^{n}\cdot \sum _{k=0}^{n}x^{k}{n \choose k}{{\frac {n+k-1}{2}} \choose n}\end{aligned}}

などのように書くことができる。後段はルジャンドル多項式を単に単項式として表して二項係数の乗法公式を使えば、漸化式から直ちに得られる。

具体的に最初のいくつかのルジャンドル多項式を挙げれば以下のようになる：

n	$P_{n}(x)$
0	$1$
1	$x$
2	$\textstyle {\frac {1}{2}}(3x^{2}-1)$
3	$\textstyle {\frac {1}{2}}(5x^{3}-3x)$
4	$\textstyle {\frac {1}{8}}(35x^{4}-30x^{2}+3)$
5	$\textstyle {\frac {1}{8}}(63x^{5}-70x^{3}+15x)$
6	$\textstyle {\frac {1}{16}}(231x^{6}-315x^{4}+105x^{2}-5)$
7	$\textstyle {\frac {1}{16}}(429x^{7}-693x^{5}+315x^{3}-35x)$
8	$\textstyle {\frac {1}{128}}(6435x^{8}-12012x^{6}+6930x^{4}-1260x^{2}+35)$
9	$\textstyle {\frac {1}{128}}(12155x^{9}-25740x^{7}+18018x^{5}-4620x^{3}+315x)$
10	$\textstyle {\frac {1}{256}}(46189x^{10}-109395x^{8}+90090x^{6}-30030x^{4}+3465x^{2}-63)$

直交性

ルジャンドル多項式の重要な性質の一つは、これらが閉区間 $[-1, 1]$ 上の L²-内積に関して直交すること、即ち以下の式を満たすことである。

\int _{-1}^{1}P_{m}(x)P_{n}(x)~\mathrm {d} x={2 \over {2n+1}}\delta _{mn}

ここで $δ mn$ はクロネッカーのデルタ、即ち $m = n$ のとき 1 で、それ以外のときは 0 である。すなわち、関数系 {1, x, x²,...} にシュミットの直交化法を適用することによってルジャンドル多項式を導出法とすることが可能である。この直交性により、ルジャンドル多項式系がエルミート微分作用素

{\mathrm {d}  \over \mathrm {d} x}\left[(1-x^{2}){\mathrm {d}  \over \mathrm {d} x}P(x)\right]=-\lambda P(x)

の固有値 $λ = n (n + 1)$ に属する固有関数系となるようなスツルム・リウヴィル理論としてルジャンドルの微分方程式を捉えることができる。

物理学における応用

ルジャンドル多項式は初め、1782年にアドリアン＝マリ・ルジャンドル^[4]により、ニュートン・ポテンシャル（英語版）

{\frac {1}{\left|{\boldsymbol {x}}-{\boldsymbol {x}}^{\prime }\right|}}={\frac {1}{\sqrt {r^{2}+r^{\prime 2}-2rr'\cos \gamma }}}=\sum _{\ell =0}^{\infty }{\frac {r^{\prime \ell }}{r^{\ell +1}}}P_{\ell }(\cos \gamma )

の展開の係数として定義された。ここに、 $r, r'$ はそれぞれベクトル $x, x'$ の長さであり、 $γ$ はそれらのベクトルのなす角である。上記の級数は $r > r'$ が満たされる場合に収束し、質点に対応する重力ポテンシャルもしくは点電荷に対応するクーロンポテンシャルを極座標表示する際に用いることができる。このルジャンドル多項式を用いた展開は、例えば連続質量や電荷分布の上でこの展開を積分するときなどに有用である。

ルジャンドル多項式は、空間の無電荷領域における電位に関するラプラス方程式

\nabla ^{2}\Phi ({\boldsymbol {x}})=0

を軸対称な（方位角に依存しない）境界条件のもとで、変数分離法を用いて解く際にも登場する。ここで、 $ˆ z$ を対称軸、 $θ$ を観測者の位置と $ˆ z$ -軸との間の角（天頂角）とするとき、電位は

\Phi (r,\theta )=\sum _{\ell =0}^{\infty }\left[A_{\ell }r^{\ell }+B_{\ell }r^{-(\ell +1)}\right]P_{\ell }(\cos \theta )

となる。 $A ℓ$ と $B ℓ$ は各問題の境界条件に従って決定される^[5]。

三次元における中心力に対するシュレーディンガー方程式を解く際にもルジャンドル多項式は現れる。

多重極展開におけるルジャンドル多項式

Figure 2

ルジャンドル多項式は、多重極展開で自然に現れる

{\frac {1}{\sqrt {1+\eta ^{2}-2\eta x}}}=\sum _{k=0}^{\infty }\eta ^{k}P_{k}(x)

なる形の関数（記号を少し変えてあるが、上で述べたものと同じ）の展開においても有用である。等式の左辺はルジャンドル多項式の母関数の閉じた形である。

例として、（球座標系での）電位 $Φ(r, θ)$ が $z$ -軸上の点 $z = a$ にある点電荷によるものとすれば、

\Phi (r,\theta )\propto {\frac {1}{R}}={\frac {1}{\sqrt {r^{2}+a^{2}-2ar\cos \theta }}}

と書くことができる。観測点 $P$ の半径 $r$ が $a$ より大きければ、電位はルジャンドル多項式を用いて

\Phi (r,\theta )\propto {\frac {1}{r}}\sum _{k=0}^{\infty }\left({\frac {a}{r}}\right)^{k}P_{k}(\cos \theta )

と展開することができる。ここでは $η = a / r < 1$ および $x = cos θ$ と置いた。この展開は通常の多重極展開を行うのに用いられる。

逆に、観測点 $P$ の半径 $r$ が $a$ より小さいならば、電位を上記のようにルジャンドル多項式展開することはできるが、 $a$ と $r$ とは入れ替わる。この展開は内部多重極展開 (interior multipole expansion) の基本となる。

その他の性質

ルジャンドル多項式は対称または反対称、即ち

P_{n}(-x)=(-1)^{n}P_{n}(x)

を満たす^[6]。

微分方程式と直交性はスケール変換に依らない性質だから、ルジャンドル多項式はその定義において適当に定数倍して

P_{k}(1)=1

を満たすように「標準化」される（「正規化」とも言うが、実際にノルムが 1 というわけではないので紛らわしい）。端点における微分係数は

P_{k}'(1)={\frac {k(k+1)}{2}}

で与えられる。既に述べたとおり、ルジャンドル多項式はボネの漸化式と呼ばれる三項間漸化式

(n+1)P_{n+1}(x)=(2n+1)xP_{n}(x)-nP_{n-1}(x)

と、公式

{\frac {x^{2}-1}{n}}{\frac {\mathrm {d} }{\mathrm {d} x}}P_{n}(x)=xP_{n}(x)-P_{n-1}(x)

に従うが、これらから得られる等式

(2n+1)P_{n}(x)={\frac {\mathrm {d} }{\mathrm {d} x}}\left[P_{n+1}(x)-P_{n-1}(x)\right]

ルジャンドル多項式の積分に有効である。これを繰り返し用いて

{\frac {\mathrm {d} }{\mathrm {d} x}}P_{n+1}(x)=(2n+1)P_{n}(x)+(2(n-2)+1)P_{n-2}(x)+(2(n-4)+1)P_{n-4}(x)+\dotsb

あるいは同じことだが、

{\mathrm {d}  \over \mathrm {d} x}P_{n+1}(x)={2P_{n}(x) \over \|P_{n}(x)\|^{2}}+{2P_{n-2}(x) \over \|P_{n-2}(x)\|^{2}}+\dotsb

が得られる。ただし、ǁP_n(x)ǁ は閉区間 [−1, 1] 上のノルム

\|P_{n}(x)\|={\sqrt {\int _{-1}^{1}(P_{n}(x))^{2}\,\mathrm {d} x}}={\sqrt {\frac {2}{2n+1}}}

である。ボネの漸化式から帰納的に陽な表現

P_{n}(x)=\sum _{k=0}^{n}(-1)^{k}{\binom {n}{k}}^{2}\left({\frac {1+x}{2}}\right)^{n-k}\left({\frac {1-x}{2}}\right)^{k}

が得られる。ルジャンドル多項式に対するアスキー-ギャスパーの不等式（英語版）は

\sum _{j=0}^{n}P_{j}(x)\geq 0\qquad (x\geq -1)

を導く。

ずらしルジャンドル多項式

ずらしルジャンドル多項式 (shifted Legendre polynomial) は

{\tilde {P}}_{n}(x)=P_{n}(2x-1)

で定義される。ここで、ずらし写像（実はアフィン変換） $x \mapsto 2 x - 1$ は、区間 $[0, 1]$ を区間 $[-1, 1]$ へ写す全単射として選ばれたもので、それゆえ多項式系 $~ P n (x)$ の区間 $[0, 1]$ 上での直交性

\int _{0}^{1}{\tilde {P_{m}}}(x){\tilde {P_{n}}}(x)\,\mathrm {d} x={1 \over {2n+1}}\delta _{mn}

が従う。ずらしルジャンドル多項式の明示式は

{\tilde {P_{n}}}(x)=(-1)^{n}\sum _{k=0}^{n}{n \choose k}{n+k \choose k}(-x)^{k}

で与えられる。ロドリゲスの公式のずらしルジャンドル多項式版は

{\tilde {P_{n}}}(x)={\frac {1}{n!}}{\mathrm {d} ^{n} \over \mathrm {d} x^{n}}\left[(x^{2}-x)^{n}\right]

となる。ずらしルジャンドル多項式の最初の方のいくつかは以下のようになる。

n	${\tilde {P_{n}}}(x)$
0	1
1	$2x-1$
2	$6x^{2}-6x+1$
3	$20x^{3}-30x^{2}+12x-1$

ルジャンドル陪多項式

→詳細は「en:Associated Legendre polynomials」を参照

非負整数 $k$ 、 $m$ で

k ≧ m

を満たすものに対し、ルジャンドル陪多項式 $P k m (t)$ を

P_{k}{}^{m}(t)={\frac {1}{2^{k}}}(1-t^{2})^{m/2}\sum _{j=0}^{\lfloor (k-m)/2\rfloor }{(-1)^{j}(2k-2j)! \over j!(k-j)!(k-2j-m)!}t^{k-2j-m}

と定義する^[7]。 $P k m (t)$ はルジャンドルの陪微分方程式

(1-t^{2})y''(t)-2ty'(t)+\left(k(k+1)-{m^{2} \over 1-t^{2}}\right)y(t)=0

の解である。なお、ルジャンドルの陪微分方程式は $k ≧ m$ を満たすときのみ解を持つことが知られている。また、 $Y k m (θ, φ)$ の定義における係数は、後述するノルムが 1 になるよう選んだものである。

$P k m (t)$ とルジャンドル多項式 $P k (t)$ は以下の関係を満たす：

P_{k}{}^{m}(t)=(1-t^{2})^{m/2}{\frac {\mathrm {d} ^{m}P_{k}(t)}{\mathrm {d} t^{m}}}

脚注

[脚注の使い方]

出典

^ 永宮健夫『応用微分方程式論』、共立出版社、1967年、pp46-52。
^ Courant & Hilbert 1953, II, §8
^ George B. Arfken, Hans J. Weber (2005), Mathematical Methods for Physicists, Elsevier Academic Press, p. 743, ISBN 0-12-059876-0
^ M. Le Gendre, “Recherches sur l'attraction des sphéroïdes homogènes”, Mémoires de Mathématiques et de Physique, présentés à l'Académie Royale des Sciences, par divers savans, et lus dans ses Assemblées, Tome X, pp. 411-435 (Paris, 1785). [注: ルジャンドルは彼の発見を1782年に科学アカデミーに提出したが、出版されたのは1785年であった。]
^ Jackson, J.D. Classical Electrodynamics, 3rd edition, Wiley & Sons, 1999. page 103
^ George B. Arfken, Hans J. Weber (2005), Mathematical Methods for Physicists, Elsevier Academic Press, p. 753, ISBN 0-12-059876-0
^ 日本測地学会 2004

参考文献

Abramowitz, Milton; Stegun, Irene A. (1965), Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover, ISBN 978-0486612720 .
Bayin, S.S. (2006), Mathematical Methods in Science and Engineering, Wiley , Chapter 2.
Belousov, S. L. (1962), Tables of normalized associated Legendre polynomials, Mathematical tables, 18, Pergamon Press .
Courant, Richard; Hilbert, David (1953), Methods of Mathematical Physics, Volume 1, New York: Interscience Publischer, Inc .
Dunster, T. M. (2010), “Legendre and Related Functions”, in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0521192255, http://dlmf.nist.gov/14
Koornwinder, Tom H.; Wong, Roderick S. C.; Koekoek, Roelof; Swarttouw, René F. (2010), “Orthogonal Polynomials”, in Olver, Frank W. J.; Lozier, Daniel M.; Boisvert, Ronald F. et al., NIST Handbook of Mathematical Functions, Cambridge University Press, ISBN 978-0521192255, http://dlmf.nist.gov/18
Refaat El Attar (2009), Legendre Polynomials and Functions, CreateSpace, ISBN 978-1-4414-9012-4
高知大学自然科学系田部井隆雄、神奈川県温泉地学研究所里村幹夫、京都大学大学院理学研究科福田洋一 (2004年). “4-4. ルジャンドルの多項式, 陪多項式”. 日本測地学会. 2017年1月4日閲覧。

外部リンク

[nag-1] 永宮健夫『応用微分方程式論』、共立出版社、1967年、pp46-52。

[2] Courant & Hilbert 1953, II, §8

[arfken-743-3] George B. Arfken, Hans J. Weber (2005), Mathematical Methods for Physicists, Elsevier Academic Press, p. 743, ISBN 0-12-059876-0

[4] M. Le Gendre, “Recherches sur l'attraction des sphéroïdes homogènes”, Mémoires de Mathématiques et de Physique, présentés à l'Académie Royale des Sciences, par divers savans, et lus dans ses Assemblées, Tome X, pp. 411-435 (Paris, 1785). [注: ルジャンドルは彼の発見を1782年に科学アカデミーに提出したが、出版されたのは1785年であった。]

[5] Jackson, J.D. Classical Electrodynamics, 3rd edition, Wiley & Sons, 1999. page 103

[arfken-753-6] George B. Arfken, Hans J. Weber (2005), Mathematical Methods for Physicists, Elsevier Academic Press, p. 753, ISBN 0-12-059876-0

[7] 日本測地学会 2004

[1]

[2]

[3]

[4]

[5]

[6]

[7]