3次エルミートスプライン

数値解析において、3次エルミートスプラインまたは3次エルミート補間は、各部分がエルミート形式で指定される3次多項式、つまり対応する区間の両端の値と1次微分によって指定されるスプラインである

。

3次エルミートスプラインは、通常、与えられる引数値で指定される数値データを補間して連続関数を得るために使用される。データは、各についての所望の関数値と微分から構成される必要がある（値のみが与えられている場合、微分はそれから推定する必要がある）。エルミート補間多項式英語版Wikipediaでは、単に "The Hermite formula"は各区間に個別に適用される。結果のスプラインは連続、かつ、一次微分連続となる。

3次多項式スプラインは他の方法でも指定でき、ベジェ形式が最も一般的である。しかし、これら2つの方法で得られるスプラインセットは同一であり、ベジェ形式とエルミート形式の間でデータの変換も容易であるため、これらの名称はしばしば同義語のように使用される。

3次多項式スプラインは、コンピュータグラフィックスやにおいて、平面または3次元空間上の指定された点を通過する曲線や運動軌跡を得るために広く用いられている。これらのアプリケーションでは、平面または空間の各座標は、個別のパラメータ t を持つ3次スプライン関数によって個別に補間される。 3次多項式スプラインは、などの構造解析アプリケーションでも広く用いられている。

3次多項式スプラインは、死亡率分析や死亡率予測にも応用されている。

3次スプラインは、いくつかの方法で2つ以上のパラメータを持つ関数に拡張できる。バイキュービックスプライン（バイキュービック補間）は、デジタル画像のピクセル値や地形の高度データなど、規則的な矩形グリッド上のデータを補間するためによく使用される。バイキュービック曲面パッチ英語版Wikipediaでは "Bicubic surface patches, defined by three bicubic splines" で、直訳すると『3つのバイキュービックスプラインで定義されるバイキュービック曲面パッチ』と理解困難な文章になってしまうため、"defined by three bicubic splines" を省いて和訳した。なお、3Dグラフィックス関連の文献で、3次（双3次ではない）ベジェ曲線を "bicubic Bézier curve" としている例が散見される。一例：は、コンピュータグラフィックスにおいて不可欠なツールである。

3次スプラインは、特にコンピュータグラフィックスではしばしばcsplineと呼ばれる。エルミートスプラインは、シャルル・エルミートにちなんで名付けられた。

単一区間の補間単位区間 [0, 1]

単位間隔について、における始点、における終点、および、における開始接ベクトル、における終了接ベクトルが与えられると、多項式は次のように定義される。

ただし t ∈ [0, 1] である。

任意区間の補間任意の区間内のでの補間は、アフィン（1次）変数変換により区間をにマッピングすることで行われる。式は次のようになる。

ただしとは、以下に定義される基底関数を参照。なお接ベクトルの長さは、単位区間の方程式に対し倍にスケーリングされる。

一意性上記の式は、指定された接ベクトルを持つ2点間の、一意の3次多項式の経路を示す。

証明は、与えられた境界条件を満たす2つの3次多項式とする。とすると、

である。とはいずれも3次多項式であるので、は最大で3次の多項式である。そのため、は以下の式でなければならない。

導関数を計算すると、

さらに、以下が明らかである。

式 () と () を合わせると、が導かれ、したがって、よってとなる。

表現

単位区間上の補間多項式は次のように書ける（任意区間については、上記再スケール版を参照）。

ただし , , , はエルミート基底関数である。これらは様々な方法で記述することができ、それぞれ異なる性質を示す。

展開因数分解バーンスタイン

「展開」列は上記の定義で使用された表現を示している。「因数分解」列は、とが端点でゼロであることを直接示している。さらに、とは0で重複度2の零点を持ち、とは1で重複度2の零点を持つため、それら端点で傾きが0になると言える。「バーンスタイン」列は、エルミート基底関数を3次のバーンスタイン多項式に分解したものを示しており、

である。

この関係を利用すると、3次エルミート補間を 3次ベジェ曲線の形式で、次の4つの値によって表現でき、ド・カステリョのアルゴリズムを用いてエルミート補間を行うことができる。これは、3次ベジェ曲線英語版Wikipediaでは "a cubic Bézier patch" だが、曲線についての記述であるため、"patch" を "curve" に置き換えて和訳した。において、補間曲線の両端点における接ベクトルが、中間の2つの制御点により決定されることを示している。

この多項式は次の標準形式で記述することもできる：

ただし、制御点と接ベクトルが係数である。この形式では係数が一定となり、一度計算して再利用できるため、t の様々な値について多項式を効率的に評価できる。

行列形式では

となる。

ベジェ曲線への変換 3次ベジェ曲線の行列形式は

であることから、3次エルミートスプラインと等しくなる3次ベジェ曲線の制御点は

である。

データセットの補間あるデータセット（）は、各区間に上記手順を適用することで補間できる。ただし、接ベクトルを適切な方法で選択する、つまり端点を共有する区間どうしで接線を等しくする。補間された曲線は区分3次エルミートスプラインで構成され、において大域的に連続して微分可能である。

接ベクトルの選択は一意ではなく、いくつかのオプションが利用可能である。

有限差分

最も単純な選択肢は3点差分であり、一定間隔のデータ点を必要としない。内部の点に対しては、

とし、データ集合の端点では片側差分とする。

カーディナルスプライン

を用いて接ベクトルを計算すると、カーディナルスプライン（英語: cardinal spline）が得られる

。なお、カーディナルスプラインは正準スプライン（英語: canonical spline）と呼ばれることもある

。パラメータはテンションパラメータであり、の範囲でなければならない。ある意味では、これは接ベクトルの「長さ」と解釈できる。を選択するとすべての接ベクトルがゼロになり、を選択すると一様パラメータ化（英語: uniform parameterization）Catmull-Romスプラインになる。

Catmull–Romスプライン

接ベクトルとして

を選択することによりカーディナルスプラインの特殊な場合のひとつであるCatmull-Romスプラインが得られる。「一様な」パラメータ間隔を前提とする。

Kochanek–Bartelsスプライン

Kochanek–Bartels スプラインは、データポイント , , が与えられた場合に接ベクトルを選択する方法をさらに一般化したものであり、テンション、バイアス、連続性の3つのパラメータを指定できる。

単調3次補間

上記のいずれかのタイプの3次エルミートスプラインを単調データセットの補間に使用する場合、補間された関数は必ずしも単調にはならないが、接ベクトルを調整することで単調性を維持できる。

端点の微分が一致する単位区間の補間点のひとつの座標成分を、整数座標 x = n − 1, n, n + 1, n + 2 で関数 f(x) が取る値とすると、

と表すことができる。

さらに、端点における接ベクトルは、隣接点の中心差分で定義されると仮定する。

実数値 x における補間値 f(x) を評価するには、まず x を整数部 n と小数部 u に分離する。

ここでは床関数を表し、x 以下の最大の整数を返すものとする。

Catmull-Romスプラインは Two hierarchies of spline interpolations. Practical algorithms for multivariate higher order splines.

ただし、は行列の転置を表す。最後の式はホーナー法を用いた計算手法を示す。

この記述はトリキュービック補間（1回の最適化につき、同じ u と異なる p を使用して CINTu を16回計算する必要がある）にも関連する。

3次エルミートスプライン

Źródło i licencja