内挿のソースを表示

'''内挿'''（ないそう、{{Lang-en-short|interpolation}}）や'''補間'''（ほかん）とは、ある既知の数値データ列を基にして、そのデータ列の各区間の範囲内を埋める数値を求めること、またはそのような関数を与えること。またその手法を'''内挿法'''（{{lang-en-short|interpolation method}}）や'''補間法'''という。対義語は[[外挿]]や補外。
==概要==
内挿するためには、各区間の範囲内で成り立つと期待される関数と境界での振舞い（境界条件）を決めることが必要である。

最も一般的で容易に適用できるものは、一次関数（直線）による内挿（直線内挿）である。ゼロ次関数（[[ヘヴィサイドの階段関数|ステップ関数]]）によってデータ列を埋めること（0次補間）を内挿と呼ぶことはあまりないが、内挿の一種である。

内挿と[[外挿]]（補外）とのアルゴリズムの類似性から、それぞれ内挿補間、外挿補間と誤って呼称されることがある。本来、補間と内挿は同義であり、内挿補間と重ねて呼ぶ必要はない。

== 内挿法の選択 ==
{{節スタブ}}
内挿のもたらす結果は、[[平滑化]]や[[最小自乗法|最小自乗近似]]と似ているが、これらは全く違ったものである。内挿は、ある区間の間に成り立つ関数モデルや境界条件を仮定し、その関数のパラメータのうちのいくつか(または全て)を決定する。このため、入力数値データ列には誤差が含まれないか、無視できると仮定している。一方、平滑化や最小自乗近似は誤差が含まれる数値データ列の関係をもっともらしく推定する数列や関数モデルを与える。

物理現象を測定したデータを入力とする内挿では、その物理現象に適用できるもっともらしい内挿法を選択することが必要である。しばしば、そうした測定値や[[コンピュータアニメーション]]におけるキャラクターの運動などで線形補間や多項式補間が好まれて適用されるのは、単にアルゴリズムのソフトウェアへの実装が容易で計算機負荷が少ないというだけでなく、多くの物理現象を表す関数が[[テイラー展開]]可能であり、その高次の項が無視できるほど小さいと仮定できるからである。

そうでない場合は、適した内挿法を選択する必要がある。
<!--ここでは、どのような内挿法がどのようなケースで用いられるのか具体例が示されるのが望ましいとおもいます。-->

== 代表的な補間法、補間関数 ==
[[Image:Interpolation Data.svg|right|thumb|230px|例として、このような点がデータとして与えられたとき、これらの点の間の値を補間することを考える。]]
[[Image:Interpolation example polynomial.svg|right|thumb|230px|多項式補間。]]
*[[多項式補間]]
*:全てのデータ点を通る[[多項式]]を用いた補間法。
**[[ニュートン補間]]
**:[[有限差分|差分]]を用いる補間公式の一種である[[アイザック・ニュートン|ニュートン]]の補間公式を使う補間法。
**[[ラグランジュ補間]]
**:[[ジョゼフ＝ルイ・ラグランジュ|ラグランジュ]]の補間公式を使う補間法。
**ガウス補間
**:ガウスの補間式を使う補間法。
*エルミート補間(Hermite)
指定した分点において，関数の値だけでなくて微分の値も一致するような多項式を用いる補間法。さらに一般化されたものとして，より高次の微分の値も一致するような多項式による補間。
[[Image:Interpolation example spline.svg|right|thumb|230px|スプライン補間。]]
*有理関数補間
指定された分点において関数と値が一致する有理関数による補間法。さらに一般化されたものとして，関数の値だけでなく微分の値（さらに高次の微分の値）も一致するような有理関数による補間法も考えることができる。
*重心形式補間法(barycentric interpolation) 
多項式や有理関数などによる関数の補間を行なう際に，重心形式と呼ばれる形式を用いて補間を行なう方法である<ref>[https://people.maths.ox.ac.uk/trefethen/barycentric.pdf Jean-Paul Berrut and Lloyd N. Trefethen: "Barycentric Lagrange Interpolation",
SIAM Review, Vol.46, No.3, pp.501-517 (2004)]</ref><ref>Claus Schneider and Wilhelm Werner: "Some New Aspects of Rational Interpolation", Math. Comp. Vol.47, No.175, (July, 1986), pp.285-299.</ref>。
*[[スプライン補間]]
*:隣り合う点に挟まれた各区間に対し、個別の多項式を用いた補間法。各区間で、境界条件として導関数の連続性を仮定する。[[CAD]]や[[グラフィックソフトウェア]]では、滑らかな曲線や曲面を与える機能として知られる。
*フーリエ級数補間
指定された分点において関数と値が一致する有限フーリエ級数による補間法。関数が周期的なものである場合には特に有用。
[[Image:Piecewise constant.svg|thumb|right|230px|最近傍補間。]]
*0次補間（最近傍補間、最近傍点補間）
[[Image:Interpolation example linear.svg|right|thumb|230px|線形補間。]]
*[[線形補間]]（直線補間、1次補間）
*放物線補間（2次補間）
*キュービック補間（3次補間）
*:2次元信号の補間の場合、たとえば直交座標では直行する二つの軸に沿った二つの関数を計算することになる。このため、線形補間はバイリニア、3次補間はバイキュービック（双三次補間、双三次関数補間）と呼ばれる。
*キュービックコンボリューション
*:字義的には3次畳み込みという意味であるが、下記の補間関数を用いる3次補間を指すことがある。aは補間関数の性質を制御するための変数（-0.5～-2程度が用いられる）
*:<math>
h(t) = \begin{cases} 
(a + 2)|t|^3 - (a + 3)|t|^2 + 1, &  0 \leq |t| < 1 \\
a|t|^3 - 5a|t|^2 + 8a|t| -4a, & 1 \leq  |t| < 2 \\
0, & 2 \leq |t| \\
\end{cases}
</math>
*[[Sinc関数]]
*[[窓関数|Lanczos]]-n補間（ランツォシュ補間）　
*:<math>n</math> は補間関数の性質を制御するための変数。<math>n=2</math> とした補間関数はLanczos-2、<math>n=3</math> とした補間関数はLanczos-3と呼ばれる。
*:<math>
h(t) = \begin{cases} 
\mathrm{sinc} (t) \cdot \mathrm{sinc} (\frac{t}{n}), &  |t| \leq n \\
0, & n < |t| \\
\end{cases}
</math>
*[[クリギング]] - 空間的な内挿を行う[[地球統計学]]の手法

== 脚注 ==
{{脚注ヘルプ}}
{{Reflist}}

== 関連図書 ==
（今後リストを拡充する予定）
*  J. F. Steffensen: ''Interpolation'', 2nd Ed., Dover pub., ISBN 978-0-48645009-4 (2006). ※ 初出は Chelsea Pub. (1965).

== 関連項目 ==
*[[外挿]]（補外）
*[[統計学]]
*[[数値解析]]
*[[曲線あてはめ]]
*[[要素内補間]]

{{デフォルトソート:ないそう}}
[[Category:統計学]]
[[Category:補間]]
[[Category:数学に関する記事]]