logo

ラゲール多項式のロドリゲスの公式 📂関数

ラゲール多項式のロドリゲスの公式

数式

ラゲール多項式の明示的な式は以下のようである。

Ln(x)=1n!exdndxn(xnex)(1) L_{n}(x) = \frac{1}{n!}e^{x}\frac{ d ^{n}}{ dx^{n} }(x^{n}e^{-x}) \tag{1}

説明

上の式はラゲール多項式に対するロドリゲスの公式と言われる。ロドリゲスの公式という用語はもともとルジャンドル多項式の明示的な形を示したが、後に多項式で表される特別な関数の明示的な形を示す公式の一般的な名称になった。最初の数個の多項式を書き出すと以下のようになる。

L0(x)=1L1(x)=x+1L2(x)=12(x24x+2)L3(x)=16(x3+9x218x+6) \begin{align*} L_{0}(x) &= 1 \\ L_{1}(x) &= -x+1 \\ L_{2}(x) &= \frac{1}{2}\left( x^{2}-4x+2 \right) \\ L_{3}(x) &=\frac{1}{6}\left( -x^{3}+9x^{2}-18x+6 \right) \\ & \vdots \end{align*}

証明1

戦略: (1)(1)ラゲール微分方程式を満たすことを示せばよい。

便宜上微分演算子 D=ddxD=\frac{d}{dx}を導入する。Df=dfdx=fDf=\frac{ d f}{ d x }=f^{\prime}で、証明で適切に両方の表現を使用する。

ロドリゲスの公式をf(x)=exn!dndxn(xnex)f(x)=\frac{e^{x}}{n!}\frac{ d ^{n}}{ d x^{n} }(x^{n}e^{-x})とし、v=xnexv=x^{n}e^{-x}とする。まず、xv=(nx)vxv^{\prime}=(n-x)vになることを示す。 v=nxn1exxnex    xv=nxnexxxnex=(nx)v \begin{align*} && v^{\prime}&=nx^{n-1}e^{-x}-x^{n}e^{-x} \\ \implies && xv^{\prime}&=nx^{n}e^{-x}-xx^{n}e^{-x} \\ && &=(n-x)v \end{align*}

今、両側をn+1n+1回微分しよう。まず左側を微分すると、ライプニッツの法則により以下のようになる。

Dn+1(xv)=k=0n(n+1)!(n+1k)!k!(Dkx)(Dn+1kv)=k=01(n+1)!(n+1k)!k!(Dkx)(Dn+1kv) \begin{align*} D^{n+1}(xv^{\prime}) &= \sum \limits _{k=0}^{n} \frac{(n+1)!}{(n+1-k)!k!}(D^{k}x)(D^{n+1-k}v^{\prime}) \\ &= \sum \limits _{k=0}^{1} \frac{(n+1)!}{(n+1-k)!k!}(D^{k}x)(D^{n+1-k}v^{\prime}) \end{align*}

2番目の等号はk2k \ge 2の時Dkx=0D^{k}x=0であるため成立する。シグマ記号を解くと、以下の結果を得る。

Dn+1(xv)=xDn+1v+(n+1)Dnv=x(Dnv)+(n+1)(Dnv) \begin{align*} D^{n+1}(xv^{\prime}) &= xD^{n+1}v^{\prime}+(n+1)D^{n}v^{\prime} \\ &= x(D^{n}v)^{\prime \prime} +(n+1)(D^{n}v)^{\prime} \end{align*}

同様に右側は以下のようになる。

Dn+1[(nx)v]=k=0n(n+1)!(n+1k)!k![Dk(nx)](Dn+1kv)=k=01(n+1)!(n+1k)!k![Dk(nx)](Dn+1kv)=(nx)Dn+1v+(n+1)(1)Dnv=(nx)(Dnv)(n+1)Dnv \begin{align*} D^{n+1}\left[ (n-x)v \right] &= \sum \limits _{k=0}^{n}\frac{(n+1)!}{(n+1-k)!k!}\left[D^{k}(n-x)\right] (D^{n+1-k}v) \\ &=\sum \limits _{k=0}^{1}\frac{(n+1)!}{(n+1-k)!k!}\left[D^{k}(n-x)\right] (D^{n+1-k}v) \\ &=(n-x)D^{n+1}v+(n+1)(-1)D^{n}v \\ &= (n-x)(D^{n}v)^{\prime}-(n+1)D^{n}v \end{align*}

したがって、以下を得る。

Dn+1(xv)=Dn+1[(nx)v]    x(Dnv)+(n+1)(Dnv)=(nx)(Dnv)(n+1)Dnv \begin{align*} && D^{n+1}(xv^{\prime}) &= D^{n+1}[(n-x)v] \\ \implies && x(D^{n}v)^{\prime \prime} +(n+1)(D^{n}v)^{\prime} &= (n-x)(D^{n}v)^{\prime}-(n+1)D^{n}v \end{align*}

これを整理すると、以下のようになる。

x(Dnv)+(1+x)(Dnv)+(n+1)Dnv=0 x(D^{n}v)^{\prime \prime} +(1+x)(D^{n}v)^{\prime}+(n+1)D^{n}v=0

しかし、仮定により以下が成り立つ。

Dnv=dndxn(xnex)=n!exf(x) D^{n}v=\frac{ d ^{n}}{ dx^{n} }(x^{n}e^{-x})=n!e^{-x}f(x)

これを上の式に代入すると、以下の結果を得る。

x[n!exf(x)]+(1+x)[n!exf(x)]+(n+1)[n!exf(x)]=0 x[n!e^{-x}f(x)]^{\prime \prime}+(1+x)[n!e^{-x}f(x)]^{\prime}+(n+1)[n!e^{-x}f(x)]=0

微分を全て展開すると、以下のようになる。計算は簡単なので、過程は省略する。

n!ex[xf(x)+(1x)f(x)+nf(x)]=0 n!e^{-x} \left[ xf^{\prime \prime}(x)+(1-x)f^{\prime}(x)+nf(x) \right]=0

n!ex0n!e^{-x} \ne 0であるため、括弧内の式は00でなければならない。

xf(x)+(1x)f(x)+nf(x)=0 xf^{\prime \prime}(x)+(1-x)f^{\prime}(x)+nf(x)=0

これはラゲール微分方程式であるため、f(x)f(x)はラゲール微分方程式の解、ラゲール多項式である。

f(x)=exn!dndxn(xnex)=Ln(x) f(x)=\frac{e^{x}}{n!}\frac{ d ^{n}}{ d x^{n} }(x^{n}e^{-x})=L_{n}(x)


  1. Mary L. Boas, 최준곤(역)(2008). 수리물리학 개정3판 p. 627. 연습문제 14 ↩︎