Regula falsi

Regula falsi (łac. fałszywa linia prosta, fałszywa reguła) — algorytm rozwiązywania równań nieliniowych jednej zmiennej.

Plik:Regula falsi pl.svg

Pierwsze dwie iteracje metody regula falsi. Czerwona krzywa pokazuje funkcje f, niebieskie linie są cięciwami.

Na funkcję $y=f(x)$ nakładane są następujące ograniczenia:

1. Jest ciągła w przedziale $[a,b]$.
2. W przedziale znajduje się dokładnie jeden pojedynczy pierwiastek.
3. Na końcach przedziału funkcja ma różne znaki: $f(a)f(b)<0$.
4. Pierwsza i druga pochodna mają na tym przedziale stałe znaki.

Algorytm przebiega następująco:

• Na początku przez punkty $A=(a,f(a))$ i $B=(b,f(b))$ przeprowadzana jest cięciwa.
• Punkt przecięcia $x_{1}$ z osią OX jest brany jako pierwsze przybliżenie pierwiastka.
• Jeśli to przybliżenie jest wystarczająco dobre, algorytm kończy się.
• Jeśli nie, to prowadzona jest cięciwa przez punkty $(x_{1},f(x_{1}))$ oraz $A$ lub $B$ – wybierany jest ten punkt, którego rzędna ma znak przeciwny do $f(x_{1})$. Jednak w praktyce, dzięki ograniczeniu nr 4 już na początku algorytmu wiadomo, który z tych punktów będzie stały, tzn. wybierany za każdym razem.
• Następnie wyznaczane jest przecięcie nowo wyznaczonej cięciwy z osią OX ($x_{i}$) i algorytm powtarza się.

Nazwa metody pochodzi od łacińskich słów: regula¹ znaczące zarówno linię prostą, jak i regułę i falsus, fałszywy — metoda bazuje na fałszywym twierdzeniu (regule), że na pewnym przedziale funkcja jest liniowa. Można więc tę nazwę przetłumaczyć zarówno jako "fałszywa linia prosta" jak i "fałszywa reguła" i odbydwa te tłumaczenia mają w tym kontekście sens.

Wzory

$x_{1}={\frac {af(b)-bf(a)}{f(b)-f(a)}}$

$x_{i+1}=\left\{{\begin{matrix}{\frac {x_{i}f(a)-af(x_{i})}{f(a)-f(x_{i})}}&gdy&f(a)f(x_{i})\leq 0\\\\{\frac {x_{i}f(b)-bf(x_{i})}{f(b)-f(x_{i})}}&gdy&f(b)f(x_{i})<0\end{matrix}}\right.$

dla $i=1,2,...$

Inne numeryczne metody wyznaczania pierwiastków równania nieliniowego:

• metoda bisekcji
• metoda stycznych
• metoda siecznych
• Metoda Newtona
• algorytm Illinois (zmodyfikowana metoda siecznych)

Linki zewnętrzne

• http://www.algorytm.org/index.php?option=com_content&task=view&id=79&Itemid=28
• http://www.hpcsoft.com/products/MathSoL/nonlinear/falseposition.html (ang.)