... Hab da noch düster etwas mit dem Newton Verfahren im Kopf, aber war da nicht was, dass ich mir mit den Newton Verfahren nur funktion bei 3 stützstellen erzeugen kann, die Parabelform also x^2 Form haben?
Rein mathematisch "passt" eine Kurve (n-1)-ter Ordnung immer genau durch n Punkte.
Nur ist das in der Praxis nicht hübsch und bei realen Aufgabenstellungen meißt ziemlich falsch, weil Polynome zwischen den Stützstellen ziemlich "wilde Werte" annehmen.
Wie bastelt man sich hieraus am Besten ne Funktion wo man z.B. x=250 einsetzt und ein y erhält, das nen Fehler von max 0,01 hat?
In der Praxis wird gerne Spline verwendet, weil der so schön brav und einfach ist.
Für eine Fehlerabschätzung musst Du allerdings entweder die Originalfunktion oder weitere Messwerte in Deinem Interpolationsbereich kennen.
Gibt es für sowas keine netten Programme?
Regressions- und Interpolations-Programme gibt's mit Sicherheit mehr als genug - Das nutzt Dir nur nix, wenn Du "den Weg" kennen willst.
Guck' doch mal nach - Der (kubische) Spline ist wirklich einfach:
Zwischen zwei Stützstellen (n und n+1) gibt's jeweils ein Polynom 3ter Ordnung (a
nx^3 + b
nx^2 + c
nx + d
n) mit folgenden Bedingungen:
Die Werte an den Stützstellen müssen stimmen (2 Gleichungen)
Die Ableitungen an den Stützstellen müssen gleich sein, wie die Ableitungen der angrenzenden Polynome (nochmal zwei Gleichungen)
Das gibt dann ein nettes Gleichungssystem für die ganzen unbekannte Koeffizienten a
n b
n b
n d
n
Die (dünn besetzte) Matrix mit einer diagonalen "Band-Struktur" erschlägt man z.B. mit Gauß-Jordan (mein persönlicher Liebling) -> fertig.
Viele Diagonal-Matrizen sind positiv definit (ich mag den Nachweis jetzt nicht nachgucken) - falls das der Fall ist, geht's mit dem Cholesky noch schneller.als mit Pivots.
mfg, Thomas
Edit:
Da Du nur 3 Punkte hast, kann man natürlich einfach eine (quadratische) Parabel hineinbasteln:
Punkte: (x
1, y
1), (x
2, y
2), (x
3, y
3)
Das gibt die Gleichungen: a * x
n2 + b * x
n + c = y
n
1. a * x
12 + b * x
1 + c = y
1
2. a * x
22 + b * x
2 + c = y
2
3. a * x
32 + b * x
3 + c = y
3
Wie man 3 Gleichungen (1, 2, 3) mit 3 Unbekannten (a, b, c) auflöst, muß ich ja nicht erklären.