ルジャンドル多項式

ルジャンドル多項式（ルジャンドルたこうしき、英: Legendre polynomial）とは、ルジャンドルの微分方程式を満たすルジャンドル関数のうち次数が非負整数のものを言う。直交多項式の一種である。 解析学においてルジャンドルの微分方程式 d d x [ ( 1 − x 2 ) d d x f ( x ) ] + λ ( λ + 1 ) f ( x ) = 0 {\displaystyle {\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 = ±1 に確定特異点（英語版）を持つから、一般には原点の周りでの級数解の収束半径は 1 である。 λ が非負整数 n = 0, 1, 2, … のときの解は x = ±1 の両点においても正則であり、かつ級数は途中で止まって多項式となる。さらに 、x = 1 において値 1 を取るという初期条件を課すと、解は一意に定まる。これを n次のルジャンドル多項式と呼び、普通は Pn(x) と記す。また、全ての非負整数についての n次のルジャンドル多項式全体が成す関数族を総称的にルジャンドル多項式と呼ぶ。ルジャンドル多項式は後述する関数空間の内積に関して直交系を成す。ただし、この内積についての各 Pn(x) の大きさは 1 ではないため (これは Pn(1) = 1 という初期条件を課したためである)、正規直交系にはなっていない点は注意を要する。各ルジャンドル多項式 Pn(x) は n次多項式で、ロドリゲスの公式 P n ( x ) = 1 2 n n ! d n d x n [ ( x 2 − 1 ) n ] {\displaystyle 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 ) d d x ( x 2 − 1 ) n = 2 n x ( x 2 − 1 ) n {\displaystyle (x^{2}-1){\frac {\mathrm {d} }{\mathrm {d} x}}(x^{2}-1)^{n}=2nx(x^{2}-1)^{n}} の両辺を n + 1 回微分して、高階微分に関する一般ライプニッツ則を適用すればわかる。各ルジャンドル多項式 Pn は以下のテイラー級数 1 1 − 2 x t + t 2 = ∑ n = 0 ∞ P n ( x ) t n {\displaystyle {\frac {1}{\sqrt {1-2xt+t^{2}}}}=\sum _{n=0}^{\infty }P_{n}(x)t^{n}} (1) の係数として定義することもできる。この母函数は物理学において多重極展開（英語版）に利用される。 上記の式 (1) で与えられたテイラー展開の最初の 2 項から、最初の 2 つのルジャンドル多項式が P 0 ( x ) = 1 , P 1 ( x ) = x {\displaystyle P_{0}(x)=1,P_{1}(x)=x} となることがわかる。残りの多項式を得るのには、上記のテイラー展開を直截に計算するよりも、ボネの漸化式 ( n + 1 ) P n + 1 ( x ) = ( 2 n + 1 ) x P n ( x ) − n P n − 1 ( x ) {\displaystyle (n+1)P_{n+1}(x)=(2n+1)xP_{n}(x)-nP_{n-1}(x)} を用いるのが適当である。この漸化式は、式 (1) の両辺を t に関して微分したものを整理して得られる等式 x − t 1 − 2 x t + t 2 = ( 1 − 2 x t + t 2 ) ∑ n = 1 ∞ n P n ( x ) t n − 1 {\displaystyle {\frac {x-t}{\sqrt {1-2xt+t^{2}}}}=(1-2xt+t^{2})\sum _{n=1}^{\infty }nP_{n}(x)t^{n-1}} の分母に現れる平方根を式 (1) で置き換えて、t の冪に対する係数比較（英語版）を行えば得られる。漸化式に初期条件としてすでに得られている P0, P1 を当てはめれば、全てのルジャンドル多項式が帰納的に生成される。 漸化式を解いて陽に表せば P n ( x ) = 1 2 n ∑ k = 0 n ( n k ) 2 ( x − 1 ) n − k ( x + 1 ) k = ∑ k = 0 n ( n k ) ( − n − 1 k ) ( 1 − x 2 ) k = 2 n ⋅ ∑ k = 0 n x k ( n k ) ( n + k − 1 2 n ) {\displaystyle {\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}}} などのように書くことができる。後段はルジャンドル多項式を単に単項式として表して二項係数の乗法公式を使えば、漸化式から直ちに得られる。 具体的に最初のいくつかのルジャンドル多項式を挙げれば以下のようになる： ルジャンドル多項式の重要な性質の一つは、これらが閉区間 [−1, 1] 上の L2-内積に関して直交すること、即ち以下の式を満たすことである。 ∫ − 1 1 P m ( x ) P n ( x )   d x = 2 2 n + 1 δ m n {\displaystyle \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, x2,...} にシュミットの直交化法を適用することによっても導出可能であることが言える。この直交性により、ルジャンドル多項式系がエルミート微分作用素 d d x [ ( 1 − x 2 ) d d x P ( x ) ] = − λ P ( x ) {\displaystyle {\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年にアドリアン＝マリ・ルジャンドルにより、ニュートン・ポテンシャル（英語版） 1 | x − x ′ | = 1 r 2 + r ′ 2 − 2 r r ′ cos ⁡ γ = ∑ ℓ = 0 ∞ r ′ ℓ r ℓ + 1 P ℓ ( cos ⁡ γ ) {\displaystyle {\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′ が満たされる場合に収束し、質点に対応する重力ポテンシャルもしくは点電荷に対応するクーロンポテンシャルを極座標表示する際に用いることができる。このルジャンドル多項式を用いた展開は、例えば連続質量や電荷分布の上でこの展開を積分するときなどに有用である。 ルジャンドル多項式は、空間の無電荷領域における電位に関するラプラス方程式 ∇ 2 Φ ( x ) = 0 {\displaystyle \nabla ^{2}\Phi ({\boldsymbol {x}})=0} を軸対称な（方位角に依存しない）境界条件のもとで、変数分離法を用いて解く際にも登場する。ここで、 ^z を対称軸、θ を観測者の位置と ^z-軸との間の角（天頂角）とするとき、電位は Φ ( r , θ ) = ∑ ℓ = 0 ∞ [ A ℓ r ℓ + B ℓ r − ( ℓ + 1 ) ] P ℓ ( cos ⁡ θ ) {\displaystyle \Phi (r,\theta )=\sum _{\ell =0}^{\infty }\left[A_{\ell }r^{\ell }+B_{\ell }r^{-(\ell +1)}\right]P_{\ell }(\cos \theta )} となる。Aℓ と Bℓ は各問題の境界条件に従って決定される。 三次元における中心力に対するシュレーディンガー方程式を解く際にもルジャンドル多項式は現れる。 ルジャンドル多項式は、多重極展開で自然に現れる 1 1 + η 2 − 2 η x = ∑ k = 0 ∞ η k P k ( x ) {\displaystyle {\frac {1}{\sqrt {1+\eta ^{2}-2\eta x}}}=\sum _{k=0}^{\infty }\eta ^{k}P_{k}(x)} なる形の関数（記号を少し変えてあるが、上で述べたものと同じ）の展開においても有用である。等式の左辺はルジャンドル多項式の母関数の閉じた形である。 例として、（球座標系での）電位 Φ(r, θ) が z-軸上の点 z = a にある点電荷によるものとすれば、 Φ ( r , θ ) ∝ 1 R = 1 r 2 + a 2 − 2 a r cos ⁡ θ {\displaystyle \Phi (r,\theta )\propto {\frac {1}{R}}={\frac {1}{\sqrt {r^{2}+a^{2}-2ar\cos \theta }}}} と書くことができる。観測点 P の半径 r が a より大きければ、電位はルジャンドル多項式を用いて Φ ( r , θ ) ∝ 1 r ∑ k = 0 ∞ ( a r ) k P k ( cos ⁡ θ ) {\displaystyle \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 ) {\displaystyle P_{n}(-x)=(-1)^{n}P_{n}(x)} を満たす。 微分方程式と直交性はスケール変換に依らない性質だから、ルジャンドル多項式はその定義において適当に定数倍して P k ( 1 ) = 1 {\displaystyle P_{k}(1)=1} を満たすように「標準化」される（「正規化」とも言うが、実際にノルムが 1 というわけではないので紛らわしい）。端点における微分係数は P k ′ ( 1 ) = k ( k + 1 ) 2 {\displaystyle P_{k}'(1)={\frac {k(k+1)}{2}}} で与えられる。既に述べたとおり、ルジャンドル多項式はボネの漸化式と呼ばれる三項間漸化式 ( n + 1 ) P n + 1 ( x ) = ( 2 n + 1 ) x P n ( x ) − n P n − 1 ( x ) {\displaystyle (n+1)P_{n+1}(x)=(2n+1)xP_{n}(x)-nP_{n-1}(x)} と、公式 x 2 − 1 n d d x P n ( x ) = x P n ( x ) − P n − 1 ( x ) {\displaystyle {\frac {x^{2}-1}{n}}{\frac {\mathrm {d} }{\mathrm {d} x}}P_{n}(x)=xP_{n}(x)-P_{n-1}(x)} に従うが、これらから得られる等式 ( 2 n + 1 ) P n ( x ) = d d x [ P n + 1 ( x ) − P n − 1 ( x ) ] {\displaystyle (2n+1)P_{n}(x)={\frac {\mathrm {d} }{\mathrm {d} x}}\left[P_{n+1}(x)-P_{n-1}(x)\right]} ルジャンドル多項式の積分に有効である。これを繰り返し用いて d d x P n + 1 ( x ) = ( 2 n + 1 ) P n ( x ) + ( 2 ( n − 2 ) + 1 ) P n − 2 ( x ) + ( 2 ( n − 4 ) + 1 ) P n − 4 ( x ) + ⋯ {\displaystyle {\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 } あるいは同じことだが、 d d x P n + 1 ( x ) = 2 P n ( x ) ‖ P n ( x ) ‖ 2 + 2 P n − 2 ( x ) ‖ P n − 2 ( x ) ‖ 2 + ⋯ {\displaystyle {\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 } が得られる。ただし、ǁPn(x)ǁ は閉区間 [−1, 1] 上のノルム ‖ P n ( x ) ‖ = ∫ − 1 1 ( P n ( x ) ) 2 d x = 2 2 n + 1 {\displaystyle \|P_{n}(x)\|={\sqrt {\int _{-1}^{1}(P_{n}(x))^{2}\,\mathrm {d} x}}={\sqrt {\frac {2}{2n+1}}}} である。ボネの漸化式から帰納的に陽な表現 P n ( x ) = ∑ k = 0 n ( − 1 ) k ( n k ) 2 ( 1 + x 2 ) n − k ( 1 − x 2 ) k {\displaystyle 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}} が得られる。ルジャンドル多項式に対するアスキー-ギャスパーの不等式（英語版）は ∑ j = 0 n P j ( x ) ≥ 0 ( x ≥ − 1 ) {\displaystyle \sum _{j=0}^{n}P_{j}(x)\geq 0\qquad (x\geq -1)} を導く。 ずらしルジャンドル多項式 (shifted Legendre polynomial) は P ~ n ( x ) = P n ( 2 x − 1 ) {\displaystyle {\tilde {P}}_{n}(x)=P_{n}(2x-1)} で定義される。ここで、ずらし写像（実はアフィン変換）x ↦ 2x − 1 は、区間 [0, 1] を区間 [−1, 1] へ写す全単射として選ばれたもので、それゆえ多項式系 ~Pn(x) の区間 [0, 1] 上での直交性 ∫ 0 1 P m ~ ( x ) P n ~ ( x ) d x = 1 2 n + 1 δ m n {\displaystyle \int _{0}^{1}{\tilde {P_{m}}}(x){\tilde {P_{n}}}(x)\,\mathrm {d} x={1 \over {2n+1}}\delta _{mn}} が従う。ずらしルジャンドル多項式の明示式は P n ~ ( x ) = ( − 1 ) n ∑ k = 0 n ( n k ) ( n + k k ) ( − x ) k {\displaystyle {\tilde {P_{n}}}(x)=(-1)^{n}\sum _{k=0}^{n}{n \choose k}{n+k \choose k}(-x)^{k}} で与えられる。ロドリゲスの公式のずらしルジャンドル多項式版は P n ~ ( x ) = 1 n ! d n d x n [ ( x 2 − x ) n ] {\displaystyle {\tilde {P_{n}}}(x)={\frac {1}{n!}}{\mathrm {d} ^{n} \over \mathrm {d} x^{n}}\left[(x^{2}-x)^{n}\right]} となる。ずらしルジャンドル多項式の最初の方のいくつかは以下のようになる。 非負整数 k、m で k ≧ m を満たすものに対し、ルジャンドル陪多項式Pkm(t) を P k m ( t ) = 1 2 k ( 1 − t 2 ) m / 2 ∑ j = 1 ⌊ ( k − m ) / 2 ⌋ ( − 1 ) j ( 2 k − 2 j ) ! j ! ( k − j ) ! ( k − 2 j − m ) ! t k − 2 j − m {\displaystyle P_{k}{}^{m}(t)={\frac {1}{2^{k}}}(1-t^{2})^{m/2}\sum _{j=1}^{\lfloor (k-m)/2\rfloor }{(-1)^{j}(2k-2j)! \over j!(k-j)!(k-2j-m)!}t^{k-2j-m}} 　　 と定義する。Pkm(t) はルジャンドルの陪微分方程式 ( 1 − t 2 ) y ″ ( t ) − 2 t y ′ ( t ) + ( k ( k + 1 ) − m 2 1 − t 2 ) y ( t ) = 0 {\displaystyle (1-t^{2})y''(t)-2ty'(t)+\left(k(k+1)-{m^{2} \over 1-t^{2}}\right)y(t)=0} の解である。なお、ルジャンドルの陪微分方程式は k ≧ m を満たすときのみ解を持つことが知られている。また、Ykm (θ, φ) の定義における係数は、後述するノルムが 1 になるよう選んだものである。 Pkm(t) とルジャンドル多項式 Pk(t) は以下の関係を満たす： P k m ( t ) = ( 1 − t 2 ) m / 2 d m P k ( t ) d t m {\displaystyle 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, https://books.google.fr/books?id=qLFo_Z-PoGIC&printsec=frontcover&hl=ja&source=gbs_ge_summary_r&cad=0#v=onepage&q&f=false  ^ 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, https://books.google.fr/books?id=qLFo_Z-PoGIC&printsec=frontcover&hl=ja&source=gbs_ge_summary_r&cad=0#v=onepage&q&f=false  ^ 日本測地学会 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日閲覧。 A quick informal derivation of the Legendre polynomial in the context of the quantum mechanics of hydrogen Hazewinkel, Michiel, ed. (2001), “Legendre polynomials”, Encyclopaedia of Mathematics, Springer, ISBN 978-1-55608-010-4, http://eom.springer.de/p/l058050.htm  Wolfram MathWorld entry on Legendre polynomials Module for Legendre Polynomials by John H. Mathews Dr James B. Calvert's article on Legendre polynomials from his personal collection of mathematics The Legendre Polynomials by Carlyle E. Moore Legendre Polynomials from Hyperphysics