Näherungsverfahren in der Mathematik - Nullstellenberechnung nach Newton, Flächenberechnung mit 2 verschidenen Rechtecksformeln, mit der Trapezformel und der Formel nach Simpson

Das NEWTONsche Näherungsverfahren

Grundlegend anders funktioniert das von Sir Isaac Newton ¹ entwickelte Näherungsverfahren zur Nullstellenberechnung. Der Grundgedanke dabei ist, daß der Schnittpunkt einer Tangente mit der x-Achse eines beliebigen (Näherungs-)Punktes der gesuchten Nullstelle einen genaueren Näherungswert liefert. Wiederholt man nun mit diesem neuen Wert diese Prozedur bis zu einem Näherungswert x n so entspricht dieser der Approximation an eine Nullstelle:

. Wählt man unterschiedliche Startwerte, so erhält man alle Lösungen einer Gleichung; entspricht der Startwert aber der x-Kordinate eines Extrema, das heißt, daß dessen Tangente parallel zur x-Achse läuft und daher den Anstieg besitzt, so führt das

Newtonsche Näherungsverfahren zu keiner Lösung. In diesem Fall muß lediglich mit einem anderen Startwert angefangen werden.

Wie für jede andere Gerade gilt auch für die Tangente t an einen beliebigen Punkt

der Funktion :

.

Anstelle von können wir auch schreiben:

daraus folgt durch Einsetzen in obige Gleichung:

Das entspricht der 1. Ableitung im Schnittpunkt der Tangente mit der x-Achse, also , und der Anstieg ist an allen Stellen der Gerade gleich:

Setzten wir nun in die Tangentengleichung ein, so folgt:

Schneiden wir diese nun mit der x-Achse, so führt dies zum Newtonschen Näherungsverfahren:

Eine andere Ableitungsmöglichkeit für das Newtonsche Näherungsverfahren basiert auf der Tatsache

. Und da ja bekanntlich

Taschenrechnerprogramm: [HP 32SII]

LBL A

INPUT A

INPUT B INPUT C INPUT D INPUT E INPUT X

LBL B

RCL X

4 _ RCL_A RCL X 3 _ RCL_B + RCL X 2 _ RCL_C + RCL X RCL_D + RCL E + STO G VIEW G PSE RCL X

3

_ RCL_A 4 _ RCL X 2 _ RCL_B 3 _ + RCL X RCL_C + RCL+D STO H VIEW H PSE RCL G RCL÷H STO I VIEW I PSE RCL X RCL-I STO N VIEW N PSE PSE RCL N STO X GTO B

Dieses Programm eignet sich für Funktionen des Ausdruckes

oder solche mit niederen Potenzen. Nach Eingabe des Terms (in unserem folgenden Beispiel:

besonders beim ,,händischen" Ermitteln der Nullstelle hilfreich ist, da man dabei mit einer Tabelle am einfachsten arbeiten kann.

Zuerst bilden wir die erste Ableitung der Funktion:

Jetzt setzten wir einfach in die Tabelle ein mit dem Startwert , da bei diesem Punkt

(unter Betrachtung der 1. Ableitung) kein Extremwert zu erwarten ist.

Dieser Wert entspricht bereits sehr genau dem mittels Solver errechneten Punkt:

Für weitere Nullstellen müssen (in der Regel) nur andere Startwerte eingegeben werden.

Flächenberechnung

Von ähnlich elementarer Bedeutung wie das Berechnen von Nullstellen einer Funktion ist die Integration einer solchen - also die Flächenberechnung zwischen einer Kurve und der x-

Achse. Dies spielt auch in der Physik - z.B. bei einem Zeit/Geschwindigkeit - Diagrammeine wichtige Rolle. In diesem Fall ist der zurückgelegte Weg nach der Formel die in

einer bestimmten Zeit von der Geschwindigkeitskurve aufgespannte Fläche. Bedauerlicher Weise läßt sich durch elementares Integrieren nicht immer beziehungsweise nur äußerst schwer die Stammfunktion von jeder beliebigen Kurve finden. Für diesen Fall gibt es sogenannte numerische Näherungsverfahren, die den Wert des in Rede stehenden Integrals im allgemeinen nicht exakt liefern, dafür aber ,,beliebig genau", wenn man nur genügend oft eine Rechenprozedur wiederholt. Dadurch versucht man, den Flächeninhalt durch möglichst einfache geometrische Figuren zu approximieren.

Die Rechtecksformel

Wie auf nebenstehender Abbildung gezeigt, soll die Fläche A der stetigen Funktion

zwischen den Grenzen berechnet werden, unter der Bedingung, daß es keine Nulllstellen innerhalb der Integrationsgrenzen gibt.

Dieser Flächeninhalt könnte nun also in beliebig viele schmale Rechtecke von derselben Breite unterteilt werden. Approximieren wir durch n Rechtecke so ist:

und für die ,,Zwischenpunkte" x i gilt allgemein:

Da jedes Rechteck die Breite und die Höhe bei Verwendung von n Rechtecken:

LBL A

0 STO S INPUT A INPUT B INPUT N RCL B RCL-A RCL÷N STO D

LBL B

T (statt X aber A !)

E R M STO+S RCL D STO+A RCL B RCL A

x

Da diese Methode bei einigen (steilen) Kurven unter Umständen erst für sehr hohe n Werte sich an das tatsächliche A nähert, verwendet man in der Praxis als Höhe eines einzelnen

Rechteckes

Dadurch hebt sich im Mittel das ,,Zukleinsein" des Rechteckes zirka auf:

LBL A

0 STO S INPUT A INPUT B INPUT N RCL B RCL-A

RCL÷N STO D 2 ÷ STO+A

LBL B

T (statt X aber A !) E R M STO+S RCL D STO+A RCL B RCL A x

Beispiel:

Die Fläche soll zwischen

Intervall!): Dazu integrieren wir einfach elementar.

Die Trapezformel

Ähnlich wie die Zerlegung des Flächenstückes in Rechtecke erfolgt das Verfahren mit der Trapezformel, wobei aber anstelle der Rechtecke Trapeze die Fläche annähern.

Unter an sonst gleichen Voraussetzungen gilt auch

für die Breite eines Trapezes.

Und für die ,,Zwischenpunkte" x i gilt ebenso allgemein:

Für den Flächeninhalt des i-ten Trapezes gilt somit:

Bei der Verwendung von n Trapezen gilt also:

LBL A

INPUT A

INPUT B INPUT N 1 x=y? GTO A RCL B RCL-A RCL÷N STO D

LBL B

T (statt X aber A !) E R M STO S RCL D STO+A

LBL C

T (statt X aber A !) E R M 2 * STO+S RCL D STO+A RCL A RCL B xy? GTO C

LBL D

T (statt X aber A !) E R M STO+S RCL D RCL*S

2

÷ STO F VIEW F GTO A

Die Fläche soll wieder zwischen berechnet werden.

Die SIMPSONsche Formel

Obwohl die Trapezformel meist recht gute Näherungswerte liefert, erhält man noch genauere Werte, wenn die zu gehörige Kurve nicht durch Geradenstücke, sondern durch

Parabelbögen approximiert wird. Wie bei den beiden bereits besprochenen Methoden wird auch hier das Grundintervall in n gleich breite Teilintervalle zerlegt. Dabei wird jedes Intervall durch einen Parabelbogen mit der allgemeinen Gleichung

ersetzt, der durch die Kurvenpunkte an den Rändern dieses Teilintervalls und durch den

Kurvenpunkt in der Mitte des Intervalls hindurchgeht. Die Kurvenpunkt a 0 , a 1 und a 2 findet

man durch einen sogenannten unbestimmten Ansatz; das heißt, daß für eine Funktion n-ten

Grades n+1 Punkte/Informationen der Kurve benötigt werden. Durch Einsetzten aller n+1

Punkte in den allgemeinen Funktionsterm, erhält man n+1 Gleichungen mit ebenso n+1

Unbekannten, was leicht durch Eliminieren gelöst werden kann. Im Fall der Parabel benötigt

man also 3 Punkte. Insgesamt hat man für die n Parabelstücke viele äquidistante, das heißt mit gleichem Abstand zu einander, Zwischenpunkte bis zu wählen: im Intervall

Die Punkte bestimmen die 1-te Parabel, und die Punkte die 2-te Parabel, und

die Punkte die k-te Parabel, und

die Punkte die n-te Parabel. und

Danach wird der Flächeninhalt unter jedem der n Parabelbögen im zugehörigen Teilintervall

durch Integration ermittelt. Die Summe dieser Flächeninhalte ist dann ein (im allgemeinen

sehr guter) Näherungswert für das gesuchte Integral.

Eine formelmäßige Darstellung dieser Vorgangsweise wurde 1743 von Thomas Simpson ² veröffentlicht. Zu ihrer Ableitung berechnen wir den Flächeninhalt unter dem k-ten

Parabelstück allgemein:

Wir setzten zur Abkürzung:

Dann gilt:

und wir erhalten:

Durch Summation von

LBL A

INPUT A

STO G INPUT B INPUT N 1 x=y? GTO A RCL B RCL-A RCL÷N 2 ÷ STO D

LBL B

T (statt X aber A !) E R M STO S LBL C RCL D STO+A T (statt X aber A !)

E R M 4 * STO+S RCL D STO+A T (statt X aber A !) E R M 2 * STO+S RCL A RCL+D RCL+D RCL-B RND x<0? GTO C

LBL D

RCL D

STO+A T (statt X aber A !) E R M 4 * STO+S RCL D

STO+A T (statt X aber A !) E R M STO+S RCL S RCL*D 3 ÷ STO F VIEW F GTO A

Die Fläche soll wieder zwischen berechnet werden.

1 englischer Mathematiker, Physiker und Astronom; * 4.1.1643, _ 13.3.1727

2 englischer Mathematiker; * 1710, _ 1761