next up previous
Next: Newton-Raphsonova (tečnová) metoda Up: Řešení jedné nelineární rovnice Previous: Metody, užívající sečnu

Brentova metoda

Metody je založena na přepínání mezi lineární metodou (metodou půlení intervalů) a superlineární metodou (inverzní kvadratická interpolace). Pokud je superlineární metoda pomalá (daleko od kořene), využívá se půlení intervalů.

Inverzní kvadratická interpolace využívá funkci $x = g(y)$, hledáme $x = g(y=0)$. Při iteraci z 3 známých bodů $a,\; b,\; c$, je funkce $y$ interpolována podle Lagrangeova vzorce

$\displaystyle x$ $\textstyle =$ $\displaystyle \frac{[y - f(a)][y - f(b)]c}{[f(c) - f(a)][f(c) - f(b)]} +$  
  $\textstyle +$ $\displaystyle \frac{[y - f(b)][y - f(c)]a}{[f(a) - f(b)][f(a) - f(c)]} +
\frac{[y - f(c)][y - f(a)]b}{[f(b) - f(c)][f(b) - f(a)]}\ \ .$  

Pro $y \equiv 0$ lze Lagrangeův vzorec napsat ve tvaru
$\displaystyle x = b + \frac{P}{Q}\ , \quad {\rm kde} \quad
P$ $\textstyle =$ $\displaystyle S[T(R-T)(c-b) - (1-R)(b-a)] \quad$  
$\displaystyle {\rm a} \quad Q$ $\textstyle =$ $\displaystyle (T-1)(R-1)(S-1)\ .$  


\begin{displaymath}
\!\!\!\!\!\!\!\!\! {\rm a\ \ kde}\qquad R \equiv \frac{f(b)}...
...quiv \frac{f(b)}{f(a)}\ ,
\qquad T \equiv \frac{f(a)}{f(c)}\ .
\end{displaymath}



Jiri Limpouch
2000-04-04