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)

0x01 graphic
, jeśli f(xk)·f(xk-1)<0

Zakończenie obliczeń: 0x01 graphic

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:

0x08 graphic

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:

  1. funkcje f(x), f'(x), f”(x) są określone i ciągłe w przedziale [a,b]

  2. f(a)·f(b)<0

  3. 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:

  1. 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.

  1. kolejne przybliżenia obliczamy ze wzoru

0x01 graphic

  1. obliczenia kończymy, jeśli spełniona jest nierówność

0x01 graphic

METODA SIECZNYCH (reguła falsi, metoda interpolacji liniowej)

Algorytm tej metody:

  1. 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.

  1. kolejne przybliżenia obliczamy ze wzoru

0x01 graphic

d = a, gdy x0=b lub d = b, gdy x0=a

  1. obliczenia kończymy, jeśli spełniona jest nierówność

0x01 graphic

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:

  1. wybór funkcji iteracyjnej 0x01 graphic

  2. obliczenia wg schematu:

0x01 graphic
,

gdzie 0x01 graphic

bk-1 jeśli f(bk-1)f(ck-1)<0

bk =

ck-1 jeśli f(bk-1)f(ck-1)≥0