logo

ルジャンドル多項式のロドリゲスの公式 📂関数

ルジャンドル多項式のロドリゲスの公式

公式

ルジャンドル多項式の明示的explicitな公式は以下の通りです。

Pl(x)=12ll!dldxl(x21)l(1) P_{l}(x)=\dfrac{1}{2^{l} l!} \dfrac{d^{l}}{dx^{l}}(x^{2}-1)^{l} \tag{1}

説明

ll番目のルジャンドル多項式を得る公式であり、これをロドリゲスの公式と言います。元々はルジャンドル多項式の明示的な形を示す言葉でしたが、その後、多項式で表される特殊関数の明示的な形を示す公式の一般的な名称となりました。

導出

ルジャンドル多項式PlP_{l}は以下のようなルジャンドルの微分方程式の解を指します。

(1x2)d2ydx22xdydx+l(l+1)y=0 (1 - x^{2}) \dfrac{d^{2} y}{d x^{2}} - 2x \dfrac{d y}{d x} + l(l+1)y = 0

したがって、(1)(1)が上記の微分方程式の解であることを示せば、証明が完了します。

まず、v=(x21)lv=(x^2-1)^lとしたとき、dlvdxl\dfrac{d^lv}{dx^l}がルジャンドル方程式の解であることを示すつもりです。その後、Pl(1)=1P_{l}(1) = 1を満たすように正規化して、(1)(1)を得ます。

dvdx=l(2x)(x21)l1 \dfrac{dv}{dx}=l(2x)(x^2-1)^{l-1}

両辺に(x21)(x^2-1)を掛けると、以下の式を得ます。

(x21)dvdx=2lx(x21)l=2lxv (x^2-1)\dfrac{dv}{dx}=2lx(x^2-1)^l=2lxv

両辺をl+1l+1回微分すると、ライプニッツの法則により以下のようになります。

k=0l+1l+1Ckdl+1kdxl+1k(dvdx)dkdxk(x21)=2lk=0l+1l+1Ckdl+1kvdxl+1kdkxdxk \sum \limits_{k=0}^{l+1} {}_{l+1}\mathrm{C}_{k} \dfrac{ d^{l+1-k}}{dx^{l+1-k} } \left( \dfrac{dv}{dx} \right) \dfrac{d^k}{dx^k} (x^2-1) = 2l\sum \limits_{k=0}^{l+1} {}_{l+1}\mathrm{C} _{k} \dfrac{d^{l+1-k} v}{dx^{l+1-k}} \dfrac{d^k x}{dx^k}

このとき、左辺はk3k \ge 3のときdkdxk(x21)=0\dfrac{d^k}{dx^k}(x^2-1)=0であるため、k=0,2,3k=0,2,3の項のみが残ります。右辺はk2k \ge 2のときdkxdxk=0\dfrac{d^kx}{dx^k}=0であるため、k=1,2k=1,2の項のみが残ります。したがって、次のように得られます。

(x21)dl+2vdxl+2+(l+1)(2x)dl+1vdxl+1+l(l+1)2!2dlvdxl=2lxdl+1vdxl+1+2l(l+1)dlvdxl (x^2-1)\dfrac{d^{l+2} v}{dx^{l+2}} + (l+1)(2x)\dfrac{d^{l+1}v}{dx^{l+1}}+\dfrac{l(l+1)}{2!}2\dfrac{d^l v}{dx^l}=2lx\dfrac{d^{l+1} v}{dx^{l+1}} + 2l(l+1)\dfrac{d^lv}{dx^l}

同じ係数項をまとめて整理すると、以下のようになります。

(1x2)(dlvdxl)2x(dlvdxl)+l(l+1)dlvdxl=0 (1-x^2)\left( \dfrac{d^l v}{dx^l} \right)^{\prime \prime} -2x\left( \dfrac{d^lv}{dx^l} \right)^{\prime} + l(l+1)\dfrac{d^lv}{dx^l}=0

これはルジャンドル方程式と同じ形です。つまり、dlvdxl\dfrac{d^l v}{dx^l}がルジャンドル方程式の解になります。

Pl(x)=dldxl(x21)l P_{l}(x)= \dfrac{d^l}{dx^l}(x^2-1)^l

Pl(1)=1P_{l}(1) = 1を満たす係数を求めてみましょう。(x21)l(x^2-1)^l(x1)l(x+1)l(x-1)^l(x+1)^lで因数分解し、ライプニッツの法則でll回微分すると、以下のようになります。

 Pl(x)=dldxl[(x1)l(x+1)l]=k=0llCkdlkdxlk(x1)ldkdxk(x+1)l=lC0l!(x+1)l+lC1l!(x1)l(x+1)l1+lC2l!2(x1)2l(l1)(x+1)l2+ \begin{align*} &\quad \ P_{l}(x) \\ &= \dfrac{d^l}{dx^l} \left[ (x-1)^l (x+1)^l \right] \\ &= \sum\limits_{k=0}^l {}_{l}\mathrm{C}_{k} \dfrac{d^{l-k}}{dx^{l-k}}(x-1)^l \dfrac{d^k}{dx^k}(x+1)^l \\ &= {}_{l}\mathrm{C}_{0} l! (x+1)^l + {}_{l}\mathrm{C}_{1} l!(x-1) l(x+1)^{l-1}+{}_{l}\mathrm{C}_2\dfrac{l!}{2}(x-1)^2l(l-1)(x+1)^{l-2}+\cdots \end{align*}

2番目の項からは因数として(x1)(x-1)を含むため、x=1x=1のとき00です。したがって、Pl(1)=l!2lP_{l}(1)=l! 2^lであり、この値が11になるためには、12ll!\dfrac{1}{2^l l!}で割ればよいです。したがって、最終的に以下のようなロドリゲスの公式を得ます。

Pl(x)=12ll!dldxl(x21)l P_{l}(x)=\dfrac{1}{2^l l!}\dfrac{d^l}{dx^l}(x^2-1)^l