0

ゼータ関数の特殊値を計算するお話

31
0
$$$$

こんにちは。前回の記事ではリーマンゼータ関数の特殊値を使っていたので、ここで証明しようと思います。また、ついでにζ(4)~ζ(8)までも求めてしまいましょう!

本題

フーリエ級数展開を用いて求めたいと思います。それが成り立つ証明は力不足でできません。ごめん。

フーリエ級数展開

$f$が周期$T$のまともな関数なら、
$\displaystyle{f(x)=\frac{a_0}{2}+\sum_{k=1}^{\infty}a_k\cos\frac{2\pi kx}{T}+b_k\sin\frac{2\pi kx}{T}}$
ただし、
$\displaystyle{a_k=\frac{2}{T}\int_{0}^{T}f(x)\cos\frac{2\pi kx}{T}dx} $
$\displaystyle{b_k=\frac{2}{T}\int_{0}^{T}f(x)\sin\frac{2\pi kx}{T}dx}$

今回は$f(x)=x^{2m}$$[0,2\pi]$でフーリエ級数展開するのですが、周期性を加味すると、$[-\pi,\pi]$でもよいことがわかります。(図1参照)

図1において、関数の区切れが$π$の倍数で起きています。この関数を$[0,2\pi]$で積分しても$[-\pi,\pi]$で積分しても値が変わらないから、今回は$[-\pi,\pi]$で計算します。

$\zeta(2)$

\begin{align*} f(x)=x^2\quad (-\pi\le x\le\pi) \end{align*}
をフーリエ級数展開する。このとき、
\begin{align*} a_k&=\frac{2}{2\pi}\int_{-\pi}^{\pi}x^2\cos\frac{2\pi kx}{2\pi}dx\\ &=\frac{1}{\pi}\int_{-\pi}^{\pi}x^2\cos(kx)dx\\ b_k&=\frac{2}{2\pi}\int_{-\pi}^{\pi}x^2\sin\frac{2\pi kx}{2\pi}dx\\ &=\frac{1}{\pi}\int_{-\pi}^{\pi}x^2\sin(kx)dx \end{align*}
ここで、$b_k$は奇関数のため、積分値は$0$とわかる。なので$a_k$のみ頑張ればよい。$a_k$は部分積分を繰り返すことで計算され、実際に計算すると
\begin{align*} a_k&=\frac{1}{\pi}\frac{4}{k^2}\pi\cos(k\pi)\\ &=\frac{4}{k^2}(-1)^{k} \end{align*}
また、$a_0$
\begin{align*} a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^2dx\\ &=\frac{2}{3}\pi^2 \end{align*}
以上より、$f(x)$は次の様にフーリエ級数展開できることがわかる。
\begin{align*} f(x)=\frac{1}{2}\frac{2}{3}\pi^2+\sum_{k=1}^{\infty}\cos(kx)\cdot\frac{4}{k^2}(-1)^{k} \end{align*}
ここで$x=\pi$とすると、$f(x)=x^2$であったため、
\begin{align*} \pi^2&=\frac{\pi^2}{3}+\sum_{k=1}^{\infty}\frac{4}{k^2}\\ \sum_{k=1}^{\infty}\frac{1}{k^2}&=\frac{1}{6}\pi^2\\ \zeta(2)&=\frac{\pi^2}{6} \end{align*}
なんと、きれいにゼータ関数の形が表れ、求まりました。この方法で$\zeta(4),\zeta(6),\zeta(8)$を計算したいと思います。また、それらの過程でも$b_k$$0$であり、残る積分も部分積分を繰り返すことで計算できるので、細かな計算は省きたいと思います。


$\zeta(4)$

\begin{align*} f(x)=x^4\quad(-\pi\le x\le\pi) \end{align*}
これをフーリエ級数展開する。このとき、$\zeta(2)$と同様に$b_k$$0$である。よって
\begin{align*} f(x)=\frac{a_0}{2}+\sum_{k=1}^{\infty}a_k\cos(kx) \end{align*}
$a_0$を計算すると、
\begin{align*} a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^4dx\\ &=\frac{1}{\pi}\frac{2}{5}\pi^5\\ &=\frac{2}{5}\pi^4 \end{align*}
また、$a_k$は部分積分を繰り返すことで次のように計算できる。
\begin{align*} a_k&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^4\cos(kx)dx\\ &=\frac{1}{\pi}\left(\frac{8}{k^2}\pi^3\cos(k\pi)-\frac{48}{k^4}\pi\cos(k\pi)\right)\\ &=\frac{8}{k^2}\pi^2(-1)^{k}-\frac{48}{k^4}(-1)^{k} \end{align*}
以上より、
\begin{align*} x^4&=\frac{1}{2}\frac{2}{5}\pi^4+\sum_{k=1}^{\infty}\left(\frac{8}{k^2}\pi^2(-1)^{k}-\frac{48}{k^4}(-1)^{k}\right)\cos(kx) \end{align*}
が言え、このとき$x=\pi$とすると、
\begin{align*} \pi^4&=\frac{\pi^4}{5}+\sum_{k=1}^{\infty}\frac{8}{k^2}\pi^2-\frac{48}{k^4}\\ \end{align*}
また、$\displaystyle\zeta(2)=\frac{\pi^2}{6}$を用いてこれを整理すると、
\begin{align*} \frac{4}{5}\pi^4&=8\pi^2\zeta(2)-48\zeta(4)\\ 48\zeta(4)&=\frac{8}{6}\pi^4-\frac{4}{5}\pi^4\\ &=\frac{8}{15}\pi^4\\ \zeta(4)&=\frac{\pi^4}{90} \end{align*}
よって値を導出することができた。


$\zeta(6)$

\begin{align*} f(x)=x^6\quad(-\pi\le x\le\pi) \end{align*}
をフーリエ級数展開する。このとき、$\zeta(2)$の場合と同様に$b_k$$0$であり、$a_k$のみ計算すればよい。よって、
\begin{align*} f(x)=\frac{a_0}{2}+\sum_{k=1}^{\infty}a_k\cos(kx) \end{align*}
$a_0$を計算すると、
\begin{align*} a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^6dx\\ &=\frac{1}{\pi}\frac{2}{7}\pi^7\\ &=\frac{2}{7}\pi^6 \end{align*}
また、$a_k$は部分積分を繰り返すことで次の様に計算できる。
\begin{align*} a_k&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^6\cos(kx)dx\\ &=\frac{1}{\pi}\left(\frac{12}{k^2}\pi^5\cos(k\pi)-\frac{240}{k^4}\pi^3\cos(k\pi)+\frac{1440}{k^6}\pi\cos(k\pi)\right)\\ &=\frac{12}{k^2}\pi^4(-1)^{k}-\frac{240}{k^4}\pi^4(-1)^{k}+\frac{1440}{k^6}(-1)^{k} \end{align*}
以上より、
\begin{align*} x^6&=\frac{1}{2}\frac{2}{7}\pi^6+\sum_{k=1}^{\infty}\left(\frac{12}{k^2}\pi^4(-1)^{k}-\frac{240}{k^4}\pi^4(-1)^{k}+\frac{1440}{k^6}(-1)^{k}\right)\cos(kx) \end{align*}
がわかる。このとき、$x=\pi$とすると
\begin{align*} \pi^6&=\frac{\pi^6}{7}+\sum_{k=1}^{\infty}\frac{12}{k^2}\pi^4-\frac{240}{k^4}\pi^4+\frac{1440}{k^6}\\ \end{align*}
$\displaystyle \zeta(2)=\frac{\pi^2}{6}\:,\:\zeta(4)=\frac{\pi^4}{90}$を踏まえてこれを整理すると、
\begin{align*} \frac{6}{7}\pi^6&=12\pi^4\zeta(2)-240\pi^2\zeta(4)+1440\zeta(6)\\ 1440\zeta(6)&=\frac{6}{7}\pi^6-\frac{12}{6}\pi^6+\frac{240}{90}\pi^6\\ &=\frac{32}{21}\pi^6\\ \zeta(6)&=\frac{\pi^6}{945} \end{align*}


$\zeta(8)$

\begin{align*} f(x)=x^8\quad(-\pi\le x\le\pi) \end{align*}
をフーリエ級数展開する。このとき、$\zeta(2)$の場合と同様に$b_k$$0$であるため$a_k$に注目すればよい。よって、
\begin{align*} f(x)&=\frac{a_0}{2}+\sum_{k=1}^{\infty}a_k\cos(kx) \end{align*}
と展開できる。ここで、$a_0$
\begin{align*} a_0&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^8dx\\ &=\frac{1}{\pi}\frac{2}{9}\pi^9\\ &=\frac{2}{9}\pi^8 \end{align*}
また、$a_k$は部分積分を繰り返すことで計算でき、実際に計算すると次の様になる。
\begin{align*} a_k&=\frac{1}{\pi}\int_{-\pi}^{\pi}x^8\cos(kx)dx\\ &=\frac{1}{\pi}\left(\frac{16}{k^2}\pi^7\cos(k\pi)-\frac{672}{k^4}\pi^5\cos(k\pi)+\frac{13440}{k^6}\pi^3\cos(k\pi)-\frac{80640}{k^8}\pi\cos(k\pi)\right)\\ &=\frac{16}{k^2}\pi^6\cos(k\pi)-\frac{672}{k^4}\pi^4\cos(k\pi)+\frac{13440}{k^6}\pi^2\cos(k\pi)-\frac{80640}{k^8}\cos(k\pi) \end{align*}
以上より、
\begin{align*} x^8&=\frac{1}{2}\frac{2}{9}\pi^8+\sum_{k=1}^{\infty}\left(\frac{16}{k^2}\pi^6\cos(k\pi)-\frac{672}{k^4}\pi^4\cos(k\pi)+\frac{13440}{k^6}\pi^2\cos(k\pi)-\frac{80640}{k^8}\cos(k\pi)\right)\cos(kx) \end{align*}
とフーリエ級数展開できる。ここで、$x=\pi$とすると、
\begin{align*} \pi^8&=\frac{\pi^8}{9}+\sum_{k=1}^{\infty}\frac{16}{k^2}\pi^6-\frac{672}{k^4}\pi^4+\frac{13440}{k^6}\pi^2-\frac{80640}{k^8} \end{align*}
この式について、$\displaystyle\zeta(2)=\frac{\pi^2}{6}\:,\:\zeta(4)=\frac{\pi^4}{90}\:,\:\zeta(6)=\frac{\pi^6}{945}$を踏まえて計算をすると、
\begin{align*} \frac{8}{9}\pi^8&=16\pi^6\zeta(2)-672\pi^4\zeta(4)+13440\pi^2\zeta(6)-80640\zeta(8)\\ 80640\zeta(8)&=-\frac{8}{9}\pi^8+\frac{16}{6}\pi^8-\frac{672}{90}\pi^8+\frac{13440}{945}\pi^8\\ &=\frac{128}{15}\pi^8\\ \zeta(8)&=\frac{\pi^8}{9450} \end{align*}


記事の趣旨とは少しずれますが、$\zeta(2)$に関して、オイラーの面白い話を紹介したいと思います。
「平方数の逆数和の収束値は何だろうか。」という問題に対し、オイラー自身も解決しようと取り組んでいました。オイラーは大まかな値をつかもうとし、級数の収束スピードを上げた級数を導き出し、正確な値を目の当たりにしました。その瞬間、オイラーは何と、「これは$\pi^2$$6$で割った値に違いない。」と確信し、このことを証明しにかかりました。オイラーはそれを$\sin{x}$のマクローリン展開と無限乗積展開の係数を比較することにより証明しています。個人的にはラマヌジャンと同等レベルに気持ち悪いことしてて好きなエピソードです^^
次回は、この方法で$\zeta(3)$をつついてみようと思います。ほな、さいなら!

投稿日:8日前
更新日:7日前
数学の力で現場を変える アルゴリズムエンジニア募集 - Mathlog served by OptHub

この記事を高評価した人

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

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

バッジはありません。

投稿者

いつの間にか高校生になった翁です。 書きたくなったことを適当に書いていきます。 注意:ミス多いです。見つけたら指摘のコメントをしていただけると助かります。自分でも努力してます_(_×-×)_

コメント

他の人のコメント

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