直交多項式と超幾何関数(1)〜チェビシェフ多項式と三角関数〜
いよいよ
この記事を書くときが来てしまった気がした。
近ごろ急に直交多項式を勉強する動機ができたので。
その動機の内容は追々語れる日が来るといいのだが。
ここから連載企画として、「直交多項式と超幾何関数」について記事を書くことにする。
本日は第1回目なので、チェビシェフ多項式の話。
この多項式に関しては本Mathlog内でも充分なほど記事があるが
改めて記事を立てることとする。
意外と(チェビシェフ多項式だけで)証明を完結させるのが大変だったり。
#前記事のベッセル多項式は第0回ってことで
#ベッセル多項式も改めて書き直しますが
古くから
「チェビシェフ多項式は数値解析において至るところ稠密である」
などと言われてきたほど重要な多項式族である。
まずは、それらが顔を出す1番よく知られた状況の説明をする。
$n$倍角の公式
さて、三角関数と言えば、高校生のときに
$\sin\theta$や$\cos\theta$に関する2倍角や3倍角の公式を覚えさせられるのが世の宿命。
\begin{align*}
\sin 2\theta &= 2\sin\theta\cos\theta, &
\sin 3\theta &= 3\sin\theta-4\sin^3\theta \\
\cos 2\theta &= 2\cos^2\theta-1, &
\cos 3\theta &= 4\cos^3\theta-3\cos\theta \\
\end{align*}
2倍角3倍角と言えば私の中学入試を思い出して頭がｳｯｯってなる。笑
($\cos$の2倍角3倍角を頑張って誘導なしに初等的に示せ(意訳)って問題が出たw)
それはさておき、
$n$倍角の公式のようなものを考えるとどうなるのであろうか。
(注意) ここで$\sin$の3倍角に対して
\begin{align*}
\sin 3\theta
&= \sin\theta (3-4\sin^2\theta) \\
&= \sin\theta (4\cos^2\theta-1)
\end{align*}
という変形をしておいた方が今後のために良いと思われる。
さて、$n$倍角の公式を求めるために、和積の公式を使う。すなわち
\begin{align*}
\sin(n+2)\theta + \sin n\theta
&=2\sin(n+1)\theta \cos\theta \\
\cos(n+2)\theta + \cos n\theta
&=2\cos(n+1)\theta \cos\theta \\
\end{align*}
これにより、$\sin n\theta$や$\cos n\theta$の値を帰納的に計算することが可能である。
例えば$n=2$を上の式に代入すると
\begin{align*}
\sin 4\theta +\sin 2\theta = 2\sin 3\theta \cos\theta
\end{align*}
がわかり、ここに上の$\sin 2\theta = 2\sin\theta \cos\theta$や$\sin 3\theta = \sin\theta (4\cos^2\theta-1) $を代入すると
\begin{align*}
\sin 4\theta
&= 2\sin 3\theta \cos\theta - \sin 2\theta \\
&= 2\sin\theta (4\cos^2\theta-1) \cos\theta
- 2\sin\theta \cos\theta \\
&= \sin\theta (8\cos^3\theta- 4\cos\theta)
\end{align*}
などと計算ができる。この計算を繰り返すと一般の$\sin n\theta$及び$\cos n\theta$が計算できる。
…と書いてもわかった気にならない。閉じた形の式が欲しい。
漸化式が出ている状況なので文字を置こう。
チェビシェフ多項式
ここで次のように第一種・第二種チェビシェフ多項式を定める。
$n$は非負整数とする。
次の式を満たす整数係数の多項式族$\{T_n(x)\}, \{U_n(x)\}$が存在する：
\begin{align*}
\cos n\theta = T_n(\cos\theta), \quad
\sin (n+1)\theta = U_n(\cos\theta)\sin\theta
\end{align*}
このとき$\{T_n(x)\}, \{U_n(x)\}$をそれぞれ第一種・第二種チェビシェフ多項式と呼ぶ。
(注)$\sin 0\theta=1$は$\sin\theta$で割れないので、初項が1つずれている。
さてこのように定めたチェビシェフ多項式が満たす漸化式を、先の和積の公式から導く。
(厳密なことを言うと、その漸化式で定義するべき)
例えば$\cos$の方は
\begin{align*}
\cos(n+2)\theta + \cos n\theta = 2 \cos(n+1)\theta\cos\theta
\end{align*}
であったので
\begin{align*}
T_{n+2}(\cos\theta)+T_n(\cos\theta)
&=2T_{n+1}(\cos\theta)\cos\theta \\
T_{n+2}(x)
&=2xT_{n+1}(x)-T_n(x)
\end{align*}
がわかる。一方で$\sin$の方は、和積の式の$n$を1増やしておくと
\begin{align*}
\sin(n+3)\theta + \sin(n+1)\theta
=2\sin(n+2)\theta \cos\theta
\end{align*}
であったので
\begin{align*}
U_{n+2}(\cos\theta)\sin\theta + U_n(\cos\theta)\sin\theta
&=2U_{n+1}(\cos\theta)\sin\theta \cdot \cos\theta \\
U_{n+2}(\cos\theta)
&=2U_{n+1}(\cos\theta)\cos\theta - U_n(\cos\theta) \\
U_{n+2}(x)
&=2xU_{n+1}(x) - U_n(x)
\end{align*}
がわかる。
以上をまとめると、次の定理が従う。
上で定めたチェビシェフ多項式$\{T_n(x)\}, \{U_n(x)\}$は次の漸化式を満たす。
\begin{align*}
\begin{cases}
T_0(x)=1, \,\, T_1(x)=x \\
T_{n+2}(x)=2xT_{n+1}(x)-T_n(x)
\end{cases} \quad
\begin{cases}
U_0(x)=1, \,\, U_1(x)=2x \\
U_{n+2}(x)=2xU_{n+1}(x)-U_n(x)
\end{cases}
\end{align*}
三項間漸化式は複数の直交多項式が満たす性質の1つである。(後の記事で記述予定)
チェビシェフ多項式の数表
第一種・第二種それぞれ最初の10項を掲載する。
＜第一種チェビシェフ多項式＞
\begin{align*}
T_0(x)&= 1 \\
T_1(x)&= x \\
T_2(x)&= 2x^2-1 \\
T_3(x)&= 4x^3-3x \\
T_4(x)&= 8x^4-8x^2+1 \\
T_5(x)&= 16x^5-20x^3+5x \\
T_6(x)&= 32x^6-48x^4+18x^2-1 \\
T_7(x)&= 64x^7-112x^5+56x^3-7x \\
T_8(x)&= 128x^8-256x^6+160x^4-32x^2+1 \\
T_9(x)&= 256x^9-576x^7+432x^5-120x^3+9x
\end{align*}
＜第二種チェビシェフ多項式＞
\begin{align*}
U_0(x)&= 1 \\
U_1(x)&= 2x \\
U_2(x)&= 4x^2-1 \\
U_3(x)&= 8x^3-4x \\
U_4(x)&= 16x^4-12x^2+1 \\
U_5(x)&= 32x^5-32x^3+6x \\
U_6(x)&= 64x^6-80x^4+24x^2-1 \\
U_7(x)&= 128x^7-192x^5+80x^3-8x \\
U_8(x)&= 256x^8-448x^6+240x^4-40x^2+1 \\
U_9(x)&= 512x^9-1024x^7+672x^5-160x^3+10x \\
\end{align*}
さてここからチェビシェフ多項式に関わる性質を紹介する。
チェビシェフ多項式の満たす諸性質
まず残り2つの和積の(と言うより差を積にする)公式を用いることで、次がわかる。
\begin{align*} U_{n+2}(x)-U_n(x)&=2T_{n+2}(x) \\ T_{n+2}(x)-T_n(x)&=2(x^2-1)U_n(x) \end{align*}
この定理が意味するのは、一方がもう一方を用いて表されることである。
(なお後の定理4の証明内で、より使いやすい関係式$\frac{d}{dx}T_n(x)=nU_{n-1}(x)$を得る)
次の定理は、直交多項式の名前で呼ばれる所以の定理である。
$i, j$を非負整数とする。このとき次が成り立つ。
\begin{align*}
\int_{-1}^1 T_i(x)T_j(x)\frac{dx}{\sqrt{1-x^2}}
&=\begin{cases}
\pi & (i=j=0) \\
\frac{\pi}{2}\delta_{ij} & (\text{それ以外})
\end{cases} \\
\int_{-1}^1 U_i(x)U_j(x)\sqrt{1-x^2}dx
&=\frac{\pi}{2}\delta_{ij}
\end{align*}
ここで$\delta_{ij}$はクロネッカーのデルタであり、$i$と$j$が等しい時に1、異なる時に0を返す関数である。
また$T_n(x)$に対する$\frac{1}{\sqrt{1-x^2}}$、$U_n(x)$に対する$\sqrt{1-x^2}$を重み関数と呼ぶ。
この定理は三角関数の積分に帰着させることで証明する。
まず最初に上の$T_n(x)$に関わる積分から示す。
$x=\cos\theta$とおくことで
\begin{align*}
\int_{-1}^1 T_i(x)T_j(x)\frac{dx}{\sqrt{1-x^2}}
=\int_\pi^0 \cos i\theta \cos j\theta \,
\frac{-\sin\theta d\theta}{\sin\theta}
=\int_0^\pi \cos i\theta \cos j\theta \, d\theta
\end{align*}
三角関数の積和の公式を用いて計算する。以下3つに分けて計算する。
・ $i=j=0$のとき
\begin{align*}
\int_0^\pi \cos i\theta \cos j\theta \, d\theta
&=\int_0^\pi d\theta = \pi
\end{align*}
・ $i=j\ge1$のとき
\begin{align*}
\int_0^\pi \cos i\theta \cos j\theta \, d\theta
=\int_0^\pi \cos^2 i\theta \, d\theta
=\int_0^\pi \frac{1+\cos 2i\theta}{2} \, d\theta
=\left[\frac{\theta+\frac{1}{2i}\sin 2i\theta}{2}\right]_0^\pi
=\frac{\pi}{2}
\end{align*}
・$i\neq j$のとき
\begin{align*}
\int_0^\pi \cos i\theta \cos j\theta \, d\theta
&=\frac{1}{2}\int_0^\pi \left\{\cos (i+j)\theta - \cos (i-j)\theta\right\} \, d\theta \\
&=\frac{1}{2}\left[\frac{\sin (i+j)\theta}{i+j}-\frac{\sin (i-j)\theta}{i-j}\right]_0^\pi=0
\end{align*}
となり計算ができる。
次に$U_n(x)$の方も同じく、$x=\cos\theta$とおくと
\begin{align*}
\int_{-1}^1 U_i(x)U_j(x)\sqrt{1-x^2}dx
&=\int_{\pi}^0
\frac{\sin (i+1)\theta}{\sin\theta}
\frac{\sin (j+1)\theta}{\sin\theta}
\sin\theta \cdot (-\sin\theta) d\theta \\
&=\int_0^{\pi}
\sin (i+1)\theta\sin (j+1)\theta \, d\theta
\end{align*}
となり、これも積和の公式を用いることで同様の計算ができる。（証明終わり）
note: チェビシェフ多項式は、これら重み関数に関する内積による多項式の空間の直交基底とみなせる
(シュミットの直交化法、これも後の記事で)
次の微分方程式、及びその次のロドリゲスの公式も直交多項式論においてはかなり重要な役目を担う。
$T_n(x)$や$U_n(x)$は次の二階微分方程式を満たしている。
\begin{align*}
(1-x^2)T_n''(x)-xT_n'(x)+n^2T_n(x)&= 0 \\
(1-x^2)U_n''(x)-3xU_n'(x)+n(n+2)U_n(x)&= 0 \\
\end{align*}
(note: それぞれの微分方程式の一般解は、前者は$T_n(x)$と$\sqrt{1-x^2}U_{n-1}(x)$の一次結合、後者は$U_n(x)$と$\frac{1}{\sqrt{1-x^2}}T_{n+1}(x)$の一次結合で書けることもわかる。証明略)
かなりゴリゴリに計算をする。と言うか定義から計算可能（もう少しいい証明があるかも）
まず$T_n(\cos\theta) = \cos n\theta$の両辺を$\theta$で微分する。
\begin{align*}
\frac{d}{d\theta}T_n(\cos\theta)
&=\frac{d}{d\theta}\cos n\theta \\
&=-n\sin n\theta \\
\end{align*}
ここで$\frac{d(\cos\theta)}{d\theta}=-\sin\theta$を使うと
\begin{align*}
\frac{d}{d(\cos\theta)}(-\sin\theta)T_n(\cos\theta)
&=-n\sin n\theta \\
\frac{d}{dx}T_n(x)&=n\frac{\sin n\theta}{\sin \theta} \\
&=nU_{n-1}(x)
\end{align*}
となり、$T_n(x)$の一階微分が求まった。
次に$U_n(x)$の方はも同様に計算する。
$U_n(\cos\theta)=\frac{\sin(n+1)\theta}{\sin\theta}$より
\begin{align*}
(-\sin\theta)\frac{d}{dx}U_n(x)
&=\frac{(n+1)\cos(n+1)\theta\sin\theta-\sin(n+1)\theta\cos\theta}{\sin^2\theta} \\
\frac{d}{dx}U_n(x)
&=-\frac{(n+1)\cos(n+1)\theta\sin\theta-\sin(n+1)\theta\cos\theta}{\sin\theta\sin^2\theta} \\
&=\frac{(n+1)T_{n+1}(x)-U_n(x)x}{x^2-1}
\end{align*}
と商の微分公式を用いることで計算ができる。
次に2階微分だが、$T$の方はすぐに求まる。
\begin{align*}
\frac{d^2}{dx^2}T_n(x)
=\frac{d}{dx}nU_{n-1}(x)
=\frac{n\{nT_n(x)-xU_{n-1}(x)\}}{x^2-1}
\end{align*}
以上より$T$の方の求める微分方程式の左辺は
\begin{align*}
&(1-x^2)\frac{d^2}{dx^2}T_n(x)-x\frac{d}{dx}T_n(x)+n^2T_n(x) \\
&\quad=
-n\{nT_n(x)-xU_{n-1}(x)\}-xnU_{n-1}(x)+n^2T_n(x)=0
\end{align*}
となり、$T_n(x)$に関する微分方程式が正しいことが示された。
次に計算の恐ろしい$U$の2階微分だが、それは計算したくないので工夫をする。
本定理の証明内で示した$\frac{d}{dx}T_n(x)=nU_{n-1}(x)$及び$T_n(x)$に関する微分方程式を使う。
すなわち、微分方程式
\begin{align*}
(1-x^2)\frac{d^2}{dx^2}T_n(x)-x\frac{d}{dx}T_n(x)+n^2T_n(x)=0
\end{align*}
を1回$x$で微分することを考える。すると
\begin{align*}
\left\{(-2x)\frac{d^2}{dx^2}T_n(x)+(1-x^2)\frac{d^3}{dx^3}T_n(x)\right\}
-\left\{\frac{d}{dx}T_n(x)+x\frac{d^2}{dx^2}T_n(x)\right\}
+n^2\frac{d}{dx}T_n(x)=0
\end{align*}
がわかる。整理すると次を得る。
\begin{align*}
(1-x^2)\frac{d^3}{dx^3}T_n(x)
-3x\frac{d^2}{dx^2}T_n(x)
+(n^2-1)\frac{d}{dx}T_n(x)=0
\end{align*}
ここで$\frac{d}{dx}T_n(x)=nU_{n-1}(x)$を用いて書き直すと
\begin{align*}
(1-x^2)\frac{d^2}{dx^2}nU_{n-1}(x)
-3x\frac{d}{dx}nU_{n-1}(x)
+(n^2-1)nU_{n-1}(x)=0
\end{align*}
$n$で割り、$n$を$1$ずらすことで
\begin{align*}
(1-x^2)\frac{d^2}{dx^2}U_n(x)
-3x\frac{d}{dx}U_n(x)
+n(n+2)U_n(x)=0
\end{align*}
が得られるので$U$に関する微分方程式も示された。（証明終わり）
続いてチェビシェフ多項式の一般項を与える公式の1つ、ロドリゲスの公式である。
$T_n(x)$や$U_n(x)$は次の$n$階微分を用いた表示を持つ：
\begin{align*}
T_n(x)
&=\frac{(-1)^n}{(2n-1)!!}\sqrt{1-x^2}
\frac{d^n}{dx^n}(1-x^2)^{n-\frac{1}{2}} \\
U_n(x)
&=\frac{(-1)^n2^n(n+1)!}{(2n+1)!}\frac{1}{\sqrt{1-x^2}}
\frac{d^n}{dx^n}(1-x^2)^{n+\frac{1}{2}}
\end{align*}
証明は上の微分方程式を利用して行う。(Koekoek pp.62-64による)
もちろんこちらを先に直接示すこともできる；その場合微分方程式の証明は簡略になるが、こちらの証明がやや大変になる。(後の記事を考えて、その方針は別記事で証明を行う。)
まず、次の式を示す。[これの$k=1$の時を前定理の証明内で示した]
\begin{align*}
(1-x^2)\frac{dT_n^{k+2}}{dx^{k+2}}(x)
-(2k+1)x\frac{dT_n^{k+1}}{dx^{k+1}}(x)
=(k^2-n^2)\frac{dT_n^k}{dx^k}(x)
\end{align*}
帰納法で簡単に示せる。$k=0$の時は
\begin{align*}
(1-x^2)T_n''(x)-xT_n'(x)=-n^2T_n(x)
\end{align*}
でありこれは$T_n$の満たす微分方程式そのもの。
さて$k-1$の時の式
\begin{align*}
(1-x^2)\frac{dT_n^{k+1}}{dx^{k+1}}(x)
-(2k-1)x\frac{dT_n^k}{dx^k}(x)
=((k-1)^2-n^2)\frac{dT_n^{k-1}}{dx^{k-1}}(x)
\end{align*}
を$x$で微分すると、左辺は
\begin{align*}
&(-2x)\frac{dT_n^{k+1}}{dx^{k+1}}(x)
+(1-x^2)\frac{dT_n^{k+2}}{dx^{k+2}}(x)
-(2k-1)\frac{dT_n^k}{dx^k}(x)
-(2k-1)x\frac{dT_n^{k+1}}{dx^{k+1}}(x) \\
&\quad=
(1-x^2)\frac{dT_n^{k+2}}{dx^{k+2}}(x)
-(2k+1)x\frac{dT_n^{k+1}}{dx^{k+1}}(x)
-(2k-1)\frac{dT_n^k}{dx^k}(x)
\end{align*}
である。最後の$k$階微分の項を右辺に回すと右辺は
\begin{align*}
&(2k-1)\frac{dT_n^k}{dx^k}(x)
+((k-1)^2-n^2)\frac{dT_n^k}{dx^k}(x)
=(k^2-n^2)\frac{dT_n^k}{dx^k}(x)
\end{align*}
となり、示された。
今示した式を使ってロドリゲスの公式を示す。
天下り的だが次の式を計算する：
\begin{align*}
\frac{d}{dx}\left\{
(1-x^2)^{k+\frac{1}{2}}
\frac{d^{k+1}}{dx^{k+1}}T_n(x)
\right\}
\end{align*}
積の微分法を用いて計算すると次のように今示した式が使えて計算ができる。
\begin{align*}
&\frac{d}{dx}\left\{
(1-x^2)^{k+\frac{1}{2}}
\frac{d^{k+1}}{dx^{k+1}}T_n(x)
\right\} \\
&\quad=
\left(k+\frac{1}{2}\right)(-2x)(1-x^2)^{k-\frac{1}{2}}
\frac{d^{k+1}}{dx^{k+1}}T_n(x)
+(1-x^2)^{k+\frac{1}{2}}\frac{d^{k+2}}{dx^{k+2}}T_n(x) \\
&\quad=
(1-x^2)^{k-\frac{1}{2}}
\left\{
-(2k+1)x\frac{d^{k+1}}{dx^{k+1}}T_n(x)
+(1-x^2)\frac{d^{k+2}}{dx^{k+2}}T_n(x)
\right\} \\
&\quad=
(k^2-n^2)(1-x^2)^{k-\frac{1}{2}}
\frac{d^k}{dx^k}T_n(x)
\end{align*}
すなわち$g_k(x)=(1-x^2)^{k-\frac{1}{2}}\frac{d^k}{dx^k}T_n(x)$とおくと$\frac{d}{dx}g_{k+1}(x)=(k^2-n^2)g_k(x)$
という簡単な式が導かれる。
$g_n(x)$を計算するために、$\frac{d^n}{dx^n}T_n(x)$を計算する必要がある。
$n$次多項式の$n$階微分なので、最高次係数の分だけが残り
\begin{align*}
\frac{d^n}{dx^n}T_n(x)
&=\frac{d^n}{dx^n}(2^{n-1}x^n+(\text{$n-1$ 次以下の項})) \\
&=2^{n-1}n!
\end{align*}
と計算ができる。以上を合わせると
\begin{align*}
T_n(x)
&=(1-x^2)^{\frac{1}{2}}g_0(x) \\
&=(1-x^2)^{\frac{1}{2}}\frac{1}{-n^2}\frac{d}{dx}g_1(x) \\
&=(1-x^2)^{\frac{1}{2}}\frac{1}{-n^2(1-n^2)}\frac{d^2}{dx^2}g_2(x) \\
&=\cdots \\
&=(1-x^2)^{\frac{1}{2}}\frac{1}{-n^2(1-n^2)(4-n^2)\cdots((n-1)^2-n^2)}\frac{d^n}{dx^n}g_n(x) \\
&=(1-x^2)^{\frac{1}{2}}\frac{2^{n-1}n!}{-n^2(1-n^2)(4-n^2)\cdots((n-1)^2-n^2)}\frac{d^n}{dx^n}(1-x^2)^{n-\frac{1}{2}} \\
&=\frac{(-1)^n}{(2n-1)!!}\sqrt{1-x^2}
\frac{d^n}{dx^n}(1-x^2)^{n-\frac{1}{2}}
\end{align*}
がわかる。
$U_n(x)$に関しても証明は全く同様なので省略する。（証明終わり）
このロドリゲスの公式もかなり良いが、具体的に$x$の何乗の係数が幾つであるかは$n$階微分をしないとわからない状況である。
今度はその係数を求めることを考える。
チェビシェフ多項式の一般項と超幾何関数
もっと一般的に書くためには、三項間漸化式を解けば良い。
\begin{align*}
T_{n+2}(x)-2xT_{n+1}(x)+T_n(x)&=0 \\
U_{n+2}(x)-2xU_{n+1}(x)+U_n(x)&=0
\end{align*}
と同じ漸化式を満たしている点に注意する。(初項・第二項の条件が違うだけ)
高校数学である通り、特性方程式を解くことで実際に解を与えることができる。
特性方程式の解は$\lambda=x\pm\sqrt{x^2-1}$となるので、この漸化式は
\begin{align*}
C_1(x)\left(x+\sqrt{x^2-1}\right)^n
+C_2(x)\left(x-\sqrt{x^2-1}\right)^n
\end{align*}
の形の一般項を持つことがわかる。初項・第二項の条件で係数を比較して
\begin{align*}
\begin{cases}
\displaystyle T_n(x)
=\frac{1}{2}\left\{\left(x+\sqrt{x^2-1}\right)^n
+\left(x-\sqrt{x^2-1}\right)^n\right\} \\
\\
\displaystyle U_n(x)
=\frac{1}{2\sqrt{x^2-1}}
\left\{\left(x+\sqrt{x^2-1}\right)^{n+1}
-\left(x-\sqrt{x^2-1}\right)^{n+1}\right\}
\end{cases}
\end{align*}
という形の一般項を得ることができる。(もう少し変形したいので定理にはしない)
しかし、このまま二項定理を用いて変形しても、係数は二項係数の積の和というあまり綺麗な形にならず、そこで行き詰まる。例えば具体的には
\begin{align*}
T_n(x)=\sum_{k=0}^{[n/2]}
\left\{
(-1)^k\sum_{j=k}^{[n/2]}\binom{n}{2j}\binom{j}{k}
\right\}x^{n-2k}
\end{align*}
とまでは変形が容易にできるが、ここから先同じ結論に辿り着くために
\begin{align*} \sum_{j=k}^{[n/2]}\binom{n}{2j}\binom{j}{k} =2^{n-2k-1}n\frac{(n-k-1)!}{k!(n-2k)!} \end{align*}
という等式を示さねばならず、ここで詰まってしまった。
（ここまで読んで頂いた方で、ご存じな方がおられましたらコメント等書いていただけると有難いです。）
note: Gaussの${}_2F_1$和公式から示せるのは確認できた。本末転倒な気しかしないが。
note:
TKSS
さんが$k$重和と望遠鏡和を用いた
すごく見た目の美しい証明記事を投稿して下さりました。
ありがとうございます。記事はこちら↓↓↓
チェビシェフ多項式の係数計算で現れる二項係数の和の計算について
($U_n(x)$に関しても同様の別の二項係数の積の和を考えねばならないが、上記事のやり方で解決する)
別の方針
さて、方針を1から変えよう。
三項間漸化式の解き方は他にもある。
特に三角関数の対称性からもわかるが、$T_n(x), U_n(x)$ともに$n$が偶数の時は偶関数、奇数の時は奇関数である。
そのため、
\begin{align*}
T_n(x)=\sum_{k=0}^{[n/2]}c_k^{(n)}x^{n-2k}
\end{align*}
とおき(先の中括弧の中を文字で置いた)これに関する漸化式として見たい。
$T_n(x)$の漸化式でも良いが、定理4の微分方程式
\begin{align*}
(1-x^2)T_n''(x)-xT_n'(x)+n^2T_n(x)=0
\end{align*}
を使った方が($n$は変化しないので)見通しが立ちやすい。
$x^{n-2k}$の係数を取り出すと
\begin{align*}
(n-2k+2)(n-2k+1)c_{k-1}^{(n)}
-(n-2k)(n-2k-1)c_k^{(n)}
-(n-2k)c_k^{(n)}+n^2c_k^{(n)}=0
\end{align*}
整理して
\begin{align*}
c_k^{(n)}=-\frac{(n-2k+2)(n-2k+1)}{4k(n-k)}c_{k-1}^{(n)}
\end{align*}
及び初項$c_0^{(n)}=2^{n-1}$であるので$c_k^{(n)}$は
\begin{align*}
c_k^{(n)}
&=-\frac{(n-2k+2)(n-2k+1)}{4k(n-k)}c_{k-1}^{(n)} \\
&=\frac{(n-2k+2)(n-2k+1)}{4k(n-k)}\frac{(n-2k+4)(n-2k+3)}{4(k-1)(n-k+1)}c_{k-2}^{(n)} \\
&=\cdots \\
&=(-1)^k\frac{(n-2k+1)(n-2k+2)\cdots(n-1)n}
{4^kk(k-1)\cdots 1(n-k)(n-k+1)\cdots(n-1)}c_0^n \\
&=(-1)^k2^{n-1-2k}n\frac{(n-k-1)!}{k!(n-2k)!}
\end{align*}
これを用いると次のように書き表せる。
\begin{align*}
T_n(x)
&=\sum_{k=0}^{[n/2]}
(-1)^k2^{n-1-2k}n\frac{(n-k-1)!}{k!(n-2k)!}x^{n-2k} \\
&=\frac{n}{2}\sum_{k=0}^{[n/2]}
(-1)^k\frac{(n-k-1)!}{k!(n-2k)!}(2x)^{n-2k}
\end{align*}
ここまでで十分綺麗な形にはなった。なお$U$の方も同様にすると
\begin{align*}
U_n(x)
=\sum_{k=0}^{[n/2]}
(-1)^k\frac{(n-k)!}{k!(n-2k)!}(2x)^{n-2k}
\end{align*}
であることが確かめられる。
しかし次の変形を考える：2倍角の公式を考えると
\begin{align*}
T_{2n}(x)
&=T_{2n}(\cos\theta) \\
&=\cos 2n\theta \\
&=T_n(\cos 2\theta)=T_n(2x^2-1)
\end{align*}
と書けるので、上の式で$T_n(2x^2-1)$を考えると
\begin{align*}
T_n(2x^2-1)&=T_{2n}(x) \\
&=n\sum_{k=0}^{n}
(-1)^k\frac{(2n-k-1)!}{k!(2n-2k)!}(2x)^{2n-2k} \\
&=n\sum_{k=0}^{n}
(-1)^k\frac{(2n-k-1)!}{k!(2n-2k)!}2^{2n-2k}(x^2)^{n-k}
\end{align*}
となる。最後に$2x^2-1$を$x$と置き直すと
\begin{align*}
T_n(x)
&=n\sum_{k=0}^{n}
(-1)^k\frac{(2n-k-1)!}{k!(2n-2k)!}2^{2n-2k}
\left(\frac{x+1}{2}\right)^{n-k} \\
&=n\sum_{k=0}^{n}
(-1)^{n-k}\frac{(n+k-1)!}{(n-k)!(2k)!}2^k(x+1)^k
\quad\text{($k$ を $n-k$ と取り替えた)} \\
&=n\sum_{k=0}^{n}
(-2)^k\frac{(n+k-1)!}{(n-k)!(2k)!}(1-x)^k
\quad\text{($x$ を $-x$ と取り替えた)}
\end{align*}
ここで最後の式では関数の偶奇性で($n$が奇数の時は奇関数で)$(-1)^n$を掛けたことに注意されたい。
$U$に関しては、同様の流れを追うよりは
今わかった$T_n(x)$の一般項に、
定理4の証明中で用いた$\frac{d}{dx}T_n(x)=nU_{n-1}(x)$を用いると計算が楽だろう。
これらを合わせて、次を得る。
チェビシェフ多項式$T_n(x)$及び$U_n(x)$は次の表示がある。
\begin{align*}
T_n(x)
&=n\sum_{k=0}^n(-2)^k\frac{(n+k-1)!}{(n-k)!(2k)!}(1-x)^k
=\sum_{k=0}^n\frac{(n)_k(-n)_k}{(1/2)_k(1)_k}\left(\frac{1-x}{2}\right)^k \\
U_n(x)
&=\sum_{k=0}^n(-2)^k\frac{(n+k+1)!}{(n-k)!(2k+1)!}(1-x)^k
=(n+1)\sum_{k=0}^n\frac{(n+2)_k(-n)_k}{(3/2)_k(1)_k}\left(\frac{1-x}{2}\right)^k
\end{align*}
ここで$(a)_k=a(a+1)\cdots(a+k-1)=\frac{(a+k-1)!}{(a-1)!}$は上昇階乗である。
note: 上の変形で用いた上昇階乗の値を幾つか説明をつける。
\begin{align*}
(1)_k
&=1\cdot2\cdots n=n! \\
(-n)_k
&=(-n)(-n+1)\cdots(-n+k-1)
=(-1)^k\frac{n!}{(n-k)!} \\
(1/2)_k
&=\frac{1}{2}\frac{3}{2}\cdots\frac{2k-1}{2}
=\frac{(2k-1)!!}{2^k}
=\frac{(2k)!}{2^{2k}k!}
=\frac{(2k-1)!}{2^{2k-1}(k-1)!} \\
(3/2)_k
&=\frac{(2k+1)!}{2^{2k}k!}
\end{align*}
以上からチェビシェフ多項式は、上昇階乗のいくつかの比を係数にもつべき級数として書き表すことができる。その形の級数を超幾何関数と呼ぶ。（後の記事で改めて定義します）
超幾何関数は様々な直交多項式や初等関数・特殊関数を包含する。さらに波動関数・場の理論をはじめとする物理学、量子力学、統計学、確率論、と応用例を挙げると枚挙にいとまがない。
まとめ
第1回目の本記事では、三角関数の$n$倍角の公式から自然に考えられるチェビシェフ多項式を例に挙げ、それらがもたらす諸性質について述べた。次の記事以降では、直交多項式の言葉を定義しそれらの性質について述べ、また様々な直交多項式の例を挙げようと思う。
# いつ超幾何関数の議論ができる日が来るのか…
参考文献
この記事を高評価した人
この記事に送られたバッジ
投稿者
