ガウス=ルジャンドルのアルゴリズム

出典: フリー百科事典『ウィキペディア(Wikipedia)』
2014年8月20日 (水) 15:24時点における36.245.196.153 (トーク)による版 (アルゴリズム)
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
移動先: 案内検索

ガウス=ルジャンドルのアルゴリズム円周率を計算する際に用いられる数学の反復計算アルゴリズムである。円周率を計算する物の中では非常に収束が速く、2009年にこの式を用いて2,576,980,370,000桁(約2兆6000億桁)の計算がされた。

このアルゴリズムはカール・フリードリヒ・ガウスアドリアン=マリ・ルジャンドルがそれぞれ別個に研究したものである。これは二つの数値の算術幾何平均を求めるために、それぞれの数値を算術平均(相加平均)と幾何平均(相乗平均)で置き換えていくものである。

アルゴリズム

これによる円周率の計算方法は以下の通りである。

初期値の設定

<math>a_0 = 1 \qquad b_0 = \frac {1}{\sqrt {2}} \qquad t_0 = \frac {1}{4} \qquad p_0 = 1\!</math>

反復式

a,bが希望する精度(桁数)になるまで以下の計算を繰り返す。 小数第n位まで求めるとき log2 n回程度の反復でよい。

<math> \begin{align} a_{n+1} & = \frac {a_n + b_n}{2} \\
                     b_{n+1} & = \sqrt {a_n b_n} \\
                     t_{n+1} & = t_n - p_n(a_n - a_{n+1})^2 \\
                     p_{n+1} & = 2p_n
       \end{align}

</math>

πの算出

円周率 π は、a、b、tを用いて以下のように近似される。

<math> \pi \approx \frac{(a+b)^2}{4t} </math>

最初の3回の反復で得られる数値(最後の桁は正確ではない)

<math>3.140\dots\!</math>
<math>3.14159264\dots\!</math>
<math>3.1415926535897932382\dots\!</math>

この反復プロセスは自然収束し、反復1回について前に正常だった桁の2倍の桁の数値までが収束する。ガウス自身、この式で4回まで反復を行い12桁まで正しいことを確認したことが知られている。

テンプレート:Math-stub