Splines

Beispiel eines Splines mit 8 Knoten

Geflecht: Die Querstreben verhalten sich theoretisch wie Splines

Ein Spline n-ten Grades ist eine Funktion, die stückweise aus Polynomen mit maximalem Grad n zusammengesetzt ist. Dabei werden an den Stellen, an denen zwei Polynomstücke zusammenstoßen (man spricht auch von Knoten) bestimmte Bedingungen gestellt, etwa dass der Spline (n-1) mal stetig differenzierbar ist.

Handelt es sich bei dem Spline um eine stückweise lineare Funktion, so nennt man den Spline linear (es handelt sich dann um einen Polygonzug), analog gibt es quadratische, kubische usw. Splines.

Allgemeines

Splines wurden zuerst in einer Veröffentlichung von Isaac Jacob Schoenberg im Jahr 1946 erwähnt.

Splines werden vor allem zur Interpolation und Approximation benutzt. Durch die stückweise Definition sind Splines flexibler als Polynome und dennoch relativ einfach und glatt. Dadurch ergeben sich bei der Spline-Interpolation nicht die Nachteile, die durch die starke Oszillation von Polynomen höheren Grades und deren Unbeschränktheit bei der Polynominterpolation entstehen. Splines lassen sich auch gut benutzen, um Kurven darzustellen. Hier finden sie Einsatz im CAD. Mathematisch analog lassen sich auf beide Weisen nicht nur Kurven, sondern auch Flächen beschreiben.

Wortherkunft: Der Begriff stammt aus dem Schiffbau: eine lange dünne Latte (Straklatte, englisch spline), die an einzelnen Punkten durch Molche fixiert wird, biegt sich genau wie ein kubischer Spline mit natürlicher Randbedingung. Die Straklatte ist dabei bestrebt, ihre durch die Biegungen hervorgerufene innere Spannung zu minimieren bzw. zu verteilen.

Kubische Splines

Siehe auch Kubischer C²-Spline und Bikubische Filterung.

Kubische Splines werden unter anderem zur Berechnung des Bahnverlaufes bei Achterbahnen verwendet, um die Beschleunigungen auf die Fahrgäste zu minimieren. Kubische Splines finden weitere Anwendung bei der exakten Verlegung der Schienen bei Hochgeschwindigkeitsstrecken der Eisenbahn. Auch beim Entwurf von Kurven und Oberflächen (sogen. "Freiformkurven und -flächen"), wie sie häufig im Schiff-, Flugzeug- und Automobilbau vorkommen, sind Splines von Bedeutung.

Die Eignung von Splines für solche Anwendungen liegt daran, dass für jeden Polynomabschnitt Randbedingungen erstens in Form von Punkten aber auch in Form von Werten für erste und zweite Ableitung (und in Abhängigkeit davon Steigung und Krümmung/Kurvenradius) vorgeben lassen. Dadurch kann erreicht werden, dass die Krümmung entlang der Kurve immer stetig ist. Das hat den Vorteil, dass Querbeschleunigungen beim Abfahren der Kurve immer allmählich aufgebaut werden bzw. an den Knotenpunkten vorgegebene Werte einhalten.

B-Splines

Wie auch der Raum der Polynome ist der Raum der stückweisen Polynome ein Vektorraum und hat eine Basis. Im Kontext numerischer Verfahren, wo Splines häufig eingesetzt werden, ist die Wahl der Basis entscheidend für eventuelle Rundungsfehler und damit für die praktische Einsetzbarkeit. Eine bestimmte Basis hat sich hier als am besten geeignet herausgestellt: sie ist numerisch stabil und erlaubt die Berechnung von Werten der Spline-Funktion mittels einer Drei-Term-Rekursion. Diese so genannten B-Spline-Basisfunktionen haben einen kompakten Träger, sie sind nur auf einem kleinen Intervall von Null verschieden. Änderungen der Koeffizienten wirken sich also nur lokal aus. Splines, die in dieser Basis dargestellt werden, nennt man B-Splines.

Siehe auch: Faltung (Mathematik)

Definition

Die B-Spline-Basisfunktionen $N_{i,p,\tau}\ (i=1,\ldots,n-p)$ der Ordnung $p$ mit Knotenvektor $\tau = (\tau_1,\ldots,\tau_n) \quad (n\ge 2\,p)$ werden durch die Rekursionsformel von de Boor/Cox/Mansfield definiert:

$N_{i,1,\tau}(u) = \begin{cases} 1, &amp;amp; u\in\left[\tau_i,\tau_{i+1}\right[ \\ 0, &amp;amp; \mbox{sonst} \end{cases}$

und

$N_{i,p+1,\tau}(u) = \frac{u-\tau_{i}}{\tau_{i+p}-\tau_{i}}\,N_{i,p,\tau}(u) \;+\; \frac{\tau_{i+p+1}-u}{\tau_{i+p+1}-\tau_{i+1}}\,N_{i+1,p,\tau}(u)$ .

Die Elemente des Knotenvektors heißen auch Knotenpunkte (im engl. breakpoints) und müssen die Bedingungen $\tau_{i}\le\tau_{i+1}$ und $τ i < τ i + p$ erfüllen. Die Ordnung $p$ einer B-Spline-Basisfunktion gibt dabei die Anzahl der Koeffizienten für jedes stückweise definierte Polynom innerhalb einer Basis an. Sie ist um eins höher als der Grad der verwendeten Polynome.

Im CAD werden fast immer die ersten $p$ und letzten $p$ Knotenpunkte auf den gleichen Wert gesetzt.

Eigenschaften:

Nicht-Negativität: $N_{i,p,\tau}(u)\ge 0$
Lokaler Träger: $N_{i,p,\tau}(u)= 0\mbox{ falls }u\not\in [\tau_{i},\tau_{i+p}[$
Zerlegung der Eins: $\sum_{i=1}^{n-p} N_{i,p,\tau}(u) = 1$ für $u\in[\tau_{p},\tau_{n-p+1}[$

Bemerkung:

Die Bedingungen an die Knotenpunkte $τ i$ erlauben es, dass in der Rekursionsformel unter Umständen 0 als Nenner auftritt (nämlich wenn $τ i + p = τ i$ bzw. $τ i + p + 1 = τ i + 1$ gilt). Allerdings ist dann die Funktion $N i, p,τ$ bzw. $N i + 1, p,τ$ automatisch die Nullfunktion. Auf die entsprechende Fallunterscheidung wird hier verzichtet, man ignoriere die entsprechenden Summanden in diesen Fällen (ersetze sie durch 0). Dies entspricht auch dem Grenzverhalten für z. B. $\tau_{i+p}\to\tau_{i}$ .

B-Spline-Kurve

Eine Spline-Kurve, deren Darstellung auf B-Splines beruht, nennt man B-Spline-Kurve. Bestimmt wird die Kurve durch so genannte De-Boor-Punkte, mit denen sich das Aussehen der Kurve leicht steuern lässt: Die Kurve liegt immer in der konvexen Hülle der De-Boor-Punkte, wird also von ihnen eingeschlossen.

Eine B-Spline-Kurve $C(u),\ u\in [\tau_{p},\tau_{n-p+1}]$ der Ordnung $p$ mit Knotenvektor $τ$ (s. o.) und Kontrollpunkten $P_i\ (i=1,\ldots,n-p)$ (auch De-Boor-Punkte genannt) wird definiert durch

$C(u) = \sum_{i=1}^{n-p}\; P_i\,N_{i,p,\tau}(u)$ .

Für Kurven in der Ebene sind die Kontrollpunkte 2-dimensional, für Kurven im Raum 3-dimensional.

Eigenschaften:

Lokalität: Der Kontrollpunkt $P i$ beeinflusst die Kurve nur im Intervall $[τ i,τ i + p]$
Endpunkt-Interpolation: Es ist $P 1 = C (τ p)$ , falls der erste Knotenpunkt $p$ -mal wiederholt wird und $P n - p = C (τ n - p + 1)$ , falls der letzte Knotenpunkt $p$ -mal wiederholt wird.

Eine ähnliche Darstellung haben Bézier-Kurven. Diese basieren nicht auf der oben genannten Basis, sondern auf den Bernsteinpolynomen. Genau wie bei B-Spline-Kurven die De-Boor-Punkte gibt es hier die Bézier-Punkte, die das so genannte Kontrollpolygon bilden und mit denen man die Kurve leicht graphisch darstellen kann.

B-Spline-Fläche

Eine B-Spline-Fläche der Ordnung $p$ und $q$ mit Knotenvektor $\tau = (\tau_{1},\ldots,\tau_{n})\ (n\ge 2p)$ und $\mu = (\mu_{1},\ldots,\mu_{m})\ (m\ge 2q)$ und Kontrollpunkten (bzw. De Boor Punkten) $P i j$ wird definiert durch

$C(u,v) = \sum_{i=1}^{n-p} \sum_{j=1}^{m-q}\; P_{ij}\,N_{i,p,\tau}(u)\,N_{j,q,\mu}(v)$

Die Fläche ist definiert über dem Rechteck $[\tau_{p},\tau_{n-p+1}] \times [\mu_{q},\mu_{m-q+1}]$ .

Eigenschaften:

Lokalität: Der Kontrollpunkt $P i j$ beeinflusst die Fläche nur im Rechteck $[\tau_{i},\tau_{i+p}]\times [\mu_{j},\mu_{j+q}]$
Endpunktinterpolation: Werden die ersten $p$ Knotenpunkte in $τ$ auf den gleichen Wert gesetzt, die letzten $p$ Knotenpunkte in $τ$ auf den gleichen Wert gesetzt, die ersten $q$ Knotenpunkte in $μ$ auf den gleichen Wert gesetzt und die letzten $q$ Knotenpunkte in $μ$ auf den gleichen Wert gesetzt, dann gilt die Endpunktinterpolation, d. h. $P 1,1 = C (τ p,μ q)$ , $P 1, m - q = C (τ p,μ m - q + 1)$ , $P n - p,1 = C (τ n - p + 1,μ q)$ und $P n - p, m - q = C (τ n - p + 1,μ m - q + 1)$

Literatur

Carl de Boor: A practical guide to splines. Springer Verlag, New York 1978 ISBN 0-387-90356-9 und ISBN 3-540-90356-9 (Rev. Aufl. von 2001 ISBN 0-387-95366-3)
Gerald Farin: Curves and Surfaces for CAGD. A practical guide. 5. Aufl. Academic Press, San Diego 2002 ISBN 1-55860-737-4
Günther Nürnberger: "Approximation by Spline Functions", Springer Verlag, 1989, ISBN 3-540-51618-2 und ISBN 0-387-51618-2
Hartmut Prautzsch, Wolfgang Böhm, Marco Paluszny: Bezier and B-Spline Techniques. Springer Verlag, Berlin 2001 ISBN 3-540-43761-4
David Salomon: Curves and Surfaces for Computer Graphics. 2006 Springer Science+Business Media, Inc.; ISBN 0-387-24196-5
Isaac Jacob Schoenberg: Contributions to the problem of approximation of equidistant data by analytic functions. Quart. Appl. Math., vol. 4, S. 45–99 und 112–141, 1946.

Weblinks

Wikimedia Foundation.

Игры ⚽ Поможем сделать НИР

Schlagen Sie auch in anderen Wörterbüchern nach:

splines — splaÉªn n. slat; long and flexible strip of wood or metal used for drawing curves; any of a series of ridges on a shaft which fit into slots on a corresponding shaft … English contemporary dictionary
B-Splines — Beispiel eines Splines mit 8 Knoten Geflecht: Die Querstreben verhalten sich theoretisch wie Splines Ein Spline n ten Grades ist eine Funktion, die stückweise aus Poly … Deutsch Wikipedia
Multivariate adaptive regression splines — (MARS) is a form of regression analysis introduced by Jerome Friedman in 1991.[1] It is a non parametric regression technique and can be seen as an extension of linear models that automatically models non linearities and interactions. The term… … Wikipedia
Nonuniform Rational B-Splines — Nonuniform Rational B Splines, NURBS … Universal-Lexikon
PDIFF — Splines are piecewise smooth, hence in PDIFF, but not globally smooth or piecewise linear, hence not in DIFF or PL. In geometric topology, PDIFF, for piecewise differentiable, is the category of piecewise smooth manifolds and piecewise smooth… … Wikipedia
Кубический сплайн — Некоторая функция f(x) задана на отрезке , разбитом на части , . Кубическим сплайном дефекта 1 называется функция , которая: на каждом отрезке является многочленом степени не выше третьей; имеет непрерывные первую и вторую производные на всём… … Википедия
B-Spline — Beispiel eines Splines mit 8 Knoten Geflecht: Die Querstreben verhalten sich theoretisch wie Splines Ein Spline n ten Grades ist eine Funktion, die stückweise aus Poly … Deutsch Wikipedia
B-Spline-Fläche — Beispiel eines Splines mit 8 Knoten Geflecht: Die Querstreben verhalten sich theoretisch wie Splines Ein Spline n ten Grades ist eine Funktion, die stückweise aus Poly … Deutsch Wikipedia
B-Spline-Kurve — Beispiel eines Splines mit 8 Knoten Geflecht: Die Querstreben verhalten sich theoretisch wie Splines Ein Spline n ten Grades ist eine Funktion, die stückweise aus Poly … Deutsch Wikipedia
Kubischer Spline — Beispiel eines Splines mit 8 Knoten Geflecht: Die Querstreben verhalten sich theoretisch wie Splines Ein Spline n ten Grades ist eine Funktion, die stückweise aus Poly … Deutsch Wikipedia

Academic dictionaries and encyclopedias

Splines

Inhaltsverzeichnis

Allgemeines

Kubische Splines