1.2. LAGRANGE-Interpolation

Bei gegebenen (n + 1) paarweise verschiedenen Stützstellen (xi,f(xi)), i = 0, 1, ..., n, lässt sich mittels (1.2) immer eindeutig ein Polynom finden, das die Interpolationsforderung erfüllt. Aber selbst mit Hilfsmitteln wie dem GAUSS-Algorithmus ist es schwierig, die Koeffizienten eines Polynom höheren Grades zu bestimmen.

Für die LAGRANGEsche Darstellung des Interpolationspolynoms wird folgender Ansatz gemacht:
pn(x) = f0L0(x) + f1L1(x) + ... + fnLn(x)(1.5)
Mit gegebenen f(xi) = fi, xi - paarweise verschieden - gilt:
Li(x) = nprod k=0 k≠i x - xk
xi - xk
Dabei heißen die Li(x) LAGRANGEsche Multiplikatoren.

Zuerst muss bewiesen werden, dass der Ansatz der LAGRANGE-Darstellung die Interpolationsforderung erfüllt.
Beweis 1.1:

Voraussetzung: x = xj
Behauptung: pn(xj) = f(xj), j = 0, 1, ..., n
Beweis:
pn(xj) = f0L0(xj) + f1L1(xj) + ... + fn-1Ln-1(xj) + fnLn(xj)
pn(xj) = nsum i=0 ( fi nprod k=0,k≠i xj - xk
xi - xk
)
pn(xj) = fj + nsum i=0,i≠j ( fi nprod k=0,k≠i xj - xk
xi - xk
)
pn(xj) = fj + nsum i=0,i≠j ( fi xj - xj
xi - xj
nprod k=0;k≠i,j xj - xk
xi - xk
)
pn(xj) = fj + nsum i=0,i≠j 0
pn(xj) = fj
QED

Satz 1.1:
Sei Li(x) = nprod k=0 k≠i x - xk
xi - xk
, dann gilt nsum j=0 Lj(x) = 1.
Beweis 1.2:
Es gibt nur ein Polynom n-ten Grades pn(x), für das pn(xj) = fj gilt. Deshalb gibt es für die Stützwerte fj = 1, j = 0, 1, ..., n, auch nur das eindeutige Polynom pn(x) = 1. Wegen fj = 1, j = 0, 1, ..., n, ist pn(x) = nsum j=0 Lj(x) = 1.
QED

Bemerkung 1.4:
Es gilt Lj(xi) =
{

0, falls i ≠ j
1, falls i = j

Beispiel 1.2:
Gegeben sind 4 Stützstellen (xi,f(xi)), i = 0, 1, 2. Stelle p3(x) als LAGRANGEsches Interpolationspolynom dar.
x -1 0 2 3
f(x) -1 3 11 27
p3(x) = -1
-12
(x-0)(x-2)(x-3) + 3
6
(x+1)(x-2)(x-3) + 11
-6
(x+1)(x-0)(x-3) + 27
12
(x+1)(x-0)(x-2)

x
f(x)
p2(x) =
(x-) (x-) +
(x-) (x-) +
(x-) (x-)

Bemerkung 1.5: Ein Nachteil der LAGRANGE-Interpolation ist, dass bei Hinzunahme einer weiteren Stützstelle die LAGRANGEschen Multiplikatoren neu berechnet werden müssen.

Hosted by www.Geocities.ws

1