Numerik FB 3 — Mathematisches Institut Prof. Dr. Thomas Götz Computerpraktikum 9 10. Januar 2017 Aufgabe 22: Schreiben Sie eine FreeMat–Funktion function [x, y] = newtR1(f, df, x0, TOL) die die Nullstelle der Funktion f ausgehend von Startwert x0 mit dem Newtonverfahren berechnet. Als Abbruchkriterium soll hierbei |f (x)| <TOL bzw. maximal 100 Iterationen verwendet werden. Die Ableitung der Funktion f wird als df übergeben. Im Vektor y speichern Sie die einzelnen Glieder der Newtonfolge. Testen Sie Ihre Funktion für die beiden Probleme 1. f (x) = x3 − 2x + 2 und x0 = 0. x 1 2. = 2 mit dem Anfangswert x0 = −2. Die exakte Lösung ist x = −1. 1 + x2 Als Toleranz können Sie jeweils TOL=1E-5 verwenden. Überprüfen Sie die Konvergenzordung. Aufgabe 23: Schreiben Sie eine FreeMat-Funktion function [x,y] = newtRn(f, df, x0, TOL) die die Nullstelle der Funktion f : Rn → Rn mit Jacobimatrix df : Rn → Rn×n ausgehend von Startpunkt x0 ∈ Rn mit dem Newtonverfahren berechnet. Als Abbruchkriterium soll hierbei kf (x)k <TOL bzw. maximal 100 Iterationen verwendet werden. Im Vektor y speichern Sie die einzelnen Glieder der Newtonfolge. Testen Sie Ihre Funktion für das Problem sin2 x − y = 0 x + y2 − 1 = 0 mit Startpunkt (x0 , y0 ) = (0, 0). Besprechung in der 3. Woche 1/1
© Copyright 2024 ExpyDoc