Metoda sečen

Metoda sečen je iterační numerická metoda užívaná v numerické matematice k hledání kořene funkce jedné reálné proměnné (tj. hledání nějakého řešení rovnice $f(x)=0$). K nalezení řešení je obvykle potřeba méně iterací než u metody půlení intervalů. Jako u jiných numerických metod, ani metoda sečen není univerzální a může nastat případ, kdy nekonverguje ke správnému řešení. V těchto případech je nutné použít jinou numerickou metodu.

Popis algoritmu

Jak funguje metoda sečen

Na začátku je zapotřebí určit dvě počáteční hodnoty $x_{0}$ a $x_{1}$, které by měly být co nejblíže řešení rovnice $f(x)=0$. Tento kořen nemusí ležet mezi hodnotami $x_{0}$ a $x_{1}$. Pokud však je funkce $f$ spojitá a její funkční hodnoty v bodech $x_{0}$ a $x_{1}$ mají opačná znaménka, nachází se v intervalu $(x_{0},x_{1})$ alespoň jeden kořen. Sečna je v tomto případě lineární interpolací funkce $f$. Nejsou-li tyto podmínky splněny, výpočet metodou sečen může selhat (a funkce ani nemusí mít reálný kořen).

Pomocí bodů $x_{0}$ a $x_{1}$ je dána rovnice sečny ke grafu funkce $f(x)=0$ v bodech $P_{0}=[x_{0},f(x_{0})]$ a $P_{1}=[x_{1},f(x_{1})]$ ve tvaru:

$y=f(x_{1})+{\frac {f(x_{1})-f(x_{0})}{x_{1}-x_{0}}}(x-x_{1})$.

V druhé iteraci musíme určit hodnotu $x_{2}$. Hodnota $x_{2}$ je průsečíkem sečny dané body $x_{0}$ a $x_{1}$ a osy x. Stačí do rovnice sečny dosadit $y=0$ a vyjádřit $x$. Dostaneme rovnici:

$x=x_{1}-{\frac {f(x_{1})(x_{1}-x_{0})}{f(x_{1})-f(x_{0})}}$,

jejímž řešením je hledaná hodnota $x_{2}$. Dále pokračujeme v algoritmu, kde pro $n>0$ se hodnota $x_{n+1}$ vypočítá z hodnot $x_{n}$ a $x_{n-1}$ pomocí vztahu:

$x_{n+1}=x_{n}-{\frac {f(x_{n})(x_{n}-x_{n-1})}{f(x_{n})-f(x_{n-1})}}$.

V algoritmu se pokračuje, dokud není splněna jedna z ukončujících podmínek.

Ukončení výpočtu

Metoda sečen je přibližná metoda. Výpočet se obvykle ukončuje, pokud absolutní hodnota rozdílu hodnot $x_{n}$ a $x_{n-1}$ je menší než požadovaná přesnost výsledku $p>0$ tedy $|x_{n}-x_{n-1}|<p.$

Pokud nejsou splněny podmínky konvergence metody, může se stát, že by výpočet nikdy neskončil. Proto by měl algoritmus obsahovat omezení počtu iterací, jehož dosažení se interpretuje jako selhání metody.

Příklad

Zadání

Najděte kořen rovnice $x^{3}-4x^{2}+6x-24=0$ v intervalu $\langle 3,5\rangle$ a přesností $p=0{,}01$.

Řešení

Za počáteční hodnoty zvolíme $x_{0}=3$ a $x_{1}=5$. Počítáme dosazováním do vzorce: 

$x_{2}=x_{1}-{\frac {f(x_{1})(x_{1}-x_{0})}{f(x_{1})-f(x_{0})}}=5-{\frac {31(5-3)}{31-(-15)}}=3{,}652174$.

Další iterace jsou znázorněny v tabulce:

$n$	$x_{n}$	$f(x_{n})$	$p$
2	3,652174	-6,726391	1,347826
3	3,892483	-2,274132	0,240309
4	4,015229	0,336894	0,122746
5	3,999391	-0,013388	0,015838
6	3,999997	-0,000074	0,000605

Z tabulky je patrné, že u šesté iterace je $p<0{,}01$ a hledaný kořen rovnice $x^{3}-4x^{2}+6x-24=0$ v intervalu $\langle 3,5\rangle$ a přesností $p=0{,}01$ je $x=4$.
(V tabulce je kořen 3,999997, ale s zaokrouhlením na setiny dostáváme 4)

Výhody a nevýhody

Mezi výhody této metody patří malý počet iterací pro nalezení kořene a snadné naprogramování algoritmu v programovacím jazyce. Naopak mezi nevýhody patří, že ne vždy dojde v této metodě ke konvergenci.

Odkazy

Související články

• Metoda tečen
• Metoda Monte Carlo
• Půlení intervalů
• Iterace

Literatura

CHAPRA, Steven C. Numerical Methods for Engineers 6 edition. 1. vyd. [s.l.]: McGraw-Hill, 2009. ISBN 9780073401065. 

VITÁSEK, Emil. Numerické metody. 1. vyd. Praha: SNTL, 1987. 512 s.