ガンマ関数,ベータ関数

== ガンマ関数 ==

(1.1)　正の数 $p\gt 0$ に対して次の定積分で定義される関数を（オイラーの第2積分）ガンマ関数 $\Gamma(p)$ という．
$\Gamma(p)=\int_{0}^{\infty}x^{p-1}e^{-x}dx$
1) ガンマ関数について，次の関係が成り立つ． $\Gamma(p+ 1)=p\Gamma(p)$
2) $n$ が正の整数であるとき
$\Gamma(n)=(n-1)!$
3) $\Gamma(\frac{1}{2})=\sqrt{\pi}$

（証明）
1) $\Gamma(p+ 1)=\int_{0}^{\infty}x^{p}e^{-x}dx$
次の部分積分を行う

$f(x)=x^{p}\rightarrow f\apos (x)=px^{p-1}$
$g(x)=-e^{-x}\leftarrow g\apos (x)=e^{-x}$

$\Gamma(p+ 1)=\Bigl[-x^{p}e^{-x}\Bigr]_0^{\infty}+ p\int_0^{\infty}x^{p-1}e^{-x}dx$
$=p\int_0^{\infty}x^{p-1}e^{-x}dx=p\Gamma(p)$
2) 1)の結果を使うと
$\Gamma(n)=(n-1)\Gamma(n-1)=(n-1)(n-2)\Gamma(n-2)$
$=(n-1)(n-2)\cdot\cdot\cdot2\cdot 1\cdot\Gamma(1)$
ここで
$\Gamma(1)=\int_{0}^{\infty}e^{-x}dx=\Bigl[-e^{-x}\Bigr]_0^{\infty}=1$
だから
$\Gamma(n)=(n-1)!$
3) 後出(1.3) 1)により示される．

== ベータ関数 ==

(1.2) $p,\hspace{2}q\gt 0$ を正の数とするとき，次の定積分で定義される関数を（オイラーの第1積分）ベータ関数という．
$B(p,\hspace{2}q)=\int_{0}^{1}x^{p-1}(1-x)^{q-1}dx$
1) $B(p,\hspace{2}q)=B(q,\hspace{2}p)$ が成り立つ
2)　置換積分により，次の式を示せる． $B(p,\hspace{2}q)=2\int_0^{\frac{\pi}{2}}\sin^{2p-1}\theta\cos^{2q-1}\theta d\theta$
3)　ベータ関数はガンマ関数と次の関係にある．
$B(p,\hspace{2}q)=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+ q)}$
4)　ベータ関数$B(m, n)$は$m, n\gt 0$で定義されるが，特に$m, n$が正（または０）の整数のとき，次のように階乗で表せる．
$B(m,\hspace{2}n)=\frac{(m-1)!(n-1)!}{(m+ n-1)!}$

（証明）
1)
$B(p,\hspace{2}q)=\int_{0}^{1}x^{p-1}(1-x)^{q-1}dx$
において，次の置換積分を行う

$1-x=t$
$dx=-dt$

$x$ $ \displaystyle x $	$0\rightarrow 1$ $ \displaystyle 0\rightarrow 1 $
$t$ $ \displaystyle t $	$1\rightarrow 0$ $ \displaystyle 1\rightarrow 0 $

$B(p,\hspace{2}q)=\int_{1}^{0}(1-t)^{p-1}t^{q-1}(-dt)$
$=\int_{0}^{1}t^{q-1}(1-t)^{p-1}dt=B(q,\hspace{2}p)$

2)
$B(p,\hspace{2}q)=\int_{0}^{1}x^{p-1}(1-x)^{q-1}dx$
に対して次の置換積分を行う

$x=\sin^2\theta$
$1-x=1-\sin^2\theta=\cos^2\theta$
$\frac{dx}{d\theta}=2\sin\theta\cos\theta$

$x$ $ \displaystyle x $	$0\rightarrow 1$ $ \displaystyle 0\rightarrow 1 $
$\theta$ $ \displaystyle \theta $	$0\rightarrow\frac{\pi}{2}$

$B(p,\hspace{2}q)=\int_{0}^{\frac{\pi}{2}}\sin^{2p-2}\theta\cos^^{2q-2}\cdot 2\sin\theta\cos\theta d\theta$
$=2\int_{0}^{\frac{\pi}{2}}\sin^{2p-1}\theta\cos^^{2q-1}d\theta$

$\Gamma(p)\Gamma(q)=\int_{0}^{\infty}x^{p-1}e^{-x}dx\int_{0}^{\infty}y^{q-1}e^{-y}dy$
のような2つの積分の積を
$B(p,\hspace{2}q)=\int_{0}^{1}x^{p-1}(1-x)^{q-1}dx$
と結び付けるには，１つの重積分として直交座標から極座標へ変換する方法が考えられる．その場合， $x^2+ y^2=r^2,\hspace{5}dxdy=rdrd\theta$
とするために，2乗の形をした式が使いやすい．

$\Gamma(p)=\int_{0}^{\infty}t^{p-1}e^{-t}dt$
において， $t=x^2$ の変換を行うと

$\frac{dt}{dx}=2x$

$t$ $ \displaystyle t $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $
$x$ $ \displaystyle x $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $

$\Gamma(p)=\int_{0}^{\infty}x^{2(p-1)}e^{-x^2}\cdot 2xdx$
$=2\int_{0}^{\infty}x^{2p-1}e^{-x^2}dx$
同様にして
$\Gamma(q)=2\int_{0}^{\infty}y^{2q-1}e^{-y^2}dy$
これらの積を求めると
$\Gamma(p)\Gamma(q)=4\int_{0}^{\infty}x^{2p-1}e^{-x^2}dx\int_{0}^{\infty}y^{2q-1}e^{-y^2}dy$
$=4\int_{0}^{\infty}\int_{0}^{\infty}x^{2p-1}y^{2q-1}e^{-(x^2+ y^2)}dxdy$

ここで， $x=r\cos\theta,\hspace{2}y=r\sin\theta$ の関係を用いて，極座標に変換する
$\frac{\partial x}{\partial r}=\cos\theta,\hspace{20}\frac{\partial y}{\partial r}=\sin\theta$
$\frac{\partial x}{\partial\theta}=-r\sin\theta,\hspace{2}\frac{\partial y}{\partial\theta}=r\cos\theta$
ヤコビアンは
$J=\frac{\partial(x,y)}{\partial(r,\theta)}=\begin{vmatrix}\frac{\partial x}{\partial r}&\frac{\partial x}{\partial\theta}\\ \frac{\partial y}{\partial r}&\frac{\partial y}{\partial\theta}\end{vmatrix}$ $=\begin{vmatrix}\cos\theta& -r\sin\theta\\ \sin\theta& r\cos\theta\end{vmatrix}$
$=r\cos^2\theta+ r\sin^2\theta=r$
だから
$dxdy=\mid J\mid drd\theta=rdrd\theta$
$x,\hspace{3}y\hspace{2}:\hspace{2}0\rightarrow\infty$ のとき $r\hspace{2}:\hspace{2}0\rightarrow\infty,\hspace{3}\theta\hspace{2}:\hspace{2}0\rightarrow\frac{\pi}{2}$

$\Gamma(p)\Gamma(q)$
$=4\int_{0}^{\infty}\int_{0}^{\infty}r^{2p-1}\cos^{2p-1}\theta r^{2q-1}\sin^{2q-1}\theta e^{-r^2}\cdot rdrd\theta$
$=2\int_{0}^{\infty}r^{2(p+ q)-1}\cdot e^{-r^2}dr\cdot 2\int_{0}^{\frac{\pi}{2}}\cos^{2p-1}\theta\cdot\sin^{2q-1}\theta d\theta$
$=\Gamma(p+ q)B(p,\hspace{2}q)$

4)
(1.1) 2)により， $m,\hspace{2}n$ が正の整数であるとき
$\Gamma(m)=(m-1)!,\hspace{2}\Gamma(n)=(n-1)!$
$\Gamma(m+ n)=(m+ n-1)!$
だから
$B(m,\hspace{2}n)=\frac{\Gamma(m)\Gamma(n)}{\Gamma(m+ n)}=\frac{(m-1)!(n-1)!}{(m+ n-1)!}$

(1.3) 誤差関数（error function→erf）
$\rm{erf}(x)=\int_0^{x}e^{-t^2}dt$
1) ガウスの積分公式
$\int_0^{\infty}e^{-x^2}dx=\frac{\sqrt{\pi}}{2}$
2) $\int_0^{\infty}\frac{e^{-x}}{\sqrt{x}}dx=\sqrt{\pi}$
3) $\int_{-\infty}^{\infty}\frac{e^{-\frac{x^2}{2}}}{\sqrt{2\pi}}dx=1$
4) $\int_{0}^{\infty}e^{-kx^2}dx=\frac{1}{2}\sqrt{\frac{\pi}{k}}\hspace{10}(k\gt 0)$

(1.3)の誤差関数は，初等的に（多項式，分数関数，無理関数，指数関数，対数関数，三角関数，およびそれらの有限回の組合せや合成）では書けない．
　右図のような正規分布（ただし定数倍の調整あり）で0～xまでの面積を表す．
1)　誤差関数で積分区間が0～無限の定積分は，ガウスの積分公式と呼ばれ，値を求めることができる．
　重積分を利用する次の方法がよく見られる．
$I=\int_0^{\infty}e^{-x^2}dx,\hspace{20}I=\int_0^{\infty}e^{-y^2}dy$
とおくと
$I^2=\int_0^{\infty}\int_0^{\infty}e^{-(x^2+ y^2)}dxdy$

$I^2=\int_0^{\infty}\int_0^{\frac{\pi}{2}}e^{-r^2}rdrd\theta=\frac{\pi}{2}\int_0^{\infty}e^{-r^2}rdr$

$z=r^2$ とおく置換積分を行うと
$\frac{dz}{dr}=2r$

$r$ $ \displaystyle r $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $
$z$ $ \displaystyle z $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $

$I^2=\frac{\pi}{2}\int_0^{\infty}e^{-z}r\frac{dz}{2r}=\frac{\pi}{4}\int_0^{\infty}e^{-z}dz$
$=\frac{\pi}{4}\Bigl[-e^{-z}\Bigr]_0^{\infty}=\frac{\pi}{4}(0+ 1)=\frac{\pi}{4}$
$I=\frac{\sqrt{\pi}}{2}$

2)←
1)において， $x^2=t$ とおく置換積分を行うと

$\frac{dt}{dx}=2x=2\sqrt{t}$

$x$ $ \displaystyle x $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $
$t$ $ \displaystyle t $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $

$\int_0^{\infty}e^{-x^2}dx=\int_0^{\infty}e^{-t}\frac{dt}{2\sqrt{t}}=\frac{1}{2}\int_0^{\infty}\frac{e^{-t}}{\sqrt{t}}dt$
したがって，1)の結果から
$\frac{1}{2}\int_0^{\infty}\frac{e^{-t}}{\sqrt{t}}dt=\frac{\sqrt{\pi}}{2}$
$\int_0^{\infty}\frac{e^{-t}}{\sqrt{t}}dt=\sqrt{\pi}$

3)←
1)において， $x=\frac{t}{\sqrt{2}}$ とおく置換積分を行うと

$\frac{dx}{dt}=\frac{1}{\sqrt{2}}$

$x$ $ \displaystyle x $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $
$t$ $ \displaystyle t $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $

$\int_0^{\infty}e^{-x^2}dx=\int_0^{\infty}e^{-\frac{t^2}{2}}\frac{dt}{\sqrt{2}}=\int_0^{\infty}e^{-\frac{t^2}{2}}\frac{dt}{\sqrt{2}}$
$=\frac{1}{\sqrt{2}}\int_0^{\infty}e^{-\frac{t^2}{2}}dt$
1)の結果から
$\frac{1}{\sqrt{2}}\int_0^{\infty}e^{-\frac{t^2}{2}}dt=\frac{\sqrt{\pi}}{2}$
$\int_0^{\infty}e^{-\frac{x^2}{2}}dx=\frac{\sqrt{\pi}}{\sqrt{2}}$
偶関数だから
$\int_{-\infty}^{0}e^{-\frac{x^2}{2}}dx=\frac{\sqrt{\pi}}{\sqrt{2}}$
結局
$\int_{-\infty}^{\infty}e^{-\frac{x^2}{2}}dx=\sqrt{2\pi}$
$\int_{-\infty}^{\infty}\frac{e^{-\frac{x^2}{2}}}{\sqrt{2\pi}}dx=1$

4)←
1)において， $x=\sqrt{k}t$ とおく置換積分を行うと

$\frac{dx}{dt}=\sqrt{k}$

$x$ $ \displaystyle x $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $
$t$ $ \displaystyle t $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $

$\int_0^{\infty}e^{-x^2}dx=\int_0^{\infty}e^{-kt^2}\sqrt{k}dt$
1)の結果から
$=\frac{\sqrt{\pi}}{2}$
したがって
$\int_{0}^{\infty}e^{-kx^2}dx=\frac{1}{\sqrt{k}}\cdot\frac{\sqrt{\pi}}{2}=\frac{1}{2}\sqrt{\frac{\pi}{k}}\hspace{10}(k\gt 0)$

Γ関数，B関数について，次の値を求めてください．

【問題1-1】
$\Gamma(6)$

（解答）
(1.1) 2)により， $n$ が正の整数であるとき
$\Gamma(n)=(n-1)!$
だから
$\Gamma(6)=5!=120$ …（答）

【問題1-2】
$\Gamma(\frac{3}{2})$

（解答）
(1.1) 1)により， $\Gamma(\frac{1}{2}+ 1)=\frac{1}{2}\Gamma(\frac{1}{2})$
さらに，(1.1) 3)により， $\Gamma(\frac{1}{2})=\sqrt{\pi}$
だから， $\Gamma(\frac{3}{2})=\frac{\sqrt{\pi}}{2}$ …（答）

【問題1-3】
$B(4,\hspace{2}3)$

（解答）
(1.2) 4)により， $B(4,\hspace{2}3)=\frac{3!2!}{6!}=\frac{1}{60}$ …（答）

【問題1-4】
$B(\frac{1}{2},\hspace{2}\frac{3}{2})$

（解答）
(1.2) 3)により，
$B(\frac{1}{2},\hspace{2}\frac{3}{2})=\frac{\Gamma(\frac{1}{2})\Gamma(\frac{3}{2})}{\Gamma(\frac{1}{2}+\frac{3}{2})}=\frac{\sqrt{\pi}\cdot\frac{\sqrt{\pi}}{2}}{1}=\frac{\pi}{2}$ …（答）

【問題1-5】
$\int_0^{\infty}e^{-x^n}dx\hspace{2}(n\gt 0)$ をΓ関数を用いて表してください．

（解答）
ガンマ関数の定義
$\Gamma(p)=\int_{0}^{\infty}t^{p-1}e^{-t}dt$
において， $t=x^n$ と変換する置換積分を行うと

$\frac{dt}{dx}=nx^{n-1}$

$t$ $ \displaystyle t $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $
$x$ $ \displaystyle x $	$0\rightarrow\infty$ $ \displaystyle 0\rightarrow\infty $

$\Gamma(p)=\int_{0}^{\infty}t^{p-1}e^{-t}dt=\int_{0}^{\infty}x^{n(p-1)}e^{-x^n}\cdot nx^{n-1}dx$
$=n\int_{0}^{\infty}x^{np-1}e^{-x^n}dx$
$p=\frac{1}{n}$ のとき
$\Gamma(\frac{1}{n})=n\int_{0}^{\infty}e^{-x^n}dx$
だから
$\int_{0}^{\infty}e^{-x^n}dx=\frac{1}{n}\Gamma(\frac{1}{n})$ …（答）

【問題1-6】
$\int_{\alpha}^{\beta}(x-\alpha)^{p-1}(\beta-x)^{q-1}dx$
をΓ関数を用いて表してください．
ただし， $p,\hspace{2}q\gt 0,\hspace{2}\alpha\lt\beta$ とする．

（解答）
ベータ関数の定義およびガンマ関数との関係式
$B(p,\hspace{2}q)=\int_{0}^{1}t^{p-1}(1-t)^{q-1}dt=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+ q)}$
において
$t=\frac{x-\alpha}{\beta-\alpha}$
$1-t=1-\frac{x-\alpha}{\beta-\alpha}=\frac{\beta-x}{\beta-\alpha}$
の変換により置換積分を行うと

$\frac{dt}{dx}=\frac{1}{\beta-\alpha}$

$t$ $ \displaystyle t $	$0\rightarrow 1$ $ \displaystyle 0\rightarrow 1 $
$x$ $ \displaystyle x $	$\alpha\rightarrow\beta$

$B(p,\hspace{2}q)=\int_{0}^{1}t^{p-1}(1-t)^{q-1}dt$
$=\int_{\alpha}^{\beta}\left(\frac{x-\alpha}{\beta-\alpha}\right)^{p-1}\left(\frac{\beta-x}{\beta-\alpha}\right)^{q-1}\frac{dx}{\beta-\alpha}$
$=\frac{1}{(\beta-\alpha)^{p+ q-1}}\int_{\alpha}^{\beta}(x-\alpha)^{p-1}(\beta-x)^{q-1}dx$
$=\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+ q)}$
だから
$\int_{\alpha}^{\beta}(x-\alpha)^{p-1}(\beta-x)^{q-1}dx$
$=(\beta-\alpha)^{p+ q-1}\frac{\Gamma(p)\Gamma(q)}{\Gamma(p+ q)}$ …（答）

【問題1-7】
$\int_{0}^{1}\frac{dx}{\sqrt{1-x^3}}dx$
をベータ関数を使って表してください．

（解答）
ベータ関数の定義式
$B(p,\hspace{2}q)=\int_{0}^{1}t^{p-1}(1-t)^{q-1}dt$
において
$t=x^3$
の変換により置換積分を行うと

$\frac{dt}{dx}=3x^2$

$t$ $ \displaystyle t $	$0\rightarrow 1$ $ \displaystyle 0\rightarrow 1 $
$x$ $ \displaystyle x $	$0\rightarrow 1$ $ \displaystyle 0\rightarrow 1 $

$B(p,\hspace{2}q)=\int_{0}^{1}x^{3p-3}(1-x^3)^{q-1}\cdot 3x^2dx$
$=3\int_{0}^{1}x^{3p-1}(1-x^3)^{q-1}dx$
ここで
$p=\frac{1}{3},\hspace{5}q=\frac{1}{2}$
を代入すると
$B(\frac{1}{3},\hspace{2}\frac{1}{2})=3\int_{0}^{1}\frac{dx}{\sqrt{1-x^3}}$
$\int_{0}^{1}\frac{dx}{\sqrt{1-x^3}}=\frac{1}{3}B(\frac{1}{3},\hspace{2}\frac{1}{2})$ …（答）

【問題1-8】
$\int_{0}^{1}\frac{dx}{\sqrt{1-x^4}}dx$
をベータ関数を使って表してください．

（解答）
前問とほぼ同様の途中経過でできます．
$\int_{0}^{1}\frac{dx}{\sqrt{1-x^4}}=\frac{1}{4}\cdot B(\frac{1}{4},\hspace{2}\frac{1}{2})$ …（答）

【問題1-9】
$\int_{0}^{\frac{\pi}{2}}\sin^{n}\theta d\theta$
をベータ関数を使って表してください．ただし， $n\gt 0$ とする．

（解答）
(1.2) 2)の式
$B(p,\hspace{2}q)=2\int_0^{\frac{\pi}{2}}\sin^{2p-1}\theta\cos^{2q-1}\theta d\theta$
において
$p=\frac{n+ 1}{2},\hspace{5}q=\frac{1}{2}$
を代入すると
$B(\frac{n+ 1}{2},\hspace{2}\frac{1}{2})=2\int_0^{\frac{\pi}{2}}\sin^{n}\theta d\theta$ …（答）
だから
$\int_0^{\frac{\pi}{2}}\sin^{n}\theta d\theta=\frac{1}{2}B(\frac{n+ 1}{2},\hspace{2}\frac{1}{2})$

【問題1-10】
　 $p,\hspace{2}q\gt 0$ を正の数とするとき，次の関係式が成り立つことを証明してください．
$B(p,\hspace{2}q)=\frac{q-1}{p}B(p+ 1,\hspace{2}q-1)$

（解答）
$B(p,\hspace{2}q)=\int_{0}^{1}x^{p-1}(1-x)^{q-1}dx$
に対して次の部分積分で変換する．

$f\apos(x)=x^{p-1}\rightarrow f(x)=\frac{x^p}{p}$

$g(x)=(1-x)^{q-1}\rightarrow g\apos(x)=-(q-1)(1-x)^{q-2}$

$B(p,\hspace{2}q)=\Bigl[\frac{x^p}{p}(1-x)^{q-1}\Bigr]_0^1+\frac{q-1}{p}\int_{0}^{1}x^{p}(1-x)^{q-2}dx$
$=\frac{q-1}{p}\int_{0}^{1}x^{p}(1-x)^{q-2}dx=\frac{q-1}{p}B(p+ 1,\hspace{2}q-1)$

●==目次に戻る