logo

エルミート多項式の直交性 📂関数

エルミート多項式の直交性

定理

エルミート多項式 {Hn}n=0\left\{ H_{n} \right\}_{n=0}^{\infty}は区間 (,)(-\infty, \infty)重み関数 w(x)=ex2w(x)=e^{-x^{2}}により直交している。

HnHmex2=ex2Hn(x)Hm(x)dx=π2nn!δnm \braket{ H_{n} | H_{m} }_{e^{-x^{2}}} =\int_{-\infty}^{\infty}e^{-x^{2}}H_{n}(x)H_{m}(x)dx=\sqrt{\pi}2^{n}n!\delta_{nm}

ここで、δnm\delta_{nm}クロネッカーのデルタである。

証明

場合 1: n=mn=m

微分作用素D=ddxD = \dfrac{d}{dx}と表記する。

ex2Hn(x)Hn(x)dx \int_{-\infty}^{\infty} e^{-x^{2}}H_{n}(x)H_{n}(x)dx

エルミート多項式

Hn(x)=(1)nex2dndxnex2=(1)nex2Dnex2 H_{n}(x) = (-1)^{n}e^{x^{2}}\frac{d^{n}}{dx^{n}}e^{-x^{2}} = (-1)^{n}e^{x^{2}}D^{n}e^{-x^{2}}

上の式の前の部分Hn(x)H_{n}(x)を解くと、以下のようになる。

ex2(1)nex2[Dnex2]Hn(x)dx=(1)n[Dnex2]Hn(x)dx \begin{align*} \int_{-\infty}^{\infty} e^{-x^{2}}(-1)^{n}e^{x^{2}}\left[D^{n}e^{-x^{2}} \right]H_{n}(x)dx &= \int_{-\infty}^{\infty} (-1)^{n}\left[D^{n}e^{-x^{2}} \right]H_{n}(x)dx \end{align*}

式を部分的に解くと、以下の式を得る。

(1)n[Dnex2]Hn(x)dx=[(1)n(Dn1ex2)Hn(x)](1)n[Dn1ex2]Hn(x)dx \begin{align*} &\int_{-\infty}^{\infty} (-1)^{n}\left[D^{n}e^{-x^{2}} \right]H_{n}(x)dx \\ &=\left[ (-1)^{n}\left(D^{n-1}e^{-x^{2}}\right)H_{n}(x) \right]_{-\infty}^{\infty}-\int_{-\infty}^{\infty}(-1)^{n}\left[ D^{n-1}e^{-x^{2}}\right]H^{\prime}_{n}(x)dx \tag{1} \end{align*}

ここで、最初の項はlimx±Dn1ex2=0\lim \limits_{x\rightarrow \pm\infty}D^{n-1}e^{-x^{2}}=0であるため、00である。この極限が00に収束する理由は、任意のnnに対して、

Dnex2=p(x)ex2 D^{n}e^{-x^{2}}=p(x)e^{-x^{2}} と表されるためである。ここで、p(x)p(x)は任意の多項式である。x±x \rightarrow \pm \inftyの時、ex2e^{-x^{2}}00に行く速度の方が多項式が発散する速度よりもずっと速いため、極限は00に収束する。

エルミート多項式の再帰関係 Hn(x)=2nHn1(x) H_{n}^{\prime}(x) =2nH_{n-1}(x)

エルミート多項式の再帰関係により、(1)(1)は次のようになる。

2n(1)n[Dn1ex2]Hn1(x)dx -2n\int_{-\infty}^{\infty}(-1)^{n}\left[ D^{n-1}e^{-x^{2}} \right]H_{n-1}(x)dx

さっきと同じ論理で部分積分をもう一度行うと、次を得る。

(1)222n(n1)(1)n[Dn2ex2]Hn2(x)dx (-1)^{2}2^{2}n(n-1)\int_{-\infty}^{\infty}(-1)^{n}\left[ D^{n-2}e^{-x^{2}} \right]H_{n-2}(x)dx

したがって、部分積分をnn回行うと、以下の式を得る。

(1)n2nn!(1)nex2H0(x)dx (-1)^{n}2^{n}n!\int_{-\infty}^{\infty}(-1)^{n}e^{-x^{2}}H_{0}(x)dx

H0(x)=1H_{0}(x)=1であるため、以下のように整理される。

2nn!ex2dx 2^{n}n!\int_{-\infty}^{\infty}e^{-x^{2}}dx

上記の積分はガウス積分で、その値はπ\sqrt{\pi}である。したがって、以下を得る。

ex2Hn(x)Hn(x)dx=π2nn! \int_{-\infty}^{\infty} e^{-x^{2}}H_{n}(x)H_{n}(x)dx=\sqrt{\pi}2^{n}n!

場合 2: nmn\ne m

一般性を失わずに、n>mn \gt mとする。

ex2Hn(x)Hm(x)dx \int_{-\infty}^{\infty}e^{-x^{2}}H_{n}(x)H_{m}(x)dx

ここで、Hn(x)H_{n}(x)だけを解くと、次のようになる。

ex2(1)nex2[Dnex2]Hm(x)dx=(1)n[Dnex2]Hm(x)dx \int_{-\infty}^{\infty}e^{-x^{2}}(-1)^{n}e^{x^{2}}\left[D^{n}e^{-x^{2}}\right]H_{m}(x)dx=\int_{-\infty}^{\infty}(-1)^{n}\left[D^{n}e^{-x^{2}}\right]H_{m}(x)dx

n=mn=mの場合の証明と同じ方法で、mm回の部分積分を行うと、以下の式を得る。

(1)n+m2mm![Dnmex2]1dx (-1)^{n+m}2^{m}m!\int_{-\infty}^{\infty}\left[ D^{n-m}e^{-x^{2}} \right]\cdot 1 dx

もう一度部分積分すると、次のようになる。

(1)n+m2mm![Dnmex2]1dx=(1)n+m2mm!([Dnm1ex2]+2(m+1)[Dnm1ex2]0dx) \begin{align*} &(-1)^{n+m}2^{m}m!\int_{-\infty}^{\infty}\left[ D^{n-m}e^{-x^{2}} \right]\cdot 1 dx \\ &= (-1)^{n+m}2^{m}m!\left(\left[D^{n-m-1}e^{-x^{2}} \right]_{-\infty}^{\infty}+2(m+1)\int_{-\infty}^{\infty}\left[ D^{n-m-1}e^{-x^{2}}\right]\cdot 0 dx \right) \end{align*}

最初の項は、上の証明で説明したように、00であり、二番目の項も00である。したがって、以下を得る。

ex2Hn(x)Hm(x)dx=0,nm \int_{-\infty}^{\infty}e^{-x^{2}}H_{n}(x)H_{m}(x)dx=0, \quad n \ne m