第一種オイラー積分 / ベータ凾数

第一種オイラー積分 / ベータ凾数

\(p,q\geq0\)の整数の元で以下の式を第一種オイラー積分という． $$ \begin{eqnarray} \int_\alpha^\beta(x-\alpha)^p(\beta-x)^q \mathrm{d}x \end{eqnarray} $$ 部分積分を一回適用すると以下のように\(p+1\),\(q-1\)の第一種オイラー積分がでてくる． $$ \begin{eqnarray} \int_\alpha^\beta(x-\alpha)^p(\beta-x)^q \mathrm{d}x &=&\int_\alpha^\beta\left[\frac{\mathrm{d}}{\mathrm{d}x}\left\{\frac{1}{p+1}(x-\alpha)^{(p+1)}\right\}\right](\beta-x)^q \mathrm{d}x \;\cdots\;\frac{\mathrm{d}}{\mathrm{d}x}\left\{\frac{1}{p+1}(x-\alpha)^{(p+1)}\right\}=(x-\alpha)^p\\ &=&\left[\frac{1}{p+1}(x-\alpha)^{(p+1)}(\beta-x)^q \right]_\alpha^\beta -\int_\alpha^\beta\frac{1}{p+1}(x-\alpha)^{(p+1)}\left\{q(\beta-x)^{(q-1)}(-1)\right\} \mathrm{d}x \;\cdots\;\int f'g\;\mathrm{d}x= fg - \int fg' \mathrm{d}x\\ &=&0-\int_\alpha^\beta\frac{-q}{p+1}(x-\alpha)^{(p+1)}(\beta-x)^{(q-1)} \mathrm{d}x \;\cdots\;(\alpha-\alpha)=0, (\beta-\beta)=0\\ &=&\frac{q}{p+1}\int_\alpha^\beta(x-\alpha)^{(p+1)}(\beta-x)^{(q-1)} \mathrm{d}x\\ \end{eqnarray} $$ 第一種オイラー積分の乗数\(p,q\)に着目し元の式を\(I_{p,q}\),部分積分によって出てきた積分を\(I_{p+1,q-1}\)と表すと以下の様になる． $$ \begin{eqnarray} I_{p,q}&=&\frac{q}{p+1}I_{p+1,q-1} \;\cdots\;I_{a,b}=\int_\alpha^\beta(x-\alpha)^a(\beta-x)^b \mathrm{d}x\\ &=&\frac{q}{p+1}\frac{q-1}{p+2}I_{p+2,q-2}\\ &=&\frac{q}{p+1}\frac{q-1}{p+2}\cdots \frac{2}{p+(q-1)}\frac{1}{p+q}I_{p+q,0}\\ I_{p+q,0} &=&\int_\alpha^\beta(x-\alpha)^{(p+q)}(\beta-x)^0 \mathrm{d}x\\ &=&\int_\alpha^\beta(x-\alpha)^{(p+q)} \mathrm{d}x\\ &=&\left[\frac{1}{p+q+1}(x-\alpha)^{\left\{(p+q)+1\right\}}\right]_\alpha^\beta\\ &=&\frac{1}{p+q+1}(\beta-\alpha)^{(p+q+1)}\\ &=&\frac{(\beta-\alpha)^{(p+q+1)}}{p+q+1}\\ I_{p,q} &=&\frac{q}{p+1}\frac{q-1}{p+2}\cdots \frac{2}{p+q-1}\frac{1}{p+q}\frac{(\beta-\alpha)^{(p+q+1)}}{p+q+1}\\ &=&\frac{q(q-1)\cdots 2\;1}{(p+1)(p+2)\cdots (p+q-1)(p+q)(p+q+1)}(\beta-\alpha)^{(p+q+1)}\\ &=&\frac{q!}{\frac{(p+q+1)!}{p!}}(\beta-\alpha)^{(p+q+1)}\\ &=&\frac{p!\;q!}{(p+q+1)!}(\beta-\alpha)^{(p+q+1)}\\ \end{eqnarray} $$ 特に\(\alpha=0,\beta=1\)の時 $$ \begin{eqnarray} \int_\alpha^\beta(x-\alpha)^p(\beta-x)^q \mathrm{d}x &=&\int_0^1(x-0)^p(1-x)^q \mathrm{d}x\\ &=&\int_0^1x^p(1-x)^q \mathrm{d}x\\ &=&\frac{p!\;q!}{(p+q+1)!}(1-0)^{(p+q+1)} \;\cdots\;\href{}{\int_\alpha^\beta(x-\alpha)^p(\beta-x)^q \mathrm{d}x=\frac{p!\;q!}{(p+q+1)!}(\beta-\alpha)^{(p+q+1)}}\\ &=&\frac{p!\;q!}{(p+q+1)!}\\ \end{eqnarray} $$ 特に\(\alpha=0,\beta=1,p=m-1,q=n-1\)の時 $$ \begin{eqnarray} \int_\alpha^\beta(x-\alpha)^p(\beta-x)^q \mathrm{d}x &=&\int_0^1x^{(m-1)}(1-x)^{(n-1)} \mathrm{d}x\\ &=&\frac{(m-1)!\;(n-1)!}{((m-1)+(n-1)+1)!} \;\cdots\;\href{}{\int_0^1x^p(1-x)^q \mathrm{d}x=\frac{p!\;q!}{(p+q+1)!}}\\ &=&\frac{(m-1)!\;(n-1)!}{(m+n-1)!}\\ &=&\frac{\Gamma(m)\;\Gamma(n)}{\Gamma(m+n)}\;\cdots\;\Gamma(a)=(a-1)!\\ &=&B(m,n)\;\cdots\;ベータ凾数\\ \end{eqnarray} $$ 逆を言えば上記積分の値を1にしたい時はベータ凾数で割れば(逆数を掛ければ)よい(確率密度凾数との関係性)． $$ \begin{eqnarray} \frac{1}{B(m,n)}\int_0^1x^{(m-1)}(1-x)^{(n-1)} \mathrm{d}x&=&1\\ f(x;m,n)&=&\frac{x^{(m-1)}(1-x)^{(n-1)}}{B(m,n)}\;\cdots\;ベータ分布\\ \end{eqnarray} $$