Spline-Interpolation

Bei der Spline-Interpolation versucht man, eine Funktion mit Hilfe von Splines zu interpolieren. Die einfachste Methode dazu ist die Verwendung von Geraden zwischen jeweils zwei benachbarten Punkten. Je mehr Punkte man nimmt und die Entfernungen zwischen diesen verkürzt, desto genauer approximiert man die gegebene Funktion. Die Extrapolation kann allerdings teilweise sehr schlecht sein. Bei der Interpolation von z. B. Messwerten muss berücksichtigt werden, dass wirklich alle Punkte einbezogen werden, d. h. so genannte Ausreißer beeinflussen das Ergebnis.

Einfacher Ansatz (Streckenzug)

Die Berechnung eines einfachen Splines als Streckenzug erfolgt auf die bekannte Weise, mit der man auch den Graphen zwischen zwei Punkten ermittelt:

$s(x)=m\cdot x+b$ , d.h.
$s(x)={\frac {y_{2}-y_{1}}{x_{2}-x_{1}}}\cdot x+y_{1}-{\frac {y_{2}-y_{1}}{x_{2}-x_{1}}}\cdot x_{1}$ (mit $m={\frac {y_{2}-y_{1}}{x_{2}-x_{1}}}$ und $b=y_{1}-m\cdot x_{1}$ )

Es ist klar, dass diese „einfachen“ Spline-Polynome – wie oben angesprochen – sehr ungenau sein können. Wesentlich bessere Ergebnisse liefern kubische Spline-Polynome.

Der kubische Spline

Es soll eine Funktion $f$ interpoliert werden.

Hierzu unterteilt man ein Intervall $[a,b]$ durch $a=x_{0}<x_{1}<\ldots <x_{n-1}<x_{n}=b$ und definiert die Stützstellen $(x_{i},f_{i}=f(x_{i})),i=0,\ldots ,n$ .

Auf jedem Teilintervall soll nun ein Polynom 3. Grades, das 2-fach stetig differenzierbar ist, angesetzt werden:

$S_{j}(x)=a_{j}+b_{j}\cdot (x-x_{j})+c_{j}\cdot (x-x_{j})^{2}+d_{j}\cdot (x-x_{j})^{3}$ für $x_{j}\leq x\leq x_{j+1}$ und $j=0,\ldots ,n-1$

Um das Gleichungssystem eindeutig zu lösen werden $4n$ Bedingungen benötigt. Für jedes der $n$ Intervalle sind zwei Interpolationsbedigungen zu erfüllen:

S_{j}(x_{j})=f_{j}\qquad \qquad j=0,1,\ldots ,n-1

S_{j}(x_{j+1})=f_{j+1}\qquad j=0,1,\ldots ,n-1

Dadurch entstehen $2n$ Bedingungen. Weitere $2n-2$ Bedingungen erhält man dadurch, dass der Spline an allen $n-1$ inneren Stützstellen zweimal stetig differenzierbar sein muss:

S'_{j}(x_{j})=S'_{j+1}(x_{j})\qquad j=1,\ldots ,n-1

S''_{j}(x_{j})=S''_{j+1}(x_{j})\qquad j=1,\ldots ,n-1

Für die weiteren 2 Bedingungen (Randbedingungen) gibt es verschiedene Möglichkeiten, so z. B.:

freier Rand oder natürlicher Spline: $S_{0}''(x_{0})=0,S_{n-1}''(x_{n})=0$
eingespannter Rand: $S_{0}'(x_{0})=f_{0}',S_{n-1}'(x_{n})=f_{n}'$ wobei $f_{0}'$ und $f_{n}'$ vorgegeben, normalerweise entweder durch die Ableitung der zu interpolierenden Funktion f oder durch eine Approximation derselben.
periodische Randbedinung: $S_{0}'(x_{0})=S_{n-1}'(x_{n})$ , $S_{0}''(x_{0})=S_{n-1}''(x_{n})$
not-a-knot (verwendet zum Beispiel vom Programm "Matlab"): Die äußeren drei Punkte werden je durch ein Polynom interpoliert.
In der Programmiersprache wird die Funktion mit dem Befehl CSpline() abgekürzt. Sie ist ebenfalls Besandteil der IEC1131.