Kurs:Numerik I/Fixpunktiteration

Einleitung

Diese Seite kann als Wiki2Reveal Folien angezeigt werden. Einzelne Abschnitte werden als Folien betrachtet und Änderungen an den Folien wirken sich sofort auf den Inhalt der Folien aus.

Einführung - mehrdimensionale Fixpunktiteration

Zunächst wird die Fixpunktiteration auf Funktionen $g:D\to D$ mit $D\subseteq {ℝ}^n$ verallgemeinern, die bereits für $n=1$, $D=[a,b]$ und hinreichend glattes $g$ diskutiert wurde.

Fixpunktiteration

Für die Bestimmung eines Fixpunktes von $g$, d. h. eines Punktes $x^{(\ast )}\in {ℝ}^n$ mit $g(x^{(\ast )})=x^{\ast }$, betrachten wir also die Iterationsvorschrift

$x^{(k+1)}:=g(x^{(k)}),k=0,1,2,\dots \text{ (FPI)}$

mit einem gegebenen Startwert $x^{(0)}\in K\subset {ℝ}^n$.

Bemerkung 1 - Fixpunktiteration

Im Folgenden werden die Abkürzung (FPI) für Fixpunktiteration $x^{(k+1)}:=g(x^{(k)})$ verwenden.

Vollständigkeit des Grundraumes

Man zeigt in dem Beweis wieder, dass die Fixpunktiteration eine Cauchyfolge liefert. Die Vollständigkeit liefert dann, dass auch ein Grenzwert der Cauchyfolge in $x^{(\ast )}\in D\subset {ℝ}^n$ existiert, d

Bemerkung 2 - Bezug zu Nullstellenverfahren

Für eine Selbstabbildung $g:D\to D$ mit $D\subseteq {ℝ}^n$

Bemerkung 3 - Stetigkeit

Die Lipschitz-Stetigkeit ist eine besondere Form der Stetigkeit, die im eindimensionalen Fall ein Beschränktheit der Sekantensteigungen liefert. Ist $g:D\to D$ differenzierbar, so ist im eindimensionalen Fall die Ableitung $|g^{\prime }(x)|\le L$. Für den mehrdimensionalen Fall wird nun die Lipschitz-Stetigkeit definiert.

Definition - Lipschitz-stetig

Sei $D\subseteq {ℝ}^n$ und $\Vert \cdot \Vert :{ℝ}^n\to {ℝ}_{+}$ eine Norm auf dem ${ℝ}^n$

• (Lipschitz-stetig) Eine Abbildung $g:D\to {ℝ}^n$ heißt Lipschitz-stetig auf $D$ mit (Lipschitz-)Konstante $L>0$, wenn gilt: $\Vert g(x)-g(y)\Vert \le L\Vert x-y\Vert ,x,y\in D(LSS).$
• (Kontraktion) Eine Lipschitz-stetige Abbildung $g:D\to D$ mit Konstante $L>0$ heißt eine Kontraktion auf $D$, wenn $L<1$ ist.

Zusammenhang - partielle Ableitungen und Lipschitz-Stetigkeit

Sei nun speziell $g\in C^1(D,{ℝ}^n)$ für $D\subseteq {ℝ}^n$, wobei wir damit eine Funktion $g:D\subseteq {ℝ}^n\to {ℝ}^n$ mit

$g(x)=(g_1(x),\dots ,g_n(x){)}^T,x\in D$

und stetigen partiellen Ableitungen $\frac{\partial g_i}{\partial x_j}(x)$ $(i,j=1,\dots ,n)$ für alle $x\in D$ meinen.

Bestimmung der Lipschitz-Konstante

Für ein solches $g$ kann man in der Praxis häufig eine Konstante $L$ aus der Definition kann mittels der partiellen Ableitungen von den Komponentenfunktionen $g_i$ bzw. der Jacobi-Matrix von $g$ bestimmt werden.

Jacobi-Matrix

Die Jacobi-Matrix ist mit $I:=\{1,\dots ,n\}$ wie folgt definiert.

${𝒥}_g(x):={(\frac{\partial g_i}{\partial x_j}(x))}_{i,j\in I}=\left(\begin{array}{cccc}\frac{\partial g_1}{\partial x_1}(x)& \frac{\partial g_1}{\partial x_2}(x)& \dots & \frac{\partial g_1}{\partial x_n}(x)\\ \frac{\partial g_2}{\partial x_1}(x)& \frac{\partial g_2}{\partial x_2}(x)& \dots & \frac{\partial g_2}{\partial x_n}(x)\\ ⋮& ⋮& \ddots & ⋮\\ \frac{\partial g_n}{\partial x_1}(x)& \frac{\partial g_n}{\partial x_2}(x)& \dots & \frac{\partial g_n}{\partial x_n}(x)\end{array}\right)\in {ℝ}^{n\times n}$

gegeben ist.

Bemerkung - Lipschitz-Stetigkeit im eindimensionalen Fall

Im eindimensionalen reellwertigen Fall bedeutet Lipschitz-Stetigkeit für differenzierbare Funktionen $f:D\to ℝ$, dass die Ableitungen $f^{\prime }(x_o)$ betragsmäßig durch $L$ für alle $x_o\in D\subseteq ℝ$ beschränkt sind, d.h.:

$|f^{\prime }(x_o)|=\underset{x\to x_o}{lim}\left|\frac{f(x)-f(x_o)}{x-x_o}\right|\le \underset{x\to x_o}{lim}\frac{L\cdot |x-x_o|}{|x-x_o|}=L$

Für die Lipschitz-Stetigkeit ist natürlich die Differnzierbarkeit keine Voraussetzung.

Definition - konvexe Menge

Eine Menge $D\subseteq {ℝ}^n$ heißt konvex, falls für je zwei Elemente $x,y\in D$ auch alle Konvexkombinationen von $x$ und $y$ zu $D$ gehört, d. h. falls

$x,y\in D,t\in [0,1]\Rightarrow (1-t)\cdot x+t\cdot y\in D.$

Bemerkung - Mittelwertsatz der Differentialrechnung

Damit lässt sich bekanntlich das folgende Lemma aus dem Mittelwertsatz der Differentialrechnung ableiten (siehe auch Konvexkombination).

Lemma - Jacobi-Matrix - Lipschitz-Stetigkeit

Es sei $g\in {𝒞}^2(D,{ℝ}^n)$ für eine offene, konvexe Menge $D\subseteq {ℝ}^n$ und für eine Konstante $L>0$ gelte

$\Vert {𝒥}_g(x)\Vert \le L,x\in D.$

Dann folgt die Abschätzung

$\Vert g(x)-g(y)\Vert \le L\Vert x-y\Vert ,x,y\in D.$

Bemerkung

Zum oben angegebenen Zusammenhang zwischen Jacobi-Matrix und Lipschitz-Stetigkeit geben wir zunächst zwei Beispiele im eindimensionalen Fall an.

Beispiel 1 - eindimensionaler Definitionsbereich

Die Funktion $f(x):=x^2$ ist auf $[0,0.4]$ eine Kontraktion, denn es ist

$f([0,0.4])=[0,0.16]\subseteq [0,0.4]$

und mit $L:=\underset{x\in [0,0.4]}{max}(2x)=0.8$

$|x^2-y^2|\le 0.8|x-y|,x,y\in [0,0.4].$

Beispiel 2 - eindimensionaler Definitionsbereich

Die Funktion $f(x):=\sqrt{x}$ ist auf $[0,a]$ mit $a>0$ nicht Lipschitz-stetig, denn für $y=0$ hat man

$|\sqrt{x}-\sqrt{0}|=\sqrt{x}=\frac{1}{\sqrt{x}}(x-0)=\frac{1}{\sqrt{x}}|x-0|,x\in (0,a]$

und $\underset{x\to 0}{lim}(1/\sqrt{x})=\mathrm{\infty }$.

Beispiel - Aufgabe 3 - mehrdimensionaler Definitionsbereich

Der Definitionsbereich $D:=[-1,+1]\times [-1,+1]\times [-1,+1]=[-1,+1{]}^3\subset {ℝ}^3$ und der Funktion $g:D\to D$ mit

$g(x_1,x_2,x_3):=(\frac{x_1^2+x_3^2}{10},\frac{x_2+x_3}{5},\frac{x_3^2}{5})$

Bestimmen Sie die Lipschitz-Konstante $L>0$. Überprüfen Sie, ob die Abbildung $g:D\to D$ eine Kontraktion mit $L<1$ ist?

Bemerkung - Banachsche Fixpunktsatz

Der folgende sog. Banachsche Fixpunktsatz gibt an, dass die Fixpunktiteration (FPI) konvergiert, und zwar für allgemeines $n$ und ohne Differenzierbarkeitsforderungen an $g$, wobei als Startvektor alle Elemente $x^{(0)}\in D\subset {ℝ}^n$ des Definitionsbereiches $D$ von $g$ zugelassen sind.

Bemerkung - Eindeutigkeit des Fixpunkt

Überdies liefert die Kontraktionseigenschaft unter den genannten Voraussetzungen die Existenz eines eindeutigen Fixpunktes. Die geforderte Kontraktionseigenschaft für die Iterationsfunktion $g$ ist allerdings eine relativ starke Voraussetzung.

Satz - Banachscher Fixpunktsatz

Sei $D\subseteq {ℝ}^n$ abgeschlossen und die Abbildung $g:D\to D$ bezüglich der Vektornorm $\Vert \cdot \Vert :{ℝ}^n\to {ℝ}_o^{+}$ eine Kontraktion mit Konstante $L\in (0,1)$. Dann gilt:

• (BF1) $g$ besitzt genau einen Fixpunkt $x^{\ast }\in D$.
• (BF2) Für jeden Startpunkt $x^{(0)}\in D$ liefert die Fixpunktiteration $x^{(k+1)}=g(x^{(k)}),k=0,1,\dots$ eine Folge $(x^{(k)}{)}_{k\in {\mathbb{N}}_o}$ in $D$, welche gegen $x^{\ast }$ konvergiert und
• (BF3) $\Vert x^{(k)}-x^{(\ast )}\Vert \le \frac{L}{1-L}\Vert x^{(k)}-x^{(k-1)}\Vert \le \frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert ,$ für $k=1,2,\dots .$

Beweis - Banachscher Fixpunktsatz

Zunächst wird die Eindeutigkeit des Fixpunktes durch Widerspruch und Anwendung der Kontraktionseigenschaft gezeigt.

Beweis (BF1) Eindeutigkeit Fixpunkt

Annahme: Es existieren zwei Fixpunkte $x^{\ast },y^{\ast }\in D$ Fixpunkte von $g$, so gilt

$\Vert x^{\ast }-y^{\ast }\Vert =\Vert g(x^{\ast })-g(y^{\ast })\Vert \le L\cdot \Vert x^{\ast }-y^{\ast }\Vert$

Beweisschritt (BF1.1) Eindeutigkeit Fixpunkt

Weil $x^{\ast }≠y^{\ast }$ nach Annahme gilt, erhält man $\Vert x^{\ast }-y^{\ast }\Vert >0$. Man erhält durch Division der obigen Gleichung durch $\Vert x^{\ast }-y^{\ast }\Vert$ einen Widerspruch mit $L\ge 1$, da nach Voraussetzung $g:D\to D$ ein Kontraktion mit $0<L<1$ ist.

Beweisschritt (BF1.2) Eindeutigkeit Fixpunkt

Also erhält man, dass $x^{\ast }=y^{\ast }$ gilt und $g$ höchstens einen Fixpunkt besitzt.

Beweis (BF2) Konvergenz für beliebige Startpunkte

Sei nun der Startvektor $x^{(0)}\in D$ beliebig gewählt und $(x^{(k)}{)}_{k\in {\mathbb{N}}_o}$ bezeichne die damit durch die Fixpunktiteration (FPI) erzeugte Folge.

Beweis (BF2.1) Konvergenz für beliebige Startpunkte

Für $x^{(k)}\in D$ ist nach (FPI) und $g:D\to D$ offenbar auch $g(x^{(k)})=x^{(k+1)}$ in $D$, so dass $x^{(k)}\in D$ für $k\in {\mathbb{N}}_0)$ und der Lipschitz-Stetigkeit wie folgt abgeschätzt werden kann.

Beweis (BF2.2) Konvergenz für beliebige Startpunkte

Man erhält somit für $g(x^{(k)})=x^{(k+1)}$ und $g(x^{(k-1)})=x^{(k)}$ folgende Abschätzung:

$\Vert x^{(k+1)}-x^{(k)}\Vert =\Vert g(x^{(k)})-g(x^{(k-1)})\Vert \le L\Vert x^{(k)}-x^{(k-1)}\Vert ,k\in {\mathbb{N}}_0,$

Beweis (BF2.3) Konvergenz für beliebige Startpunkte

Wiederholt man diese Abschätzung mit $g(x^{(k-2)})=x^{(k-1)}$ ergibt sich: :$L\cdot \Vert x^{(k)}-x^{(k)}\Vert =L\Vert g(x^{(k-1)})-g(x^{(k-2)})\Vert \le L^2\Vert x^{(k-1)}-x^{(k-2)}\Vert$.

Beweis (BF2.4) Konvergenz für beliebige Startpunkte

Durch rekursive Anwendung bis zu den Folgengliedern $x^{(0)},x^{(1)}\in D$ erhält man

$\Vert x^{(k+1)}-x^{(k)}\Vert \le L^k\Vert x^{(1)}-x^{(0)}\Vert$.

Damit kann man nun den Abstand zweier aufeinanderfolgender Folgenglieder $\Vert x^{(k+1)}-x^{(k)}\Vert$ gegen eine Potenz der Lipschitz-Konstante und dem Abstand der ersten beiden Folgenglieder abschätzen.

Beweis (BF2.5) Geometrische Reihe

Zur Vorbereitung der Cauchy-Folgeneigenschaft wird nun der Abstand von beliebigen Folgengliedern $\Vert x^{(k+n)}-x^{(k)}\Vert$ mit $n\in \mathbb{N}$ nach oben abgeschätzt. Da Potenzen der Form $L^k$ nutzt man den Grenzwert der geometrischen Reihe für diese Zweck mit:

$\sum _{k=0}^{\mathrm{\infty }}{L}^k=\frac{1}{1-L.}$

Beweis (BF2.6) Dreiecksungleichung und telekopierende Summe

Man ergänzt Nullvektoren $0_V=x^{(k+i)}-x^{(k+i)}$ und erzeugt damit eine teleskopierende Summe von Differenzen und schätzt diese mit Dreieckungleichung nach oben ab. Dies erfolgt, um die Abschätzung (BF2.4) für den Abstand von zwei aufeinander folgenden Folgenglieder $\Vert x^{(k+i+1)}-x^{(k+i)}\Vert$ nach oben abschätzen zu können.

Beweis (BF2.7) Dreiecksungleichung und telekopierende Summe

Die vorhergehenden Überlegungen liefern die folgenden Abschätzungen für $k,n\in {\mathbb{N}}_0$:

$\begin{array}{c}\Vert x^{(k+n)}-x^{(k)}\Vert =\Vert x^{k+n}+\underset{=0_V}{\underbrace{x^{(k+1)}-x^{(k+1)}}}-x^{(k)}\Vert \\ =\Vert (x^{(k+1)}-x^{(k)})+(x^{(k+2)}-x^{(k+1)})+\dots +(x^{(k+n)}-x^{(k+n-1)})\Vert \\ \le \Vert x^{(k+1)}-x^{(k)}\Vert +\Vert x^{(k+2)}-x^{(k+1)}\Vert +\dots +\Vert x^{(k+n)}-x^{(k+n-1)}\Vert \\ \le \Vert x^{(k+1)}-x^{(k)}\Vert +L\Vert x^{(k+1)}-x^{(k)}\Vert +\dots +L^{n-1}\Vert x^{(k+1)}-x^{(k)}\Vert \\ =\Vert x^{(k+1)}-x^{(k)}\Vert \cdot \sum _{k=0}^{n-1}{L}^k\le \Vert x^{(k+1)}-x^{(k)}\Vert \cdot \underset{=\frac{1}{1-L}}{\underbrace{\sum _{k=0}^{\mathrm{\infty }}{L}^k}}\end{array}$

Beweis (BF2.8) Abschätzung von aufeinander folgenden Folgengliedern

Mit (BF2.4) kann man nun zwei aufeinander folgende Folgenglieder gegen den Abstand von den ersten beiden Folgengliedern abschätzen und man erhält mit (BF2.7) für $k,n\in {\mathbb{N}}_0$:

$\begin{array}{ccc}\Vert x^{(k+n)}-x^{(k)}\Vert & \le & \Vert x^{(k+1)}-x^{(k)}\Vert \cdot \frac{1}{1-L}\\ & \le & \underset{(BF2.4)}{\underbrace{L^k\cdot \Vert x^{(1)}-x^{(0)}\Vert }}\cdot \frac{1}{1-L}\end{array}$

Beweis (BF2.9) Abschätzung für Cauchy-Folge

Also gilt für alle $k,n\in {\mathbb{N}}_0$

$\Vert x^{(k+n)}-x^{(k)}\Vert \le \frac{L}{1-L}\Vert x^{(k)}-x^{(k-1)}\Vert \le \frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert$

Demnach ist die Folge ${\left(x^{(k)}\right)}_{k\in {\mathbb{N}}_0}$ eine Cauchy-Folge.

Beweis (BF2.10) Begründung Cauchy-Folge

Da $(L^k{)}_{k\in \mathbb{N}}$ mit $0<L<1$ eine Nullfolge ist wählt man für ein beliebiges $\epsilon >0$ das $k_{\epsilon }\in \mathbb{N}$ so, dass folgende Ungleichung gilt:

$L^{k_{\epsilon }}<\frac{\epsilon \cdot (1-L)}{\Vert x^{(1)}-x^{(0)}\Vert }$

Beweis (BF2.11) Begründung Cauchy-Folge

Dies liefert für alle $k,n\in {\mathbb{N}}_0$ mit $k\ge k_{\epsilon }$ die Eigenschaft:

$\Vert x^{(k+n)}-x^{(k)}\Vert \le \frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert \le \underset{<\epsilon }{\underbrace{\frac{L^{k_{\epsilon }}}{1-L}\Vert x^{(1)}-x^{(0)}\Vert }}<\epsilon$

Beweis (BF2.12) Vollständigkeit und Cauchy-Folge

Da der Vektorraum ${ℝ}^n$ vollständig ist und $D\subset {ℝ}^n$ eine abgeschlossene Teilmenge von ${ℝ}^n$ ist, existiert ein Grenzwert $x^{\ast }\in {ℝ}^n$ und dieser liegt mit der Abgeschlossenheit von $D$ auch in $D$. Liegt der Grenz

Beweis (BF3) - Konvergenz

Mit der Abschätzung (BF2.9) erhält man für alle $k,n\in {\mathbb{N}}_0$

$\Vert x^{(k+n)}-x^{(k)}\Vert \le \frac{L}{1-L}\Vert x^{(k)}-x^{(k-1)}\Vert \le \frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert$

Damit gilt die Aussage für alle $n\in \mathbb{N}$ un die Ungleichung "$\le$" bleibt erhalten beim Grenzübergang „$n\to \mathrm{\infty }$“.

Beweis (BF3.1) - Konvergenz

Mit der Abschätzung (BF2.9) erhält man für alle $k,n\in {\mathbb{N}}_0$

$\Vert x^{\ast }-x^{(k)}\Vert \le \frac{L}{1-L}\Vert x^{(k)}-x^{(k-1)}\Vert \le \frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert$

Bemerkung (BF3.2)

Wenn nun einer solcher Grenzwert $x^{\ast }\in D$ und eindeutig ist für eine Kontraktion $g$ ist und nach Voraussetzung (Lipschitz-)stetig ist, kommte man die Abschätzung in (BF3) mit (BF3.1) für die Fixpunktiteration auch auf den Grenzwert $x^{\ast }$ der Fixpunktinteration übertragen. Der Grenzübergang „$n\to \mathrm{\infty }$“ in (BF3.1) liefert schließlich die Abschätzung (BF3).

q.e.d.

Konvergenzgeschwindigkeit

Man beachte, dass man unter den Voraussetzungen des Banachschen Fixpunktsatzes mindestens lineare Konvergenz für die Fixpunktiteration hat, denn dann gilt

$\Vert x^{k+1}-x^{\ast }\Vert =\Vert g(x^k)-g(x^{\ast })\Vert \le L\Vert x^k-x^{\ast }\Vert ,k\in {\mathbb{N}}_0,$

wobei $L\in (0,1)$ ist.

Konvergenzgeschwindigkeit

Der Ausdruck

$\frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert$

in (BF2) kann, nachdem $x^1$ berechnet wurde, für jedes $k\in \mathbb{N}$ vor Beginn der Iteration bestimmt werden.

a-priori-Fehlerabschätzung

Er ermöglicht eine a priori Fehlerabschätzung für den Approximationsfehler $\Vert x^k-x^{\ast }\Vert$. Weiter hat man wegen $L\in (0,1)$ für eine Abbruchschranke $\epsilon >0$

$\frac{L^k}{1-L}\Vert x^{(1)}-x^{(0)}\Vert \le \epsilon \iff k\log (L)\le \log \left(\frac{(1-L)\epsilon }{\Vert x^{(1)}-x^{(0)}\Vert }\right)\iff k\ge a$

mit

$a:=\log \left(\frac{(1-L)\epsilon }{\Vert x^{(1)}-x^{(0)}\Vert }\right)/\log (L),$

so dass mit (5.21) folgt:

$k\ge aℝightarrow\Vert x^{(k+1)}-x^{\ast }\Vert \le \epsilon .$

Praktisch ist also spätestens in Schritt $k:=k(\epsilon )$ mit

(5.22) $k(\epsilon )=\lceil a\rceil$

die Fehlerabschätzung $\Vert x^k-x^{\ast }\Vert \le \epsilon$ erfüllt, wobei $\lceil a\rceil$ die kleinste ganze Zahl $\ge a$ bezeichnet.

Der mittlere Ausdruck in (5.20) kann im $k$-ten Iterationsschritt bestimmt werden und erlaubt eine a posteriori Fehlerabschätzung für den Approximationsfehler $\Vert x^k-x^{\ast }\Vert$. Praktisch wird für eine gegebene Schranke $\epsilon >0$ in Schritt $k:=k(\epsilon )$ abgebrochen, wenn erstmalig

$\frac{L}{1-L}\Vert x^k-x^{k-1}\Vert \le \epsilon$

erfüllt ist. In diesem Fall genügt $x^k$ der Fehlerabschätzung $\Vert x^k-x^{\ast }\Vert \le \epsilon$.

Wir geben dazu ein Beispiel.

Beispiel 5.16

Es sei

$f(x):=x-e^{-x},x\in ℝ.$

Dann hat man

$f(x^{\ast })=0$ für $x^{\ast }\approx 0.56714329.$

Diese Nullstelle $x^{\ast }$ soll nun approximativ berechnet werden. Mit

$g(x):=e^{-x},x\in ℝ$

gilt offenbar

$g(x^{\ast })=x^{\ast }\iff f(x^{\ast })=0,$

so dass wir dazu die Fixpunktiteration $x_{k+1}=g(x_k)$ mit der Iterationsfunktion $g$ anwenden wollen.

Dafür müssen wir zunächst die Voraussetzungen von Satz 5.15 überprüfen. Auf dem Intervall $D:=[0.5,0.69]$ ist $g$ monoton fallend und somit

$g(D)=[e^{-0.69},e^{-0.5}]\subseteq [0.5,0.61]\subseteq D$

sowie

$L:=\underset{x\in [0.5,0.69]}{max}|g^{\prime }(x)|=\underset{x\in [0.5,0.69]}{max}{e}^{-x}=e^{-1/2}\approx 0.606531.$

Also ist $g$ eine Kontraktion. Die folgende Tabelle liefert für den Startwert $x_0:=0.55$ ausgewählte Iterierte des Verfahrens:

$\begin{array}{cccccc}k& x_k& k& x_k& k& x_k\\ & & ⋮& ⋮& ⋮& ⋮\\ 0& 0.55000000& 10& 0.56708394& 20& 0.56714309\\ 1& 0.57694981& 11& 0.56717695& 21& 0.56714340\\ 2& 0.56160877& 12& 0.56712420& 22& 0.56714323\\ 3& 0.57029086& 13& 0.56715412& 23& 0.56714332\\ 4& 0.56536097& 14& 0.56713715& 24& 0.56714327\\ ⋮& ⋮& ⋮& ⋮& ⋮& ⋮\end{array}$

Die Situation soll nun für $k=12$ genauer betrachtet werden. Die Fehlerabschätzung (5.20) liefert in diesem Fall

$\frac{L^{12}}{1-L}|x_1-x_0|=\frac{0.606531^{12}}{1-0.606531}0.02694981=1.70\cdot 10^{-4},$

$\frac{L}{1-L}|x_{12}-x_{11}|=\frac{0.606531}{1-0.606531}0.00005275=8.13\cdot 10^{-5}.$

Der tatsächliche Fehler

$|x_{12}-x^{\ast }|\approx 1.91\cdot 10^{-5}$

wird also durch die a posteriori Abschätzung um etwa den Faktor 4 und durch die a priori Abschätzung um etwa den Faktor 9 überschätzt.

Praktisches Vorgehen

Das praktische Vorgehen soll nun für die spezielle Fehlerschranke $\epsilon =0.0076$ illustriert werden. Der (üblicherweise unbekannte) Approximationsfehler unterschreitet diese Schranke bereits für $k=2$, denn man hat $|x_2-x^{\ast }|\approx 0.0055\le \epsilon$. Die a posteriori Abschätzung liefert dagegen für $k=4$

$|x_4-x^{\ast }|\le \frac{e^{-1/2}}{1-e^{-1/2}}|0.56536097-0.57029086|=0.007599\le \epsilon ,$

also $k(\epsilon ):=4$ als Stoppzahl, während die a priori Abschätzung in (5.22) mit

$a=\log \left(\frac{(1-L)\epsilon }{\Vert x^1-x^0\Vert }\right)/\log (L)=\log \left(\frac{(1-e^{-1/2})0.0076}{0.57694981-0.55}\right)/\log (e^{-1/2})\approx 4.397$

die Vorhersage $k(\epsilon ):=5$ macht.

Siehe auch

• Konvexkombination

Seiteninformation

Diese Lernresource können Sie als Wiki2Reveal-Foliensatz darstellen.

Wiki2Reveal

Dieser Wiki2Reveal Foliensatz wurde für den Lerneinheit Kurs:Numerik I' erstellt der Link für die Wiki2Reveal-Folien wurde mit dem Wiki2Reveal-Linkgenerator erstellt.

• Die Seite wurde als Dokumententyp PanDocElectron-SLIDE erstellt.
• Link zur Quelle in Wikiversity: https://de.wikiversity.org/wiki/Kurs:Numerik%20I/Fixpunktiteration
• siehe auch weitere Informationen zu Wiki2Reveal und unter Wiki2Reveal-Linkgenerator.