ガンマ関数

提供: miniwiki
2018/8/19/ (日) 17:38時点におけるAdmin (トーク | 投稿記録)による版 (1版 をインポートしました)
移動先:案内検索
ファイル:Gamma.png
y = Γ(x) のグラフ
ファイル:Gamma abs.png
Γ(x + iy) の絶対値
(グラフ中「Re」は x に相当、「Im」は y に相当)

数学においてガンマ関数(ガンマかんすう、: Gamma function)とは、階乗の概念を複素数全体に拡張した特殊関数である。互いに同値となるいくつかの定義が存在するが、1729年、数学者レオンハルト・オイラーが階乗の一般化として、最初に導入した。

定義

実部が正となる複素数 z について、次の積分で定義される関数

[math]\Gamma(z)=\int^{\infty}_{0}t^{z-1}e^{-t}\,dt\qquad(\real{z}\gt 0)[/math]

をガンマ関数と呼ぶ[1]。この積分は、ルジャンドルの定義にしたがって、第二種オイラー積分とも呼ばれる。元は階乗の一般化としてオイラーが得たもので、Γ という記号は、ルジャンドルが用いたものである。それ以前は Π(x) などと表記していた(ただし Π(x) = Γ(x + 1))。

一般の複素数 z については、解析接続もしくは次の無限乗積で定義される。

[math] \Gamma(z)=\lim_{n\to\infty}\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}. [/math]

基本的性質

0 と負の整数を除く任意の複素数 z に対して

[math]\Gamma(z+1)=z\Gamma(z) \,[/math]

が成り立つ。実際、Re(z) > 0 に対してはオイラー積分による定義から

[math]\begin{align} \Gamma(z+1)&=\int_{t=0}^{\infty}{t^z(-e^{-t})'}\,dt \\ &=\left[-t^z e^{-t}\right]_{t=0}^{\infty} +z\int_{t=0}^{\infty}{t^{z-1}e^{-t}}\,dt \\ &=z\Gamma(z) \end{align}[/math]

となる。また、

[math]\begin{align} \Gamma(1)&=\int_{t=0}^{\infty}{e^{-t}}\,dt=\left[-e^{-t}\right]_{t=0}^{\infty} \\ &=1 \end{align}[/math]

である。従って、自然数 n について

[math]\begin{align} \Gamma (n+1)=n! \end{align}[/math]

が成り立ち、その意味でガンマ関数は階乗の定義域を複素平面に拡張したものとなっている。そもそもガンマ関数は「階乗の複素数への拡張となるもの」の実例としてオイラーが考案したのである。実際には、そのような関数は無数に存在するが、正の実軸上で対数凸である解析関数という条件を付ければ、それは一意に定まりガンマ関数に他ならない(→ボーア・モレルップの定理)。右半平面においてオイラー積分で定義されたガンマ関数は全平面に有理型解析接続する。ガンマ関数は零点を持たず、原点と負の整数に一位の極を持つ。その留数は、

[math]\operatorname{Res}(\Gamma , -n) = \frac{(-1)^n}{n!} [/math]

である。また、非整数でのガンマ関数の値のうちでおそらく最も有名なのは、ガウス積分になる以下の場合であろう。

[math]\Gamma\left(\frac{1}{2}\right) = \sqrt{\pi}[/math]

これより、自然数 n について

[math] \Gamma\left(\frac{1}{2}+n\right) = \frac{(2n-1)!!}{2^{n}}\sqrt{\pi} [/math]

が成立することがわかる。ここで !!二重階乗を表す。この性質を利用して高次元の体積と表面積を求めることができる。また、

[math] \Gamma\left(\frac{1}{2}-n\right) = \frac{ (-2)^n}{(2n-1)!!} \sqrt{\pi} [/math]

定義の整合性

定義の積分表示と乗積表示が一致することを示す。

[math]G_n(z)=\int_{0}^{n}{t^{z-1}\left(1-\frac{t}{n}\right)^{n}}dt[/math]

とすれば

[math]\lim_{n\to\infty}{(1-t/n)^n}=e^{-t}[/math]

であるから

[math]\lim_{n\to\infty}{G_n(z)}=\int^{\infin}_{0}t^{z-1}e^{-t}dt[/math]

である。t = nu の置換により

[math]G_n(z)=n^{z}\int_{0}^{1}{u^{z-1}(1-u)^{n}}du[/math]

となる.nz を除く部分を gn(z) として

[math]g_0(z)=\int_{0}^{1}{u^{z-1}}du=\left[\frac{u^z}{z}\right]_{u=0}^{1}=\frac{1}{z}[/math]
[math]g_n(z)=\int_{0}^{1}{\left(\frac{u^{z}}{z}\right)'(1-u)^{n}}du=\frac{n}{z}\int_{u=0}^{1}{u^{z}(1-u)^{n-1}}du=\frac{n}{z}g_{n-1}(z+1)[/math]

これにより

[math]G_n(z)=\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}[/math]

を得る。故に

[math]\int^{\infin}_{0}t^{z-1}e^{-t}dt=\lim_{n\to\infty}G_n(z)=\lim_{n\to\infty}\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}[/math]

である。

ワイエルシュトラスの乗積表示

オイラーの乗積表示からオイラーの定数

[math]\gamma=\lim_{n\to\infty}\left(\sum_{k=1}^{n}\frac{1}{k}-\log{n}\right)[/math]

を括り出すとワイエルシュトラスの乗積表示が得られる。ワイエルシュトラスはガンマ関数が負の整数にを持つことを嫌って逆数を用いた。ガンマ関数の逆数は複素平面全体で正則である。

[math]\frac{1}{\Gamma(z)}=\lim_{n\to\infty}\frac{\prod_{k=0}^{n}{(z+k)}}{n^zn!}=\lim_{n\to\infty}zn^{-z}\left(\prod_{k=1}^{n}{e^{z/k}}\right)\left(\prod_{m=1}^{n}{\frac{z+m}{m}}e^{-z/m}\right)=ze^{{\gamma}z}\prod_{m=1}^{\infty}\left(1+\frac{z}{m}\right)e^{-z/m}[/math]

ハンケルの積分表示

ガンマ関数は次の周回積分で表される[2]。積分経路は正の無限大から実軸の上側に沿って原点に至り、原点を正の向きに回り、実軸の下側に沿って無限大に戻るものとする。但し、その偏角は[math]-\pi\le\arg(-t)\le\pi,0\le\arg(s)\le2\pi[/math]とする。

[math]\begin{align} &\Gamma(z)=\frac{i}{2\sin{\pi}z}\int_C(-t)^{z-1}e^{-t}dt\qquad(z\in\mathbb{C}\setminus\mathbb{Z})\\ &\Gamma(z)=\frac{1}{e^{2{\pi}iz}-1}\int_Cs^{z-1}e^{-s}ds\qquad(z\in\mathbb{C}\setminus\mathbb{Z})\\ &\frac{1}{\Gamma(z)}=\frac{i}{2\pi}\int_C(-t)^{-z}e^{-t}dt\qquad(z\in\mathbb{C})\\ \end{align}[/math]

これをハンケルの積分表示と呼ぶ。このハンケルの積分表示は、積分経路を適当に変形し、数値積分でガンマ関数の値を求めるために使われることがある[3]

ハンケルの積分表示の導出

極座標表示[math](-t)=re^{i\theta}[/math]を用いると、実軸の上側に沿う部分は[math]\theta=-\pi[/math][math]r=\infty[/math]から[math]r=\delta[/math]まで、原点を回る部分は[math]r=\delta[/math][math]\theta=-\pi[/math]から[math]\theta=\pi[/math]まで、実軸の下側に沿う部分は[math]\theta=\pi[/math][math]r=\delta[/math]から[math]r=\infty[/math]までとなる。

[math]\begin{align}\int_C(-t)^{z-1}e^{-t}dt &=\int_{\infty}^{\delta}(re^{-{\pi}i})^{z-1}e^{-r}dr+\int_{-\pi}^{\pi}({\delta}e^{i\theta})^{z-1}e^{{\delta}e^{i\theta}}(-i{\delta}e^{i\theta})d\theta+\int_{\delta}^{\infty}(re^{{\pi}i})^{z-1}e^{-r}dr\\ &=\int_{\infty}^{\delta}r^{z-1}e^{-{\pi}i(z-1)}e^{-r}dr-\int_{-\pi}^{\pi}i\delta^ze^{i{\theta}z}e^{{\delta}e^{i\theta}}d\theta+\int_{\delta}^{\infty}r^{z-1}e^{{\pi}i(z-1)}e^{-r}dr\\ &=\left(-e^{-{\pi}i(z-1)}+e^{{\pi}i(z-1)}\right)\int_{\delta}^{\infty}r^{z-1}e^{-r}dr-\int_{-\pi}^{\pi}i\delta^ze^{i{\theta}z}e^{{\delta}e^{i\theta}}d\theta\\ &=-2i\sin{\pi}z\int_{\delta}^{\infty}r^{z-1}e^{-r}dr-\int_{-\pi}^{\pi}i\delta^ze^{i{\theta}z}e^{{\delta}e^{i\theta}}d\theta\\ \end{align}[/math]

[math]\real{z}\gt 0[/math]とすると[math]\delta\to0[/math][math]\delta^z\to0[/math]であるから

[math]\begin{align}\int_C(-t)^{z-1}e^{-t}dt &=-2i\sin{\pi}z\int_{0}^{\infty}r^{z-1}e^{-r}dr\\ &=-2i\sin{\pi}z\Gamma(z)\qquad(\real{z}\gt 0)\\ \end{align}[/math]

である。しかし、左辺の被積分関数は[math]z[/math]が有界であるかぎり正則であるから、左辺は複素平面全体に解析接続する。従って、

[math]\Gamma(z)=\frac{i}{2\sin{\pi}z}\int_C(-t)^{z-1}e^{-t}dt\qquad(z\in\mathbb{C}\setminus\mathbb{Z})[/math]

である。[math]s=re^{i\theta}[/math]とすれば、同様にして

[math]\Gamma(z)=\frac{1}{e^{2{\pi}iz}-1}\int_Cs^{z-1}e^{-t}ds\qquad(z\in\mathbb{C}\setminus\mathbb{Z})[/math]

を得る。また、相反公式により、

[math]\frac{1}{\Gamma(z)}=\frac{\sin{\pi}z}{\pi}\Gamma(1-z)=\frac{i}{2\pi}\int_C(-t)^{-z}e^{-t}dt\qquad(z\in\mathbb{C})[/math]

を得る。

スターリングの公式

[math]z \to \infty[/math]での漸近展開として、ガンマ関数はスターリングの公式で近似される。この漸近近似は複素平面全体(負の実数を除く)で成立するが、[math]|{\arg z}|={\pi}[/math]に近づくにつれ近似の誤差が大きくなる(極限の収束が遅くなる)ため、応用上は相反公式などを用いて[math]|{\arg z}|\le{\pi}/2[/math]程度に制限することもある。

[math]\Gamma(z+1)\approx\sqrt{2{\pi}z}\left(\frac{z}{e}\right)^z\qquad(|{\arg z}|\lt {\pi},|z|\gg0)[/math]
[math]\lim_{z\to\infty}\frac{\Gamma(z+1)}{\sqrt{2{\pi}z}\left(\frac{z}{e}\right)^z}=1\qquad(|{\arg z}|\lt {\pi})[/math]

相反公式

次の恒等式をオイラーの相反公式(reflection formula)という。

[math]\Gamma(z)\Gamma(1-z)=-z\Gamma(z)\Gamma(-z)=\frac{\pi}{\sin{{\pi}z}}[/math]

相補公式とも呼ばれる[4]。 この恒等式はオイラーの乗積表示から得られる。

[math]\begin{align} -z\Gamma(z)\Gamma(-z) &=-z\left(\lim_{n\to\infty}\frac{n^zn!}{\prod_{k=0}^{n}{(z+k)}}\right)\left(\lim_{n\to\infty}\frac{n^{-z}n!}{\prod_{k=0}^{n}{(-z+k)}}\right)\\ &=\frac{1}{z}\prod_{k=1}^{\infty}\frac{k^2}{k^2-z^2}\\ &=\frac{\pi}{{\pi}z\displaystyle\prod_{k=1}^{\infty}\displaystyle\frac{k^2-z^2}{k^2}}\\ \end{align}[/math]

この分母は正弦関数の無限乗積展開であるから、

[math]\Gamma(z)\Gamma(1-z)=-z\Gamma(z)\Gamma(-z)=\frac{\pi}{\sin{{\pi}z}}[/math]

である。相反公式に[math]z=\frac{1}{2}[/math]を代入すれば

[math]\Gamma\left(\frac{1}{2}\right)\Gamma\left(1-\frac{1}{2}\right)=\frac{\pi}{\sin{\frac{\pi}{2}}}=\pi[/math]

となり

[math]\Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}[/math]

を得る。

乗法公式

次の恒等式をガウスの乗法公式(multiplication formula)という。

[math]\Gamma(nz)=\frac{n^{nz-1/2}}{(2\pi)^{(n-1)/2}}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}[/math]

この証明を示す。両辺の比を[math]f(z)[/math]とすると

[math]\begin{align}f(z)= &\frac{n^{nz-1/2}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}}{(2\pi)^{(n-1)/2}\Gamma(nz)}\\ \end{align}[/math]
[math]\begin{align}f(z+1) &=\frac{n^{nz-1/2}n^n\left[\prod_{k=0}^{n-1}\left(z+\frac{k}{n}\right)\Gamma{\left(z+\frac{k}{n}\right)}\right]}{(2\pi)^{(n-1)/2}\left[\prod_{k=0}^{n-1}(nz+k)\right]\Gamma(nz)}\\ &=\frac{n^{nz-1/2}\left[\prod_{k=0}^{n-1}\left(nz+k\right)\right]\prod_{k=0}^{n-1}\Gamma{\left(z+\frac{k}{n}\right)}}{(2\pi)^{(n-1)/2}\left[\prod_{k=0}^{n-1}(nz+k)\right]\Gamma(nz)}\\ &=f(z)\\ \end{align}[/math]

故に、任意に大きな自然数[math]m[/math]について[math]f(z+m)=f(z)[/math]が成立する。スターリングの公式により

[math]\begin{align}\lim_{\real{z}\to+\infty}f(z) &=\lim_{\real{z}\to+\infty}\frac{n^{nz-1/2}\left[\prod_{k=0}^{n-1}{\sqrt{\frac{2{\pi}}{z+k/n}}\left(\frac{z+k/n}{e}\right)^{z+k/n}}\right]}{(2\pi)^{(n-1)/2}\sqrt{\frac{2{\pi}}{nz}}\left(\frac{nz}{e}\right)^{nz}}\\ &=\lim_{\real{z}\to+\infty}z^{1/2}\left[\prod_{k=0}^{n-1}z^{k/n-1/2}(1+k/nz)^{z+k/n-1/2}e^{-k/n}\right]\\ &=\lim_{\real{z}\to+\infty}z^{1/2}\left[\prod_{k=0}^{n-1}z^{k/n-1/2}e^{k/n}e^{-k/n}\right]\\ &=1 \end{align}[/math]

途中で

[math]\lim_{\real{z}\to+\infty}(1+k/nz)^{z+k/n-1/2}=\lim_{\real{z}\to+\infty}(1+k/nz)^{z}=e^{n/k}[/math]

を適用した。

[math]f(z)=\lim_{n\to\infty}f(z+n)=1[/math]

であり、故に

[math]\Gamma(nz)=\frac{n^{nz-1/2}}{(2\pi)^{(n-1)/2}}\prod_{k=0}^{n-1}{\Gamma{\left(z+\frac{k}{n}\right)}}[/math]

が成立する。

微分方程式

[math](x,\ y,\ y_1,\ \ldots ,\ y_n)[/math]を変数とする多項式[math]F(x,\ y,\ y_1,\ \ldots ,\ y_n)[/math]に対し、

[math]F(x, y, y_1, \cdots, y_n)=0 , \quad y_i=\frac{d^i y}{dx^i} \quad (i=1, \cdots,n)[/math]

の形で表される微分方程式を代数的微分方程式という。ガンマ関数はいかなる代数的微分方程式も満たさないことが知られている。ヘルダーが1887年に最初に証明を与えた後 [5]E. H. ムーアEnglish版[6]A. オストロフスキEnglish版[7] [8]E. バーンズEnglish版[9]ハウスドルフ[10]により、別証明や一般化がなされた。

いくつかの具体的な値

[math]\Gamma\left(-\frac{3}{2}\right)\,= \frac {4\sqrt{\pi}} {3} \approx 2.363\,[/math]
[math]\Gamma\left(-\frac{1}{2}\right)\,= -2\sqrt{\pi} \approx -3.545\,[/math]
[math]\Gamma\left(\frac{1}{2}\right)\,= \sqrt{\pi} \approx 1.772\,[/math]
[math]\Gamma(1)\,=0!=1 \,[/math]
[math]\Gamma\left(\frac{3}{2}\right)\,= \frac {\sqrt{\pi}} {2} \approx 0.886\,[/math]
[math]\Gamma(2)\,=1!=1 \,[/math]
[math]\Gamma\left(\frac{5}{2}\right)\,= \frac {3 \sqrt{\pi}} {4} \approx 1.329\,[/math]
[math]\Gamma(3)\,=2!=2 \,[/math]
[math]\Gamma\left(\frac{7}{2}\right)\,= \frac {15\sqrt{\pi}} {8} \approx 3.323\,[/math]
[math]\Gamma(4)\,=3!=6 \,[/math]

ポリガンマ関数

ガンマ関数の対数微分

[math]\psi(z)=\frac{d}{dz}\log \Gamma(z)[/math]

ディガンマ関数(Digamma function)と呼ぶ。同様の対数微分を繰り返した関数

[math]\psi^{(n)}(z)=\frac{d^n}{dz^n}\log \Gamma(z)[/math]

を、ポリガンマ関数(Polygamma function)と呼ぶ。

脚注

  1. Wolfram mathworld: Gamma Function
  2. Springer Online Reference Works: Gamma-function
  3. Schmelzer & Trefethen (2007), Computing the Gamma function using contour integrals and rational approximations
  4. 神保 2003, 定理 5.15.
  5. Otto Ludwig Hölder, "Über die Eigenschaft der Gammafunction keiner algebraischen Differentialgleichung zu genügen," Math. Ann., 28, (1887) pp. 1–13. doi:10.1007/BF02430507
  6. Eliakim Hastings Moore, "Concerning transcendentally transcendental functions," Math. Ann., 48 (1897), pp. 49–74. doi:10.1007/BF01446334
  7. A. Ostrowski, "Neuer Beweis der Hölderschen Satzes, dass die Gammafunktion keiner algebraischen Differntialgleichung genügt." Math. Ann. 79 (1919), pp. 286–288. doi:10.1007/BF01458212
  8. A. Ostrowski, "Zum Hölderschen Satz über Γ(x). Math. Ann. 94 (1925), pp. 248–251. doi:10.1007/BF01208657
  9. E. W. Barnes, "The theory of the Gamma function," Messenger of Math. 29 (1900), pp. 64–128.
  10. F. Hausdorff, "Zum Hölderschen Satz über Γ(x)," Math. Ann. 94 (1925), pp. 244–247. doi:10.1007/BF01208656

参考文献

関連項目

外部リンク