Newtonsches und lagrangesches Interpolationsverfahren

Bei der sogenannten linearen Interpolation wird zum Berechnen von Funktionswerten das Bild einer Funktion f partiell (d.h. zwischen zwei Punkten P1undP2) durch eine Gerade ersetzt. Eine bessere Annäherung an das Bild von f und damit einer größere Genauigkeit des interpolierten Wertes erreicht man, wenn man mehr Punkte heranzieht und eine Funktion ermittelt, deren Bild durch alle diese Punkte geht.

Newtonsches Interpolationsverfahren

In einem Verfahren, das auf ISAAC NEWTON (1643 bis 1727; Bild 1) zurückgeht, nähert man die Funktion durch eine ganzrationale Funktion, eine sogenannte Polynomfunktion , von möglichst kleinem Grad an. Dieses newtonsche Interpolationsverfahren sei im Folgenden dargestellt.
Von der Funktion f mit f(x) seien gegeben (s. Bild 1 bzw. folgendes Textbild)

  • die Argumente x0;x1;x2;...;xn, die sogenannten Stützstellen ,
  • die zugehörigen Funktionswerte y0;y1;y2;...;yn, die sogenannten Stützwerte.

Es wird vorausgesetzt, dass die Stützstellen (alle) voneinander verschieden sind, für die Stützwerte ist dies nicht erforderlich.BildDas Interpolationspolynom p(x) erhält man durch die folgende newtonsche Interpolationsformel :
p(x)=a0+a1(xx0)+a2(xx0)(xx1)+...+an(xx0)(xx1)(xx2)...(xxn1)
Die Koeffizienten a0;a1;a2;...;an kann man dann schrittweise aus folgenden Gleichungen berechnen:
y0=a0y1=a0+a1(x1x0)y2=a0+a1(x2x0)+a2(x2x0)(x2x1)...yn=a0+a1(xnx0)+a2(xnx0)(xnx1)+...+an(xnx0)(xnx1)(xnx2)...(xnxn1)
Wir betrachten dazu ein Beispiel.

Gegeben seien die (im untenstehenden Bild dargestellten) Punkte p(x)=a0+a1(x1)+a2(x1)(x2)+a3(x1)(x2)(x3)P0(1;2),P1(2;3),P2(3;1)undP3(4;3).
Dann gilt (wegen n=3):

Bild


Die Koeffizienten ai werden nach den oben genannten Verfahren wie folgt berechnet:
2=a0a0=23=2+a1(21)a1=11=2+1(31)+a2(31)(32)a2=32=1,53=2+1(41)1,5(41)(42)+a3(41)(42)(43)a3=76
Setzt man diese Werte in die Ausgangsformel ein, erhält man nach einigen Umformungen als Polynomfunktion:
p(x)=16(7x351x2+110x54)
Die Probe zeigt, dass das Bild dieser ganzrationalen Funktion 3. Grades durch die vier gegebenen Punkte geht.

Isaac Newton (1643 bis 1727)

Lagrangesches Interpolationsverfahren

Ein anderes Verfahren stammt vom französischen Mathematiker JOSEPH LOUIS LAGRANGE (1736 bis 1813, Bild 2).
Das k-te lagrangesche Polynom n-ten Grades (mitkn) ist folgendermaßen definiert:
Lk(x)=(xxo)(xx1)...(xxk1)(xxk+1)...(xxn)(xkxo)(xkx1)...(xkxk1)(xkxk+1)...(xkxn)
Als lagrangesches Interpolationspolynom ergibt sich dann:
L(x)=y0L0+y1L1+y2L2+...+ynLn

Beispiel: Gegeben seien (wiederum) die Punkte P0(1;2),P1(2;3),P2(3;1)undP3(4;3).

Dann sind zunächst die lagrangeschen Polynome 3. Grades L0;L1;L2undL3 zu berechnen, und es ist:
L0=(xx1)(xx2)(xx3)(x0x1)(x0x2)(x0x3)=(x2)(x3)(x4)(12)(13)(14)=16(x39x2+26x24)

L1=(xx0)(xx2)(xx3)(x1x0)(x1x2)(x1x3)=(x1)(x3)(x4)(21)(23)(24)=12(x38x2+19x12)

L2=(xx0)(xx1)(xx3)(x2x0)(x2x1)(x2x3)=(x1)(x2)(x4)(31)(32)(34)=12(x37x2+14x8)

L3=(xx0)(xx1)(xx2)(x3x0)(x3x1)(x3x2)=(x1)(x2)(x3)(41)(42)(43)=16(x36x2+11x+6)

Dann ergibt sich das lagrangesche Interpolationspolynom in der Form
L(x)=2L0+3L1+1L2+3L3,
woraus nach Einsetzen und Umformen
L(x)=16(7x351x2+110x54)
folgt. Erwartungsgemäß ist dieses Ergebnis identisch mit dem des vorangehenden Beispiels.

Ein Vergleich beider Interpolationsverfahren zeigt, dass beim newtonschen Verfahren durch die Hinzunahme einer neuen Stützstelle die bisherigen Rechnungen Gültigkeit behalten während beim lagrangeschen Verfahren die gesamte Rechnung neu begonnen werden muss.

Weitere Interpolationsformeln wurden von den Mathematikern JAMES STIRLING (1692 bis 1770), PIERRE SIMON LAPLACE (1749 bis 1829), CARL FRIEDRICH GAUSS (1777 bis 1855) und FRIEDRICH WILHELM BESSEL (1784 bis 1846) entwickelt.
Der Interpolation, bei der es darum geht, eine Funktion zu finden, deren Bild durch eine vorgegebene Menge von Punkten geht, ist die Ausgleichsrechnung verwandt. Bei dieser ist ebenfalls eine Menge von Punkten (z.B. als Ergebnis von Messungen) vorgegeben, und es wird eine Funktion (lineare, quadratische, trigonometrische Funktion bzw. Potenz-, Wurzel, Exponential-, Logarithmusfunktion) gesucht, die einen den gegebenen Punkten zugrunde liegenden Zusammenhang möglichst gut widerspiegelt. Im Gegensatz zur Interpolation wird bei der Ausgleichsrechnung nicht gefordert, dass das Bild der Funktion durch alle gegebenen Punkte geht. Im Gegenteil, da die Existenz von Messfehlern angenommen werden muss, geht es darum, die Abweichungen zu minimieren und auszugleichen sowie Aussagen über die erreichte Genauigkeit zu treffen.

Joseph Louis Lagrange (1736 bis 1813)
Lernhelfer-App für dein Smartphone oder Tablet

Learnattack

Gemeinsam zu besseren Noten: Kooperation mit Duden Learnattack

Lernvideos, interaktive Übungen und WhatsApp-Nachhilfe – jetzt Duden Learnattack 48 Stunden kostenlos testen.

Du wirst automatisch zu Learnattack weitergeleitet.
Lexikon Share
Beliebte Artikel
alle anzeigen

Einloggen