0 Daumen
411 Aufrufe

Aufgabe:

Zur nichtlinearen Funktion \( f(x)=\frac{1}{x}-a \) mit Konstante \( a>0 \) soll die Nullstelle näherungsweise bestimmt werden.

a) Stellen Sie die Formel der gewöhnlichen Newton-Iteration zur Funktion \( f \) auf. Zur Kontrolle: Es folgt \( x^{k+1}=2 x^{k}-a\left(x^{k}\right)^{2} \).

b) Begründen Sie durch eine Taylor-Entwicklung der Iterationsfunktion im gewöhnlichen Newton-Verfahren, dass die Iteration hei obiger Funkion f lokal konvergent von Ordnung genau \( p=2 \) ist.


Problem/Ansatz:

Halo zusammen, könnte mir jemand bitte bei b) helfen? Ich habe was gemacht aber bin mir nicht sicher ob es richtig ist.

Meine Lösung:

Wir betrachten die Taylor-Entwicklung der Iterationsfunktion \( \Phi(x) = 2x - a \cdot x^2 \) um die Nullstelle \( \hat{x} \).

Die Taylor-Entwicklung für \( \Phi(x) \) um \( \hat{x} \) ist gegeben durch:

\( \Phi(x) = \Phi(\hat{x}) + \Phi'(\hat{x})(x - \hat{x}) + \frac{\Phi''(\hat{x})}{2!}(x - \hat{x})^2 + \ldots \)

Da wir wissen, dass \( \Phi(\hat{x}) = \hat{x} \) (da \( \hat{x} \) die Nullstelle von \( f(x) \) ist), betrachten wir nur den linearen und quadratischen Term:

\( \Phi(x) = \Phi(\hat{x}) + \Phi'(\hat{x})(x - \hat{x}) + \frac{\Phi''(\hat{x})}{2!}(x - \hat{x})^2 \)

Nun können wir den linearen Term ignorieren, da er verschwindet, da \( \hat{x} - \hat{x} = 0 \). Wir betrachten nur den quadratischen Term:

\( \Phi(x) = \hat{x} + \Phi''(\hat{x}) \cdot \frac{(x - \hat{x})^2}{2} \)

Da die Iteration des gewöhnlichen Newton-Verfahrens \( \Phi(x) \) linear gegen \( \hat{x} \) konvergiert, bedeutet dies, dass der Fehler nach einem Iterationsschritt proportional zum Quadrat des vorherigen Fehlers ist. Dies zeigt eine Konvergenzordnung von genau \( p = 2 \).

Avatar von

1 Antwort

+1 Daumen
 
Beste Antwort

Du kannst nicht einfach \(\hat x\) in Teiltermen einsetzen, um sie loszuwerden.

Schau dir nochmal die Definition der Konvergenzordnung an.


In deinem Fall ist \(\Phi\) ein Polynom 2. Grades. Damit bricht die Taylorentwicklung nach dem Glied 2. Grades ab, denn es gilt \(\boxed{\Phi''(x) = -2a}\) und beachte, dass \(\boxed{\hat x = \frac 1a}\) ist:

$$x_{n+1} = \Phi(x_n) =\underbrace{\Phi(\hat x)}_{=\hat x} + \underbrace{\Phi'(\hat x)}_{=0}(x_n-\hat x) +\frac 12 (-2a)(x_n - \hat x)^2$$

$$\Rightarrow |x_{n+1}-\hat x | = a|x_{n}-\hat x|^2 \text{ für } n\in \mathbb N$$

Das ist exakte Konvergenz der Ordnung 2.

Avatar von 10 k

Danke für deine Antwort:)


Wieso ist \( \hat{x}=\frac{1}{a} \)  ? Wie hast du das gerechnet? Und ist das wichtig zu rechnen? Wir haben es nirgendwo verwendet oder?


Und nich eine Frage: 1/2 * (-2a) = -a. Dann sollte die Lösung so aussehen :

\(\Rightarrow\left|x_{n+1}-\hat{x}\right|=-a\left|x_{n}-\hat{x}\right|^{2} \text { für } n \in \mathbb{N} \)

Oder?

Zu deiner 1. Frage:
Der interessierende Fixpunkt des Newton-Verfahrens ist die positive Nullstelle von \(f(x) = \frac 1x - a\). Das ist offenbar \(x= \frac 1a\).

(Das Newton-Verfahren dient der numerischen Bestimmung von Nullstellen in Form eines Fixpunktproblems. D.h., der Fixpunkt der Iterationsfunktion ist dann eine Nullstelle der Ausgansfunktion.)


Zu deiner 2. Frage:
Es ist \(a>0\) laut Voraussetzung. Durch die Beträge verschwindet das Minus:

$$x_{n+1}-\hat x = -a(x_{n+1}-\hat x)^2 \Rightarrow |x_{n+1}-\hat x| = a|x_{n+1}-\hat x|^2$$

\( \underbrace{\Phi^{\prime}(\hat{x})}_{=0}\left(x_{n}-\hat{x}\right) \)

Ist Φ′( x^) = 0, weil Φ′(1/a)  = 0 ist oder gibt es einen anderen Grund?

Müssen wir So schreiben oder ist es nicht wichtig?

\( \left|x_{n+1}-1/a\right|=a\left|x_{n}-1/a\right|^{2} \)

Es gibt einen allgemeineren Grund im folgenden Sinne:

\(\Phi'(x) = \frac{f(x)\cdot f''(x)}{[f(x)]^2}\) (einfach mal selber nachrechnen und \(\Phi\) differenzieren und dann vereinfachen)

Wenn \(\hat x\) eine einfache Nullstelle von f ist (also \(f'(\hat x) \neq 0\)), dann ist \(\Phi'(\hat x) = 0\).

Zu deiner letzten Frage:

Ob du \(\hat x\) oder \(\frac 1a\) benutzt, ist hier unwesentlich, solang du irgendwo hinschreibst, dass \(\hat x = \frac 1a\) ist.

Vielen Dank! :)

Ich habe noch eine Frage die auch darum geht:

Die Funktion \( f(x) \) besitze genau eine Nullstelle \( \hat{x} \in\left(0, \frac{1}{4}\right) \). Ein Iterationsverfahren zur Nullstellenbestimmung sei quadratisch konvergent im Intervall \( \left[0, \frac{1}{4}\right] \) mit der zugehörigen Konstante \( C=3 \). Gegeben sei ein Startwert \( x^{0} \in\left[0, \frac{1}{4}\right] \) und eine Genauigkeitsforderung \( \varepsilon<\frac{1}{4} \). Leiten Sie eine untere Schranke für die Schrittzahl \( k \) in Abhängigkeit von \( \varepsilon \) her, womit ein Fehler \( \left|x^{k}-\hat{x}\right| \leq \varepsilon \) garantiert werden kann.


Könntest du mir bitte dabei helfen?

Ich sehe gerade, ich hab bei \(\Phi'\) einen Schreibfehler im Nenner:
$$\Phi'(x) = \frac{f(x)\cdot f''(x)}{[f'(x)]^2}$$
Im Nenner musss die Ableitung stehen.

Deine weitere Frage solltest du als neue Frage stellen.

Danke :)

Deine weitere Frage solltest du als neue Frage stellen.

Das habe ich schon

https://www.mathelounge.de/1027849/nichtlineare-gleichungen-quadratische-konvergenz

Ich habe bi) hinbekommen. Könntest du mir bitte da bei a) und bii) helfen?

Ein anderes Problem?

Stell deine Frage

Willkommen bei der Mathelounge! Stell deine Frage einfach und kostenlos

x
Made by a lovely community