ルジャンドル多項式

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

定義

解析学においてルジャンドルの微分方程式

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

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

ファイル:Legendrepolynomials6.svg

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

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

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

で表すことができる。

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

[math](x^2-1)\frac{\mathrm d}{\mathrm dx}(x^2-1)^n = 2nx(x^2-1)^n[/math]

の両辺を n + 1 回微分して、高階微分に関する一般ライプニッツ則を適用すればわかる。各ルジャンドル多項式 P_n は以下のテイラー級数 テンプレート:NumBlk の係数として定義することもできる。この母函数は物理学において多重極展開に利用される。

帰納的定義

上記のテンプレート:EquationNoteで与えられたテイラー展開の最初の 2 項から、最初の 2 つのルジャンドル多項式が

[math]P_0(x) = 1, P_1(x) = x[/math]

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

[math] (n+1) P_{n+1}(x) = (2n+1) x P_n(x) - n P_{n-1}(x)[/math]

を用いるのが適当である。この漸化式は、テンプレート:EquationNoteの両辺を t に関して微分したものを整理して得られる等式

[math]\frac{x-t}{\sqrt{1-2xt+t^2}} = (1-2xt+t^2) \sum_{n=1}^\infty n P_n(x) t^{n-1}[/math]

の分母に現れる平方根をテンプレート:EquationNoteで置き換えて、t の冪に対する係数比較を行えば得られる。漸化式に初期条件としてすでに得られている P₀, P₁ を当てはめれば、全てのルジャンドル多項式が帰納的に生成される。

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

[math]\begin{align} 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{align}[/math]

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

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

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

直交性

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

[math]\int_{-1}^{1} P_m(x) P_n(x) ~ \mathrm{d}x = {2 \over {2n + 1}} \delta_{mn}[/math]

ここで δ_mn はクロネッカーのデルタ、即ち m = n のとき 1 で、それ以外のときは 0 を意味する。ここから、ルジャンドル多項式のもう一つの導出法として、同内積に関する直交関数系 {1, x, x²,...} にシュミットの直交化法を適用することによっても導出可能であることが言える。この直交性により、ルジャンドル多項式系がエルミート微分作用素

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

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

物理学における応用

ルジャンドル多項式は初め、1782年にアドリアン＝マリ・ルジャンドルにより、ニュートン・ポテンシャル

[math] \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) [/math]

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

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

[math]\nabla^2 \Phi(\boldsymbol{x})=0[/math]

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

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

となる。A_ℓ と B_ℓ は各問題の境界条件に従って決定される。

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

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

ファイル:Point axial multipole.svg

Figure 2

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

[math]\frac{1}{\sqrt{1 + \eta^{2} - 2\eta x}} = \sum_{k=0}^{\infty} \eta^{k} P_{k}(x)[/math]

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

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

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

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

[math]\Phi(r, \theta) \propto \frac{1}{r} \sum_{k=0}^{\infty} \left( \frac{a}{r} \right)^{k} P_{k}(\cos \theta)[/math]

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

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

その他の性質

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

[math]P_n(-x) = (-1)^n P_n(x)[/math]

を満たす。

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

[math]P_k(1) = 1[/math]

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

[math]P_k'(1) = \frac{k(k+1)}{2}[/math]

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

[math](n+1) P_{n+1}(x) = (2n+1) x P_n(x) - n P_{n-1}(x)[/math]

と、公式

[math]\frac{x^2-1}{n} \frac{\mathrm d}{\mathrm dx} P_n(x) = xP_n(x) - P_{n-1}(x)[/math]

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

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

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

[math]\frac{\mathrm d}{\mathrm dx} 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[/math]

あるいは同じことだが、

[math]{\mathrm d \over \mathrm dx} P_{n+1}(x) = {2 P_n(x) \over \| P_n(x) \|^2} + {2 P_{n-2}(x) \over \| P_{n-2}(x) \|^2} + \dotsb[/math]

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

[math]\| P_n(x) \| = \sqrt{\int _{- 1}^{1}(P_n(x))^2 \,\mathrm dx} = \sqrt{\frac{2}{2 n + 1}}[/math]

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

[math]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[/math]

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

[math]\sum_{j=0}^n P_j(x)\ge 0\qquad (x\ge -1)[/math]

を導く。

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

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

[math]\tilde{P}_n(x) = P_n(2x-1)[/math]

で定義される。ここで、ずらし写像（実はアフィン変換）x ↦ 2x − 1 は、区間 [0, 1] を区間 [−1, 1] へ写す全単射として選ばれたもので、それゆえ多項式系 テンプレート:Tilde_n(x) の区間 [0, 1] 上での直交性

[math]\int_{0}^{1} \tilde{P_m}(x) \tilde{P_n}(x)\,\mathrm dx = {1 \over {2n + 1}} \delta_{mn}[/math]

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

[math]\tilde{P_n}(x) = (-1)^n \sum_{k=0}^n {n \choose k} {n+k \choose k} (-x)^k[/math]

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

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

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

n	[math]\tilde{P_n}(x)[/math]
0	1
1	[math]2x-1[/math]
2	[math]6x^2-6x+1[/math]
3	[math]20x^3-30x^2+12x-1[/math]

{{safesubst:#invoke:Anchor|main}}ルジャンドル陪多項式

"「en:Associated Legendre polynomials」"

非負整数 k、m で

k ≧ m

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

[math]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}[/math]　　

と定義する。P_k^m(t) はルジャンドルの陪微分方程式

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

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

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

[math]P_k{}^m (t) = (1-t^2)^{m/2} \frac{ \mathrm{d}^m P_k(t)}{\mathrm{d}t^m}[/math]

関連項目

• ガウス求積
• ゲーゲンバウアー多項式
• ルジャンドル有理関数
• トゥラーンの不等式
• ルジャンドルのウェーブレット
• ヤコビ多項式
• 球面調和関数

参考文献

• 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 .
• テンプレート:Dlmf
• テンプレート:Dlmf
• 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
• テンプレート:Springer
• 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