Interpolationsfehler bei der Interpolation mit linearen Splines
# Interpolationsfehler bei der Interpolation mit linearen Splines

Der Interpolationsfehler bei der Interpolation mit linearen Splines ist fundamental bei der Gewinnung von Fehlerabschätzungen für die Finite-Elemente-Methode.

## Eindimensionaler Fall

Auf dem Intervall $[a,b]$ werde eine stetige Funktion $u$ stückweise linear interpoliert. D.h.: in den Gitterpunkten $a=x_{1}<x_{2}<\cdots <x_{i}<\cdots <x_{N}=b$ wird interpoliert, und auf dem Intervall $[x_{i-1},x_{i}]$ der Länge $h_{i}$ ist die Interpolierende $u^{I}$ linear mit

Es sei $h=\max _{i}h_{i}$. Für eine zweimal stetig differenzierbare Funktion $u$ gilt die Interpolationsfehlerabschätzung

In Hinsicht auf partielle Differentialgleichungen und den zweidimensionalen Fall ist die Voraussetzung der zweimaligen stetigen Differenzierbarkeit unbefriedigend. Anstrebenswert sind Interpolationsfehlerabschätzungen für Funktionen aus gewissen Sobolev-Räumen. Im Fall linearer Splines ist dies der Raum $H^{2}(a,b)$. Funktionen aus diesem Raum sind stetig, so dass die Interpolierende definiert ist. Aus den Darstellungen

und

erhält man die Interpolationsfehlerabschätzungen in der $L^{2}(a,b)$-Norm und der Semi-Norm im $H^{1}(a,b)$:

## Zweidimensionaler Fall

Im zweidimensionalen Fall hat sich für die Interpolationsfehlerabschätzung eine Technik durchgesetzt, bei der Integrale über die Elemente der Zerlegung auf Integrale über ein Referenzelement transformiert werden, dann das Bramble-Hilbert-Lemma angewandt wird und im dritten Schritt zurücktransformiert wird.

Betrachtet wird ein polygonales Gebiet $\Omega$ und eine zulässige Zerlegung in Dreiecke. Für $u\in H^{2}(\Omega )$ sei $u^{I}$ die stückweise lineare Interpolierende von $u$, die in allen Ecken der Dreiecke mit $u$ übereinstimmt. Es sei $K$ ein Dreieck bzw. Element der Zerlegung in der $(x,y)-$ Ebene und $E$ mit den Ecken $(0,0)$, $(1,0)$ und $(0,1)$ das Referenzelement in der $(\xi ,\eta )-$ Ebene. Besitzt $K$ die Ecken $(x_{i},y_{i})$, so wird die Abbildung $K\leftrightarrow E$ vermittelt durch

Abzuschätzen sei nun zunächst $\|u-u^{I}\|_{0,K}=[\int _{K}(u-u^{I})^{2}dK]^{1/2}$. Bei der Transformation auf $E$ kommt die (konstante) Funktionaldeterminante $D_{K}$ ins Spiel, $({\hat {u}},{\hat {u}}^{I})$ seien die transformierten Größen auf $E$. Schritt 1 liefert

Nun wird das Bramble-Hilbert-Lemma angewandt. $[\int _{E}({\hat {u}}-{\hat {u}}^{I})^{2}dE]^{1/2}$ ist ein auf $H^{2}(E)$ beschränktes, sublineares Funktional, das für lineare ${\hat {u}}$ verschwindet. Also gilt

Im dritten Schritt wird zurücktransformiert. Es müssen dazu die zweiten Ableitungen bezüglich $(\xi ,\eta )$ in Ableitungen bezüglich $(x,y)$ umgerechnet werden. Die Kettenregel liefert dies. Die Größen $x_{\xi },y_{\xi },x_{\eta },y_{\eta }$ sind leicht berechenbar, sie können alle durch $h_{K}$ (Durchmesser von $K$) abgeschätzt werden. Das ergibt

Zusammenfügen der Teilergebnisse impliziert mit $h=max_{K}h_{K}$ die gewünschte Abschätzung in der $L^{2}-$norm:

Wichtig für die Fehlerabschätzung für die Finiten-Element-Methode ist der Interpolationsfehler in der Semi-Norm des $H^{1}$. Schritt 2 und 3 verlaufen analog wie bei der Abschätzung des $L^{2}-$Fehlers. Im ersten Schritt muss man jetzt aber die Ableitungen nach $x,y$ in Ableitungen nach $\xi ,\eta$ im Integral

umrechnen. Die Ableitungen $\xi _{x},\eta _{x}$ erhält man z. B., indem man die $x,y$ definierenden Gleichungen nach $x$ differenziert und das entsprechende Gleichungssystem löst. Die Koeffizientendeterminante dieses Systems ist $D_{K}$. $D_{K}$ ist gleich dem zweifachen Flächeninhalt von $D_{K}$, und dieser ist das Produkt vom Umfang und dem Inkreisradius $\rho _{K}$ von $K$. Daraus folgt

Schritt 2 und 3 liefern dann

Ist der Quotient ${\frac {h_{K}}{\rho _{K}}}$ für alle $K$ gleichmäßig beschränkt, liegt also eine quasiuniforme Triangulierung vor, so folgt die gewünschte Abschätzung

Die Quasiuniformität (Minimalwinkelbedingung) ist hinreichend für diese Abschätzung, aber nicht notwendig. Hinreichend ist auch die Maximalwinkelbedingung.

## Stückweise polynomiale Approximation

Interpoliert man mit einer stetigen Interpolierenden, die stückweise polynomial ist vom Grad $k\geq 2$, und ist die Triangulierung quasiuniform, so verbessert sich die Appoximationsordnung, falls $u$ glatter ist mit $u\in H^{k+1}(\Omega )$:

Der Beweis erfolgt analog zu dem für $k=1$.

