klassische Interpolationsmethode, bei der eine stetige Funktion durch eine endliche Menge von Werten eindeutig festgelegt wird.

Lagrange-Interpolation wird in der Numerischen Mathematik und der Approximationstheorie behandelt.

Es sei G ={g₀, g₁, …, g_N} ein System von N + 1 linear unabhängigen, stetigen, reell- oder komplexwertigen Funktionen, definiert auf einem Intervall [a, b], einem Gebiet G ⊂ ℂ, oder einem Kreis T. Weiter sei X ={x₀, …, x_N} eine Menge von N + 1 Punkten aus [a, b], bzw. T, mit der Eigenschaft x₀<… < x_N. Das Problem der Lagrange-Interpolation hinsichtlich G und X besteht nun darin, für reelle Werte c_i, i = 0, …, N, eine eindeutige Funktion \begin{eqnarray}g=\displaystyle \sum _{j=0}^{N}{a}_{j}{g}_{j}\end{eqnarray}

mit der Eigenschaft \begin{eqnarray}g({x}_{i})={c}_{i},\space \text{}i=0,\ldots, N,\end{eqnarray}

zu finden. Falls für beliebige c_i stets eine solche Funktion g existiert, so ist das Problem der Lagrange-Interpolation hinsichtlich G und X lösbar. X heißt in diesem Fall Lagrange-Interpolationsmenge für G.

Verallgemeinerungen der Lagrange-Interpolation, bei denen neben den Funktionswerten von g auch gewisse Ableitungen vorgeschrieben werden, nennt man Hermite-Interpolation bzw. Birk-hoff-Interpolation.

Bei der Lagrange-Interpolation spielen strukturelle Eigenschaften des zugrundeliegenden Systems G ein Rolle. Es ist bekannt, daß das Problem der Lagrange-Interpolation genau dann für jede beliebige Wahl von X lösbar ist, wenn G ein Tschebyschew-System bildet. Ein solches System G bilden beispielsweise die Polynome vom Grad N, \begin{eqnarray}G=\{1,x,\ldots, {x}^{N}\}.\end{eqnarray}

In diesem Fall läßt sich für jedes vorgegebene X und beliebig vorgegebene Werte c_i das Interpolationspolynom g in der Lagrange-Darstellung \begin{eqnarray}g(x)=\displaystyle \sum _{i=0}^{N}{c}_{i}{l}_{N,i}(x),x\in [a,b]\end{eqnarray}

angeben. Hierbei sind die Lagrange-Polynome (auch Lagrange-Fundamentalpolynome oder Lagrangesche Grundpolynome genannt) l_N,i definiert durch \begin{eqnarray}\begin{array}{cc}{l}_{N,i}(x)=\displaystyle \prod _{\begin{array}{c}j=0\\ j\ne i\end{array}}^{N}\frac{x-{x}_{j}}{{x}_{i}-{x}_{j}},i=0,\ldots, N. \end{array}\end{eqnarray}

Es gilt offenbar \begin{eqnarray}{l}_{N,i}({x}_{k})=\left\{\begin{array}{cc}1, & \text{falls}\space i=k\\ 0, & \text{falls}\space i\ne k.\end{array}\right.\end{eqnarray}

Die Lagrange-Polynome sind daher linear unabhängig. Sie hängen nur von den x_i, nicht von den c_i ab.

Man faßt diese Tatsache im klassischen Interpolationssatz von Lagrange zusammen:

Es seien N + 1 ∈ ℕ verschiedene Punkte x₀,…,x_N aus ℝ oder ℂ und N + 1 beliebige Zahlen c₀,…,c_N gegeben. Dann existiert genau ein Polynom p vom Grad höchstens N derart, daß p(x_i) = c_i für i = 1,…,N. Dieses Polynom besitzt die Darstellung\begin{eqnarray}p(z)=\displaystyle \sum _{i=0}^{N}{c}_{i}{l}_{N,i}(x)\end{eqnarray}

mit den durch (1) definierten Lagrange-Polynomen l_N,i.

Numerisch deutlich stabiler als die Lagrange-Darstellung ist die Newtonsche Interpolationsformel für g.

Falls G kein Tschebyschew-System bildet, so ist das Problem der Lagrange-Interpolation hinsichtlich G und X nur unter gewissen Zusatzvoraussetzungen an X lösbar. Betrachtet man beispielsweise ein System G von Splinefunktionen\begin{eqnarray}G=\{1,x,\ldots, {x}^{N},{(x-{t}_{j})}_{+}^{N},j=0,\ldots, k-1\},\end{eqnarray}

wobei a ≤ t₀< t₁<… < t_k−1 ≤ b und (z)₊ = max{0, z}, so gilt der folgende Satz von Schoenberg und Whitney aus dem Jahr 1953.

Das Problem der Lagrange-Interpolation hinsichtlich G und X ={x₀<… < x_N+k} ist genau dann lösbar, wenn X wie folgt über [a, b] verteilt ist:\begin{eqnarray}{x}_{j}\lt {t}_{j}\lt {x}_{j+N+1},\space j=0,\ldots, k-1.\end{eqnarray}

Einfache Charakterisierungen von Lagrange-Interpolationsmengen X dieser Art sind nicht für jedes System G möglich. So weiß man beispielsweise, daß für Splines definiert auf einem Kreis T, sogenannte periodische Splines, eine solche Verteilungsbedingung im allgemeinen nur notwendig, jedoch nicht hinreichend ist.

Für Systeme G von multivariaten Funktionen führt das Problem der Lagrange-Interpolation auf moderne und komplexe mathematische Fragestellungen, die derzeit von Approximationstheoretikern untersucht werden. Hierbei sind vor allem Systeme G multivariater Polynome und multivariater Splines von großer Bedeutung.

[1] Hämmerlin, G.; Hoffmann, K.-H.: Numerische Mathematik. Springer-Verlag Berlin, 1989.
[2] Schönhage, A.: Approximationstheorie. de Gruyter & Co. Berlin, 1971.