算術幾何平均とレムニスケートの長さ

112\sqrt{2} の算術幾何平均 M(1,2)M(1,\sqrt{2}) は,単位円周の長さ l=2πl=2\pi とレムニスケートの長さ LL の比に等しい: M(1,2)=lL M(1,\sqrt{2})=\dfrac{l}{L}

「算術幾何平均」や「レムニスケートの長さ」については後述します。非常に美しい定理です。

レムニスケートの長さ LL

r2=cos2θr^2=\cos 2\theta という式で表される曲線をレムニスケートと言います(より一般には r2=2a2cos2θr^2=2a^2\cos 2\theta )。

lemniscate

この曲線の長さを LL とすると, L4=01dx1x4\dfrac{L}{4}=\displaystyle\int_0^1\dfrac{dx}{\sqrt{1-x^4}} が成立します。

証明

弧長積分の公式の極座標版を使うと,

L4=011+r2(dθdr)2dr \dfrac{L}{4}=\int_0^1\sqrt{1+r^2\left(\frac{d\theta}{dr}\right)^2}dr となる。

r2=cos2θr^2=\cos 2\theta の両辺を rr で微分すると, 2r=2sin2θdθdr2r=-2\sin 2\theta\dfrac{d\theta}{dr} である。

よって, L4=011+r2r2sin22θdr=011+r41r4dr=01dr1r4\begin{aligned} \dfrac{L}{4}&= \int_0^1\sqrt{1+r^2\dfrac{r^2}{\sin^2 2\theta}}dr\\ &= \int_0^1\sqrt{1+\dfrac{r^4}{1-r^4}}dr\\ &= \int_0^1\dfrac{dr}{\sqrt{1-r^4}} \end{aligned}

算術幾何平均

算術幾何平均 M(a,b)M(a,b) とは,2つの正の実数 a,ba,b から定まる1つの数です。

具体的には,正の実数 a,ba,b が与えられたとき,a0=a,b0=ba_0=a,b_0=b を初期値として,以下の漸化式で数列 an,bna_n,b_n を定めます: {an+1=an+bn2bn+1=anbn \begin{cases} a_{n+1}=\dfrac{a_n+b_n}{2}\\ b_{n+1}=\sqrt{a_nb_n} \end{cases}

それぞれ相加平均(算術平均)と相乗平均(幾何平均)です。

このとき,2つの数列 ana_nbnb_nnn\to\infty で同じ値に収束します(*)。この値を aabb算術幾何平均と呼びます。この記事では aabb の算術幾何平均を M(a,b)M(a,b) と書くことにします。

(*)の証明の概略

bab\leq a の場合を考える(aba\leq b の場合も同様)。

相加相乗平均の不等式と帰納法を使うと,

bb1b2a2a1a b\leq b_1\leq b_2\leq \cdots \leq a_2\leq a_1\leq a

が分かる。

単調で有界な数列は収束するので,2つの数列は収束する。その収束先を α,β\alpha,\beta とすると,漸化式より α=α+β2\alpha=\dfrac{\alpha+\beta}{2} となり,α=β\alpha=\beta が分かる。

ちなみに,同様に算術調和平均や調和幾何平均も考えることができます(→調和平均について)。

算術幾何平均の積分表示

算術幾何平均は積分によって記述されます。

補題

M(a,b)=π2÷0π2dθa2cos2θ+b2sin2θ M(a,b) = \dfrac{\pi}{2} \div \int_0^{\frac{\pi}{2}} \dfrac{d\theta}{\sqrt{a^2\cos^2 \theta + b^2 \sin^2 \theta}}

これはガウスによって発見されました。

証明

F(a,b)=2π0π2dθa2cos2θ+b2sin2θ F(a,b) = \dfrac{2}{\pi} \int_0^{\frac{\pi}{2}} \dfrac{d\theta}{\sqrt{a^2\cos^2 \theta + b^2 \sin^2 \theta}} とおく。これが M(a,b)M(a,b) の逆数であることを示す。

x=btanθx = b \tan \theta と置換すると,dx=bcos2θdθdx = \dfrac{b}{\cos^2 \theta}d\theta であるため, 2π0π2dθa2cos2θ+b2sin2θ=2π0dx(a2+x2)(b2+x2)=1πdx(a2+x2)(b2+x2)\begin{aligned} &\dfrac{2}{\pi}\int_0^{\frac{\pi}{2}} \dfrac{d\theta}{\sqrt{a^2\cos^2 \theta + b^2 \sin^2 \theta}}\\ &= \dfrac{2}{\pi}\int_0^{\infty} \dfrac{dx}{\sqrt{(a^2+x^2)(b^2+x^2)}}\\ &= \dfrac{1}{\pi}\int_{-\infty}^{\infty} \dfrac{dx}{\sqrt{(a^2+x^2)(b^2+x^2)}} \end{aligned} となる。さらに y=12(xabx)y = \dfrac{1}{2} \left( x-\dfrac{ab}{x} \right) とおくと, 1πdx(a2+x2)(b2+x2)=1πdy((a+b2)2+y2)(ab+y2)\begin{aligned} &\dfrac{1}{\pi}\int_{-\infty}^{\infty} \dfrac{dx}{\sqrt{(a^2+x^2)(b^2+x^2)}} \\ &= \dfrac{1}{\pi} \int_{-\infty}^{\infty} \dfrac{dy}{\sqrt{\left( \left( \dfrac{a+b}{2} \right)^2 + y^2 \right) (ab+y^2)}} \end{aligned} と変形される。

これより F(a,b)=F(a+b2,ab)F(a,b) = F\left( \dfrac{a+b}{2} , \sqrt{ab} \right) となる。

ゆえに a0=aa_0 =ab0=bb_0 = b として {an+1=an+bn2bn+1=anbn \begin{cases} a_{n+1} = \dfrac{a_n + b_n}{2}\\ b_{n+1} = \sqrt{a_n b_n} \end{cases} と定まる数列について F(a,b)=F(a0,b0)=F(a1,b1)==F(an,bn)==F(M(a,b),M(a,b))\begin{aligned} F(a,b) &= F(a_0 , b_0)\\ &= F(a_1 , b_1)\\ &= \cdots \\ &= F(a_n , b_n ) \\ &= \cdots\\ &= F(M(a,b) , M(a,b)) \end{aligned} である。こうして F(a,b)=F(M(a,b),M(a,b))=2π0π2dθα2cos2θ+α2sin2θ=1M(a,b)\begin{aligned} F(a,b) &= F(M(a,b) , M(a,b)) \\ &= \dfrac{2}{\pi} \int_0^{\frac{\pi}{2}} \dfrac{d\theta}{\sqrt{\alpha^2\cos^2 \theta + \alpha^2 \sin^2 \theta}} \\ &= \dfrac{1}{M(a,b)} \end{aligned} となる。

冒頭の定理の証明

証明の概略

補題に a=1,b=2a=1,b=\sqrt{2} を代入したとき,定積分は,

0π2dθ1+cos2θ \int_0^{\frac{\pi}{2}}\dfrac{d\theta}{\sqrt{1+\cos^2\theta}}

となるが,cosθ=x\cos\theta=x とおくと,dxdθ=sinθ\dfrac{dx}{d\theta}=-\sin\theta なので上式は,

1011+x2dx1x2=01dx1x4\begin{aligned} &\int_1^{0}\dfrac{1}{\sqrt{1+x^2}}\cdot\dfrac{dx}{-\sqrt{1-x^2}}\\ &= \int_0^1\dfrac{dx}{\sqrt{1-x^4}} \end{aligned}

となる。

これはさきほど証明したように,L4\dfrac{L}{4} なので,

M(1,2)=π2÷L4=lL\begin{aligned} M(1,\sqrt{2})&=\dfrac{\pi}{2}\div\dfrac{L}{4}\\ &=\dfrac{l}{L} \end{aligned}

@nagomi_osaka さんのツイートに inspire されて書いた記事です。楽しかったです!