logo

ルジャンドル微分方程式の直列解法:ルジャンドル多項式 📂微分方程式

ルジャンドル微分方程式の直列解法:ルジャンドル多項式

定義1

以下の微分方程式をルジャンドルLegendre微分方程式と言う。

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

ルジャンドル微分方程式の解をルジャンドル多項式と言い、通常Pl(x)P_{l}(x)で示される。最初のいくつかのllによるルジャンドル多項式は次のようである。

P0(x)= 1P1(x)= xP2(x)= 12(3x21)P3(x)= 12(5x33x)P4(x)= 18(35x430x2+3)P5(x)= 18(63x570x3+15x) \begin{align*} P_{0}(x) =&\ 1 \\ P_{1}(x) =&\ x \\ P_2(x) =&\ \dfrac{1}{2}(3x^2-1) \\ P_{3}(x) =&\ \dfrac{1}{2}(5x^3-3x) \\ P_{4}(x) =&\ \dfrac{1}{8}(35x^4-30x^2+3) \\ P_{5}(x) =&\ \dfrac{1}{8}(63x^5-70x^3+15x) \\ \vdots& \end{align*}

説明

ルジャンドル微分方程式は、次のような形で紹介されることもある。

ddx[(1x)2dydx]+l(l+1)y=0 \dfrac{d}{dx}\left[ (1-x)^2 \dfrac{dy}{dx} \right] +l(l+1)y=0

これはシュツルム-リウヴィル理論Sturm-Liouville theoryで表されるものである。第一項を展開して整理すると、同じ式が得られる。ルジャンドル微分方程式を以下のように一般化したものを関連ルジャンドル微分方程式associated Legendre differential equationと言う。

(1x2)d2ydx22xdydx+(m21x2+l(l+1))y=0 (1-x^2)\dfrac{d^2 y}{dx^2} -2x\dfrac{dy}{dx}+\left( \dfrac{-m^2}{1-x^2} +l(l+1) \right) y=0

ここでm=0m=0の場合、ルジャンドル微分方程式となる。

ルジャンドル方程式は物理学や工学などで登場し、特に球面座標系でのラプラス方程式を解く時に見ることができる。物理学科ならば、電磁気学で球面座標系での電位を計算する時、量子力学で球面座標系でのシュレディンガー方程式を解く時に出会うことがある。解法の過程が長いため、教科書では通常、ロドリゲス公式で表される解答のみを記載することが多い。実際、物理学の学生は解法が非常に非常に気になるわけではなければ、知らなくても問題はない。

解法

係数に独立変数xxが含まれた形で、解が冪級数の形であると仮定すれば解くことができる


(1x2)y2xy+l(l+1)y=0 \begin{equation} (1-x^2)y^{\prime \prime} -2xy^{\prime}+l(l+1) y=0 \label{1} \end{equation}

ルジャンドル微分方程式の解を次のように仮定しよう。

y=a0+a1(xx0)+a2(xx0)2+=n=0an(xx0)n y=a_{0}+a_{1}(x-x_{0})+a_2(x-x_{0})^2+\cdots=\sum \limits_{n=0}^\infty a_{n}(x-x_{0})^n

この時x=0x=0の時、yy^{\prime \prime}の係数が(1x2)x=0=10(1-x^2)|_{x=0}=1\ne 0であるため、x0=0x_{0}=0と置く。すると級数解は

y=a0+a1x+a2x2+=n=0anxn \begin{equation} y=a_{0}+a_{1}x+a_2x^2+\cdots=\sum \limits_{n=0}^\infty a_{n}x^n \label{2} \end{equation}

解を級数と仮定したが、解法の最後に実際にはyyの項が有限であることがわかる。これで(1)\eqref{1}に代入するためにyy^{\prime}yy^{\prime \prime}を求めよう。

y=a1+2a2x+3a3x2+=n=1nanxn1 y^{\prime}=a_{1}+2a_2x+3a_{3}x^2+\cdots=\sum \limits_{n=1}^\infty na_{n}x^{n-1}

y=2a2+32a3x+43a4x2+=n=2n(n1)anxn2 y^{\prime \prime}=2a_2+3\cdot 2a_{3}x+4\cdot 3 a_{4}x^2 +\cdots = \sum \limits_{n=2} n(n-1)a_{n}x^{n-2}

これで(1)\eqref{1}y,y,yy, y^{\prime}, y^{\prime \prime}を代入すると

(1x2)n=2n(n1)anxn22xn=1nanxn1+l(l+1)n=0anxn=0 (1-x^2)\sum \limits_{n=2}^\infty n(n-1)a_{n}x^{n-2} -2x\sum \limits_{n=1}^\infty na_{n}x^{n-1}+l(l+1) \sum \limits_{n=0}^\infty a_{n}x^n=0

第一項の係数(1x2)(1-x^2)の括弧を展開して整理すると

n=2n(n1)anxn2x2n=2n(n1)anxn22xn=1nanxn1+l(l+1)n=0anxn=0 \sum \limits_{n=2}^\infty n(n-1)a_{n}x^{n-2} -x^2\sum \limits_{n=2}^\infty n(n-1)a_{n}x^{n-2} -2x\sum \limits_{n=1}^\infty na_{n}x^{n-1}+l(l+1) \sum \limits_{n=0}^\infty a_{n}x^n=0

    n=2n(n1)anxn2n=2n(n1)anxn2n=1nanxn+l(l+1)n=0anxn=0 \implies \sum \limits_{n=2} ^\infty n(n-1)a_{n}x^{n-2} -\sum \limits_{n=2}^\infty n(n-1)a_{n}x^{n} -2\sum \limits_{n=1}^\infty na_{n}x^{n}+l(l+1) \sum \limits_{n=0}^\infty a_{n}x^n=0

ここでのポイントは**xxの次数を合わせること**である。他は全てxnx^nで表されるのに対し、最初の級数だけがxn2x^{n-2}で表されているため、nnの代わりにn+2n+2を代入すると

n=0(n+2)(n+1)an+2xnn=2n(n1)anxn2n=1nanxn+l(l+1)n=0anxn=0 \sum \limits_{n=0} ^\infty (n+2)(n+1)a_{n+2}x^{n} -\sum \limits_{n=2}^\infty n(n-1)a_{n}x^{n} -2\sum \limits_{n=1}^\infty na_{n}x^{n}+l(l+1) \sum \limits_{n=0}^\infty a_{n}x^n=0

二番目の級数がx2x^2項から始まるので、他の級数からn=0,1n=0,1の項を外して、定数項は定数項同士、1次項は1次項同士をまとめると

[21a2+l(l+1)a0]+[32a32a1+l(l+1)a1]x+n=2[(n+2)(n+1)an+2n(n+1)an2nan+l(l+1)an]xn=0 \left[ 2\cdot 1 a_2+l(l+1)a_{0} \right]+\left[ 3\cdot 2 a_{3}-2a_{1}+l(l+1)a_{1} \right]x \\ + \sum \limits_{n=2}^\infty \left[ (n+2)(n+1)a_{n+2}-n(n+1)a_{n}-2na_{n}+l(l+1)a_{n} \right] x^n=0

上の式が成り立つためには全ての係数が00でなければならない

21a2+l(l+1)a0=0 2\cdot 1 a_2+l(l+1)a_{0} =0

32a32a1+l(l+1)a1=0 3\cdot 2 a_{3}-2a_{1}+l(l+1)a_{1} =0

(n+2)(n+1)an+2n(n+1)an2nan+l(l+1)an=0 (n+2)(n+1)a_{n+2}-n(n+1)a_{n}-2na_{n}+l(l+1)a_{n}=0

それぞれを整理すると

a2=l(l+1)21a0 \begin{equation} a_2=-\dfrac{l(l+1)}{2 \cdot 1}a_{0} \label{3} \end{equation}

a3=(l+2)(l1)32a1 \begin{equation} a_{3}=-\dfrac{(l+2)(l-1)}{3\cdot 2} a_{1} \label{4} \end{equation}

an+2=(l+n+1)(ln)(n+2)(n+1)an \begin{equation} a_{n+2}=-\dfrac{(l+n+1)(l-n)}{(n+2)(n+1)}a_{n} \label{5} \end{equation}

(3),(4),(5)\eqref{3}, \eqref{4}, \eqref{5}を利用すると、a0a_{0}a1a_{1}の値だけを知っていれば全ての係数を知ることができる。(3)\eqref{3}(5)\eqref{5}で偶数次項の係数を求めると

a4= (l+3)(l2)43a2=l(l2)(l+1)(l+3)4!a0a6= (l+5)(l4)65a4=l(l2)(l4)(l+1)(l+3)(l+5)6!a0 \begin{align*} a_{4} =&\ - \dfrac{(l+3)(l-2)}{ 4 \cdots 3}a_2 = \dfrac{l(l-2)(l+1)(l+3)}{4!}a_{0} \\ a_{6} =&\ -\dfrac{(l+5)(l-4)}{6\cdot5} a_{4} = -\dfrac{ l(l-2)(l-4)(l+1)(l+3)(l+5)}{6!} a_{0} \\ \vdots& \end{align*}

n=2m (m=1,2,3,)n=2m\ (m=1,2,3,\cdots)とすると

an=a2m=(1)ml(l2)(l2m+4)(l2m+2)(l+1)(l+3)(l+2m3)(l+2m1)(2m)!a0 a_{n}=a_{2m}=(-1)^m \dfrac{l(l-2)\cdots (l-2m+4)(l-2m+2)(l+1)(l+3)\cdots(l+2m-3)(l+2m-1)}{(2m)!}a_{0}

同様に(4)\eqref{4}(5)\eqref{5}で奇数次項の係数を求めると

a5= (l+4)(l3)54a3=(l+2)(l+4)(l1)(l3)5!a1a7= (l+6)(l5)76a5=(l+2)(l+4)(l+6)(l1)(l3)(l5)7!a1 \begin{align*} a_{5} =&\ -\dfrac{(l+4)(l-3)}{5\cdot 4}a_{3} = \dfrac{(l+2)(l+4)(l-1)(l-3)}{5!}a_{1} \\ a_{7} =&\ -\dfrac{(l+6)(l-5)}{7\cdot 6}a_{5} = -\dfrac{(l+2)(l+4)(l+6)(l-1)(l-3)(l-5)}{7!}a_{1} \\ \vdots& \end{align*}

n=2m+1 (m=1,2,3,)n=2m+1\ (m=1,2,3,\cdots)とすると

an=a2m+1=(1)m(l+2)(l+4)(l+2m2)(l+2m)(l1)(l3)(l2m+3)(l2m+1)(2m+1)!a1 a_{n}=a_{2m+1}=(-1)^m\dfrac{(l+2)(l+4)\cdots(l+2m-2)(l+2m)(l-1)(l-3)\cdots(l-2m+3)(l-2m+1)}{(2m+1)!}a_{1}

これで求めた係数を(2)\eqref{2}に代入して解を求めると

y=\a0+a1xl(l+1)2!a0x2(l+2)(l1)3!a1x3+l(l2)(l+1)(l+3)4!a0x4+(l+2)(l+4)(l1)(l3)5!a1x5++(1)ml(l2)(l2m+4)(l2m+2)(l+1)(l+3)(l+2m3)(l+2m1)(2m)!a0x2m+(1)m(l+2)(l+4)(l+2m2)(l+2m)(l1)(l3)(l2m+3)(l2m+1)(2m+1)!a1x2m+1+ \begin{align*} y =&\a_{0}+a_{1}x -\dfrac{l(l+1)}{2!}a_{0}x^2-\dfrac{(l+2)(l-1)}{3!}a_{1}x^3 + \dfrac{l(l-2)(l+1)(l+3)}{4!}a_{0}x^4+\dfrac{(l+2)(l+4)(l-1)(l-3)}{5!}a_{1}x^5 \\ &+ \cdots +(-1)^m \dfrac{l(l-2)\cdots (l-2m+4)(l-2m+2)(l+1)(l+3)\cdots(l+2m-3)(l+2m-1)}{(2m)!}a_{0}x^{2m} \\ &+ (-1)^m\dfrac{(l+2)(l+4)\cdots(l+2m-2)(l+2m)(l-1)(l-3)\cdots(l-2m+3)(l-2m+1)}{(2m+1)!}a_{1}x^{2m+1} +\cdots \end{align*}

(m=1,2,3,)(m=1,2,3,\cdots)偶数次項はa0a_{0}で、奇数次項はa1a_{1}でまとめると

y=\a0[1l(l+1)2!x2+l(l2)(l+1)(l+3)4!x4+m=3(1)ml(l2)(l2m+4)(l2m+2)(l+1)(l+3)(l+2m3)(l+2m1)(2m)!x2m]+a1[x(l+2)(l1)3!x3+(l+2)(l+4)(l1)(l3)5!x5+m=3(1)m(l+2)(l+4)(l+2m2)(l+2m)(l1)(l3)(l2m+3)(l2m+1)(2m+1)!x2m+1] \begin{align*} y =&\a_{0}\left[1-\dfrac{l(l+1)}{2!}x^2+\dfrac{l(l-2)(l+1)(l+3)}{4!}x^4 \right. \\ &\left.+\sum \limits_{m=3}^\infty (-1)^m \dfrac{l(l-2)\cdots (l-2m+4)(l-2m+2)(l+1)(l+3)\cdots(l+2m-3)(l+2m-1)}{(2m)!} x^{2m} \right] \\ &+ a_{1}\left[x- \dfrac{(l+2)(l-1)}{3!}x^3+\dfrac{(l+2)(l+4)(l-1)(l-3)}{5!}x^5 \right. \\ & \left. +\sum \limits_{m=3}^\infty (-1)^m\dfrac{(l+2)(l+4)\cdots(l+2m-2)(l+2m)(l-1)(l-3)\cdots(l-2m+3)(l-2m+1)}{(2m+1)!} x^{2m+1} \right] \end{align*}

最初の括弧をy0y_{0}、二番目の括弧をy1y_{1}とすると、ルジャンドル方程式の一般解は次のようになる。

y=a0y0+a1y1 y=a_{0}y_{0}+a_{1}y_{1}

二つの級数y0y_{0}y1y_{1}比率判定法により、x<1|x|<1の範囲で収束

することがわかる。(5)\eqref{5}によりan+2an=(l+n+1)(ln)(n+2)(n+1)=(n+l+1)(nl)(n+2)(n+1)\dfrac{a_{n+2}}{a_{n}}=-\dfrac{(l+n+1)(l-n)}{(n+2)(n+1)}=\dfrac{(n+l+1)(n-l)}{(n+2)(n+1)}であるため、比率判定法を使うと

limn(n+l+1)(nl)(n+2)(n+1)x2=x2<1 \lim \limits_{n \rightarrow \infty} \dfrac{(n+l+1)(n-l)}{(n+2)(n+1)}x^2=x^2<1

    1<x<1 \implies -1<x<1

しかし、多くの問題でx=cosθx=\cos \thetallは非負の整数の形で式が現れ、全てのθ\thetaに対して収束する解を得たい。つまり、x=±1x=\pm 1でも収束する解を見つけることが目標である。幸いにもllが整数の時は、欲しい解が存在し、その時llの値によって必ずy0,y1y_{0}, y_{1}のどちらかの解のみが存在する。ll00か偶数の時はy1y_{1}が発散し、y0y_{0}は偶数次項のみを持つ有限項の多項式となる。llが奇数ならばy0y_{0}が発散し、y1y_{1}は奇数次項のみを持つ有限項の多項式となる。表にまとめると以下のようになる。

llの値y0y_{0}y1y_{1}方程式の解
00か偶数有限項の多項式発散y=a0y0y=a_{0}y_{0}
奇数発散有限項の多項式y=a1y1y=a_{1}y_{1}
  • ケース1. ll00か偶数

    • l=0l=0の時、2次項からllを因数に持ち、全て00になるので、y0=1y_{0}=1

    • l=2l=2の時、4次項から(l2)(l-2)を因数に持ち、全て00になるので、y0=13x2y_{0}=1-3x^2

    • l=4l=4の時、6次項から(l4)(l-4)を因数に持ち、全て00になるので、y0=110x2+353x4y_{0}= 1-10x^2+\dfrac{35}{3}x^4

    そしてl=0l=0の時、x2=1x^2=1からy1=1+13+15+y_{1}=1+\frac{1}{3}+\frac{1}{5}+\cdotsであるが、これは積分判定法により発散する。他の偶数の時も同様である。したがって、ll00か偶数の時は、解が偶数次項のみを持つ有限項の多項式となる。つまり、級数y0y_{0}の特定の項までのみ残る形の解を得る。

  • ケース2. llが奇数

    偶数の時と反対の結果が現れる。

    • l=1l=1の時、3次項から(l1)(l-1)を因数に持ち、全て00になるので、y1=xy_{1}=x

    • l=3l=3の時、5次項から(l3)(l-3)を因数に持ち、全て00になるので、y1=x53x3y_{1}=x-\dfrac{5}{3}x^3

    • l=5l=5の時、7次項から(l5)(l-5)を因数に持ち、全て00になるので、y1=x143x3+215x5y_{1}=x-\dfrac{14}{3}x^3+\dfrac{21}{5}x^5

    l=1l=1の時、x2=1x^2=1からy0y_{0}は発散し、他の奇数の時も同様である。したがって、llが奇数の時は、解が奇数次項のみを持つ有限項の多項式となる。つまり、級数y1y_{1}の特定の項までのみ残る形の解を得る。

そして、llが負の場合は、llが0ではない整数の場合と同じであることがy0y_{0}y1y_{1}を見ればわかる。例えば、l=2l=2の場合とl=3l=-3の場合が同じであり、l=1l=1の場合とl=2l=-2の場合が同じである。したがって、llが非負の整数についてのみ考えれば良い。a0a_{0}a1a_{1}の値を上手く選んでx=1x=1の時の解がy(x)=1y(x)=1になるようにすると、これをルジャンドル多項式Legendre polynomialと言い、Pl(x)P_{l}(x)と書く。最初のいくつかのルジャンドル多項式は以下の通りである。

P0(x)= 1P1(x)= xP2(x)= 12(3x21)P3(x)= 12(5x33x)P4(x)= 18(35x430x2+3)P5(x)= 18(63x570x3+15x) \begin{align*} P_{0}(x) =&\ 1 \\ P_{1}(x) =&\ x \\ P_2(x) =&\ \dfrac{1}{2}(3x^2-1) \\ P_{3}(x) =&\ \dfrac{1}{2}(5x^3-3x) \\ P_{4}(x) =&\ \dfrac{1}{8}(35x^4-30x^2+3) \\ P_{5}(x) =&\ \dfrac{1}{8}(63x^5-70x^3+15x) \end{align*}

この結果はロドリゲス公式Rodrigues’ formulaで直接得ることもできる。


  1. Mary L. Boas, 数理物理学(Mathematical Methods in the Physical Sciences, 최준곤 訳) (3rd Edition, 2008), p577-580 ↩︎