楕円積分(elliptic integral)とは、重根をもたない 3 次または 4 次の多項式 $P$ と有理関数 $R$ による積分 $\int R(x,\sqrt{P(x)})\,dx$ の総称であり、一般に原始関数は初等関数で書けない。代表は Legendre の標準形で、第 1 種・第 2 種の完全楕円積分 $K(k)=\int_0^{\pi/2}(1-k^2\sin^2\theta)^{-1/2}d\theta$、$E(k)=\int_0^{\pi/2}(1-k^2\sin^2\theta)^{1/2}d\theta$ は冪級数展開、Legendre の関係式 $EK'+E'K-KK'=\pi/2$、算術幾何平均による表示 $K(k)=\pi/(2M(1,k'))$ をもつ。楕円の周、振り子の周期、レムニスケートの弧長がこれらで表され、第 1 種の積分の逆関数が楕円関数である。
楕円積分は、重根をもたない 3 次または 4 次の多項式 $P$ の平方根 $\sqrt{P(x)}$ を含む積分の総称である。名前は楕円の周の長さを表す積分に由来する。根号の中が 2 次式なら $\int\frac{dx}{\sqrt{1-x^2}}=\arcsin x$ のように置換で計算できるが、3 次・4 次になると一般に原始関数は初等関数で書けない(この記事では証明しない。初等関数で表せない原始関数の例「根号の中が 2 次式と 4 次式の場合」)。
そこで、代表的な積分を新しい関数として名前を付ける。Legendre の標準形
$$
K(k)=\int_0^{\pi/2}\frac{d\theta}{\sqrt{1-k^2\sin^2\theta}},\qquad E(k)=\int_0^{\pi/2}\sqrt{1-k^2\sin^2\theta}\,d\theta
$$
(第 1 種・第 2 種の完全楕円積分)がその中心である。楕円の周、振り子の周期、レムニスケートの弧長はこれらで表され、$K$ と $E$ には冪級数展開、Legendre の関係式、算術幾何平均による表示がある。第 1 種の不完全楕円積分を上端の関数とみたときの逆関数が楕円関数である。
$P$ を重根をもたない 3 次または 4 次の多項式、$R(x,y)$ を 2 変数の有理関数とする。$R(x,y)$ が $y$ の奇数次の冪を含むとき、積分
$$
\int R\bigl(x,\sqrt{P(x)}\bigr)\,dx
$$
を楕円積分(elliptic integral)という。
$y^2=P(x)$ を使うと $R(x,y)$ は $\rho(x)/y+\sigma(x)$($\rho,\sigma$ は有理関数)の形に直せるので、本質的な部分は $\int\rho(x)\,dx/\sqrt{P(x)}$ である(DLMF 19.2.2–19.2.3)。Legendre は、すべての楕円積分が次の 3 種類の積分と初等関数(代数関数・対数・三角関数)で表せることを示した。この記事ではこの還元を証明しない(DLMF §19.14(ii))。
$0\le k<1$ とし、$k':=\sqrt{1-k^2}$ とおく。$0\le\varphi\le\frac\pi2$ と実数 $n<1$ に対し
$$
F(\varphi,k):=\int_0^\varphi\frac{d\theta}{\sqrt{1-k^2\sin^2\theta}},\qquad
E(\varphi,k):=\int_0^\varphi\sqrt{1-k^2\sin^2\theta}\,d\theta,
$$
$$
\Pi(\varphi,n,k):=\int_0^\varphi\frac{d\theta}{(1-n\sin^2\theta)\sqrt{1-k^2\sin^2\theta}}
$$
をそれぞれ第 1 種・第 2 種・第 3 種の(不完全)楕円積分という。$k$ を母数、$k'$ を補母数という。$\varphi=\frac\pi2$ のものを完全楕円積分といい、$K(k):=F(\frac\pi2,k)$, $E(k):=E(\frac\pi2,k)$ と書く。さらに $0< k<1$ に対し $K'(k):=K(k')$, $E'(k):=E(k')$ とおく($k=0$ では $k'=1$ で、$K(1)$ は定まらない)。
置換 $t=\sin\theta$ により $F(\varphi,k)=\int_0^{\sin\varphi}\frac{dt}{\sqrt{(1-t^2)(1-k^2t^2)}}$ である(DLMF 19.2.4)。$0< k<1$ なら根号の中は重根のない 4 次式で、def-ell-general の形になっている。文献によっては $k^2$ を $m$ と書いて $K(m)$ とし、第 3 種の $n\sin^2\theta$ の符号を逆にとる。$E(k)$ は $\theta\mapsto\frac\pi2-\theta$ により $\int_0^{\pi/2}\sqrt{1-k^2\cos^2\theta}\,d\theta$ とも書ける。
$\Delta(k,\theta):=\sqrt{1-k^2\sin^2\theta}$ とおく。1:$0<\theta\le\frac\pi2$ では $\Delta$ は $k$ について狭義単調減少で、$0<\Delta\le1$ である。よって $1/\Delta$ の積分 $K$ は増加し、$\Delta$ の積分 $E$ は減少する。$k=0$ では $\Delta=1$ である。
2:$k_0<1$ を固定すると、$\Delta$ と $1/\Delta$ は $[0,k_0]\times[0,\frac\pi2]$ 上で連続、したがって一様連続なので、$K$ と $E$ は $[0,k_0]$ 上連続である。$\Delta$ は $[0,1]\times[0,\frac\pi2]$ 上でも連続なので $E$ は $[0,1]$ 上連続で、$E(1)=\int_0^{\pi/2}\cos\theta\,d\theta=1$ である。
3:$0<\delta<\frac\pi2$ とする。$[0,\frac\pi2-\delta]$ では $\cos\theta>0$ なので、$k\to1$ で $1/\Delta\to1/\cos\theta$ が一様に成り立ち、$\liminf_{k\to1}K(k)\ge\int_0^{\pi/2-\delta}\frac{d\theta}{\cos\theta}=\log\frac{1+\cos\delta}{\sin\delta}$ である。$\delta\to0$ で右辺は $+\infty$ に近づく。$\square$
$a\ge b>0$ とし、楕円 $\frac{x^2}{a^2}+\frac{y^2}{b^2}=1$ の離心率を $e:=\sqrt{1-b^2/a^2}$ とする。媒介変数表示 $(x,y)=(a\sin\theta,b\cos\theta)$ では速さが $\sqrt{a^2\cos^2\theta+b^2\sin^2\theta}=a\sqrt{1-e^2\sin^2\theta}$ なので、周の長さは
$$
L=4a\int_0^{\pi/2}\sqrt{1-e^2\sin^2\theta}\,d\theta=4a\,E(e)
$$
である。$e=0$(円)では $L=2\pi a$ に戻る。高校数学の範囲での評価と数値計算は 楕円の周の長さ にある。
長さ $\ell$ の単振り子が角度 $\theta_0$($0<\theta_0<\pi$)から静かに放されると、運動方程式 $\ddot\theta=-\frac g\ell\sin\theta$ とエネルギー保存から $\dot\theta^2=\frac{2g}\ell(\cos\theta-\cos\theta_0)$ である。運動が周期的で、$\theta_0$ から $0$ までの時間が周期の $\frac14$ であることを認めると、周期は
$$
T=4\int_0^{\theta_0}\frac{d\theta}{\sqrt{\frac{2g}\ell(\cos\theta-\cos\theta_0)}}=4\sqrt{\frac\ell g}\,K\Bigl(\sin\frac{\theta_0}2\Bigr)
$$
となる。実際、$k:=\sin\frac{\theta_0}2$ とおき $\sin\frac\theta2=k\sin\varphi$ と置換すると、$\cos\theta-\cos\theta_0=2(k^2-\sin^2\frac\theta2)=2k^2\cos^2\varphi$、$d\theta=\frac{2k\cos\varphi}{\sqrt{1-k^2\sin^2\varphi}}\,d\varphi$ となり、被積分関数は $\sqrt{\ell/g}\,(1-k^2\sin^2\varphi)^{-1/2}$ に変わる。振れ幅が小さい($k\to0$)とき $T\to2\pi\sqrt{\ell/g}$ であり、振れ幅が大きいほど周期は長い(prop-ell-basic の 1)。
極座標で $r^2=\cos2\theta$ と表される曲線(レムニスケート)の全長の $\frac14$ は $\int_0^1\frac{dt}{\sqrt{1-t^4}}$ である(この積分の 2 倍 $\varpi$ をレムニスケート定数という)。置換 $t=\cos\varphi$ により $1-t^4=\sin^2\varphi\,(1+\cos^2\varphi)=\sin^2\varphi\,(2-\sin^2\varphi)$ なので
$$
\int_0^1\frac{dt}{\sqrt{1-t^4}}=\int_0^{\pi/2}\frac{d\varphi}{\sqrt{2-\sin^2\varphi}}=\frac1{\sqrt2}K\Bigl(\frac1{\sqrt2}\Bigr)
$$
となる。一方 ベータ関数の系「相反公式との組み合わせ」によりこの積分は $\frac{\Gamma(1/4)^2}{4\sqrt{2\pi}}=1.3110\ldots$ に等しい(DLMF 19.20.2)ので、$K(1/\sqrt2)=\frac{\Gamma(1/4)^2}{4\sqrt\pi}=1.8540\ldots$ である。
以下、$c_n:=\dfrac{(2n)!}{4^n(n!)^2}=\dfrac{1\cdot3\cdots(2n-1)}{2\cdot4\cdots(2n)}$($c_0=1$)とおく。
$0\le k<1$ に対し
$$
K(k)=\frac\pi2\sum_{n=0}^\infty c_n^2k^{2n}=\frac\pi2\Bigl(1+\frac14k^2+\frac9{64}k^4+\cdots\Bigr),
$$
$$
E(k)=\frac\pi2\sum_{n=0}^\infty\frac{c_n^2}{1-2n}k^{2n}=\frac\pi2\Bigl(1-\frac14k^2-\frac3{64}k^4-\cdots\Bigr)
$$
が成り立つ(DLMF 19.5.1–19.5.2)。
一般二項定理(二項定理と組合せの恒等式の定理「一般二項定理」)を $\alpha=\mp\frac12$ と $-x$ に使う。$\binom{-1/2}n(-1)^n=\frac{\frac12\cdot\frac32\cdots\frac{2n-1}2}{n!}=c_n$ であり、$n\ge1$ では $\binom{1/2}n(-1)^n=-\frac{\frac12\cdot\frac12\cdot\frac32\cdots\frac{2n-3}2}{n!}=-\frac{c_n}{2n-1}$ である。よって $0\le x<1$ で
$$
\frac1{\sqrt{1-x}}=\sum_{n=0}^\infty c_nx^n,\qquad\sqrt{1-x}=\sum_{n=0}^\infty\frac{c_n}{1-2n}x^n
$$
である。$x=k^2\sin^2\theta$ とおくと、$\theta\in[0,\frac\pi2]$ で各項の絶対値は $c_nk^{2n}$ 以下であり、$\sum_nc_nk^{2n}=(1-k^2)^{-1/2}<\infty$ である。一様収束の定理「Weierstrassの優級数判定法」により 2 つの級数は $[0,\frac\pi2]$ 上一様収束し、定理「一様収束と積分の交換」により項別に積分できる。Wallis積分(高校数学)の命題「Wallis 積分の値」により $\int_0^{\pi/2}\sin^{2n}\theta\,d\theta=\frac\pi2c_n$ なので、主張を得る。$\square$
$K$ と $E$ は母数について簡単な微分方程式を満たし、そこから $k$ と $k'$ を入れ替えた 4 つの積分のあいだの恒等式が得られる。以下 $\Delta:=\sqrt{1-k^2\sin^2\theta}$ と書く。
$0\le k<1$ に対し $\displaystyle\int_0^{\pi/2}\frac{d\theta}{\Delta^3}=\frac{E(k)}{k'^2}$ である。
$k=0$ では両辺とも $\frac\pi2$ である。$0< k<1$ とし、$f(\theta):=\sin\theta\cos\theta/\Delta$ とおく。$s:=\sin^2\theta$ と書くと $\Delta^2=1-k^2s$ であり、
$$
f'(\theta)=\frac{(\cos^2\theta-\sin^2\theta)\Delta^2+k^2\sin^2\theta\cos^2\theta}{\Delta^3}=\frac{(1-2s)(1-k^2s)+k^2s(1-s)}{\Delta^3}=\frac{1-2s+k^2s^2}{\Delta^3}
$$
である。$s=(1-\Delta^2)/k^2$ を代入すると分子は $\frac{\Delta^4-k'^2}{k^2}$ となるので、$f'=\frac1{k^2}\bigl(\Delta-\frac{k'^2}{\Delta^3}\bigr)$ である。$f(0)=f(\frac\pi2)=0$ なので、$0$ から $\frac\pi2$ まで積分して $E(k)=k'^2\int_0^{\pi/2}\Delta^{-3}\,d\theta$ を得る。$\square$
$0< k<1$ で
$$
\frac{dE}{dk}=\frac{E-K}{k},\qquad\frac{dK}{dk}=\frac{E-k'^2K}{k\,k'^2}
$$
が成り立つ(DLMF 19.4.1–19.4.2)。
$k_0<1$ を固定すると、$\Delta$ と $1/\Delta$ の $k$ についての偏導関数は $[0,k_0]\times[0,\frac\pi2]$ 上連続なので、積分記号下で微分できる(Rud76 Theorem 9.42, p. 236)。$k^2\sin^2\theta=1-\Delta^2$ を使うと
$$
\frac{dE}{dk}=\int_0^{\pi/2}\frac{-k\sin^2\theta}{\Delta}\,d\theta=\frac1k\int_0^{\pi/2}\frac{\Delta^2-1}{\Delta}\,d\theta=\frac{E-K}k
$$
であり、
$$
\frac{dK}{dk}=\int_0^{\pi/2}\frac{k\sin^2\theta}{\Delta^3}\,d\theta=\frac1k\int_0^{\pi/2}\Bigl(\frac1{\Delta^3}-\frac1\Delta\Bigr)d\theta=\frac1k\Bigl(\frac{E}{k'^2}-K\Bigr)
$$
である(最後は lem-ell-delta-cube)。$\square$
$0<\kappa<1$ に対し $K(\kappa)\le\dfrac2{\sqrt\kappa}\log\dfrac{1+\sqrt\kappa}{\sqrt{1-\kappa}}$ である。
置換 $t=\sin\theta$ により $K(\kappa)=\int_0^1\frac{dt}{\sqrt{(1-t^2)(1-\kappa^2t^2)}}$ である。$0\le t<1$ で $1-t^2\ge1-t$、$1-\kappa^2t^2\ge1-\kappa t$ なので、$s=1-t$ と置換して
$$
K(\kappa)\le\int_0^1\frac{dt}{\sqrt{(1-t)(1-\kappa t)}}=\int_0^1\frac{ds}{\sqrt{s\,(1-\kappa+\kappa s)}}
$$
である。$G(s):=\frac2{\sqrt\kappa}\log\bigl(\sqrt{\kappa s}+\sqrt{1-\kappa+\kappa s}\bigr)$ を微分すると $G'(s)=\frac1{\sqrt s\sqrt{1-\kappa+\kappa s}}$ となるので、右辺は $G(1)-G(0)=\frac2{\sqrt\kappa}\log\frac{\sqrt\kappa+1}{\sqrt{1-\kappa}}$ である。$\square$
$0< k<1$ に対し
$$
E(k)K'(k)+E'(k)K(k)-K(k)K'(k)=\frac\pi2
$$
が成り立つ(DLMF 19.7.1)。
左辺を $L(k)$ とおく。$k'=\sqrt{1-k^2}$ は $\frac{dk'}{dk}=-\frac k{k'}$ を満たすので、prop-ell-derivatives を母数 $k'$ で使い、$k'^2$ と $k^2$ の役割を入れ替えると
$$
\frac{dE'}{dk}=-\frac{k(E'-K')}{k'^2},\qquad\frac{dK'}{dk}=-\frac{E'-k^2K'}{k\,k'^2}
$$
である。$\frac{dE}{dk}-\frac{dK}{dk}=\frac{k'^2(E-K)-E+k'^2K}{kk'^2}=-\frac{kE}{k'^2}$ に注意して $L=(E-K)K'+E'K$ を微分し、$kk'^2$ を掛けると
$$
kk'^2\frac{dL}{dk}=-k^2EK'-(E-K)(E'-k^2K')-k^2(E'-K')K+E'(E-k'^2K)
$$
である。$k'^2=1-k^2$ を使って右辺を展開すると
$$
-k^2EK'+\bigl(-EE'+k^2EK'+KE'-k^2KK'\bigr)+\bigl(-k^2E'K+k^2KK'\bigr)+\bigl(EE'-E'K+k^2E'K\bigr)=0
$$
となる。よって $L$ は $(0,1)$ 上で定数である。
$k\to0$ の極限を求める。$K-E=\int_0^{\pi/2}\frac{1-\Delta^2}{\Delta}\,d\theta=\int_0^{\pi/2}\frac{k^2\sin^2\theta}{\Delta}\,d\theta$ なので $0\le K-E\le\frac{\pi k^2}{2k'}$ である。lem-ell-log-bound を $\kappa=k'$ に使い、$1-k'=\frac{k^2}{1+k'}\ge\frac{k^2}2$ を用いると $K'(k)\le\frac2{\sqrt{k'}}\log\frac{2\sqrt2}{k}$ である。よって $k\to0$ で $(E-K)K'\to0$ である。また $E'(k)=E(k')\to E(1)=1$、$K(k)\to K(0)=\frac\pi2$ である(prop-ell-basic)。したがって $L=(E-K)K'+E'K\to\frac\pi2$ であり、定数 $L$ は $\frac\pi2$ に等しい。$\square$
$k=\frac1{\sqrt2}$ では $k'=k$ なので、Legendre の関係式は $2EK-K^2=\frac\pi2$ となる。ex-ell-lemniscate の $K(\frac1{\sqrt2})=\frac{\Gamma(1/4)^2}{4\sqrt\pi}$ を代入すると
$$
E\Bigl(\frac1{\sqrt2}\Bigr)=\frac K2+\frac\pi{4K}=\frac{\Gamma(1/4)^2}{8\sqrt\pi}+\frac{\pi^{3/2}}{\Gamma(1/4)^2}=1.3506\ldots
$$
となり、第 2 種の完全楕円積分の値も Γ関数で求まる。
$0\le k<1$ に対し $K(k)=\dfrac{\pi}{2M(1,k')}$ である。ここで $M(a,b)$ は $a,b>0$ の算術幾何平均である(DLMF 19.8.5)。
$\cos^2\theta+k'^2\sin^2\theta=1-k^2\sin^2\theta$ なので、算術幾何平均と積分の定理「算術幾何平均の積分表示」を $a=1$, $b=k'$ に使えばよい。$\square$
算術幾何平均の反復は 2 次で収束するので、cor-ell-agm は $K(k)$ の数値計算の標準的な方法である。たとえば $k=\frac1{\sqrt2}$ では $M(1,\frac1{\sqrt2})=0.8472\ldots$ から $K=1.8540\ldots$ を得て、ex-ell-lemniscate と一致する。
第 1 種の積分 $u=F(\varphi,k)$ を $\varphi$ について解いた関数を $\varphi=\operatorname{am}u$ と書き、$\operatorname{sn}u:=\sin(\operatorname{am}u)$ とおく。すなわち $\operatorname{sn}$ は $u=\int_0^{x}\frac{dt}{\sqrt{(1-t^2)(1-k^2t^2)}}$ の逆関数 $x=\operatorname{sn}u$ である。$0< k<1$ のとき、これを複素数の $u$ へ延ばすと、$4K$ と $2iK'$ を周期にもつ有理型関数、すなわち楕円関数になる(この記事では証明しない。Ahl66 Chapter 6 §2.3、DLMF22 Table 22.4.1)。$k\to0$ で $\operatorname{sn}u\to\sin u$ となるので、楕円関数は三角関数を 2 重周期に広げたものとみなせる。Weierstrass の $\wp$ 関数も、楕円積分 $\int_{\wp(z)}^\infty\frac{dx}{\sqrt{4x^3-g_2x-g_3}}$ の逆関数として現れる(Weierstrassの楕円関数)。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| $P$ は重根をもたない | $P(t)=(1-t^2)^2$($F$ で $k=1$ としたもの) | $\int_0^x\frac{dt}{1-t^2}=\operatorname{artanh}x$ と初等関数で書ける |
| $P$ の次数は 3 以上 | $P(t)=1-t^2$ | $\int_0^x\frac{dt}{\sqrt{1-t^2}}=\arcsin x$ と初等関数で書ける |
| 完全楕円積分の $k<1$ | $k=1$ | $K(1)=\int_0^{\pi/2}\frac{d\theta}{\cos\theta}=+\infty$ で、$K$ は定まらない($E(1)=1$ は有限) |
| 冪級数展開の $k<1$ | $k=1$ | $K$ の級数 $\sum_nc_n^2$ は発散する |
1・2 行目は $-1< x<1$ での計算である。1 行目では根号がはずれて被積分関数が有理関数になり、2 行目では $t=\sin\theta$ で被積分関数が $1$ になる。3 行目は prf-ell-basic の 3 と同じ計算による。4 行目:$j\ge2$ で $(2j-1)^2\ge(2j-2)\cdot2j$ なので $\bigl(\frac{2j-1}{2j}\bigr)^2\ge\frac{j-1}j$ であり、$c_n^2=\frac14\prod_{j=2}^n\bigl(\frac{2j-1}{2j}\bigr)^2\ge\frac14\prod_{j=2}^n\frac{j-1}j=\frac1{4n}$ である。調和級数は発散するので $\sum_nc_n^2=+\infty$ である。一方、$(2j-1)(2j+1)\le(2j)^2$ から同様に $c_n^2\le\prod_{j=1}^n\frac{2j-1}{2j+1}=\frac1{2n+1}$ なので、$E$ の級数は $k=1$ でも $\sum_{n\ge1}\frac1{(2n-1)(2n+1)}$ と比べて収束する。$K(1)$ が無限大で $E(1)$ が有限であることと合っている。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する