Laguerre多項式

同義語:ラゲール多項式Laguerre polynomialsLaguerre polynomial一般化Laguerre多項式generalized Laguerre polynomials

概要

Laguerre多項式(Laguerre polynomials)とは、$\alpha>-1$ について $L_n^{(\alpha)}(x)=\frac{x^{-\alpha}e^x}{n!}\frac{d^n}{dx^n}(e^{-x}x^{n+\alpha})$ で定まる多項式の列で、$\alpha\ne0$ の場合は一般化 Laguerre 多項式ともいう。半直線 $(0,\infty)$ 上の重み $x^\alpha e^{-x}$ に関する直交多項式である。3 項漸化式、微分方程式 $xy''+(\alpha+1-x)y'+ny=0$、母関数 $(1-t)^{-\alpha-1}e^{-xt/(1-t)}$ をもち、Gauss–Laguerre 公式や水素原子の固有関数の動径部分に現れる。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 直交多項式, Γ関数, 積の微分とLeibnizの公式

Laguerre 多項式(Laguerre polynomials)は、半直線 $(0,\infty)$ の上で重み $x^\alpha e^{-x}$($\alpha>-1$)を付けた積分に関して互いに直交する多項式の列である。$\alpha=0$ の場合を単に Laguerre 多項式、一般の $\alpha$ の場合を 一般化 Laguerre 多項式 ということが多い。Legendre多項式 が有界な区間の直交多項式の基本例であるのに対し、Laguerre 多項式は無限区間の基本例であり、Rodrigues の公式・3 項漸化式・2 階の微分方程式・指数関数を含む母関数をもつ。指数関数で減衰する重みの積分を近似する Gauss–Laguerre 公式や、水素原子のエネルギー固有関数の動径部分に現れる。
この記事では、以下つねに $\alpha>-1$ とする。

定義

Laguerre 多項式

$n=0,1,2,\dots$ と $x>0$ について
$$ L_n^{(\alpha)}(x):=\frac{x^{-\alpha}e^x}{n!}\,\frac{d^n}{dx^n}\bigl(e^{-x}x^{n+\alpha}\bigr) $$
とおき、$L_n^{(\alpha)}$ を $n$ 次の(一般化)Laguerre 多項式 という。この式を Rodrigues の公式 という。$\alpha=0$ のとき $L_n:=L_n^{(0)}$ と書く。

右辺が多項式であることは次の命題で示す。以下、実数 $\gamma$ と整数 $m\ge0$ について $\binom\gamma m:=\frac{\gamma(\gamma-1)\cdots(\gamma-m+1)}{m!}$($m=0$ では $1$)とする。

係数の公式

$$ L_n^{(\alpha)}(x)=\sum_{k=0}^n\binom{n+\alpha}{n-k}\frac{(-x)^k}{k!} $$
である。とくに $L_n^{(\alpha)}$ は次数ちょうど $n$ の多項式で、最高次の係数は $\frac{(-1)^n}{n!}$、定数項は $L_n^{(\alpha)}(0)=\binom{n+\alpha}n=\frac{(\alpha+1)(\alpha+2)\cdots(\alpha+n)}{n!}>0$ である。この式で $L_n^{(\alpha)}$ を $x\le0$ にも延ばして考える。

Leibniz の公式で展開する

積の微分とLeibnizの公式 の定理「Leibniz の公式」により
$$ \frac{d^n}{dx^n}\bigl(e^{-x}x^{n+\alpha}\bigr)=\sum_{j=0}^n\binom nj(-1)^{n-j}e^{-x}\,(n+\alpha)(n+\alpha-1)\cdots(n+\alpha-j+1)\,x^{n+\alpha-j} $$
である($e^{-x}$ の $n-j$ 階導関数は $(-1)^{n-j}e^{-x}$)。$\frac{x^{-\alpha}e^x}{n!}$ を掛けて $k:=n-j$ とおくと、$x^k$ の係数は
$$ \frac{(-1)^k}{n!}\binom nk(n+\alpha)\cdots(k+\alpha+1)=\frac{(-1)^k}{k!}\cdot\frac{(n+\alpha)\cdots(k+\alpha+1)}{(n-k)!}=\frac{(-1)^k}{k!}\binom{n+\alpha}{n-k} $$
である(積 $(n+\alpha)\cdots(k+\alpha+1)$ は $n-k$ 個の因子からなる)。$k=n$ の項が最高次で係数 $\frac{(-1)^n}{n!}$、$k=0$ の項が定数項である。$\alpha>-1$ なので $\alpha+1,\dots,\alpha+n$ はすべて正である。$\square$

最初のいくつか

$$ L_0^{(\alpha)}=1,\qquad L_1^{(\alpha)}=1+\alpha-x,\qquad L_2^{(\alpha)}=\frac{(\alpha+1)(\alpha+2)}2-(\alpha+2)x+\frac{x^2}2. $$
$\alpha=0$ では
$$ L_1=1-x,\qquad L_2=1-2x+\frac{x^2}2,\qquad L_3=1-3x+\frac{3x^2}2-\frac{x^3}6 $$
である。たとえば $\alpha=0$、$n=2$ の Rodrigues の公式では $\frac{d^2}{dx^2}(e^{-x}x^2)=e^{-x}(x^2-4x+2)$ で、$\frac{e^x}{2}$ を掛けると $L_2$ になる。

主定理 1:直交性

端点で消える導関数

$u(x):=e^{-x}x^{n+\alpha}$ とする。$0\le j\le n$ について、$u^{(j)}(x)=e^{-x}x^{n+\alpha-j}r_j(x)$ と書ける。ここで $r_j$ は次数 $j$ 以下の多項式である。とくに $0\le j\le n-1$ と任意の多項式 $q$ について、$x\to+0$ と $x\to\infty$ のどちらでも $q(x)u^{(j)}(x)\to0$ である。

証明を開く

$j=0$ では $r_0=1$ である。$u^{(j)}=e^{-x}x^{\beta}r_j$($\beta:=n+\alpha-j$)なら、微分して $u^{(j+1)}=e^{-x}x^{\beta-1}\bigl(-xr_j+\beta r_j+xr_j'\bigr)$ で、括弧の中は次数 $j+1$ 以下の多項式である。帰納法で前半を得る。$j\le n-1$ なら $\beta\ge1+\alpha>0$ なので、$x\to+0$ で $x^\beta\to0$ であり、$q\,r_j$ は $0$ の近くで有界だから $q\,u^{(j)}\to0$ である。$x\to\infty$ では、$q\,r_j\,x^{\beta}$ は多項式と冪の積で、Γ関数 の補題「指数減衰による優関数」により $x^{a}e^{-x/2}$ は $x\ge1$ で有界なので、$|q\,r_j|x^\beta e^{-x}\le Cx^{N}e^{-x}\le C'e^{-x/2}\to0$ である($N$ は十分大きい整数)。$\square$

直交性と大きさ

$m,n\ge0$ について
$$ \int_0^\infty L_m^{(\alpha)}(x)L_n^{(\alpha)}(x)\,x^\alpha e^{-x}\,dx=\begin{cases}0&(m\ne n),\\[4pt]\dfrac{\Gamma(n+\alpha+1)}{n!}&(m=n)\end{cases} $$
である。さらに、次数 $n-1$ 以下のすべての多項式 $q$ について $\int_0^\infty q\,L_n^{(\alpha)}\,x^\alpha e^{-x}\,dx=0$ である。

部分積分で微分を $q$ に移す

$u:=e^{-x}x^{n+\alpha}$ とおくと、定義から $L_n^{(\alpha)}x^\alpha e^{-x}=\frac1{n!}u^{(n)}$ である。$q$ を多項式とする。$0\le i\le n$ について、$q^{(i)}u^{(n-i)}$ は lem-lag-vanish により $e^{-x}x^{\alpha+i}\times(\text{多項式})$ の形で、$\alpha+i>-1$ なので $(0,\infty)$ で絶対可積分である(多項式は冪 $x^k$ の和で、$s>-1$ なら $\int_0^\infty x^se^{-x}\,dx=\Gamma(s+1)$ は絶対収束する。Γ関数 の定義「$\Gamma$ 関数」)。$0< \varepsilon< R$ で部分積分すると
$$ \int_\varepsilon^Rq^{(i)}u^{(n-i)}\,dx=\Bigl[q^{(i)}u^{(n-i-1)}\Bigr]_\varepsilon^R-\int_\varepsilon^Rq^{(i+1)}u^{(n-i-1)}\,dx\qquad(0\le i\le n-1) $$
であり、$\varepsilon\to+0$、$R\to\infty$ とすると lem-lag-vanish($j=n-i-1\le n-1$)により境界の項は $0$ になる。これを $i=0,1,\dots$ と繰り返して
$$ \int_0^\infty q\,u^{(n)}\,dx=(-1)^i\int_0^\infty q^{(i)}u^{(n-i)}\,dx\qquad(0\le i\le n) $$
を得る。
低い次数との直交性。 $q$ の次数 $m$ が $n-1$ 以下なら、$i=m+1\le n$ として $q^{(m+1)}=0$ なので積分は $0$ である。$m\ne n$ のときは、次数の小さいほうを $q$ にとれば第 1 式の $0$ を得る。
大きさ。 $q=L_n^{(\alpha)}$、$i=n$ とする。prop-lag-explicit により $q^{(n)}=n!\cdot\frac{(-1)^n}{n!}=(-1)^n$ なので
$$ \int_0^\infty\bigl(L_n^{(\alpha)}\bigr)^2x^\alpha e^{-x}\,dx=\frac1{n!}(-1)^n(-1)^n\int_0^\infty e^{-x}x^{n+\alpha}\,dx=\frac{\Gamma(n+\alpha+1)}{n!} $$
である(Γ関数 の定義「$\Gamma$ 関数」)。$\square$

$w(x)=x^\alpha e^{-x}$ は $(0,\infty)$ 上の正の連続関数で、モーメント $\int_0^\infty x^kw\,dx=\Gamma(k+\alpha+1)$ はすべて有限なので、直交多項式 の定義「重み関数と内積」の重み関数である。最高次の係数で割った $(-1)^nn!\,L_n^{(\alpha)}$ はこの重みのモニック直交多項式であり(直交多項式 の命題「モニック直交多項式の存在と一意性」の 3)、同記事の一般論がすべて当てはまる。

  • $L_n^{(\alpha)}$ の根は $(0,\infty)$ の中の相異なる $n$ 個の実数で、$L_n^{(\alpha)}$ と $L_{n+1}^{(\alpha)}$ の根は交互に並ぶ(同記事の定理「零点の位置」「零点の交互性」)。
  • $L_n^{(\alpha)}$ の根を節点にとる $n$ 点の公式 $\int_0^\infty f(x)x^\alpha e^{-x}\,dx\approx\sum_k\lambda_kf(x_k)$ は、次数 $2n-1$ 以下の多項式で正確である(同記事の定理「Gauss 求積」)。これを Gauss–Laguerre 公式 という。
Gauss–Laguerre 公式($\alpha=0$、$n=2$)

$L_2=1-2x+\frac{x^2}2$ の根は $x_{1,2}=2\mp\sqrt2$ である。重み $\lambda_k=\int_0^\infty\ell_k(x)e^{-x}\,dx$ は、$\ell_1(x)=\frac{x-x_2}{x_1-x_2}$ と $\int_0^\infty xe^{-x}\,dx=1$ から $\lambda_1=\frac{1-x_2}{x_1-x_2}=\frac{2+\sqrt2}4$、同様に $\lambda_2=\frac{2-\sqrt2}4$ である。公式
$$ \int_0^\infty f(x)e^{-x}\,dx\approx\frac{2+\sqrt2}4f(2-\sqrt2)+\frac{2-\sqrt2}4f(2+\sqrt2) $$
は $f=1,x,x^2,x^3$ で正しい値 $0!,1!,2!,3!=1,1,2,6$ を与える(たとえば $x_1^3=20-14\sqrt2$、$x_2^3=20+14\sqrt2$ から $\frac{(2+\sqrt2)(20-14\sqrt2)+(2-\sqrt2)(20+14\sqrt2)}4=\frac{24}4=6$)。$f=x^4$ では $\frac{80}4=20$ となり、正しい値 $4!=24$ と一致しない。

母関数

母関数

実数 $x$ と $|t|<1$ について
$$ \frac1{(1-t)^{\alpha+1}}\exp\Bigl(-\frac{xt}{1-t}\Bigr)=\sum_{n=0}^\infty L_n^{(\alpha)}(x)\,t^n $$
であり、右辺は絶対収束する。

指数関数と一般二項定理の 2 重級数をまとめ直す

指数関数の級数から
$$ \exp\Bigl(-\frac{xt}{1-t}\Bigr)=\sum_{k=0}^\infty\frac{(-xt)^k}{k!}\,(1-t)^{-k} $$
なので、左辺は $\sum_k\frac{(-xt)^k}{k!}(1-t)^{-(\alpha+k+1)}$ である。$\gamma:=\alpha+k+1>0$ について、二項定理と組合せの恒等式 の定理「一般二項定理」から $(1-t)^{-\gamma}=\sum_{j\ge0}\binom{-\gamma}j(-t)^j=\sum_{j\ge0}\binom{\gamma+j-1}jt^j$ で、係数 $\binom{\gamma+j-1}j=\binom{\alpha+k+j}j$ はすべて正である。よって左辺は 2 重級数
$$ \sum_{k=0}^\infty\sum_{j=0}^\infty\frac{(-x)^k}{k!}\binom{\alpha+k+j}j\,t^{k+j} $$
に等しい。各項を絶対値に替えた和は、同じ計算を $x,t$ の代わりに $-|x|,|t|$ で行ったもので、$(1-|t|)^{-\alpha-1}\exp\bigl(\frac{|x||t|}{1-|t|}\bigr)<\infty$ である。したがって項の順序を変えてよく、$n=k+j$ ごとにまとめると $t^n$ の係数は
$$ \sum_{k=0}^n\frac{(-x)^k}{k!}\binom{\alpha+n}{n-k}=L_n^{(\alpha)}(x) $$
である(prop-lag-explicit)。$\square$

$x=0$ では左辺は $(1-t)^{-\alpha-1}$ で、$t^n$ の係数は $\binom{n+\alpha}n=L_n^{(\alpha)}(0)$ である。同じ等式は DLMF (18.12.13) にある。

主定理 2:漸化式と微分方程式

3 項漸化式

$n\ge1$ について
$$ (n+1)L_{n+1}^{(\alpha)}(x)=(2n+\alpha+1-x)L_n^{(\alpha)}(x)-(n+\alpha)L_{n-1}^{(\alpha)}(x) $$
である。$L_0^{(\alpha)}=1$、$L_1^{(\alpha)}=1+\alpha-x$ とこの漸化式から、すべての $L_n^{(\alpha)}$ が順に決まる。

母関数が満たす微分方程式の係数を比べる

$x$ を固定し、$G(t):=(1-t)^{-\alpha-1}\exp\bigl(-\frac{xt}{1-t}\bigr)$ とおく。$\frac{d}{dt}\frac{t}{1-t}=\frac1{(1-t)^2}$ から
$$ G'(t)=\Bigl(\frac{\alpha+1}{1-t}-\frac x{(1-t)^2}\Bigr)G(t),\qquad\text{すなわち}\qquad(1-t)^2G'(t)=\bigl((\alpha+1)(1-t)-x\bigr)G(t) $$
である。prop-lag-generating により $G(t)=\sum_na_nt^n$($a_n:=L_n^{(\alpha)}(x)$)は $|t|<1$ で収束する冪級数なので、項別に微分でき、和の関数が等しい冪級数の係数は等しい(冪級数 の定理「項別微分」と系「何回でも微分できることと係数の公式」)。両辺の $t^n$($n\ge1$)の係数を比べると
$$ (n+1)a_{n+1}-2na_n+(n-1)a_{n-1}=(\alpha+1-x)a_n-(\alpha+1)a_{n-1} $$
で、整理すると主張の式になる。$\square$

漸化式を最高次の係数で割ると、モニック直交多項式 $p_n=(-1)^nn!\,L_n^{(\alpha)}$ について $p_{n+1}=(x-(2n+\alpha+1))p_n-n(n+\alpha)p_{n-1}$ となる。$\alpha=0$ では 直交多項式 の例「古典直交多項式」の表の $\beta_n=2n+1$、$\gamma_n=n^2$ である。

微分と添字の関係
  1. $n\ge1$ について $\dfrac{d}{dx}L_n^{(\alpha)}=-L_{n-1}^{(\alpha+1)}$。
  2. $n\ge1$ について $L_n^{(\alpha)}=L_n^{(\alpha+1)}-L_{n-1}^{(\alpha+1)}$。

1:prop-lag-explicit の式を微分して $j:=k-1$ とおくと
$$ \frac{d}{dx}L_n^{(\alpha)}=\sum_{k=1}^n\binom{n+\alpha}{n-k}\frac{(-1)^kx^{k-1}}{(k-1)!}=-\sum_{j=0}^{n-1}\binom{(n-1)+(\alpha+1)}{(n-1)-j}\frac{(-x)^j}{j!}=-L_{n-1}^{(\alpha+1)} $$
である。2:母関数は $\alpha$ を $\alpha+1$ に替えたものの $(1-t)$ 倍なので、$\sum_nL_n^{(\alpha)}t^n=(1-t)\sum_nL_n^{(\alpha+1)}t^n$ の $t^n$ の係数を比べればよい(prop-lag-generating の右辺は絶対収束する冪級数なので、積の係数は各項の積の和である)。$\square$

Laguerre の微分方程式

$y=L_n^{(\alpha)}(x)$ は
$$ xy''+(\alpha+1-x)y'+ny=0,\qquad\text{すなわち}\qquad\bigl(x^{\alpha+1}e^{-x}y'\bigr)'+n\,x^\alpha e^{-x}y=0\quad(x>0) $$
を満たす。

係数の比を使う

$y=\sum_{k=0}^nc_kx^k$、$c_k:=\binom{n+\alpha}{n-k}\frac{(-1)^k}{k!}$ とする($c_{n+1}:=0$)。$xy''+(\alpha+1)y'-xy'+ny$ の $x^k$ の係数は
$$ (k+1)k\,c_{k+1}+(\alpha+1)(k+1)c_{k+1}-kc_k+nc_k=(k+1)(k+\alpha+1)c_{k+1}+(n-k)c_k $$
である。$0\le k\le n-1$ では、$\binom{n+\alpha}{n-k-1}\big/\binom{n+\alpha}{n-k}=\frac{n-k}{k+\alpha+1}$(分子は分母より因子 $k+\alpha+1$ が 1 つ少なく、階乗は $(n-k-1)!$ と $(n-k)!$)から $c_{k+1}=-\frac{n-k}{(k+1)(k+\alpha+1)}c_k$ なので、この係数は $0$ である。$k=n$ でも $c_{n+1}=0$、$n-k=0$ で $0$ である。後半の形は $\bigl(x^{\alpha+1}e^{-x}y'\bigr)'=x^\alpha e^{-x}\bigl(xy''+(\alpha+1-x)y'\bigr)$ による。$\square$

後半の形は Sturm–Liouville 型の方程式 $-(py')'=\lambda wy$($p=x^{\alpha+1}e^{-x}$、$w=x^\alpha e^{-x}$、$\lambda=n$)である。区間が有界でなく、端点 $0$ で $p$ が $0$ になるので、Sturm–Liouville問題 の正則な問題の枠には入らない。それでも thm-lag-orthogonality の $m\ne n$ の場合は、$\int_0^\infty(pL_m')'L_n\,dx=-\int_0^\infty pL_m'L_n'\,dx$(境界の項は $x^{\alpha+1}e^{-x}$ が両端で $0$ に近づくので消える)が $m,n$ について対称であることからも分かる。

応用

水素原子の動径関数

$\gamma>0$ とし、$\mathbb{R}^3$ の Schrödinger 作用素 $H=-\Delta-\frac\gamma{|x|}$ を考える。$H$ の固有値は $E_n=-\frac{\gamma^2}{4(n+1)^2}$($n=0,1,2,\dots$)で、固有値 $E_n$ の固有空間の正規直交基底として、$0\le l\le n$、$|m|\le l$ について
$$ \psi_{n,l,m}(x)=R_{n,l}(r)\,Y_l^m(x/r),\qquad R_{n,l}(r)=C_{n,l}\,r^le^{-\frac{\gamma r}{2(n+1)}}L_{n-l}^{(2l+1)}\Bigl(\frac{\gamma r}{n+1}\Bigr)\qquad(r=|x|) $$
がとれる。ここで $C_{n,l}>0$ は正規化の定数、$Y_l^m$ は 球面調和関数 である(Tes09 Theorem 10.9、(10.64)–(10.65)、pp. 231–232。この記事では証明しない)。動径の方程式を、変数 $\frac{\gamma r}{n+1}$ と未知関数の置き換えで書き直すと、thm-lag-ode で $\alpha=2l+1$ とした方程式(Kummer の方程式の特別な場合)が現れ、その多項式解が Laguerre 多項式である。固有値 $E_n$ の重複度 $\sum_{l=0}^n(2l+1)=(n+1)^2$ は、量子力学の関数解析的構造 の定理「水素型Hamilton作用素の離散スペクトル」の重複度と同じものである(同記事は番号を $1$ ずらして書いている)。

Hermite 多項式(実数直線上の重み $e^{-x^2}$ の直交多項式)は、$\alpha=\mp\frac12$ の Laguerre 多項式で $H_{2n}(x)=(-1)^n2^{2n}n!\,L_n^{(-1/2)}(x^2)$、$H_{2n+1}(x)=(-1)^n2^{2n+1}n!\,x\,L_n^{(1/2)}(x^2)$ と表される(DLMF (18.7.19)–(18.7.20)。この記事では証明しない)。

反例:条件を外すと成り立たない

外す条件反例成り立たなくなること
$\alpha>-1$$\alpha=-1$重み $x^{-1}e^{-x}$ が $0$ の近くで積分できず、$\int_0^\infty\bigl(L_0^{(-1)}\bigr)^2x^{-1}e^{-x}\,dx=\infty$ となって大きさの式が意味をもたない
積分区間が $(0,\infty)$ であること区間 $(0,1)$、$\alpha=0$$\int_0^1L_0L_1e^{-x}\,dx=e^{-1}\ne0$ で、直交しない
重みの $x^\alpha$ が多項式の $\alpha$ と一致すること$L_n^{(1)}$ を重み $e^{-x}$ で積分$\int_0^\infty L_0^{(1)}L_1^{(1)}e^{-x}\,dx=1\ne0$ で、直交しない
母関数の $\lvert t\rvert<1$$x=0$、$\alpha=0$、$t=1$左辺は定義されず、右辺 $\sum_nL_n(0)=\sum_n1$ は発散する
Gauss–Laguerre 公式の次数 $2n-1$ 以下$n=2$、$f=x^4$公式の値 $20$ が積分 $24$ と一致しない
反例の確かめ
  1. $L_0^{(-1)}=1$ で、$\int_0^1x^{-1}e^{-x}\,dx\ge e^{-1}\int_0^1x^{-1}\,dx=\infty$ である。
  2. $\frac{d}{dx}(xe^{-x})=(1-x)e^{-x}$ なので $\int_0^1(1-x)e^{-x}\,dx=e^{-1}$ である。
  3. $L_0^{(1)}=1$、$L_1^{(1)}=2-x$ で、$\int_0^\infty(2-x)e^{-x}\,dx=2-1=1$ である。重み $xe^{-x}$ なら $\int_0^\infty(2-x)xe^{-x}\,dx=2-2=0$ で直交する。
  4. $x=0$、$\alpha=0$ では左辺は $\frac1{1-t}$ で、$t\to1$ で限りなく大きくなる。
  5. ex-lag-gauss の計算である($x_1^4=68-48\sqrt2$、$x_2^4=68+48\sqrt2$)。

関連項目

参考文献

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する