Interpolation mit Splinefunktionen

[Dipl Phys Ernst Winter]

Interpolation mit Splinefunktionen
Die Interpolation einer Funktion y = f(x), die mit n + 1 Stützstellen x0, x1,... xn definiert ist, beginnt mit wachsendem n zwischen den Stützstellen zu oszillieren (Satz von Faber).
Aus dieser Erkenntnis entwickelten sich die Verfahren zur Interpolation mit Splinefunktionen. Der Name ist aus der englischen Bezeichnung eines Zeichengerätes mit einem biegsamen Stahlband hergeleitet.
Kubische Splinefunktionen
Am häufigsten werden kubische Splinefunktionen verwendet, die mit stetiger Steigung und Krümmung durch die Stützpunkte gehen.
Durch die Stützpunkte und die an den inneren Stützpunkten zu fordernde Gleichheit von Funktionswert, 1. und 2. Ableitung stehen 4n-2 Gleichungen zur Berechnung der 4n Koeffizienten bereit. Für die 2 fehlende sind 2 Randbedingungen hinzuzufügen.

Ich habe alles notwendige in der Unit Spline zusammengefasst, sie exportiert
- eine Funktionen zur Erzeugung von Splinefunktionen und
- Funktionen zur Berechnung der Funktionswerte, der Werten der 1. und 2. Ableitung und einer Stammfunktion..

Delphi-Quellcode:

function MakeSpline(Flag, n: integer; var P: array of TP; 
                    c0, cn: extended; var K: array of TK):boolean;

erzeugt das Koeffizientenarray K einer Splinefunktion aus n Stützstellen und zwei Randbedingungen. Flag, c0 und cn legen die Randbedingungen fest:
Flag Art der Randbedingungen
0 Natürliche Spline: c0, yn'' = cn
1 Spline mit s'(0) = c0, s''(n) = cn
2 Spline mit s'(0) = c0, s'(n) = cn
3 Spline mit s'(n) = cn, s''(0) = c(0)
4 Periodische Spline
In K gibt die Funktion die Koeffizienten des Ansatzes (1) zurück. Mit MakeSpline=true wird angezeigt, dass die Koeffizienten in K gültig sind. Bei MakeSpline=false ist ein Fehler aufgetreten.

Delphi-Quellcode:

function Spl(n: integer; x: extended;  var P: array of TP;
                                       var K: array of TK): extended; 
function SplAbl1(n: integer; x: extended; var P: array of TP; 
                                          var K: array of TK): extended;
function SplAbl2(n: integer; x: extended; var P: array of TP;  
                                          var K: array of TK): extended;
function SplStamm(n: integer; x: extended; var P: array of TP; 
                          var K: array of TK; var Stamm: array of extended): extended;

berechnen Werte von Splinefunktionen bzw. ihren 1. und 2. Ableitungen. Die Spline wird mit den Arrays P der Stützstellen und K der Koeffizienten übergeben. Das aufrufenden Programms hat dafür zu sorgen, dass x im Definitionsbereich [x0, xn] liegt, wobei geringe Bereichsüberschreitungen durch Rundungsfehler zulässig sind.
Man kann eine Splinefunktion analytisch Integrieren und gewinnt so die zugehörige Stammfunktion zur Berechnung bestimmter Integrale. Dies ist ein effektives Verfahren zur numerischen Integration.

[brechi]

Und was willst du mit dem Beitrag jetzt sagen? Wenn das in die Code-LIB soll dann solltest du schon den Quellcode anhängen, bzw. den Quellcode im Post mit beifügen. Ansonsten hilt der Beitrag nix und ist nur Werbung.

Zu Splines: Der Name kommt aus dem Schiffsbau.
Von welchen Splines sprichst du? Anscheinend nicht von Beziére-Splines (die interpolieren nicht, sondern approximieren nur). Also B-Splines oder NURBS etc.? (solltest du eventl. erwähnen).
Deine Programme hab ich jetzt aber nicht getestet.

[Mithrandir]

Zitat:

Zitat von brechi

Deine Programme hab ich jetzt aber nicht getestet.

Wäre vielleicht besser gewesen, dann hättest du vielleicht gemerkt, dass es sich um selbstextrahierende Archive handelt. Erst gucken, dann pupen. ;)

[Daniel]

Zitat:

Zitat von brechi

Und was willst du mit dem Beitrag jetzt sagen? Wenn das in die Code-LIB soll dann solltest du schon den Quellcode anhängen, bzw. den Quellcode im Post mit beifügen. Ansonsten hilt der Beitrag nix und ist nur Werbung.

Offenbar bist Du mit dem falschen Fuß aufgestanden. Die Sourcen sind dabei und für diesen ätzenden Tonfalls gibt's derzeit echt keinen Grund. :roll:

[jbaertle]

Wer sich für weitere (allerdings etwas angestaubte, aber funktionsfähige) numerische Routinen interessiert, findet unter efg's Delphi Math Functions doch tatsächlich einen Link zu einer frühen Version der Numerical Recipes in Pascal (nicht Delphi).

[mschaefer]

Sowas kann man immer mal wieder gebrauchen. Hatte mal Temperaturwerte zu interpolieren.

Ein praktisches Problem habe ich bei interpolationsroutinen aber meist gehabt. Die Beispiele waren immer mit festgelegten Arrays und nie mit dynamischen. Leider weiss man oft einfach noch nicht wie lange die Routinen Messwerte sammeln soll. Da mathematische Routinen Geduld und Knobelgeist erforden ist das sicher etwas für die CodeLib !

// Grüße Martin

[brechi]

Tja lag wohl an der EXE-warum ich des nicht geladen habe :)
Wie wärs mit nem ZIP-Archiv?

[Dipl Phys Ernst Winter]

"brechi"

Zitat:

Wie wärs mit nem ZIP-Archiv?

Ich halte ein selbstentpackendes Archiv für allgemeiner.
Für ein Zip-Archiv wird ein Entpacker gebraucht. Ob den wirklich jeder hat?

[Dipl Phys Ernst Winter]

"mschaefer"

Zitat:

Leider weiss man oft einfach noch nicht wie lange die Routinen Messwerte sammeln soll.

Immer erst schön sammeln und dann Interpolieren.

NP: Messwerte sollte man meistens Approximieren statt sie zu Interpolieren!

[brechi]

Du kannst es ja zusätzlich noch als ZIP-Anhängen, dann würd ich mir das auch noch anschauen, sobla ich zu Hause bin. Interessant ist es auf jeden Fall und wenn der Source dabei sein sollte auch passend für die CodeLib.