logo

チェビシェフの微分方程式とチェビシェフ多項式 📂微分方程式

チェビシェフの微分方程式とチェビシェフ多項式

定義

次の微分方程式をチェビシェフChebyshev の微分方程式と呼ぶ。

(1x2)d2ydx2xdydx+n2y=0 \begin{equation} (1-x^2)\dfrac{d^2 y}{dx^2} -x\dfrac{dy}{dx}+n^2 y=0 \label{def1} \end{equation}

チェビシェフ微分方程式の解をチェビシェフ多項式と呼び、通常、Tn(x)T_{n}(x)で表記する。Tn(x)T_{n}(x)の一般項は以下のとおりである。

  • nnが偶数の場合

    1λ22!x2+λ2(λ222)4!x4+m=3(1)mλ2(λ222)(λ2(2m2)2)(2m)!x2m 1-\dfrac{\lambda^2}{2!}x^2+\dfrac{\lambda^2(\lambda^2-2^2)}{4!}x^4+\sum \limits_{m=3}^\infty (-1)^m \dfrac{\lambda^2(\lambda^2-2^2)\cdots(\lambda^2-(2m-2)^2)}{(2m)!} x^{2m}

  • nnが奇数の場合

xλ2123!x3+(λ212)(λ232)5!x5+m=3(1)m(λ212)(λ232)(λ2(2m1)2)(2m+1)!x2m+1 x-\dfrac{\lambda^2-1^2}{3!}x^3+\dfrac{(\lambda^2-1^2)(\lambda^2-3^2)}{5!}x^5+\sum \limits_{m=3}^\infty (-1)^m\dfrac{(\lambda^2-1^2)(\lambda^2-3^2) \cdots (\lambda^2-(2m-1)^2)}{(2m+1)!} x^{2m+1}

特に、最初の数個の多項式は以下の通りである。

T0(x)=1T1(x)=xT2(x)=2x21T3(x)=4x33x \begin{align*} T_{0}(x) &= 1 \\ T_{1}(x) &= x \\ T_2(x) &= 2x^2-1 \\ T_{3}(x) &= 4x^3-3x \\ \vdots & \end{align*}

定理

チェビシェフ多項式 TnT_{n}に対して、次の等式が成り立つ。

Tn(cost)=cos(nt) T_{n}(\cos t)= \cos (nt)

説明

Tn(x)T_{n}(x)xxに関するnn次の多項式であったので、Tn(cost)T_{n}(\cos t)cost\cos tに関する多項式である。したがってチェビシェフ多項式はcos(nt)\cos (nt)cost\cos tに関するnn次の多項式として展開したものとも理解できる。

n=2, 3n=2,\ 3の時に、うまく合うか確認すると

cos2t=cos2t1=T2(cost)    T2(x)=x21 \cos 2t=\cos ^2 t-1=T_2(\cos t)\iff T_2(x)=x^2-1

cos3t=4cos3t3cost=T3(cost)    T3(x)=4x33x \cos 3t=4\cos ^3 t-3\cos t=T_{3}(\cos t) \iff T_{3}(x)=4x^3-3x

また、x=costx=\cos tので、arccosx=t\arccos x=tになり、上の式に代入すると

Tn(x)=cos(narccosx)orTn(x)=cos(ncos1x) T_{n}(x)=\cos(n\arccos x) \quad \text{or} \quad T_{n}(x) = \cos (n\cos^{-1}x)

証明

戦略: x=costx=\cos tで置換した時に、y=cos(nt)y=\cos (nt)がチェビシェフの微分方程式の解になることを示す。


x=costx=\cos tとすると

dx=sintdt    dtdx=1sint dx=-\sin t dt \quad \implies \quad \dfrac{dt}{dx}=-\dfrac{1}{\sin t}

したがって、yy^{\prime}は次のようになる。

dydx=dydtdtdx=1sintdydt \dfrac{dy}{dx}=\dfrac{dy}{dt} \dfrac{dt}{dx}=-\dfrac{1}{\sin t}\dfrac{dy}{dt}

yy^{\prime \prime}は以下の通りである。

d2ydx2=ddx(dydx)=ddt(dydx)dtdx=ddt(1sintdydt)1sint=1sint(costsin2tdydt+1sintd2ydt2)=1sin2t(costsintdydt+d2ydt2) \begin{align*} \dfrac{d^2 y}{dx^2} &= \dfrac{d}{dx} \left(\dfrac{dy}{dx}\right) \\ &= \dfrac{d}{dt}\left(\dfrac{dy}{dx} \right) \dfrac{dt}{dx} \\ &=\dfrac{d}{dt}\left(\dfrac{1}{\sin t}\dfrac{dy}{dt}\right) \dfrac{1}{\sin t} \\ &= \dfrac{1}{\sin t} \left( \dfrac{-\cos t}{\sin^2 t}\dfrac{dy}{dt}+\dfrac{1}{\sin t}\dfrac{d^2y}{dt^2} \right) \\ &= \dfrac{1}{\sin ^2 t} \left( \dfrac{-\cos t}{\sin t}\dfrac{dy}{dt}+ \dfrac{d^2y}{dt^2} \right) \end{align*}

上記の式を(def1)\eqref{def1}に代入すると、次のようになる。

(1cos2t)1sin2t(costsintdydt+d2ydt2)cost(1sint)dydt+n2y=0 (1-\cos ^2t)\dfrac{1}{\sin ^2 t} \left( \dfrac{-\cos t}{\sin t}\dfrac{dy}{dt}+ \dfrac{d^2y}{dt^2} \right)-\cos t \left( -\dfrac{1}{\sin t} \right)\dfrac{dy}{dt} +n^2y=0

整理すると

y+n2y=0 y^{\prime \prime}+n^2y=0

よってTn(cost)T_{n}(\cos t)は上記の微分方程式の解になる。だが、その式はとても簡単な2次微分方程式であり、一般解はy=C1cos(nt)+C2sin(nt)y=C_{1}\cos (nt) + C_2\sin (nt)である。したがって、

Tn(cost)=cos(nt) T_{n}(\cos t)=\cos (nt)

参照