METODY ROZWIĄZYWANIA RÓWNAŃ NIELINIOWYCH
Rozwiązanie równania f(x)=0 jest równoważne lokalizacji miejsca zerowego, x*, funkcji f(x);
Jeśli f(x) jest ciągła i f(a)·f(b)<0 to x* ∈ [a,b].
METODA BISEKCJI (połowienia przedziału)
, jeśli f(xk)·f(xk-1)<0
Zakończenie obliczeń:
iteracja |
a |
b |
f(a) |
f(b) |
c=(a+b)/2 |
f(c) |
|
koniec obliczeń |
1 2 . k |
|
|
|
|
xk |
f(xk) |
|
nie . . tak |
Przedział (a,b) w k-tej iteracji wyznaczamy za pomocą wzorów:
ak-1 jeśli f(ak-1)f(ck-1)<0
ak =
ck-1 jeśli f(ak-1)f(ck-1)≥0
METODA NEWTONA (stycznych)
Warunki, które musi spełniać funkcja:
funkcje f(x), f'(x), f”(x) są określone i ciągłe w przedziale [a,b]
f(a)·f(b)<0
f'(x) I f”(x) nie są równe zeru I nie zmieniają znaku w domkniętym przedziale [a,b] , tzn. w rozważanym przedziale funkcja f(x) jest ściśle monotoniczna oraz wklęsła lub wypukła.
Algorytm metody stycznych:
w charakterze przybliżenia początkowego, x0, przyjmujemy ten koniec a lub b przedziału, który spełnia warunek: f(x0)·f”(x0)>0.
kolejne przybliżenia obliczamy ze wzoru
obliczenia kończymy, jeśli spełniona jest nierówność
METODA SIECZNYCH (reguła falsi, metoda interpolacji liniowej)
Algorytm tej metody:
w charakterze przybliżenia początkowego x0 przyjmujemy ten koniec a lub b przedziału, który spełnia warunek: f(x0)·f”(x0)<0.
kolejne przybliżenia obliczamy ze wzoru
d = a, gdy x0=b lub d = b, gdy x0=a
obliczenia kończymy, jeśli spełniona jest nierówność
METODA ITERACJI PROSTEJ
Równanie f(x) = 0
należy przekształcić do równoważnego równania x = g(x)
Znajdowanie pierwiastka jest realizowane w dwóch etapach:
wybór funkcji iteracyjnej
obliczenia wg schematu:
przybliżenie zerowe x0 = (a+b)/2
kolejne przybliżenia xk+1 = g(xk)
zakończenie obliczeń:
,
gdzie
bk-1 jeśli f(bk-1)f(ck-1)<0
bk =
ck-1 jeśli f(bk-1)f(ck-1)≥0