Sstrhsrtj

エイトケンのΔ2乗加速法（エイトケンのデルタじじょうかそくほう、Aitken's Δ² process）とは、少ない計算量で数列の収束を速めるためのアルゴリズムの一つである。

概要

数列 s_n がある極限値に収束するとき、以下の定義によって新たな数列 t_n を作ると、後者の収束の速さが大きく加速されることがある。

tn:=sn−(sn+1−sn)2sn+2−2sn+1+sn.{displaystyle t_{n}:=s_{n}-{frac {(s_{n+1}-s_{n})^{2}}{s_{n+2}-2s_{n+1}+s_{n}}}.}

この新たな数列 t_n によって、極限値の近似値の精度を上げる方法をエイトケンの Δ² 加速法と呼ぶ。

加速される理由

以下に見るように、エイトケンの Δ² 加速法は一種のはさみうち法である。今、s_n+1 は s_n によって決まり、数列はある極限値 α へ収束すると仮定する。

sn+1=g(sn)⟶α,(n⟶+∞),{displaystyle s_{n+1}=g(s_{n})longrightarrow alpha ,quad (nlongrightarrow +infty ),}

したがって、 α は g(x) の不動点 (fixed point) である。α を求めることは、方程式 x = g(x) の解を求めることと等価であり、図形的には、直線 y = x と曲線 y = g(x) の交点を求めることに等しい。ここで図を参照すると、2点 Pn:=(sn,g(sn))=(sn,sn+1),{displaystyle P_{n}:=(s_{n},g(s_{n}))=(s_{n},s_{n+1}),;} Pn+1:=(sn+1,g(sn+1))=(sn+1,sn+2){displaystyle P_{n+1}:=(s_{n+1},g(s_{n+1}))=(s_{n+1},s_{n+2})} を通る直線 L の方程式は、

y−sn+2=sn+1−sn+2sn−sn+1(x−sn+1),{displaystyle y-s_{n+2}={frac {s_{n+1}-s_{n+2}}{s_{n}-s_{n+1}}}(x-s_{n+1}),}

で与えられる。 L と直線 y = x との交点を Q = (x₀, x₀) とすると、L の方程式に y = x = x₀ を代入して、

x0=sn−(sn+1−sn)2sn+2−2sn+1+sn,{displaystyle ;x_{0}=s_{n}-{frac {(s_{n+1}-s_{n})^{2}}{s_{n+2}-2s_{n+1}+s_{n}}},;}

を得る。図では、交点 Q は P_n+2 よりも不動点 F := (α, α) に近づいている。これがエイトケンの Δ² 加速法の定義式の理由である。

注意点

収束を加速できるのは、出発値が収束値に十分近いときであり、離れているときは多くのステップを要する、または収束しないことも起こりうる。また、本来は収束しない数列に対してエイトケンの Δ² 加速法を適用した場合は、あたかも極限値が存在してそれに収束するように振る舞うことがある。エイトケンの Δ² 加速法により数列の極限値への収束を加速できるか否かは、元の数列の性質に依存する。対数収束のような収束が遅い数列に対しては殆ど効果が無いので別の加速法が必要である。

エイトケンの Δ² 加速法は、上の説明からわかるように方程式 x = g(x) の数値計算にも応用可能である。事実、エイトケンは代数方程式の近似計算にこの方法を適用した^[1][2]。

歴史

現時点に於いて判明しているところによると、エイトケンの Δ² 加速法は和算家の関孝和によって1681年頃に導出されたのが世界で最初である。関孝和は暦の作成で必要になった円周率の近似計算でこの加速法を用いて、小数点以下第16位までを正確に求めた^[2][3]。関孝和の業績は日本でも長らく忘却されていたが、フランス人Brezinski 1991に指摘されて^[2]再発見された。西洋に於いてエイトケンの Δ² 加速法が導出されたのは、それから約200年後の1876年であり、H. von. Nägelsbach 1876による。エイトケンの Δ² 加速法という名前で今日呼ばれるようになったのはさらに後の、アレクサンダー・エイトケン（英語版、ドイツ語版）の論文(Aitken 1926)にちなむ。

参考文献

脚注