こんにちは,たい焼きです.最近は問題の解説記事ばかり書いていたので,1つのテーマに絞って記事を書くのが久々な気がします.(というかこれが2つ目)今回はディガンマ関数の積分表示の導出に焦点を当てていきたいと思います
ディガンマ関数は,ガンマ関数の対数微分で定義されます.
基本性質についてですが,神鳥奈紗さんのPDFに証明とセットで書かれているので省略します.以下を参照してください.(おい)
https://twitter.com/sounansya_29/status/1320736570027110400?s=21
........................さすがに雑すぎる気がしたのでとりあえず列挙しておきます.
上の2つはガンマ関数の性質で対数微分を用いることで容易に導出でき,$\psi(1)$はガンマ関数の無限乗積表示を用いることで以下のように導出できます.
$$
\begin{eqnarray}
\psi(z)&=&\frac{d}{dz}\ln\Gamma(z) \\&=&\frac{d}{dz}\ln\lim_{n\rightarrow\infty}\frac{n^zn!}{\prod_{k=0}^{n}(z+k)} \\
&=&\frac{d}{dz}\lim_{n\rightarrow\infty}\left(z\ln n+\ln n!-\sum_{k=0}^n\ln (z+k)\right) \\
&=&\lim_{n\rightarrow\infty}\left(\ln n-\sum_{k=0}^n\frac{1}{z+k}\right)
\end{eqnarray}
$$
$z=1$として,
$$
\psi(1)=\lim_{n\rightarrow\infty}\left(\ln n-\sum_{k=1}^n\frac{1}{k}\right)=-\gamma
$$
を得ます.$\displaystyle \psi\left(\frac{1}{2}\right)$の方は,後ほど示す積分表示から得られます.
結局導出を書きました.自分偉い!!!!!
また,ディガンマ関数の積分表示の導出にポリガンマ関数(ディガンマ関数の高階微分)も用いるので,そちらについても軽く載せておきます.
ポリガンマ関数にも漸化式や相反公式は存在しますが,簡単に導けるので省略します.
この記事の本題ですね.wikipediaのディガンマ関数の解説ページに掲載されている積分表示の上から3つを全て証明していきます.
https://ja.wikipedia.org/wiki/%E3%83%87%E3%82%A3%E3%82%AC%E3%83%B3%E3%83%9E%E9%96%A2%E6%95%B0
Wikipediaより
ただし,これらを導く上でディガンマ関数の差の積分表示とオイラーの定数$\gamma$の積分表示を用います.
$$ \begin{eqnarray} \int_0^1 \frac{t^{x-1} - t^{y-1}}{1-t} dt &=&\lim_{a\rightarrow 0}\left(\int_0^1t^{x-1}(1-t)^{a-1}dt-\int_0^1t^{y-1}(1-t)^{a-1}dt \right)\\ &=&\lim_{a\rightarrow 0}\left(B(x,a)-B(y,a) \right)\\ &=&\lim_{a\rightarrow 0}\left(\frac{\Gamma(x)\Gamma(a)}{\Gamma(x+a)}-\frac{\Gamma(y)\Gamma(a)}{\Gamma(y+a)} \right)\\ &=&\lim_{a\rightarrow 0} \frac{\frac{\Gamma(x)}{\Gamma(x+a)}-\frac{\Gamma(y)}{\Gamma(y+a)}}{\frac1{\Gamma(a)}}\\ &=&\lim_{a\rightarrow 0} \frac{\frac{\Gamma(y)\psi(y+a)}{\Gamma(y+a)}-\frac{\Gamma(x)\psi(x+a)}{\Gamma(x+a)}}{-\frac{\psi(a)}{\Gamma(a)}}\\ &=&-\lim_{a\rightarrow 0}\left(\frac{\Gamma(y)\psi(y+a)}{\Gamma(y+a)}-\frac{\Gamma(x)\psi(x+a)}{\Gamma(x+a)}\right)\frac{\Gamma(a+1)}{a\psi(a+1)-1}\\ &=&\psi(y)-\psi(x) \quad\square \end{eqnarray} $$
(神鳥奈紗さんのPDFから引用)
これを用いて,ポリガンマ関数の積分表示を導出します.
ディガンマ関数の差の積分表示において$x=1$,$y=z$とすると,
$$\begin{eqnarray}
\psi(z)&=&\psi(1)+\int_0^1 \frac{1-t^{z-1}}{1-t}dt \\
&=&\psi(1)+\int_0^\infty \frac{e^{-s}-e^{-zs}}{1-e^{-s}}ds\qquad〈t=e^{-s}〉
\end{eqnarray}$$
これの両辺を$m$回微分することで,積分表示が得られる.$\square$
導出する1つ目の積分表示は以下の式です.(都合上,wikipediaに並んでいる順番とは異なります)
$$ \begin{eqnarray} \int_1^\infty \frac{s^{z-1}-1}{s^z(s-1)}ds&=&\int_1^\infty \frac{1}{s^z(s-1)}\int_0^{z-1}\frac{d}{dt}s^tdtds \\ &=&\int_1^\infty \frac{\ln s}{s^z(s-1)}\int_0^{z-1}s^tdtds \\ &=&\int_0^{z-1}\int_1^\infty \frac{s^{t-z}\ln s}{s-1}dsdt \\ &=&\int_0^{z-1}\int_0^\infty \frac{xe^{x(t-z)}}{e^x-1}e^xdxdt\qquad 〈x=\ln s〉\\ &=&\int_0^{z-1}\int_0^\infty \frac{xe^{x(t-z)}}{1-e^{-x}}dxdt \\ &=&\int_0^{z-1}\psi^{(1)}(z-t)dt \\ &=&\left[-\psi(z-t)\right]_0^{z-1} \\ &=&-\psi(1)+\psi(z) \end{eqnarray} $$
$\psi(1)=-\gamma$より,$\displaystyle \psi(z)=-\gamma+\int_1^\infty \frac{s^{z-1}-1}{s^z(s-1)}ds$が示された.$\square$
5行目から6行目の変形で,ポリガンマ関数の積分表示を用いたことに注意してください.
2つ目の証明に移る前に,オイラーの定数の積分表示を確認しておきましょう.
これの証明は,以下のページを参照してください.
https://ja.wikipedia.org/wiki/%E3%82%AA%E3%82%A4%E3%83%A9%E3%83%BC%E3%81%AE%E5%AE%9A%E6%95%B0
では,2つ目の積分表示の導出をしていきましょう.
$$ \begin{eqnarray} \int_0^\infty \left(\frac{e^{-t}}{t}-\frac{e^{-zt}}{1-e^{-t}}\right)dt &=&\int_0^\infty \left(\frac{e^{-t}}{t}-\frac{1}{e^t-1}\right)dt+\int_0^\infty \left(\frac{1}{e^t-1}-\frac{e^{-zt}}{1-e^{-t}}\right)dt \\ &=&-\gamma+\int_0^\infty \frac{e^{-t}-e^{-zt}}{1-e^{-t}}dt \end{eqnarray} $$
これは,ポリガンマ関数の積分表示の導出で出てきた形に他ならない.したがって,$\displaystyle \psi(z)=\int_0^\infty \left(\frac{e^{-t}}{t}-\frac{e^{-zt}}{1-e^{-t}}\right)dt$が示された.$\square$
積分の形がオイラーの定数の積分表示に似ていることから変形していくと上手く示せましたね.1つ目のよりはかなり易しかったと思います.
では,3つ目の積分表示の導出をしていきましょう.
先ほど導出したディガンマ関数の積分表示➁を変形して,これを示す.
$$
\begin{eqnarray}
\psi(z)&=&\int_0^\infty\left(\frac{e^{-t}}{t}-\frac{e^{-zt}}{1-e^{-t}}\right)dt \\
&=&\lim_{\delta\rightarrow+0}\left(\int_\delta^\infty\frac{e^{-t}}{t}dt-\int_\delta^\infty \frac{e^{-zt}}{1-e^{-t}}dt\right) \\
&=&\lim_{\delta\rightarrow+0}\left(\int_\delta^\infty\frac{e^{-t}}{t}dt-\int_{e^\delta-1}^\infty\frac{(1+s)^{-z}}{1-(1+s)^{-1}}\frac{ds}{1+s}\right) \qquad〈s=e^t-1〉 \\
&=&\lim_{\delta\rightarrow+0}\left(\int_\delta^\infty\frac{e^{-t}}{t}dt-\int_{e^\delta-1}^\infty\frac{1}{t(1+t)^z}dt\right)
\end{eqnarray}
$$
区間$[\delta,e^\delta-1]$で$\displaystyle\frac{1}{t(1+t)^z}$を評価すると,${\rm Re}(z)\gt0$より
$$
0\leq\frac{1}{t(1+t)^z} \leq \frac{1}{t} \leq \frac{1}{\delta}
$$
したがって,定積分を以下のように評価できる.
$$
\left|\int_\delta^{e^\delta-1}\frac{1}{t(1+t)^z}dt\right|\leq\left|\int_\delta^{e^\delta-1}\frac{1}{\delta}dt\right|\leq O(\delta)
$$
これより,
$$
\begin{eqnarray}
\psi(z)&=&\lim_{\delta\rightarrow+0}\left(\int_\delta^\infty \frac{e^{-t}}{t}dt-\int_\delta^\infty \frac{1}{t(1+t)^z}dt\right) \\
&=&\lim_{\delta\rightarrow+0}\int_\delta^\infty\left(e^{-t}-\frac{1}{(1+t)^z}\right)\frac{dt}{t} \\
&=&\int_0^\infty\left(e^{-t}-\frac{1}{(1+t)^z}\right)\frac{dt}{t}
\end{eqnarray}
$$
となるので,示された.$\square$
導出する式の形が2つ目の積分表示に似ていたので方針は立てやすかったと思いますが,発散する積分を項別に分けて計算することはできない(修正前は分けた状態で投稿してました...)ので,極限を考える必要があります.変数変換によって生まれてしまう積分区間の誤差は$\delta$の1次と同位の微少量なので,極限$\delta\rightarrow+0$を取ると0に収束しますね.
最後まで目を通していただきありがとうございます.数ヶ月前にこれらの導出にチャレンジしたときはどれも導くことが出来なかったので,自分の成長を感じられました.