🟠

ガンマ関数の積分公式

2023/10/26に公開

\int_0^{\infty} \cdots \int_0^{\infty} x_1^{\alpha_1-1} \cdots x_n^{\alpha_n-1} dx_1 \cdots dx_{n-1} = \frac{\prod_{i=1}^n\Gamma(\alpha_i)}{\Gamma(\sum_{i=1}^n\alpha_i)}

ただし、左辺の定積分の範囲は $x_i \geq 0 \hspace{3pt} (i = 1,...,n)$ および $\sum_{i=1}^n x_i =1$ とする。

注意：
$\sum_{i=1}^n x_i =1$ の制約によって $x_n= 1-\sum_{i=1}^{n-1} x_i$ となるので独立な変数は $x_1,...,x_{n-1}$
そのため積分変数も $x_{n-1}$ までとなっている。

証明

ガンマ関数

$\alpha > 0$ に対して、ガンマ関数 $\Gamma(\alpha)$ は次のように定義されます。

\Gamma(\alpha) = \int_0^{\infty} x^{\alpha-1} e^{-x} dx

ガンマ関数の定義より、

\prod_{i=1}^n\Gamma(\alpha_i) = \int_0^{\infty} \cdots \int_0^{\infty} e^{-\sum_{i=1}^nt_i} t_1^{\alpha_1-1} \cdots t_n^{\alpha_n-1} dt_1 \cdots dt_n

ここで、変数変換

t_i = x_i z \hspace{3pt} ,\hspace{3pt} x_i\ge 0 \hspace{3pt}(i=1,...,n), \hspace{3pt} \sum_{i=1}^n x_i =1

とする。

ヤコビアンを考える。

\begin{align*} J &= \left|\begin{matrix}\partial t_1 / \partial x_1 & \cdots &\partial t_1/ \partial x_{n-1} & \partial t_1 / \partial z\\ & \vdots\\ \partial t_{n-1} / \partial x_1 & \cdots &\partial t_{n-1}/ \partial x_{n-1} & \partial t_{n-1} / \partial z\\ \partial t_n / \partial x_1 & \cdots &\partial t_n/ \partial x_{n-1} & \partial t_n / \partial z \end{matrix} \right| \\ &=\left|\begin{matrix} z & 　 & 　 & & x_1\\ & \ddots & \text{\huge{0}} \\ & \text{\huge{0}} & \ddots & & \vdots \\ & & & z & x_{n-1}\\ -z & \cdots & \cdots & -z & x_n \end{matrix} \right| \\ &= \left|\begin{matrix} z & & & & x_1\\ & \ddots & \text{\huge{0}}& \\ & \text{\huge{0}} & \ddots & & \vdots\\ & & & z & x_{n-1}\\ 0 & &\cdots & 0 & 1 \end{matrix} \right| \\ & = z^{n-1} \end{align*}

ヤコビアンは $z^{n-1}$ であることが分かる。

したがって

\begin{align*} \begin{split} \prod_{i=1}^n\Gamma(\alpha_i) = &\int_0^{\infty} e^{-z\sum_{i=1}^nx_i} z^{\sum_{i=1}^n\alpha_i - n} z^{n-1}dz \\ &\times \int_0^\infty \cdots \int_0^\infty x_1^{\alpha_1-1} \cdots x_{n-1}^{\alpha_{n-1}-1} x_n^{\alpha_n-1} dx_1 \cdots dx_{n-1} \end{split} \end{align*}

1つ目の積分は $\sum_{i=1}^n x_i =1$ とガンマ関数の定義より

\begin{align*} \int_0^{\infty} e^{-z\sum_{i=1}^nx_i} z^{\sum_{i=1}^n\alpha_i - n} z^{n-1}dz &=\int_0^{\infty} e^{-z} z^{\sum_{i=1}^n\alpha_i - 1} dz\\ &=\Gamma\left(\sum_{i=1}^n\alpha_i\right) \end{align*}

2つ目の積分は「積分公式」の左辺と等しい。

以上より

\int_0^{\infty} \cdots \int_0^{\infty} x_1^{\alpha_1-1} \cdots x_n^{\alpha_n-1} dx_1 \cdots dx_{n-1} = \frac{\prod_{i=1}^n\Gamma(\alpha_i)}{\Gamma(\sum_{i=1}^n\alpha_i)}

となる。

ディリクレ分布の期待値や分散を求めるのに必要な公式を導出しました。

結構複雑で、大変でした。

最後に、どなたかMarkdownで枠を作ったりする方法ご存知でしたら、ご教示いただけると幸いです。Notionでいうコールアウトみたいにして定義などをみやすくしたいです。