Newton-Verfahren zur Nullstellenberechnung
Motivation: Warum numerische Verfahren?
Viele Gleichungen der Form \( f(x) = 0 \) lassen sich nicht analytisch lösen. Bereits für Polynome fünften Grades oder höher existiert keine allgemeine Lösungsformel, und bei transzendenten Gleichungen wie \( e^x = 3x \) oder \( \cos(x) = x \) versagen die üblichen algebraischen Methoden vollständig. In solchen Fällen greift man auf numerische Näherungsverfahren zurück, die Schritt für Schritt immer genauere Approximationen der Nullstelle liefern. Das bekannteste und effizienteste dieser Verfahren ist das Newton-Verfahren (auch Newton-Raphson-Verfahren genannt).
Die Grundidee
Die Idee des Newton-Verfahrens basiert auf der Tangente als linearer Näherung. Man startet mit einem Schätzwert \( x_0 \) in der Nähe der gesuchten Nullstelle. An der Stelle \( x_0 \) legt man die Tangente an den Graphen von \( f \) und bestimmt deren Nullstelle. Diese Nullstelle der Tangente dient als verbesserter Schätzwert \( x_1 \). Das Verfahren wird iteriert, bis die gewünschte Genauigkeit erreicht ist.
Die Tangente an \( f \) im Punkt \( (x_n, f(x_n)) \) hat die Gleichung \( t(x) = f(x_n) + f'(x_n)(x – x_n) \). Ihre Nullstelle liegt bei \( t(x) = 0 \), was auf die Iterationsformel führt:
$$x_{n+1} = x_n – \frac{f(x_n)}{f'(x_n)}$$
Voraussetzung ist, dass \( f'(x_n) \neq 0 \), also die Tangente nicht waagerecht verläuft.
Durchführung am Beispiel
Beispiel 1: Bestimme eine Nullstelle von \( f(x) = x^2 – 2 \) (also näherungsweise \( \sqrt{2} \)).
Es gilt \( f'(x) = 2x \). Startwert: \( x_0 = 1 \).
$$x_1 = 1 – \frac{1^2 – 2}{2 \cdot 1} = 1 – \frac{-1}{2} = 1{,}5$$
$$x_2 = 1{,}5 – \frac{1{,}5^2 – 2}{2 \cdot 1{,}5} = 1{,}5 – \frac{0{,}25}{3} = 1{,}5 – 0{,}0833… \approx 1{,}4167$$
$$x_3 \approx 1{,}4167 – \frac{1{,}4167^2 – 2}{2 \cdot 1{,}4167} \approx 1{,}4142$$
Bereits nach drei Iterationen ist das Ergebnis auf vier Dezimalstellen genau: \( \sqrt{2} \approx 1{,}4142 \). Diese schnelle Konvergenz ist typisch für das Newton-Verfahren.
Beispiel 2: Löse \( \cos(x) = x \), also \( f(x) = \cos(x) – x = 0 \).
Mit \( f'(x) = -\sin(x) – 1 \) und Startwert \( x_0 = 0{,}5 \):
$$x_1 = 0{,}5 – \frac{\cos(0{,}5) – 0{,}5}{-\sin(0{,}5) – 1} \approx 0{,}5 – \frac{0{,}3776}{-1{,}4794} \approx 0{,}7553$$
Nach wenigen weiteren Iterationen konvergiert das Verfahren gegen \( x \approx 0{,}7391 \), den sogenannten Dottie-Fixpunkt.
Konvergenzverhalten
Das Newton-Verfahren konvergiert unter günstigen Voraussetzungen quadratisch: Die Anzahl der korrekten Dezimalstellen verdoppelt sich in jedem Schritt. Das macht es wesentlich schneller als einfache Bisektion (die nur linear konvergiert). Allerdings ist die Konvergenz nicht garantiert. Probleme können auftreten, wenn der Startwert zu weit von der Nullstelle entfernt liegt, wenn \( f'(x_n) \) nahe null ist (flache Tangente), oder wenn die Funktion stark oszilliert.
Wann kann das Verfahren scheitern?
Ein klassisches Gegenbeispiel ist \( f(x) = x^3 – 2x + 2 \) mit Startwert \( x_0 = 0 \). Hier springt die Iteration zwischen \( 0 \) und \( 1 \) hin und her, ohne zu konvergieren. Auch bei Funktionen mit Sattelpunkten oder mehreren nahe beieinanderliegenden Nullstellen kann das Verfahren zu einer unerwarteten Nullstelle konvergieren. Deshalb ist die Wahl eines guten Startwerts entscheidend – eine grobe Skizze oder Wertetabelle hilft dabei.
Geometrische Veranschaulichung
Geometrisch betrachtet legt man in jedem Schritt eine Tangente an den Graphen und folgt dieser bis zur x-Achse. Dort liest man den neuen Näherungswert ab und wiederholt den Vorgang. Liegt der Startwert in der Nähe der Nullstelle und ist die Krümmung des Graphen moderat, „zielt“ jede Tangente näher an die tatsächliche Nullstelle. Man kann sich das wie ein Trichter vorstellen, der die Iterationswerte zur Lösung führt.
Vergleich mit anderen Verfahren
Im Vergleich zur Bisektion (Intervallhalbierung), die nur das Vorzeichenwechselprinzip nutzt und langsam, aber sicher konvergiert, ist das Newton-Verfahren deutlich schneller, aber empfindlicher gegenüber der Startwahl. Die Regula Falsi liegt dazwischen. In der Praxis kombiniert man oft Methoden: zunächst Bisektion zur groben Eingrenzung, dann Newton zur schnellen Verfeinerung.
Zusammenfassung
Das Newton-Verfahren berechnet Nullstellen iterativ über die Formel \( x_{n+1} = x_n – \frac{f(x_n)}{f'(x_n)} \), indem es in jedem Schritt die Tangente als lineare Näherung nutzt. Es konvergiert unter guten Bedingungen quadratisch schnell, erfordert aber einen geeigneten Startwert und eine differenzierbare Funktion mit nicht-verschwindender Ableitung. Das Verfahren ist eines der wichtigsten numerischen Werkzeuge der Mathematik und findet in vielen Bereichen praktische Anwendung.