21

ディガンマ関数の基本性質と積分表示

3034
2
$$$$

こんにちは,たい焼きです.最近は問題の解説記事ばかり書いていたので,1つのテーマに絞って記事を書くのが久々な気がします.(というかこれが2つ目)今回はディガンマ関数の積分表示の導出に焦点を当てていきたいと思います

ディガンマ関数の定義・基本性質

ディガンマ関数は,ガンマ関数の対数微分で定義されます.

・ディガンマ関数の定義
ガンマ関数$\Gamma(z)$対し,
$$ \psi(z)=\frac{d}{dz}\ln\Gamma(z)=\frac{\Gamma'(z)}{\Gamma(z)} $$

基本性質についてですが,神鳥奈紗さんのPDFに証明とセットで書かれているので省略します.以下を参照してください.(おい)
https://twitter.com/sounansya_29/status/1320736570027110400?s=21
........................さすがに雑すぎる気がしたのでとりあえず列挙しておきます.

・漸化式
$$\psi(z+1)=\psi(z)+\frac{1}{z}$$
・相反公式
$$\psi(1-z)-\psi(z)=\pi\cot \pi z$$
・特殊値
$\psi(1)=-\gamma$$\displaystyle \psi\left(\frac{1}{2}\right)=-\gamma-2\ln2$

上の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)$の方は,後ほど示す積分表示から得られます.
 結局導出を書きました.自分偉い!!!!!
 また,ディガンマ関数の積分表示の導出にポリガンマ関数(ディガンマ関数の高階微分)も用いるので,そちらについても軽く載せておきます.

・ポリガンマ関数の定義
ディガンマ関数$\psi(z)$対し,
$$ \psi^{(m)}(z)=\frac{d^m}{{dz}^m}\psi(z) $$

ポリガンマ関数にも漸化式や相反公式は存在しますが,簡単に導けるので省略します.

ディガンマ関数の積分表示

この記事の本題ですね.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より Wikipediaより  
 ただし,これらを導く上でディガンマ関数の差の積分表示とオイラーの定数$\gamma$の積分表示を用います.

・ディガンマ関数の差の積分表示
$$ \psi(y)-\psi(x)=\int_0^1 \frac{t^{x-1}-t^{y-1}}{1-t}dt$$

$$ \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から引用)
 これを用いて,ポリガンマ関数の積分表示を導出します.

・ポリガンマ関数の積分表示
$$\psi^{(m)}(z)=(-1)^{m+1}\int_0^\infty \frac{s^me^{-zs}}{1-e^{-s}}ds\qquad(m>0)$$

ディガンマ関数の差の積分表示において$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に並んでいる順番とは異なります)

・ディガンマ関数の積分表示➀
$$\psi(z)=-\gamma+\int_1^\infty \frac{s^{z-1}-1}{s^z(s-1)}ds$$

$$ \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つ目の証明に移る前に,オイラーの定数の積分表示を確認しておきましょう.

・オイラーの定数$\gamma$の積分表示
$$ \gamma=\int_0^\infty \left(\frac{1}{e^t-1}-\frac{e^{-t}}{t}\right)dt$$

これの証明は,以下のページを参照してください.
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つ目の積分表示の導出をしていきましょう.

・ディガンマ関数の積分表示➁
$$\psi(z)=\int_0^\infty \left(\frac{e^{-t}}{t}-\frac{e^{-zt}}{1-e^{-t}}\right)dt$$

$$ \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つ目の積分表示の導出をしていきましょう.

・ディガンマ関数の積分表示➂
$$\psi(z)=\int_0^\infty \left(e^{-t}-\frac{1}{(1+t)^z}\right)\frac{dt}{t}$$

先ほど導出したディガンマ関数の積分表示➁を変形して,これを示す.
$$ \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に収束しますね.

最後に

最後まで目を通していただきありがとうございます.数ヶ月前にこれらの導出にチャレンジしたときはどれも導くことが出来なかったので,自分の成長を感じられました.

投稿日:20201115

この記事を高評価した人

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

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

バッジはありません。

投稿者

物理学と数学を嗜んでいる高2です 食わず嫌いせずにいろんな分野に触れたいと思ってます

コメント

他の人のコメント

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