リチャードソンの補外

リチャードソンの補外（リチャードソンのほがい）とは、外挿法の一種である。パラメータ テンプレート:Math を持つ量 テンプレート:Math について、テンプレート:Math における テンプレート:Math の極限値を近似的に求めるときに用いられる。

応用例として、台形公式を用いた数値積分にリチャードソンの補外を用いることで、テンプレート:仮リンクを導くことができる。また、CAEで計算格子を限りなく小さくしていく極限での解を予想することにも使われる。

手法

既知の2つのデータ テンプレート:Math と テンプレート:Mathから、テンプレート:Math の極限値 テンプレート:Math の近似値 ${\overline{f}}_1^{(1)}$ を求めるアルゴリズムは以下である。

${\overline{f}}_1^{(1)}=\frac{f(\lambda x)-{\lambda }^{p_1}f(x)}{1-{\lambda }^{p_1}}$

ただし テンプレート:Math は、テンプレート:Math を テンプレート:Math の多項式として漸近展開した以下の式に現れる指数である。

$f(x)=F+\sum _{i\ge 1}{a}_i{x}^{p_i}=F+a_1{x}^{p_1}+a_2{x}^{p_2}+\cdots ,0<p_1<p_2<\cdots$

また、追加データとして $f({\lambda }^{2}x),\dots ,f({\lambda }^{M}x)$ を求めることができたなら、

$\begin{array}{cc}{\overline{f}}_i^{(0)}& =f({\lambda }^{i}x),i=1,\dots ,M\\ {\overline{f}}_i^{(j)}& =\frac{{\overline{f}}_i^{(j-1)}-{\lambda }^{p_j}{\overline{f}}_{i-1}^{(j-1)}}{1-{\lambda }^{p_j}},j=1,\dots ,i\end{array}$

を順次求めていくことで、より精度の高い近似値 ${\overline{f}}_M^{(M)}$ を求めていくことができる。

精度

${\overline{f}}_1^{(1)}$ は テンプレート:Math に比例する誤差を持つ近似値である。すなわち

${\overline{f}}_1^{(1)}=F+O(x^{p_2})$

同様に、${\overline{f}}_M^{(M)}$の誤差評価は次式となる。

${\overline{f}}_M^{(M)}=F+O(x^{p_{M+1}})$