Metoda Newtona

Metoda Newtona

Rodzaj	algorytm iteracyjny, wyznaczanie przybliżonej wartości pierwiastka funkcji

Metoda Newtona (zwana również metodą Newtona-Raphsona lub metodą stycznych) – algorytm iteracyjny wyznaczania przybliżonej wartości pierwiastka funkcji^[1].

Rozwiązywanie równania nieliniowego

Zadanie

Zadaniem jest znalezienie pierwiastka równania zadanej funkcji ciągłej f:

${\displaystyle \mathbb {R} \supset [a,b]\ni x\mapsto f(x)\in \mathbb {R} }$

w przedziale ${\displaystyle [a,b].}$ A zatem znalezienie takiego ${\displaystyle x^{*}\in [a,b],}$ które spełnia następujące równanie:

${\displaystyle f(x^{*})=0.}$

Opis metody

(c) Wojciech mula z polskiej Wikipedii, CC BY-SA 2.5

Ilustracja działania metody Newtona, pokazane zostały 4 pierwsze kroki.

W metodzie Newtona przyjmuje się następujące założenia dla funkcji f:

1. W przedziale [a,b] znajduje się dokładnie jeden pierwiastek.
2. Funkcja ma różne znaki na krańcach przedziału, tj. ${\displaystyle f(a)\cdot f(b)<0.}$
3. Pierwsza i druga pochodna funkcji mają stały znak w tym przedziale.

W pierwszym kroku metody wybierany jest punkt startowy x₁ (zazwyczaj jest to wartość a, b, 0 lub 1), z którego następnie wyprowadzana jest styczna w f(x₁). Odcięta punktu przecięcia stycznej z osią OX jest pierwszym przybliżeniem rozwiązania (ozn. x₂).

Jeśli to przybliżenie nie jest satysfakcjonujące, wówczas punkt x₂ jest wybierany jako nowy punkt startowy i wszystkie czynności są powtarzane. Proces jest kontynuowany, aż zostanie uzyskane wystarczająco dobre przybliżenie pierwiastka

Kolejne przybliżenia są dane rekurencyjnym wzorem:

${\displaystyle x_{k+1}=x_{k}-{\frac {f(x_{k})}{f'(x_{k})}}.}$

Szacowanie błędu

Błąd k-tego przybliżenia można oszacować poprzez nierówności (x* to dokładna wartość pierwiastka):

${\displaystyle |x^{*}-x_{k}|\leqslant {\frac {f(x_{k})}{m}}}$

lub

${\displaystyle |x^{*}-x_{k}|\leqslant {\frac {M}{2m}}(x^{*}-x_{k-1})^{2},}$

gdzie stałe wyznacza się ze wzorów:

${\displaystyle m=\min _{x\in [a,b]}|f'(x)|}$

oraz

${\displaystyle M=\max _{x\in [a,b]}|f''(x)|.}$

Warunek zakończenia obliczeń

W metodzie Newtona wykonuje się iteracyjne obliczenia, aż do momentu gdy ich wyniki będą satysfakcjonujące. W praktyce stosowanych jest kilka kryteriów warunków zakończenia obliczeń dla algorytmu (${\displaystyle \epsilon }$ to przyjęta dokładność obliczeń):

1. wartość funkcji w wyznaczonym punkcie jest bliska 0:

${\displaystyle \left|f(x_{k})\right|\leqslant \epsilon }$

2. odległość pomiędzy kolejnymi przybliżeniami jest dość mała:

${\displaystyle \left|x_{k+1}-x_{k}\right|\leqslant \epsilon }$

3. szacowany błąd jest dostatecznie mały:

${\displaystyle {\frac {M}{2m}}(x_{k}-x_{k-1})^{2}\leqslant \epsilon }$

4. kryterium mieszane (punkty 1 i 2 jednocześnie)

Zbieżność

Metoda Newtona-Raphsona jest metodą o zbieżności kwadratowej – rząd zbieżności wynosi 2 (wyjątkiem są zera wielokrotne, dla których zbieżność jest liniowa i wynosi 1), zaś współczynnik zbieżności ${\displaystyle {\frac {M}{2m}}.}$ Oznacza to, iż przy spełnionych założeniach błąd maleje kwadratowo wraz z ilością iteracji.

Metoda Newtona jest metodą rozwiązywania równań często używaną w solverach, ze względu na jej szybką zbieżność (w algorytmie liczba cyfr znaczących w kolejnych przybliżeniach podwaja się). Wadą jej jest fakt, iż zbieżność nie musi zawsze zachodzić. W wielu przypadkach metoda bywa rozbieżna, kiedy punkt startowy jest zbyt daleko od szukanego pierwiastka równania.

Profesjonalne solvery wykorzystują stabilniejsze, lecz mniej wydajne metody (jak np. metoda bisekcji) do znalezienia obszarów zbieżności w dziedzinie funkcji, a następnie używają metody Newtona-Raphsona do szybkiego i dokładniejszego obliczenia lokalnego pierwiastka równania. Dodatkowo solvery posiadają zabezpieczenia przed nadmierną ilością wykonywanych iteracji (przekroczenie ustalonej liczby iteracji jest równoznaczne z niepowodzeniem algorytmu w zadanym przedziale).

Przykład

Za pomocą metody Newtona można obliczyć pierwiastek ${\displaystyle {\sqrt {a}}}$ dla każdej liczby ${\displaystyle a\in \mathbb {R} ^{+}{:}}$

${\displaystyle {\sqrt {a}}=x\iff a=x^{2}\iff x^{2}-a=0.}$

Funkcja f(x) ma postać:

${\displaystyle f(x)=x^{2}-a,}$

${\displaystyle f\,'(x)=2x.}$

Rekurencyjny wzór wynosi:

${\displaystyle x_{k+1}=x_{k}-{\frac {x_{k}^{2}-a}{2x_{k}}},}$

${\displaystyle x_{k+1}={\frac {1}{2}}\left(x_{k}+{\frac {a}{x_{k}}}\right).}$

Dla danych ${\displaystyle a=2}$ i ${\displaystyle x_{0}=1}$ algorytm przebiega następująco:

${\displaystyle x_{0}=1,}$

${\displaystyle x_{1}={\frac {1}{2}}(1+{\frac {2}{1}})=1{,}5,}$

${\displaystyle x_{2}={\frac {1}{2}}(1{,}5+{\frac {2}{1{,}5}})\approx 1{,}416666,}$

${\displaystyle x_{3}={\frac {1}{2}}(1{,}416666+{\frac {2}{1{,}416666}})\approx 1{,}414214.}$

Rozwiązywanie układu równań nieliniowych

Przykład użycia metody Newtona do rozwiązywania układu równań nieliniowych

Metodę Newtona można zgeneralizować do przypadku wielowymiarowego i użyć jej do rozwiązywania układów równań nieliniowych.

Zadanie

Niech U będzie otwartym podzbiorem przestrzeni ${\displaystyle \mathbb {R} ^{n}}$ oraz ${\displaystyle F\colon U\to \mathbb {R} ^{n}}$ będzie funkcją różniczkowalną.

Zadaniem uogólnionej metody Newtona jest znalezienie takiej wartości x*, dla której:

${\displaystyle F(\mathbf {x^{*}} )=\mathbf {0} .}$

Opis metody

Algorytm, podobnie jak dla przypadku jednowymiarowego, polega na wyborze punktu startowego ${\displaystyle \mathbf {x_{0}} }$ (często wybiera się wektor zerowy lub wektor jedynek), a następnie rekurencyjnym przekształcaniu tego wektora aż do momentu, gdy kolejne przybliżenia będą satysfakcjonujące. Wektory przekształcane są zgodnie z równaniem macierzowym:

${\displaystyle \mathbf {x_{k+1}} =\mathbf {x_{k}} -\left(F'(\mathbf {x_{k}} )\right)^{-1}F(\mathbf {x_{k}} ),}$

gdzie ${\displaystyle F'}$ jest pochodną (Frécheta) – jest to de facto macierz wielkości ${\displaystyle n\times n.}$

Przy implementacji metody, zamiast odwracania macierzy ${\displaystyle F',}$ efektywniej jest rozwiązać układ równań (tożsamy z powyższym równaniem):

${\displaystyle F'(\mathbf {x_{k}} )\mathbf {d} =-F(\mathbf {x_{k}} ),}$

a następnie na podstawie obliczonego d wyznaczyć kolejne przybliżenie:

${\displaystyle \mathbf {x_{k+1}} =\mathbf {x_{k}} +\mathbf {d} .}$

Warunek zakończenia obliczeń

Kryterium zakończenia obliczeń podobnie jak w metodzie jednowymiarowej może być (w zadanej normie ${\displaystyle \|\cdot \|}$ oraz dokładności ${\displaystyle \epsilon }$):

1. wartość funkcji dostatecznie bliska wektorowi zerowemu:

   ${\displaystyle \|F(\mathbf {x_{k}} )\|\leqslant \epsilon }$

2. dostatecznie mała odległość pomiędzy kolejnymi punktami w iteracji:

   ${\displaystyle \|\mathbf {x_{k+1}} -\mathbf {x_{k}} \|\leqslant \epsilon }$

3. kryterium mieszane (punkty 1 oraz 2 jednocześnie)

Zbieżność

Jeśli funkcja F:

• ${\displaystyle F(\mathbf {x^{*}} )=\mathbf {0} }$ dla pewnego ${\displaystyle \mathbf {x^{*}} \in U,}$
• pochodna ${\displaystyle F'}$ jest odwzorowaniem zwężającym i nieosobliwym,

to dla punktu startowego ${\displaystyle \mathbf {x^{0}} }$ będącego dostatecznie blisko x*, wielowymiarowa metoda Newtona jest zbieżna oraz zbieżność ta jest kwadratowa.

Pierwiastki wielokrotne

Przy rozwiązywaniu równań nieliniowych kłopotliwymi dla metody Newtona mogą być pierwiastki wielokrotne, dla których zbieżność algorytmu staje się liniowa. Dla takich przypadków metoda Newtona może być dużo wolniejsza niż inne metody rozwiązywania równań o zbieżności liniowej.

Aby zaradzić tego typu problemom, w praktyce stosuje się następujące podejścia:

• Dla układu równań – przeprowadzenie optymalizacji funkcji G (znalezienie minimum zadanej funkcji celu):

${\displaystyle \mathbf {G} (\mathbf {x} )=\langle F(\mathbf {x} ),F(\mathbf {x} )\rangle \in \mathbb {R} ,}$

gdzie ${\displaystyle \langle \cdot ,\cdot \rangle }$ oznacza iloczyn skalarny dwóch wektorów.

• Dla równania nieliniowego – znalezienie pierwiastka odpowiedniej pochodnej ${\displaystyle f}$ lub przeprowadzenie minimalizacji funkcji ${\displaystyle \left(f(x)\right)^{2}}$

Zobacz też

Inne metody rozwiązywania równań nieliniowych:

• metoda bisekcji
• metoda Halleya
• metoda siecznych
• odwrotna interpolacja kwadratowa
• Joseph Raphson

Przypisy

Bibliografia

• Zenon Fortuna, Bohdan Macukow, Janusz Wąsowski, Metody numeryczne, Wydawnictwa Naukowo-Techniczne, Warszawa 2006, ISBN 978-83-7926-281-6.

Linki zewnętrzne

• Eric W.E.W. Weisstein Eric W.E.W., Newton's Method, [w:] MathWorld [online], Wolfram Research [dostęp 2020-12-12]  (ang.).
• Równania nieliniowe na wazniak.mimuw.edu.pl