Metoda sečen

Z Wikipedie, otevřené encyklopedie
Skočit na: Navigace, Hledání

Metoda sečen je iterační numerická metoda užívaná v numerické matematice k numerickému řešení soustav nelineárních rovnic jedné reálné proměnné. K nalezení řešení je obvykle potřeba méně iterací, než u numerické metody půlení intervalů. Jako u jiných numerických metod, tak i 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[editovat | editovat zdroj]

Jak funguje metoda sečen

Na začátku je zapotřebí určit dvě počáteční hodnoty a , které by měly být co nejblíže řešení rovnice . Tento kořen nemusí ležet mezi hodnotami a .
Pomocí bodů a je dána rovnice sečny ke grafu funkce v bodech a ve tvaru:

.

V druhé iteraci musíme určit hodnotu . Hodnota je průsečíkem sečny dané body a a osy x. Stačí do rovnice sečny dosadit a vyjádřit . Dostaneme rovnici:

,

jejímž řešením je hledaná hodnota . Dále pokračujeme v algoritmu, kde pro se hodnota vypočítá z hodnot a pomocí vztahu:

.

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

Podmínky pro ukončení algoritmu[editovat | editovat zdroj]

Pro ukončení algoritmu musí dojít k naplnění k jedné z těchto podmínek:

  • Rozdíl hodnot a je menší než zvolená přesnost výsledku tedy
  • Proběhne maximální počet iterací.

Počet iterací je na uvážení každého, neexistuje žádná metoda na určení optimálního počtu iterací. Obecně se volí větší počet iterací, ale musíme si dát pozor při použití algoritmu na počítači, aby nedošlo k jeho přehlcení. Naopak při malém počtu iterací metoda sečen nedokonverguje.

Příklad[editovat | editovat zdroj]

Zadání[editovat | editovat zdroj]

Najděte kořen rovnice  v intervalu  a přesností .

Řešení[editovat | editovat zdroj]

Za počáteční hodnoty zvolíme  a . Počítáme dosazováním do vzorce: 

.

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

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 šešté iterace je a hledaný kořen rovnice v intervalu a přesností je .
(V tabulce je kořen 3.999997, ale s zaokrouhlením na setiny dostáváme 4)

Výhody a nevýhody[editovat | editovat zdroj]

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.

Literatura[editovat | editovat zdroj]

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.  

Související články[editovat | editovat zdroj]