0

個人的に大変だった積分2

29
0
$$$$

こんにちは。目薬の回数を勘違いしてて、より多くしていいことに気づいて革命を感じてる今日この頃。
気分がとてもいいので前回の 記事 の最後に言った積分をやります。

本日の主役

$\displaystyle{I=\int_{0}^{\frac{\pi}{4}}x^2\tan{x}dx=\frac{\pi}{4}\beta(2)-\frac{\pi^2}{32}\ln{2}-\frac{21}{64}\zeta(3)}$

ただし、$\beta(s)$はディリクレのベータ関数、$\zeta(s)$はリーマンゼータ関数である。

一応それらの定義を載せます。

ディリクレのベータ関数

実部が$0$より大きい複素数$s$について、
$\displaystyle{\beta(s)=\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{(2k-1)^s}}$

リーマンゼータ関数

実部が実部が$1$より大きい複素数$s$について、
$\displaystyle{\zeta(s)=\sum_{k=1}^{\infty}\frac{1}{k^s}}$

準備

前回同様、補題をいくつか証明しておきましょう。

$\ln(1+x)$のマクローリン展開

$\displaystyle{\ln(1+x)=\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}x^k}$

\begin{align*} f(t)=t-\frac{1}{2}t^2+\frac{1}{3}t^3-\frac{1}{4}t^4+\cdot\:\cdot\:\cdot \end{align*}
とする。ここで、両辺を$t$で微分すると
\begin{align*} f'(t)=1-t+t^2-t^3+\cdot\:\cdot\:\cdot=\frac{1}{1+t} \end{align*}
両辺を$t $で不定積分すると、
\begin{align*} \int_{0}^{x}f'(t)dt&=\int_{0}^{x}\frac{1}{1+t}dt\\ &=\ln(1+x)\\ \end{align*}
また、
\begin{align*} \int_{0}^{x}f'(t)dt&=\left[f(t)\right]_{0}^{x}\\ &=f(x) \end{align*}
よって

$\displaystyle{f(x)=\ln(1+x)}$

一方、
\begin{align*} f(t)=\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}x^k \end{align*}
と表せるため、
\begin{align*} \ln(1+x)=\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}x^k \end{align*}
が言える。($|x|\le1$)

一応補足すると、この記事ではオイラーの公式も利用する。が、複素関数において三角関数や指数関数はマクローリン展開で定義されるものと思っているため、それの証明はしない。

本日の主役

では、本日の主役をいじり倒していきましょう!

本日の主役

\begin{align*} I&=\int_{0}^{\frac{\pi}{4}}x^2\tan{x}dx\\ &=\left[-x^2\ln(\cos{x})\right]_{0}^{\frac{\pi}{4}}+2\int_{0}^{\frac{\pi}{4}}x\ln(\cos{x})dx\\ &=\frac{\pi^2}{32}\ln{2}+2I' \end{align*}
さて、前回同様部分積分を適用しました。今回の問題の大戦犯は$I'$です。気合を入れてしばきましょう。
\begin{align*} I'&=\int_{0}^{\frac{\pi}{4}}x\ln(\cos{x})dx\\ &=\Re\int_{0}^{\frac{\pi}{4}} x\ln\left(\frac{e^{ix}+e^{-ix}}{2}\right)dx \end{align*}
この変形がこの積分のカギです。が、実はフーリエ級数展開をすればもっと楽に解くことができます。でもここでは自分の使った解法で解こうと思います^^
\begin{align*} I'&=\Re\int_{0}^{\frac{\pi}{4}}x\ln\left(\frac{e^{ix}+e^{-ix}}{2}\right)dx\\ &=\Re\int_{0}^{\frac{\pi}{4}}x\ln\left(e^{ix}+e^{-ix}\right)dx-\ln{2}\int_{0}^{\frac{\pi}{4}}xdx\\ &=-\frac{\pi^2}{32}\ln{2}+\Re\int_{0}^{\frac{\pi}{4}}x\ln\left\{e^{ix}\left(1+e^{-2ix}\right)\right\}dx\\ &=-\frac{\pi^2}{32}\ln{2}+\Re\int_{0}^{\frac{\pi}{4}}x\ln\left(1+e^{-2ix}\right)dx\\ &=-\frac{\pi^2}{32}\ln{2}+\Re\int_{0}^{\frac{\pi}{4}}x\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}\left(e^{-2ix}\right)^kdx\\ &\overset{!}{=}-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}\int_{0}^{\frac{\pi}{4}}x\cos(2kx)dx \end{align*}
最後の変形について少し補足。$\left(e^{-2ix}\right)^k=\cos(2kx)-i\sin(2kx)$と変形できる。またここでは実部を取るため$\sin$の項は捨てて、$\cos$だけとっている。
そしてもう一つ注意。この記事を書いている人は一様収束などを理解していないため、シグマとインテグラルを無言で入れ替えてます。本当はそれについても説明するべきですが、ごめんなさい。(一応、$\overset{!}{=}$で知らせています。)
\begin{align*} I'&=-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}\int_{0}^{\frac{\pi}{4}}x\cos(2kx)dx\\ &=-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}I_1\\ I_1&=\int_{0}^{\frac{\pi}{4}}x\cos(2kx)dx\\ &=\left[\frac{1}{2x}x\sin(2kx)\right]_{0}^{\frac{\pi}{4}}-\frac{1}{2k}\int_{0}^{\frac{\pi}{4}}\sin(2kx)dx\\ &=\frac{\pi}{8}\frac{1}{k}\sin\frac{k\pi}{2}+\frac{1}{4k^2}\left[\cos(2kx)\right]_{0}^{\frac{\pi}{4}}\\ &=\frac{\pi}{8}\frac{1}{k}\sin\frac{k\pi}{2}+\frac{1}{4k^2}\cos\frac{k\pi}{2}-\frac{1}{4k^2} \end{align*}
よって、$I'$
\begin{align*} I'&=-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}I_1\\ &=-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}\left(\frac{\pi}{8}\frac{1}{k}\sin\frac{k\pi}{2}+\frac{1}{4k^2}\cos\frac{k\pi}{2}-\frac{1}{4k^2}\right) \end{align*}
級数を一つ一つ簡単にしていきましょう。
\begin{align*} \frac{\pi}{8}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^2}\sin\frac{k\pi}{2} \end{align*}
から。$\displaystyle{\sin\frac{k\pi}{2}}$$k$が偶数で$0$になるため、$k$が奇数のときのみを考える。よって、
\begin{align*} \frac{\pi}{8}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^2}\sin\frac{k\pi}{2}&=\frac{\pi}{8}\sum_{k=1}^{\infty}\frac{(-1)^{(2k-1)-1}}{(2k-1)^2}\sin\frac{(2k-1)\pi}{2}\\ &=\frac{\pi}{8}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{(2k-1)^2}\\ &=\frac{\pi}{8}\beta(2) \end{align*}
最後にディリクレのベータ関数の定義がもろに出てくれましたね。これで一つ目の級数は終わり。

\begin{align*} \frac{1}{4}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^3}\cos\frac{k\pi}{2} \end{align*}
これも先ほどと同じように考える。$\displaystyle{\cos\frac{k\pi}{2}}$$k$が奇数のときに$0$となるため、偶数のときのみを考えればよい。よって、
\begin{align*} \frac{1}{4}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^3}\cos\frac{k\pi}{2}&=\frac{1}{4}\sum_{k=1}^{\infty}\frac{(-1)^{2k-1}}{(2k)^3}\cos(k\pi)\\ &=\frac{1}{32}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^3} \end{align*}
一応補足。二つ目のイコールでは外にあったマイナスをわざと中に入れて、$(-1)$の指数を一つずらした。(これのほうが個人的に楽だから)
\begin{align*} \frac{1}{32}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^3}&=\frac{1}{32}\left(\sum_{k=1}^{\infty}\frac{1}{k^3}-2\sum_{k=1}^{\infty}\frac{1}{(2k)^3}\right)\\ &=\frac{1}{32}\frac{3}{4}\sum_{k=1}^{\infty}\frac{1}{k^3}\\ &=\frac{3}{128}\zeta(3) \end{align*}
これで二つ目の級数は終わり。三つ目の級数は二つ目のやつの途中とほとんど同じため、結果のみ載せる。
\begin{align*} -\frac{1}{4}\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k^3}&=-\frac{3}{16}\zeta(3) \end{align*}
以上より、
\begin{align*} I'&=-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}I_1\\ &=-\frac{\pi^2}{32}\ln{2}+\sum_{k=1}^{\infty}\frac{(-1)^{k-1}}{k}\left(\frac{\pi}{8}\frac{1}{k}\sin\frac{k\pi}{2}+\frac{1}{4k^2}\cos\frac{k\pi}{2}-\frac{1}{4k^2}\right)\\ &=-\frac{\pi^2}{32}\ln{2}+\frac{\pi}{8}\beta(2)+\frac{3}{128}\zeta(3)-\frac{3}{16}\zeta(3)\\ &=-\frac{\pi^2}{32}\ln{2}+\frac{\pi}{8}\beta(2)-\frac{21}{128}\zeta(3) \end{align*}
そして、
\begin{align*} I&=\frac{\pi^2}{32}\ln{2}+2I'\\ &=\frac{\pi^2}{32}\ln{2}-\frac{\pi^2}{16}\ln{2}+\frac{\pi}{4}\beta(2)-\frac{21}{64}\zeta(3)\\ &=\frac{\pi}{4}\beta(2)-\frac{\pi^2}{32}\ln{2}-\frac{21}{64}\zeta(3) \end{align*}
よって、命題は示された。

やったね!

お疲れさまでした。自分が解いているときはディリクレのベータ関数が積分表示で現れるのかな?とか予想しながら解いていましたが、級数の形で出てきましたね。個人的には$\zeta(3)$がいてなんか好きな積分でした。では次回、$\displaystyle{\int_{0}^{\frac{\pi}{4}}x^3\tan{x}dx}$を解きます!
とはなりません。$\text{WolframAlpha}$に入れましたが、まだ多重対数関数に疎いのでできませんでした。代わりに自作の積分を解説しようかなと思ってます。ほな、さいなら!

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

この記事を高評価した人

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

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

バッジはありません。

投稿者

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

コメント

他の人のコメント

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