6

チェビシェフ多項式のいろいろ【工事中】

963
1
$$$$

あいさつ。

おひさーです。急にチェビシェフ多項式が気になったので色々調べて書いてみます。
急ピッチで建設してるんですが、なかなか進まないので一旦投稿してしまえ!ということで投稿しました。随時、目次や、内容を追加します。

更新履歴

1.2023/05/18 10:38 投稿
2.2023/05/23 14:21 チェビシェフ多項式の閉じた形について式の証明を追加
3.2023/05/23 21:37 チェビシェフ多項式一覧(0<=n<=10)を追加
4.2023/05/23 23:27 直交性の部分を若干追加+etc
5.2023/05/26 15:39 第二種チェビシェフ多項式についての公式などを追加
6.2023/06/19 11:28 直行性、特殊値の部分を若干追加

目次

1.チェビシェフ多項式とは
2.チェビシェフ多項式の閉じた形
3.母関数について
4.チェビシェフ多項式一覧(0<=n<=10)
5.チェビシェフ多項式の直交性
6.チェビシェフ多項式の特殊値
7.チェビシェフ多項式が満たす微分方程式
8.にょー

チェビシェフ多項式とは

第一種チェビシェフ多項式(Chebyshev polynomials of the first kind)

$\cos n\theta$$\cos \theta$$n$次多項式で表せる。

この$n$次多項式を第一種チェビシェフ多項式といい、$T_{n}(x)$と表す。つまり$T_{n}(\cos\theta)=\cos n\theta$である。

加法定理で帰納的に示す。

加法定理より、
\begin{align} \cos(k+1)\theta &= \cos k\theta\cos\theta-\sin k\theta\sin\theta\\ \cos(k+2)\theta &= \cos k\theta\cos 2\theta-\sin k\theta\sin2\theta\\ 第一式の両辺に2\cos \thetaをかけて\\ 2\cos \theta\cos (k+1)\theta &=2\cos k\theta\cos^{2}\theta-2\sin k\theta\sin \theta\cos \theta\\ &= 2\cos k\theta\cos^{2}\theta-\sin k\theta\sin 2\theta\\ 第二式と差をとって\\ \cos (k+2)\theta-2\cos\theta \cos (k+1)\theta &=\cos k\theta \cos 2\theta -2\cos k\theta \cos ^{2}\theta\\ &= (2\cos^{2}\theta - 1)\cos k\theta - 2\cos k\theta \cos^{2}\theta\\ &=-\cos k\theta\\ \\ \therefore\cos (k+2)\theta &= 2\cos \theta \cos (k+1)\theta - \cos k\theta \end{align}

を得る。$\cos 0\theta = 1,\cos 1\theta = \cos \theta$より$k$について帰納的に、$\cos k\theta$$\cos \theta$の多項式であることが従う。

第一種チェビシェフ多項式の漸化式

上で得た漸化式を$T_n$を用いて書くと、$T_{n+2}(\cos \theta)=2\cos \theta T_{n+1}(\cos \theta)-T_{n}(\cos \theta)$
$x\in \mathbb{R}$についての式にすることで第一種チェビシェフ多項式の漸化式$T_{n+2}(x)=2xT_{n+1}(x)-T_{n}(x)$を得る。

第二種チェビシェフ多項式(Chebyshev polynomials of the second kind)

$\sin n\theta$$\sin \theta$$\cos \theta$$n-1$次多項式の積で表せる。

$U_{n-1}(\cos \theta) =\displaystyle \frac{\sin n\theta}{\sin \theta}$を満たす$U_{n}$を第二種チェビシェフ多項式という。

定理1の証明と同様である。

\begin{align} \sin (k+2)\theta - 2\cos \theta \sin (k+1)\theta &= \sin k\theta \cos 2\theta -2\sin k\theta \cos^{2} \theta \\ &=-\sin k\theta\\ \therefore \sin (k+2)\theta &= 2\cos \theta \sin (k+1)\theta -\sin k\theta \end{align}
を得る。$\sin 0\theta = 0,\sin 1\theta = \sin \theta$より$k$について帰納的に、$\sin k\theta$$\sin \theta$$\cos \theta$$n-1$次多項式であることが従う。

第二種チェビシェフ多項式の漸化式

上で得た漸化式を$U_n$を用いて書くと、$U_{n+2}(\cos \theta)=2\cos \theta U_{n+1}(\cos \theta)-U_{n}(\cos \theta)$
$x\in \mathbb{R}$についての式にすることで第二種チェビシェフ多項式の漸化式$U_{n+2}(x)=2xU_{n+1}(x)-U_{n}(x)$を得る。

チェビシェフ多項式の閉じた形

第一種チェビシェフ多項式の閉じた形

\begin{align} T_n(x)&=\left\{ \begin{alignedat}{2} \frac{1}{2}& \left[\left( x+ i\sqrt{1-x^2} \right)^{n} + \left( x- i\sqrt{1-x^2} \right)^{n}\right] \qquad \qquad &\mathrm{for}\quad |x|\leq 1\\ \\ \frac{1}{2}& \left[ \left( x+ \sqrt{x^2-1} \right)^{n} + \left( x- \sqrt{x^2-1} \right)^{n} \right] \qquad \qquad &\mathrm{for}\quad |x|\geq 1 \end{alignedat} \right. \\ \\ &=\left\{ \begin{alignedat}{3} &\cos (n \arccos x)\qquad \qquad &\mathrm{for}&\quad -1\leq x \leq 1 \\ \\ &\cosh (n\ \mathrm{arcosh}\ x) \qquad \qquad & \mathrm{for}&\quad 1\leq x \\ \\ &(-1)^{n}\cosh (n\ \mathrm{arcosh}(-x))\qquad \qquad &\mathrm{for}&\quad x\leq -1 \end{alignedat} \right. \end{align}

第一式を示す。
ド・モアブルの定理を用いる。

$x=\cos \theta$と置くと
\begin{align} T_n(x)+i\sin n\theta&=(\cos \theta +i\sin \theta)^n=\left(x + i\sqrt{1 -x^2}\right)^n\\ T_n(x)-i\sin n\theta &= (\cos \theta - i\sin \theta)^n=\left(x-i\sqrt{1-x^2}\right)^n \end{align}
より第一式
$\displaystyle T_n(x)=\frac{1}{2}\left[\left(x+i\sqrt{1-x^2}\right)^n+\left(x-i\sqrt{1-x^2}\right)^n\right]\qquad \qquad \mathrm{for}\quad |x|\leq 1$
を得る。

第二式を示す。

第一種チェビシェフ多項式の漸化式
$T_{n+2}(x)=2xT_{n+1}(x)-T_{n}(x)$を用いて示す。

この三項間漸化式の特性方程式は
$t^2 - 2xt + 1 = 0$となりその解は$t=x\pm \sqrt{x^2 - 1}$
よってこの漸化式は次のように書ける。

\begin{align} T_{n+2}(x)-AT_{n+1}(x) &= B(T_{n+1}(x)-AT_{n}(x))\\ T_{n+2}(x)-BT_{n+1}(x) &= A(T_{n+1}(x)-BT_{n}(x))\\ \end{align}
ただし、$A=x+\sqrt{x^2 - 1},B=x-\sqrt{x^2 - 1}$
よって上の式から$T_{n+1}(x)-AT_n(x)$は公比$B$の等比数列となる。
したがって、
$T_{n+1}(x)-AT_n(x)=B^n(T_1(x)-AT_0(x))=B^n(x-A)=-B^n\sqrt{x^2-1}$
同様に下の式から
$T_{n+1}(x)-BT_n(x)=A^n\sqrt{x^2-1}$
この二式の差をとって
$AT_n(x)-BT_n(x)=A^n\sqrt{x^2-1}+B^n\sqrt{x^2-1}$
$2\sqrt{x^2-1}T_n(x)=\sqrt{x^2-1}(A^n+B^n)$

$\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]$
を得る。

第三式$T_n(x)=\cos (n\arccos x)\qquad \mathrm{for}\quad -1 \leq x \leq 1 $
を示す。
逆関数の定義から明らかであるが、私自身の練習のためちょっとした証明を与える。

定義より、$T_n(\cos \theta)=\cos(n\theta)$である。
$\cos\theta$$[0,\pi]$から$ [-1,1]$へ全単射である。なので、
$\forall x \in [-1,1],\exists \theta \in [0,\pi](x=\cos\theta)$となる。
よって、この$x$に対応する$\theta$$\arccos x$とかけば、
第三式$T_n(x)=\cos (n\arccos x)\qquad \mathrm{for}\quad -1 \leq x \leq 1 $
を得る。

第四式$T_n(x)=\cosh (n\ \mathrm{arcosh}\ x) \qquad \mathrm{for}\quad 1\leq x$

を示す。

$\cos(i\theta)=\cosh(\theta)$より
$T_n(\cosh \theta)=\cos (ni\theta)=\cosh(n\theta)$
また、$\cosh\theta$はその定義
$\displaystyle \cosh \theta = \frac{e^\theta + e^{-\theta}}{2}$より
$[0,\infty)\to[1,\infty)$で全単射である。
なので$\forall x \in [1,\infty),\exists \theta \in [0,\infty)(x=\cosh \theta)$である。
この$x$に対応する$\theta$$\mathrm{arcosh}\ x$とかくと、

第四式$T_n(x)=\cosh (n\ \mathrm{arcosh}\ x) \qquad \mathrm{for}\quad 1\leq x$
を得る。

第五式$T_n(x)=(-1)^{n}\cosh (n\ \mathrm{arcosh}(-x))\qquad \mathrm{for} \quad x\leq -1$
を示す。

$T_n(-x)=(-1)^nT_n(x)$(チェビシェフ多項式の特殊値を参照)と第四式より

第五式$T_n(x)=(-1)^{n}\cosh (n\ \mathrm{arcosh}(-x))\qquad \mathrm{for} \quad x\leq -1$
を得る。

第二種チェビシェフ多項式の閉じた形

\begin{align} U_n(x)=\left\{ \begin{alignedat}{2} \frac{1}{2i\sqrt{1-x^2}}& \left[\left( x+ i\sqrt{1-x^2} \right)^{n+1} - \left( x- i\sqrt{1-x^2} \right)^{n+1}\right] \qquad \qquad &\mathrm{for}\quad |x|\leq 1\\ \\ \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] \qquad \qquad &\mathrm{for}\quad |x|\geq 1 \end{alignedat} \right. \end{align}

公式3の証明と同様である。

$\displaystyle T_n(x)=\sum_{k=0}^{\lfloor n/2\rfloor}\dbinom{n}{2k}(x^2-1)^k x^{n-2k}=n\sum_{k=0}^{n}(-2)^k\frac{(n-k-1)!}{(n-k)!(2k)!}(1-x)^k=\, _2F_1(-n,n;1/2;(1-x)/2)$

$_2F_1$はガウスの超幾何関数でPochhammer記号$(・)_n$を用いて


$\displaystyle _2F_1 (a,b;c;z)=\sum_{n=0}^{\infty}\frac{(a)_n (b)_n}{(c)_n}\frac{z^n}{n!}$

と定義され、$a,b$どちらかが負の時$_2F_1$は有限和で

$\displaystyle _2F_1 (-m,b;c;z)=\sum_{n=0}^{m}(-1)^{n}\dbinom{m}{n}\frac{(b)_n}{(c)_n}z^n$

$T_n(\cos \theta)+i\sin(n\theta) = (\cos \theta + i\sin \theta)^n$
二項定理より、
$\displaystyle \sum_{k=0}^{n}\dbinom{n}{k}(i\sin \theta)^k(\cos \theta)^{n-k}$
$k:\mathrm{odd}$ならば虚数となるので$k:\mathrm{even}$の時だけ考えればよい
従って、
$\displaystyle \sum_{k=0}^{\lfloor n/2 \rfloor}\dbinom{n}{2k}(-\sin^2 \theta)^{k} (\cos \theta)^{n-2k}$
$x=\cos\theta$とおけば
$=\displaystyle \sum_{k=0}^{\lfloor n/2 \rfloor}\dbinom{n}{2k}(x^2-1)^{k} x^{n-2k}$

\begin{align*} T_n(x) &= \sum_{k=0}^{\lfloor n/2 \rfloor}\dbinom{n}{2k}(x^2 - 1)^{k}x^{n-2k}\\ &= x^n \sum_{k=0}^{\lfloor n/2 \rfloor}\dbinom{n}{2k}(1-x^{-2})^k\\ &= \frac{n}{2}\sum_{k=0}^{\lfloor n/2 \rfloor}(-1)^{k}\frac{(n-k-1)!}{k!(n-2k)!}(2x)^{n-2k}\qquad &\mathrm{for}\ n > 0\\ &= n\sum_{k=0}^{n}(-2)^k \frac{(n+k-1)!}{(n-k)!(2k)!}(1-x)^k \qquad &\mathrm{for}\ n > 0\\ &=\, _2F_1 (-n,n;1/2;(1-x)/2) \end{align*}

$\displaystyle U_n(x)=\sum_{k=0}^{\lfloor n/2\rfloor}(-1)^k \dbinom{n-k}{k}(2x)^{n-2k}\qquad \qquad \mathrm{for}\quad n > 0$

母関数について

通常型母関数

\begin{equation*} F(t):=\sum_{n=0}^{\infty}T_{n}(x)t^{n}=T_{0}(x)+T_{1}(x)t+T_{2}(x)t^{2}+\cdots \end{equation*}
とする。
\begin{array}{cc} F(t)&=&T_{0}(x)&+&T_1(x)t &+& T_2(x)t^2 &+& T_3(x)t^3 &+& \cdots\\ -2xtF(t) &=& & &-2xT_0(x)t &+& -2xT_1(x)t^2 &+& -2xT_2(x)t^3 &+& \cdots\\ t^2F(t) &=& & & & & T_0(x)t^2 &+& T_1(x)t^3 &+& \cdots \end{array}
三つの式を足して
$(1-2xt+t^2)F(t)=T_0(x)+(T_1(x) - 2xT_0(x))t + (T_2(x) -2xT_1(x) +T_0(x))t^2 + (T_3(x) -2xT_2(x) + T_1(x))t^3+\cdots$
となる。ここで$T_n$の漸化式を用いると、$t^2$以上の項が$0$となり、
\begin{align} (1-2xt+t^2)F(t)&=T_0(x)+(T_1(x)-2xT_0(x))t\\ &=1-tx \end{align}
よって、次の公式を得る。

第一種チェビシェフ多項式の通常型母関数

$\displaystyle \sum_{n=0}^{\infty} T_n(x)t^{n}=\frac{1-tx}{1-2xt+t^2}$

$U_n$も同様の漸化式なので
$(1-2xt+t^2)F(t)=U_0(x)+(U_1(x)-2xU_0(x))t$である。また
$U_0(x)=1,U_1(x)=2x$より
$(1-2xt+t^2)F(t)=1$
よって、次の公式を得る。

第二種チェビシェフ多項式の通常型母関数

$\displaystyle \sum_{n=0}^{\infty} U_n(x)t^n=\frac{1}{1-2xt+t^2}$

指数型母関数

チェビシェフ多項式の閉じた形の式を用いると、
\begin{align} \displaystyle \sum_{n=0}^{\infty} T_n(x)\frac{t^n}{n!} &= \sum_{n=0}^{\infty}\frac{1}{2}\left[ \left( x+\sqrt{x^2 - 1} \right)^n +\left( x-\sqrt{x^2 -1} \right)^n \right] \frac{t^n}{n!} \\ &=\frac{1}{2}\left( \sum_{n=0}^{\infty}\frac{1}{n!}\left( t\left( x+\sqrt{x^2 -1}\right)\right)^n+\sum_{n=0}^{\infty}\frac{1}{n!}\left( t\left( x-\sqrt{x^2 -1} \right)\right)^n \right) \\ &= \frac{1}{2}\left( \mathrm{exp}\left(t\left(x+\sqrt{x^2 - 1}\right )\right)+\mathrm{exp}\left(t\left(x-\sqrt{x^2 -1}\right)\right)\right)\\ &= e ^{xt}\frac{\mathrm{exp}\left( t\sqrt{x^2 -1}\right)+\exp \left(-t\sqrt{x^2 -1}\right)}{2}\\ &= e ^{xt}\cosh \left(t\sqrt{x^2 -1}\right) \end{align}
よって、次の公式を得る。

第一種チェビシェフ多項式の指数型母関数

\begin{align} \displaystyle \sum_{n=0}^{\infty} T_n(x)\frac{t^n}{n!} &= \frac{1}{2}\left( \mathrm{exp}\left(t\left(x+\sqrt{x^2 - 1}\right )\right)+\mathrm{exp}\left(t\left(x-\sqrt{x^2 -1}\right)\right)\right) \\ &= e^{xt}\cosh \left(t\sqrt{x^2-1} \right) \end{align}

第二種チェビシェフ多項式の指数型母関数

$\displaystyle \sum_{n=0}^{\infty} U_n(x)\frac{t^n}{n!}=e^{xt}\left(\cosh\left(t\sqrt{x^2-1}\right)+\frac{x}{\sqrt{x^2-1}}\sinh\left(t\sqrt{x^2-1}\right)\right)$

証明は公式9と同様である。

チェビシェフ多項式一覧(0<=n<=10)

第一種チェビシェフ多項式

\begin{align} T_{0}(x)&=1\\ T_{1}(x)&=x\\ T_{2}(x)&=2 x^{2} - 1\\ T_{3}(x)&=4 x^{3} - 3 x\\ T_{4}(x)&=8 x^{4} - 8 x^{2} + 1\\ T_{5}(x)&=16 x^{5} - 20 x^{3} + 5 x\\ T_{6}(x)&=32 x^{6} - 48 x^{4} + 18 x^{2} - 1\\ T_{7}(x)&=64 x^{7} - 112 x^{5} + 56 x^{3} - 7 x\\ T_{8}(x)&=128 x^{8} - 256 x^{6} + 160 x^{4} - 32 x^{2} + 1\\ T_{9}(x)&=256 x^{9} - 576 x^{7} + 432 x^{5} - 120 x^{3} + 9 x\\ T_{10}(x)&=512 x^{10} - 1280 x^{8} + 1120 x^{6} - 400 x^{4} + 50 x^{2} - 1\\ \end{align}

第二種チェビシェフ多項式

\begin{align} U_{0}(x)&=1\\ U_{1}(x)&=2 x\\ U_{2}(x)&=4 x^{2} - 1\\ U_{3}(x)&=8 x^{3} - 4 x\\ U_{4}(x)&=16 x^{4} - 12 x^{2} + 1\\ U_{5}(x)&=32 x^{5} - 32 x^{3} + 6 x\\ U_{6}(x)&=64 x^{6} - 80 x^{4} + 24 x^{2} - 1\\ U_{7}(x)&=128 x^{7} - 192 x^{5} + 80 x^{3} - 8 x\\ U_{8}(x)&=256 x^{8} - 448 x^{6} + 240 x^{4} - 40 x^{2} + 1\\ U_{9}(x)&=512 x^{9} - 1024 x^{7} + 672 x^{5} - 160 x^{3} + 10 x\\ U_{10}(x)&=1024 x^{10} - 2304 x^{8} + 1792 x^{6} - 560 x^{4} + 60 x^{2} - 1\\ \end{align}

チェビシェフ多項式の直交性

$T_n$$U_n$は直交多項式である。
$T_n$の重みは
$\displaystyle \frac{1}{\sqrt{1-x^2}}$で区間$[-1,1]$つまり、
\begin{align} \int_{-1}^{1}T_n(x)T_m(x)\frac{\mathrm{d}x}{\sqrt{1-x^2}}=\left\{ \begin{alignedat}{3} &0\qquad &\mathrm{if}&\quad &n \neq m \\ &\pi \qquad &\mathrm{if}&\quad &n = m = 0\\ &\pi/2 \qquad &\mathrm{if}&\quad &n=m\neq 0 \end{alignedat} \right. \end{align}

$x=\cos \theta$とすれば、
$$\int_{-1}^{1}T_n(x)T_m(x)\frac{\mathrm{d}x}{\sqrt{1-x^2}}$$
$$=\int_{0}^{\pi}\cos(n\theta)\cos(m\theta)\mathrm{d}\theta$$
$$=\frac{1}{2}\int_{0}^{\pi}\left(\cos(n+m)\theta+\cos(n-m)\theta\right)\mathrm{d}\theta$$
$n\neq m$のとき
$$=\frac{1}{2}\left[\frac{\sin(n+m)\theta}{n+m}+\frac{\sin(n-m)\theta}{n-m}\right]_{0}^{\pi}$$
$$=0$$
$n=m=0$のとき
$$=\int_{0}^{\pi}\mathrm{d}\theta$$
$$=\pi$$
$n=m\neq0$のとき
$$=\frac{1}{2}\int_{0}^{\pi}\left(\cos(2n\theta)+1\right)\mathrm{d}\theta$$
$$=\frac{1}{2}\left[\frac{\sin(2n\theta)}{2n}\right]_{0}^{\pi}+\frac{\pi}{2}$$
$$=\frac{\pi}{2}$$

$U_n$の重みは
$\displaystyle \sqrt{1-x^2}$で区間は$[-1,1]$つまり、
\begin{align} \int_{-1}^{1}U_n(x)U_m(x)\sqrt{1-x^2}\mathrm{d}x=\left\{ \begin{alignedat}{2} &0\qquad &\mathrm{if}\quad n\neq m\\ &\pi/2\qquad &\mathrm{if}\quad n = m \end{alignedat} \right. \end{align}

上の証明と同様に$x = \cos \theta$で置換すればよい。

チェビシェフ多項式の特殊値

\begin{align} &T_n(-x)=(-1)^n T_n(x)\\ &U_n(-x) = (-1)^n U_n (x)\\ &T_n(1) = 1\\ &T_n(-1) = (-1)^n\\ &U_n(1) = n+1\\ &U_n(-1) = (-1)^n (n+1) \end{align}

公式5や公式6を用いればよい。

チェビシェフ多項式が満たす微分方程式

投稿日:2023518

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

Fefr
Fefr
22
2308
好きなこと書きます

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中