1 ニュートンの補間の概要
1.1 目的と基本概念
ニュートンの補間は、既知の複数のデータ点からそれらを通る多項式を構成し、その多項式を用いて未知の値を求める手法である。特に、点集合が与えられた状況で新たな点が追加されても、多項式の係数を再利用しながら更新しやすい点が特徴として挙げられる。
基本となる考え方は、「補間多項式を、差の情報（等間隔では有限差分、一般には分割差）を用いた“階層的な係数”として組み立てる」ことである。これにより、ある点での値を計算するだけでなく、追加点を反映する場合の計算手順を整理できる。
1.2 補間多項式の考え方
1.2.1 点を通る多項式
補間の要請は、与えられた点 \((x_0,y_0),(x_1,y_1),\dots,(x_n,y_n)\) に対し、多項式 \(P(x)\) が各点で \(P(x_i)=y_i\) を満たすことにある。次数 \(n\) 以下の多項式は一般に \(n+1\) 個の条件を満たし得るため、補間多項式は通常は一意に定まる。
ニュートンの補間は、この補間多項式を「単に存在を主張する」のではなく、実際に計算できる形へ落とし込む。多項式は \((x-x_0)\), \((x-x_0)(x-x_1)\) といった因子を段階的に増やしながら構築されるため、係数が階層化された表現となる。
1.2.2 内挿と外挿の違い
補間多項式の利用は、目的の \(x\) がデータ点の範囲内にあるかどうかで意味が変わる。内挿（補間）は既知範囲の内部で値を推定するのに対し、外挿は範囲外での推定となる。
内挿では多項式が条件を満たす前提が活きやすい一方、外挿では高次項が支配的になりやすく、入力誤差や数値誤差が予測値へ強く影響しやすい。そのため、ニュートンの補間を用いる際も内挿を中心に考えることが実務上の基本となる。
1.3 手法の特徴
ニュートンの補間は、差分の計算結果を係数列として保持し、評価式を逐次的に適用できる。等間隔データでは前進差分や後退差分により差分表を作成し、一般データでは分割差による標準形が用いられる。
さらに、逐次追加の観点では、既存の計算を活用しながら新しい点に対応した項を追加できる。これは、計算の手戻りを減らす意味で、実装やデータ更新の運用に適した性質といえる。
2 差分に基づくニュートン補間（等間隔の場合）
2.1 前進差分と後退差分
等間隔のデータ点では、隣接する \(x\) の間隔が一定となるため、有限差分が自然に定義される。ここで \(h\) を \((x_{i+1}-x_i)\) とすると、差分表は離散列 \(\{y_i\}\) に対して計算できる。
2.1.1 前進差分表
前進差分表は、最下段に元の値列を置き、上へ向かうほど差の段数が増えるように構成する。一次差は隣り合う値の差、二次差は一次差同士の差、といった規則で積み上げる。
計算の意義は、補間多項式の係数がこれらの差分表の要素と対応する点にある。したがって表の作成は、補間計算の前処理として位置づけられる。
2.1.1.1 離散データへの適用手順
手順の流れは次のように整理できる。
- 等間隔のため \(h=x_{i+1}-x_i\) を確定する。
- 与えられた \(y_0,y_1,\dots,y_n\) を表の最下段に配置する。
- 一次差 \( \Delta y_i = y_{i+1}-y_i \) を計算し、上の段に並べる。
- 同様に二次差 \( \Delta^2 y_i\)、三次差…と繰り返し、必要段数まで作る。
- 最初の行（または特定の列）に現れる差分の要素を用いて、評価式（前進補間）を組み立てる。
これにより、目標の \(x\) に対して補間値を順次評価できる形が得られる。
2.2 ニュートン前進補間公式
前進補間は、基準点として始点 \(x_0\) を用い、増分の方向に沿って構築する。等間隔では、\((x-x_0)/h\) を基準にした無次元量を用いると表現が簡潔になる。
前進補間の形は、差分表から得られる係数を用いて \[ P(x)= y_0 + t\,\Delta y_0 + \frac{t(t-1)}{2!}\Delta^2 y_0 + \cdots \] のように逐次項を加える形になる（ここで \(t=(x-x_0)/h\)）。係数として登場する差分は、差分表の特定位置に対応しているため、計算の機械化がしやすい。
2.3 ニュートン後退補間公式
後退補間は、終点 \(x_n\) を基準に取り、減少方向の差分を用いて組み立てる。前進補間と対になる形であり、構造は類似している。
典型的には、\((x-x_n)/h\) を用いた無次元量により \[ P(x)= y_n + s\,\nabla y_n + \frac{s(s+1)}{2!}\nabla^2 y_n + \cdots \] のように表される（ここで \(\nabla\) は後退差分、\(s=(x-x_n)/h\)）。差分表からは、末尾側の要素が係数として選ばれる。
2.4 端点近傍での使い分け
どちらの公式を使うべきかは、評価する \(x\) がどの端点に近いかで判断される。前進補間は始点 \(x_0\) 付近で相対的に係数の寄与が整理されやすく、後退補間は終点 \(x_n\) 付近で同様の利点がある。
一般に、端点に近い側の形式を選ぶと、評価式の逐次計算で扱う有効な項の挙動が安定しやすい場合がある。したがって実務では、内挿領域の中心付近ではどちらでもよい一方、端に寄るほど適合する形式を優先する方針が採られることが多い。
3 分割差（一般の場合）
3.1 分割差の定義
データ点が等間隔でない場合、有限差分よりも分割差が適する。分割差は「差を一般化した量」であり、点の \(x\) 値の配置に依存して定義される。
分割差は再帰的に定義され、一次分割差は隣接する点間の傾きに対応する形を持つ。二次以上では、その差をさらに \(x\) の間隔で割った形へ拡張されるため、計算結果は点の位置関係を反映する。
3.1.1 再帰的計算
再帰的計算とは、低次の分割差を高次へ積み上げる計算方式である。一般に、同じ次数の分割差の集合を表形式で作ると、既に計算済みの値を参照しながら次の段を求められる。
この段階的構造により、ニュートンの分割差形では多項式の係数が「分割差表のある成分」として得られる。したがって、分割差表の作成は補間のコア手続きになる。
3.2 ニュートンの分割差形（標準形）
3.2.1 多項式の段階的構成
ニュートンの分割差形では、補間多項式を \[ P(x)= a_0 + a_1(x-x_0) + a_2(x-x_0)(x-x_1)+ \cdots + a_n\prod_{k=0}^{n-1}(x-x_k) \] のように、因子を段階的に増やす形で書く。ここで \(a_i\) は分割差により与えられる係数である。
この表現の利点は、評価の際に逐次的に計算できる点にある。多項式を展開して全項を明示する必要がなく、因子の積と係数の積和を順に処理する設計が可能になる。
3.2.2 係数の意味づけ
係数 \(a_i\) は、分割差表の特定の要素に一致する。一次の場合は与えられた点の間隔でスケールされた傾きの情報、二次以降は曲率に相当する変化量とみなせる。
この意味づけにより、係数は単なる数値ではなく、点列が持つ幾何学的な変化の度合いを反映する。係数がゼロに近い場合は、その次数以降の変化が小さいことを示唆し、近似の挙動を理解する手がかりにもなる。
3.3 点の追加・更新の考え方
分割差形は、点を追加したときの更新が比較的整理しやすい。新しい点を末尾側に加える運用では、既存の分割差表から再計算範囲を絞れる場合がある。
ただし、一般の点配置では新点の挿入位置によって更新手順が変わり得る。実装上は「どの順序で点を管理するか」「更新時に必要な分割差だけを再計算する設計にするか」を決めておくと、計算コストの見通しが立つ。
4 理論と実用上の注意
4.1 補間の誤差（概観）
補間多項式による推定には誤差が伴う。誤差は一般に、真の関数が多項式で表せない場合に生じ、データ点の配置や補間次数に依存する。
理論的には、誤差は高次導関数と関係する形で表されることが多い。次数を上げるほど条件を満たす多項式は柔軟になるが、データの揺れや数値誤差の影響も増えやすく、常に誤差が単調に改善するとは限らない。
4.2 数値計算における安定性
実装では、浮動小数点演算による丸め誤差が問題になることがある。特に高次数の多項式は係数の大きさや積の計算が増え、数値が不安定になり得る。
分割差計算では、差を \(x\) の差で割る操作が含まれるため、点が近すぎると桁落ちや過大な丸め誤差が起こりやすい。安定性を確保するには、点の重複回避、スケーリング、評価手順の工夫などが重要になる。
4.3 多項式の次数と過剰適合
補間は「与えた点をちょうど通る」ため、データに雑音が含まれる場合は過剰適合の形をとり得る。次数を上げることで見かけの一致は向上しても、他の点での振る舞いは悪化することがある。
このため実用では、補間次数の上限を決める、あるいは雑音除去の観点で別の手法（回帰や平滑化）を選ぶ判断が必要になる。ニュートン補間は構成が明快な一方で、入力データ品質の影響を強く受ける側面がある。
4.4 実装の指針
4.4.1 計算量とデータ構造
分割差表の作成は、一般に \(O(n^2)\) 程度の計算量を要する。評価は多項式の逐次評価により \(O(n)\) に抑えられる設計が可能である。
データ構造としては、分割差表を二次元配列で保持する方法のほか、必要成分のみを一次元配列へ畳み込んで管理する方法もある。メモリ節約の観点から、更新に必要な範囲だけを格納する構成が好ましい場合がある。
4.4.2 具体的な計算例の作り方
具体例を作るときは、まず点列を準備し、分割差（または等間隔なら有限差分）を段階的に計算する。次に、係数として利用する分割差成分を抽出し、評価点 \(x\) に対する逐次評価を実行する。
例題では、次数を小さくして手計算に近い段階から始め、結果が期待される挙動（点を通ること、評価が一貫すること）を確認する。最後に、点追加や評価位置の移動（端点近傍、中心付近）による変化を見ると、手法の特性が理解しやすい。
5 関連手法との比較
5.1 ラグランジュ補間との違い
ラグランジュ補間は、与えられた点に対して各点の寄与を基底多項式の線形結合として表す方式である。ニュートン補間は、分割差に基づく段階的な係数表現を用いる点で異なる。
計算面では、ラグランジュ形式では評価のたびに基底多項式を構成する必要が出やすい。一方、ニュートン形式は分割差を一度計算しておけば、評価式を逐次的に適用しやすい。点追加の運用でもニュートン形式が有利になりやすい。
5.2 スプライン補間との違い
スプライン補間は、多項式を区間ごとに分割し、端点や結合条件で滑らかさを保証する方式である。単一の高次数多項式で全範囲を近似するのではなく、局所性を導入する点が大きい。
ニュートン補間は基本的に全点を一つの多項式で結ぶため、高次数化の影響が出やすい場合がある。スプラインは一般に滑らかな曲線としての振る舞いが安定しやすいが、構成の自由度や実装要件は増える。用途に応じた選択が必要である。
5.3 フィッティング（回帰）との違い
フィッティング（回帰）は、データが必ずしも関数グラフ上に乗らない場合でも、傾向を捉えることを目的とする手法である。対して補間は、与えた点を厳密に通すことを前提とする。
その結果、雑音を含む観測データでは回帰の方が一般に頑健になりやすい。ニュートン補間は、入力が理想的な（あるいは測定誤差が小さい）場合に適しやすいが、データ品質が不確かなときは過剰適合が顕在化し得る。
5.4 計算目的別の選択指針
補間値の計算だけでなく、点追加や複数の評価点での再利用を重視するなら、分割差形のニュートン補間が適することが多い。等間隔データで差分表を管理できる環境では、前進・後退の公式が扱いやすい。
一方、範囲全体を滑らかに近似したい、あるいはデータが雑音を含む可能性が高い場合には、スプラインや回帰の導入が検討対象となる。最終的には、データの性質（等間隔か、雑音の有無）、要求される滑らかさ、更新頻度、計算資源などの条件で選択するのが実務的である。