ベータ関数

同義語:Beta関数B関数beta function第一種Euler積分

概要

ベータ関数(beta function)とは、$\operatorname{Re}x>0$, $\operatorname{Re}y>0$ に対し $B(x,y)=\int_0^1t^{x-1}(1-t)^{y-1}\,dt$ で定まる関数であり、正の整数では $B(m,n)=(m-1)!\,(n-1)!/(m+n-1)!$ と階乗で表される。Γ関数によって $B(x,y)=\Gamma(x)\Gamma(y)/\Gamma(x+y)$ と書け、三角関数の冪の積分や $\int_0^\infty w^{x-1}(1+w)^{-x-y}\,dw$ の形の積分はこれで計算できる。置換から $\Gamma(1/2)=\sqrt\pi$ や Legendre の倍数公式が得られ、相反公式と組み合わせるとレムニスケートの弧長などの値も求まる。ベータ分布や $t$ 分布の正規化定数にも現れる。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: Γ関数, 広義積分, 置換積分, Fubiniの定理

ベータ関数は、区間 $[0,1]$ 上の積分 $B(x,y)=\int_0^1t^{x-1}(1-t)^{y-1}\,dt$ で定まる 2 変数の関数である。正の整数では $B(m,n)=\frac{(m-1)!\,(n-1)!}{(m+n-1)!}$ と階乗で書け、放物線と直線で囲まれた面積の「$\frac16$ 公式」はその特別な場合 $B(2,2)=\frac16$ にあたる。
一般の $x,y$ では、ベータ関数は Γ関数によって
$$ B(x,y)=\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} $$
と表される。この関係により、ベータ関数に帰着できる積分(三角関数の冪の積分、$\int_0^\infty\frac{w^{x-1}}{(1+w)^{x+y}}dw$ の形の積分など)はすべて Γ関数の値で計算でき、逆にベータ関数の置換から $\Gamma(\frac12)=\sqrt\pi$ や Γ関数の倍数公式が得られる。

定義

複素数の冪は、$t>0$ と複素数 $x$ に対し $t^{x}:=e^{x\log t}$ と定める。$|t^{x}|=t^{\operatorname{Re}x}$ である。

ベータ関数

$\operatorname{Re}x>0$, $\operatorname{Re}y>0$ を満たす複素数 $x,y$ に対し
$$ B(x,y):=\int_0^1t^{x-1}(1-t)^{y-1}\,dt $$
と定める。$B$ をベータ関数(beta function)といい、右辺の積分を第 1 種 Euler 積分またはベータ積分という。積分は prop-bt-convergence により絶対収束する。

積分の収束域

複素数 $x,y$ について、$\int_0^1|t^{x-1}(1-t)^{y-1}|\,dt<\infty$ となるのは $\operatorname{Re}x>0$ かつ $\operatorname{Re}y>0$ のときに限る。

区間を半分に分ける

$a:=\operatorname{Re}x$, $b:=\operatorname{Re}y$ とおくと、被積分関数の絶対値は $t^{a-1}(1-t)^{b-1}$ で、$(0,1)$ 上連続である。$0< t\le\frac12$ では $(1-t)^{b-1}$ は $\frac12\le1-t<1$ の範囲の連続関数なので、正の定数 $c_1\le c_2$ で $c_1\le(1-t)^{b-1}\le c_2$ と挟める。$\int_0^{1/2}t^{a-1}\,dt$ は $a>0$ のとき $\frac{2^{-a}}a$ に収束し、$a\le0$ のとき $+\infty$ に発散するので、広義積分の比較判定法により $\int_0^{1/2}$ の部分が有限であることは $a>0$ と同値である。置換 $t\mapsto1-t$ により、$\int_{1/2}^1$ の部分が有限であることは $b>0$ と同値である。$\square$

基本的な性質

対称性と漸化式

$\operatorname{Re}x>0$, $\operatorname{Re}y>0$ とする。

  1. $B(x,y)=B(y,x)$。
  2. $B(x,1)=\dfrac1x$。
  3. $B(x,y)=B(x+1,y)+B(x,y+1)$。
  4. $B(x+1,y)=\dfrac{x}{x+y}\,B(x,y)$。
置換と部分積分

1 は置換 $t\mapsto1-t$ による。2 は $\int_0^1t^{x-1}\,dt=\bigl[t^x/x\bigr]_0^1=\frac1x$($|t^x|=t^{\operatorname{Re}x}\to0$)による。3 は被積分関数の等式 $t^{x-1}(1-t)^{y-1}=t^{x}(1-t)^{y-1}+t^{x-1}(1-t)^{y}$ を積分すればよい。
4:$0<\varepsilon<\frac12$ として $[\varepsilon,1-\varepsilon]$ 上で部分積分すると
$$ \int_\varepsilon^{1-\varepsilon}t^{x}(1-t)^{y-1}\,dt=\Bigl[-\frac{t^{x}(1-t)^{y}}{y}\Bigr]_\varepsilon^{1-\varepsilon}+\frac xy\int_\varepsilon^{1-\varepsilon}t^{x-1}(1-t)^{y}\,dt $$
である。境界項の絶対値は $t^{\operatorname{Re}x}(1-t)^{\operatorname{Re}y}/|y|$ で、$t\to0$ でも $t\to1$ でも $0$ に近づく。$\varepsilon\to0$ として $B(x+1,y)=\frac xy\,B(x,y+1)$ を得る。3 により $B(x,y+1)=B(x,y)-B(x+1,y)$ なので、$B(x+1,y)\bigl(1+\frac xy\bigr)=\frac xyB(x,y)$、すなわち 4 を得る。$\square$

正の整数での値

正の整数 $m,n$ に対し
$$ B(m,n)=\frac{(m-1)!\,(n-1)!}{(m+n-1)!} $$
である。

$m$ についての帰納法

$m=1$ では prop-bt-basic の 1 と 2 により $B(1,n)=\frac1n=\frac{0!\,(n-1)!}{n!}$ である。$m$ で成り立てば、prop-bt-basic の 4 により $B(m+1,n)=\frac m{m+n}\cdot\frac{(m-1)!\,(n-1)!}{(m+n-1)!}=\frac{m!\,(n-1)!}{(m+n)!}$ である。$\square$

$\frac16$ 公式

実数 $\alpha<\beta$ と 0 以上の整数 $m,n$ に対し、置換 $x=\alpha+(\beta-\alpha)t$ により $x-\alpha=(\beta-\alpha)t$, $\beta-x=(\beta-\alpha)(1-t)$ なので
$$ \int_\alpha^\beta(x-\alpha)^m(\beta-x)^n\,dx=(\beta-\alpha)^{m+n+1}B(m+1,n+1)=\frac{m!\,n!}{(m+n+1)!}(\beta-\alpha)^{m+n+1} $$
である。$m=n=1$ とすると $\int_\alpha^\beta(x-\alpha)(x-\beta)\,dx=-\frac{(\beta-\alpha)^3}6$ となり、放物線と直線で囲まれた面積の公式を得る(高校数学での扱いは 1/6公式とベータ積分)。

ほかの積分表示

$\operatorname{Re}x>0$, $\operatorname{Re}y>0$ に対し
$$ B(x,y)=2\int_0^{\pi/2}(\sin\theta)^{2x-1}(\cos\theta)^{2y-1}\,d\theta=\int_0^\infty\frac{w^{x-1}}{(1+w)^{x+y}}\,dw $$
が成り立つ($0<\theta<\frac\pi2$ で $\sin\theta,\cos\theta>0$ なので、冪は上の約束で定まる)。

2 つの置換

1 つ目:$t=\sin^2\theta$ は $(0,\frac\pi2)$ から $(0,1)$ への狭義単調増加な $C^1$ 級の全単射で、$dt=2\sin\theta\cos\theta\,d\theta$, $1-t=\cos^2\theta$ である。よって $t^{x-1}(1-t)^{y-1}\,dt=2(\sin\theta)^{2x-1}(\cos\theta)^{2y-1}\,d\theta$ となる。
2 つ目:$t=\frac w{1+w}$ は $(0,\infty)$ から $(0,1)$ への狭義単調増加な $C^1$ 級の全単射で、$1-t=\frac1{1+w}$, $dt=\frac{dw}{(1+w)^2}$ である。よって $t^{x-1}(1-t)^{y-1}\,dt=w^{x-1}(1+w)^{-(x-1)-(y-1)-2}\,dw=w^{x-1}(1+w)^{-x-y}\,dw$ となる。
どちらも、両端を除いた有界閉区間で置換積分を行い、端点の極限をとればよい(積分は prop-bt-convergence により絶対収束する)。$\square$

主定理:Γ関数による表示

Γ関数による表示

$\operatorname{Re}x>0$, $\operatorname{Re}y>0$ に対し
$$ B(x,y)=\frac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} $$
が成り立つ(DLMF 5.12.1、Rud76 Theorem 8.20)。

二重積分の順序を入れ替える

$a:=\operatorname{Re}x$, $b:=\operatorname{Re}y$ とおく。$(0,\infty)$ 上の連続関数で絶対値の広義積分が有限なものについては、広義積分は Lebesgue 積分に一致する(単調収束定理と優収束定理による)。以下この同一視のもとで Fubiniの定理 を使う。
$c>0$ と $\operatorname{Re}z>0$ に対し、置換 $s=cw$ により
$$ \int_0^\infty w^{z-1}e^{-cw}\,dw=\frac{\Gamma(z)}{c^{z}} $$
である(Γ関数の例「尺度変換と Laplace 型積分」)。$t>0$ を固定してこの式を $z=x$, $c=t$ に使うと $\Gamma(x)=t^{x}\int_0^\infty w^{x-1}e^{-tw}\,dw$ なので、
$$ \Gamma(x)\Gamma(y)=\int_0^\infty t^{y-1}e^{-t}\,\Gamma(x)\,dt=\int_0^\infty\Bigl(\int_0^\infty t^{x+y-1}w^{x-1}e^{-t(1+w)}\,dw\Bigr)dt $$
である。被積分関数 $H(t,w):=t^{x+y-1}w^{x-1}e^{-t(1+w)}$ は $(0,\infty)^2$ 上連続で、$|H(t,w)|=t^{a+b-1}w^{a-1}e^{-t(1+w)}$ である。同じ計算を実数 $a,b$ で行うと、Tonelli の定理(Fubiniの定理の記事の定理「Tonelliの定理」)により $\iint|H|\,dt\,dw=\Gamma(a)\Gamma(b)<\infty$ となるので、$H$ は $(0,\infty)^2$ 上可積分である。よって積分の順序を入れ替えてよく、同じ式を $z=x+y$, $c=1+w$ に使うと
$$ \Gamma(x)\Gamma(y)=\int_0^\infty w^{x-1}\Bigl(\int_0^\infty t^{x+y-1}e^{-(1+w)t}\,dt\Bigr)dw=\Gamma(x+y)\int_0^\infty\frac{w^{x-1}}{(1+w)^{x+y}}\,dw=\Gamma(x+y)\,B(x,y) $$
を得る(最後は prop-bt-other-forms)。$\Gamma$ は零点をもたない(Γ関数の系「$\Gamma$ 関数の零点」)ので、$\Gamma(x+y)$ で割って主張を得る。$\square$

別の証明

$x,y$ が正の実数の場合は、$y$ を固定して $f(x):=\frac{\Gamma(x+y)}{\Gamma(y)}B(x,y)$ が $f(1)=1$、$f(x+1)=xf(x)$(prop-bt-basic の 4)、対数凸性(Hölder の不等式による)を満たすことを確かめ、Γ関数の定理「Bohr–Mollerup の定理」から $f=\Gamma$ を結論する証明もある(Rud76 Theorem 8.20, p. 193)。

右辺 $\Gamma(x)\Gamma(y)/\Gamma(x+y)$ は $x,y\notin\{0,-1,-2,\ldots\}$ で意味をもつので、$B$ の定義域を広げる式として使われる。ただし ex-bt-counterexamples のとおり、広げた範囲では値はもとの積分と一致しない。

系と応用

$\Gamma(\frac12)$ と三角関数の冪の積分
  1. $\Gamma(\frac12)=\sqrt\pi$ である。
  2. 実数 $p,q>-1$ に対し
    $$ \int_0^{\pi/2}\sin^p\theta\cos^q\theta\,d\theta=\frac12B\Bigl(\frac{p+1}2,\frac{q+1}2\Bigr)=\frac{\Gamma(\frac{p+1}2)\Gamma(\frac{q+1}2)}{2\,\Gamma(\frac{p+q}2+1)} $$
    であり、とくに $\int_0^{\pi/2}\sin^p\theta\,d\theta=\frac{\sqrt\pi}2\cdot\frac{\Gamma(\frac{p+1}2)}{\Gamma(\frac p2+1)}$ である。
積分表示に代入する

1:prop-bt-other-forms で $x=y=\frac12$ とすると $B(\frac12,\frac12)=2\int_0^{\pi/2}d\theta=\pi$ である。thm-bt-gamma と $\Gamma(1)=1$ から $\Gamma(\frac12)^2=\pi$ であり、$\Gamma(\frac12)=\int_0^\infty t^{-1/2}e^{-t}\,dt>0$ なので $\Gamma(\frac12)=\sqrt\pi$ である。
2:prop-bt-other-forms で $x=\frac{p+1}2$, $y=\frac{q+1}2$(実部は正)とし、thm-bt-gamma を使う。$q=0$ では 1 により $\Gamma(\frac12)=\sqrt\pi$ である。$\square$

この証明は Gauss 積分 $\int_{-\infty}^\infty e^{-u^2}du=\sqrt\pi$ を使っていない。逆に、置換 $t=u^2$ により $\Gamma(\frac12)=\int_{-\infty}^\infty e^{-u^2}du$ なので、cor-bt-half の 1 から Gauss積分の値が従う。$p$ が整数のときの 2 は Wallis 積分の値で、高校数学での扱いは Wallis積分とガンマ関数 にある。

Legendre の倍数公式

$2z\notin\{0,-1,-2,\ldots\}$ を満たす複素数 $z$ に対し
$$ \Gamma(z)\,\Gamma\Bigl(z+\frac12\Bigr)=2^{1-2z}\sqrt\pi\,\Gamma(2z) $$
が成り立つ(DLMF 5.5.5、Rud76 8.21 の (102))。

$B(z,z)$ を 2 通りに計算する

まず $\operatorname{Re}z>0$ とする。被積分関数 $(t(1-t))^{z-1}$ は $t\mapsto1-t$ で変わらないので $B(z,z)=2\int_0^{1/2}(t(1-t))^{z-1}\,dt$ である。$u=4t(1-t)$ は $(0,\frac12)$ から $(0,1)$ への狭義単調増加な $C^1$ 級の全単射で、$du=4(1-2t)\,dt$, $(1-2t)^2=1-u$ より $dt=\frac{du}{4\sqrt{1-u}}$ である。$t(1-t)=\frac u4$ なので
$$ B(z,z)=2\int_0^1\Bigl(\frac u4\Bigr)^{z-1}\frac{du}{4\sqrt{1-u}}=2^{1-2z}\int_0^1u^{z-1}(1-u)^{-1/2}\,du=2^{1-2z}B\Bigl(z,\frac12\Bigr) $$
である。両辺を thm-bt-gamma で書くと $\frac{\Gamma(z)^2}{\Gamma(2z)}=2^{1-2z}\frac{\Gamma(z)\Gamma(\frac12)}{\Gamma(z+\frac12)}$ であり、$\Gamma(z)\ne0$ で割り、cor-bt-half の $\Gamma(\frac12)=\sqrt\pi$ を使うと主張を得る。
一般の場合:$\Omega:=\mathbb{C}\setminus\{0,-\frac12,-1,-\frac32,\ldots\}$ は連結な開集合で、主張の両辺は $\Omega$ 上正則である(Γ関数の極は $0,-1,-2,\ldots$ だけ)。両辺は $\Omega$ の空でない開部分集合 $\{\operatorname{Re}z>0\}$ で一致するので、一致の定理により $\Omega$ 全体で一致する。$\square$

相反公式との組み合わせ

$0<\operatorname{Re}x<1$ に対し
$$ B(x,1-x)=\int_0^\infty\frac{w^{x-1}}{1+w}\,dw=\frac{\pi}{\sin\pi x} $$
である。とくに $\Gamma(\frac14)\Gamma(\frac34)=\sqrt2\,\pi$ であり、
$$ \int_0^1\frac{dt}{\sqrt{1-t^4}}=\frac14B\Bigl(\frac14,\frac12\Bigr)=\frac{\Gamma(\frac14)^2}{4\sqrt{2\pi}} $$
である。

主定理と相反公式から

prop-bt-other-forms と thm-bt-gamma により $B(x,1-x)=\Gamma(x)\Gamma(1-x)/\Gamma(1)$ であり、これは 無限積の定理「Γ関数の相反公式」により $\pi/\sin\pi x$ に等しい。$x=\frac14$ とすると $\Gamma(\frac14)\Gamma(\frac34)=\pi/\sin\frac\pi4=\sqrt2\,\pi$ である。最後の積分は置換 $u=t^4$($dt=\frac14u^{-3/4}du$)により $\frac14\int_0^1u^{-3/4}(1-u)^{-1/2}\,du=\frac14B(\frac14,\frac12)$ となり、thm-bt-gamma と $\Gamma(\frac34)=\sqrt2\pi/\Gamma(\frac14)$ により
$$ \frac14\cdot\frac{\Gamma(\frac14)\sqrt\pi}{\Gamma(\frac34)}=\frac{\Gamma(\frac14)^2\sqrt\pi}{4\sqrt2\,\pi}=\frac{\Gamma(\frac14)^2}{4\sqrt{2\pi}} $$
である。$\square$

最後の積分はレムニスケートの弧長の $\frac14$ で(レムニスケート)、第 1 種の楕円積分の値の一つである。

ベータ分布と $t$ 分布の正規化定数

$p,q>0$ に対し、$(0,1)$ 上の関数 $\frac{t^{p-1}(1-t)^{q-1}}{B(p,q)}$ は積分が $1$ の正の関数で、これを密度とする確率分布をベータ分布という。その平均は $\frac{B(p+1,q)}{B(p,q)}=\frac p{p+q}$ である(prop-bt-basic の 4)。
$\nu>0$ に対し、偶関数であることと置換 $x=\sqrt{\nu w}$($dx=\frac{\sqrt\nu}{2\sqrt w}dw$)により
$$ \int_{-\infty}^\infty\Bigl(1+\frac{x^2}\nu\Bigr)^{-\frac{\nu+1}2}dx=\sqrt\nu\int_0^\infty\frac{w^{-1/2}}{(1+w)^{\frac12+\frac\nu2}}\,dw=\sqrt\nu\,B\Bigl(\frac12,\frac\nu2\Bigr) $$
である(prop-bt-other-forms)。これは自由度 $\nu$ の $t$ 分布の密度の正規化定数である。

反例:条件を外すと

外す条件反例成り立たなくなること
定義の $\operatorname{Re}x>0$$x=0$, $y=1$$\int_0^1t^{-1}\,dt=+\infty$ で、積分が定まらない
thm-bt-gamma の $\operatorname{Re}x>0$$x=-\frac12$, $y=1$右辺は $-2$ と有限に定まるが、左辺の積分は $+\infty$ に発散する
cor-bt-integers の「$m,n$ は整数」$m=n=\frac12$階乗の式は整数でない $m,n$ では定義されず、thm-bt-gamma の Γ関数の式に替える必要がある($B(\frac12,\frac12)=\pi$)
反例の確かめ

1 行目は prop-bt-convergence の $a=0$ の場合である。2 行目:Γ関数の例「半整数での値」により $\Gamma(-\frac12)=-2\sqrt\pi$, $\Gamma(\frac12)=\sqrt\pi$ なので $\frac{\Gamma(-1/2)\Gamma(1)}{\Gamma(1/2)}=-2$ である。一方 $\int_0^1t^{-3/2}\,dt=+\infty$ である。被積分関数が正なのに右辺が負であることからも、右辺を積分として読めないことがわかる。3 行目は cor-bt-half による。

関連項目

参考文献

[1]
Walter Rudin, Principles of Mathematical Analysis, McGraw-Hill, 1976, Chapter 8, 8.19 Theorem(Bohr–Mollerup)・8.20 Theorem(ベータ関数と Γ関数), p. 193;8.21(三角関数の冪の積分 (98)、Γ(1/2)=√π (99)、倍数公式 (102)), p. 194

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する