線形補間の計算式と近似誤差

更新 2021/03/06

関数 $y=f(x)$ が2点 $(x_1,y_1)$ ， $(x_2,y_2)$ を通ることがわかっているとする。このとき，関数 $f(x)$ を，2点を通る線分：

$y=y_1+\dfrac{y_2-y_1}{x_2-x_1}(x-x_1)$

で近似する手法を線形補間と言う。

線形補間の例

例えば，(1,4)(1,4)(1,4)
 と
(3,8)(3,8)(3,8)
 を通る関数の，x=2x=2x=2
 での値を線形補間すると，
y=4+8−43−1(2−1)=6y=4+\dfrac{8-4}{3-1}(2-1)=6y=4+3−18−4​(2−1)=6
になります。
なお，補間とは与えられた点の間を別の関数で近似することです。2点の外側の関数値を同様に近似する場合は，補外や外挿などと言います。

線形補間の拡張

折れ線グラフは，線形補間をつなげたもの（区分線形補間）です。
高次元バージョンも考えられます。例えば，2変数関数 $f(x,y)$ について，通る3点が与えられた場合，その3点を通る平面の方程式を使うことで線形補間できます。平面の方程式とその３通りの求め方
１次関数ではなく，より一般的に多項式を使った補間であるラグランジュ補間という手法があります。

線形補間の誤差

以下，x1<x2x_1 < x_2x1​<x2​
 とします。
y=f(x)y=f(x)y=f(x)
 は
C2C^2C2
 級（二階微分可能で，二階導関数が連続）とする。このとき，線形補間の誤差の絶対値は，
(x2−x1)2M8\dfrac{(x_2-x_1)^2M}{8}8(x2​−x1​)2M​
以下である。ただし，M=max⁡x1≤x≤x2∣f′′(x)∣M=\displaystyle\max_{x_1\leq x\leq x_2}|f''(x)|M=x1​≤x≤x2​max​∣f′′(x)∣
ここで，線形補間の誤差
E(x)E(x)E(x)
 とは，関数
f(x)f(x)f(x)
 と近似値の差です。つまり，
f(x)−{y1+y2−y1x2−x1(x−x1)}f(x)-\left\{y_1+\dfrac{y_2-y_1}{x_2-x_1}(x-x_1)\right\}f(x)−{y1​+x2​−x1​y2​−y1​​(x−x1​)}
です。この誤差（の絶対値）が上から定数で抑えられるという素敵な定理です。
誤差の上界
(x2−x1)2M8\dfrac{(x_2-x_1)^2M}{8}8(x2​−x1​)2M​
 は
x1x_1x1​ と x2x_2x2​ が近いとき，つまり補間の範囲が狭いときに小さくなります。
二階微分が全体的に小さいとき，つまり「そんなに曲がっていない」ときに小さくなります。

誤差の上界の証明

証明

$x_1\leq x\leq x_2$ において近似誤差 $E(x)$ を最大にするような $x$ の1つを $x^*$ とする（最大値の原理より，そのような $x^*$ は存在する）。

このとき，テイラーの定理（平均値の定理の一般形）より

$E(x_1)\\ =E(x^*)+E'(x^*)(x_1-x^*)+\dfrac{E''(c)}{2}(x_1-x^*)^2$

となる $c$ が $x_1$ と $x^*$ の間に存在する。

ここで，線分の端点での誤差は $0$ なので， $E(x_1)=0$ である。また， $E(x^*)$ は最大値なので， $E'(x^*)=0$ である：

$0=E(x^*)+\dfrac{E''(c)}{2}(x^*-x_1)^2$

ここで， $x^*$ が $x_2$ よりも $x_1$ に近い場合，

つまり $x^*-x_1\leq\dfrac{x_2-x_1}{2}$ の場合は，上式より

$|E(x^*)|\leq \dfrac{(x_2-x_1)^2M}{8}$

が導かれる。一方， $x^*$ が $x_1$ よりも $x_2$ に近い場合は，テイラーの定理の式で $x_1$ のかわりに $x_2$ を使えば同じ上界が導出できる。

より一般的な定理

より一般に，ラグランジュ補間による誤差は，以下の式で評価できます：
E(x)=1(n+1)!f(n+1)(cx)∏i=0n(x−xi)E(x)=\dfrac{1}{(n+1)!}f^{(n+1)}(c_x)\displaystyle\prod_{i=0}^n(x-x_i)E(x)=(n+1)!1​f(n+1)(cx​)i=0∏n​(x−xi​)
ただし，xix_ixi​
 たちは与えられた点の
xxx
 座標で，cxc_xcx​
 は「xix_ixi​
 たちの最小値以上」で「xix_ixi​
 たちの最大値以下」の実数です。
n=1n=1n=1
 とすると，さきほどの定理（よりも少し強い形）になります。
参考文献：Errors in Polynomial Interpolation

この記事の監修者

マスオ

東京大学大学院情報理工学系研究科修了／2014年にWebサイト『高校数学の美しい物語』を立ち上げ／著書累計 50,000部突破／「わかりやすいこと」と「ごまかさないこと」の両立を意識している。 →著者情報・書籍一覧を見る