Metoda tečen

Metoda tečen je iterační numerická metoda užívaná v numerické matematice k nalezení kořenů funkce nebo k řešení soustavy nelineárních algebraických rovnic. Nazývá se také Newtonova metoda (nebo Newton-Raphsonova metoda) a metodou tečen je označována, protože přesnější řešení rovnice f(x) = 0 je hledáno ve směru tečny funkce f(x).

Popis algoritmu

Newtonova metoda tečen slouží k nalezení řešení rovnice $f\left(x\right)=0$ za předpokladu, že známe derivaci funkce $f^{\prime }\left(x\right)$, tedy směrnici tečny. Pro jednoduchost dále předpokládejme, že $x$ i $f\left(x\right)$ jsou skaláry.

Dalším nezbytným předpokladem je znalost počáteční hodnoty $x_0$, v jejíž blízkosti hledáme řešení. Pokud se funkce $f\left(x\right)$ chová rozumně (je spojitá, hladká a monotónní v intervalu, ve kterém hledáme řešení), lze očekávat řešení v místě, kde tečna sestrojená z bodu $f\left(x_0\right)$ protíná osu $x$. (Směrnice této tečny je $f^{\prime }\left(x_0\right)$.) Tento průsečík označíme $x_1$ a vypočteme jej podle následujícího vztahu.

$x_1=x_0-<!--\; -\; -->\frac{f\left(x_0\right)}{f^{\prime }\left(x_0\right)}$

Za splnění výše uvedených předpokladů by měla hodnota $f\left(x_1\right)$ být blíže nule než původní $f\left(x_0\right)$. Stejný postup můžeme opakovat a najít tak ještě přesnější hodnotu $x_k$.

$x_{k+1}=x_k-<!--\; -\; -->\frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}$

Iteraci provádíme tak dlouho, dokud hodnota $f\left(x_0\right)$ neleží dostatečně blízko nuly. Pokud metoda konverguje a pokud kořen není násobný, je konvergence velice rychlá. S každou iterací se přibližně zdvojnásobí počet desetinných míst, která jsou správně (přesněji, konvergence je kvadratická). Častým kriteriem pro ukončení výpočtu proto bývá okamžik, kdy vzdálenost dvou po sobě jdoucích iterací je dostatečně malá.^[zdroj?]

Odvození iteračního vzorce

Víme, že rovnice tečny t_k funkce f(x) v libovolném bodě x_k je dána vzorcem, jelikož derivace je směrnicí tečny (koeficient lineárního členu):

$t_k:y=f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x+b$

Zároveň také víme, že tečna t_k prochází bodem $\left[x_k,f\left(x_k\right)\right]$. Po dosazení:

$f\left(x_k\right)=f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x_k+b$

$b=f\left(x_k\right)-<!--\; -\; -->f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x_k$

Hledáme bod $\left[x_{k+1},0\right]$, což je průnik tečny t_k v bodě x_k s osou x, a zároveň dosadíme vyjádřené b:

$0=f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x_{k+1}+f\left(x_k\right)-<!--\; -\; -->f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x_k$

$f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x_{k+1}=f^{\prime }\left(x_k\right)\cdot <!--\; \cdot \; -->x_k-<!--\; -\; -->f\left(x_k\right)$

Odtud již plyne:

$x_{k+1}=x_k-<!--\; -\; -->\frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}$

Příklad: Výpočet druhé odmocniny

Úkolem je vypočítat druhou odmocninu kladného reálného čísla a.

$x=\sqrt{a}$

Problém lze definovat také jako nalezení kořenu funkce $f\left(x\right)=x^2-<!--\; -\; -->a$, neboli řešení rovnice $f\left(x\right)=0$.

Vypočteme derivaci $f^{\prime }\left(x\right)$.

$f^{\prime }\left(x\right)=2x$

Dosadíme do obecného vzorce a upravíme.

$x_{k+1}=x_k-<!--\; -\; -->\frac{f\left(x_k\right)}{f^{\prime }\left(x_k\right)}=x_k-<!--\; -\; -->\frac{x_k^2-<!--\; -\; -->a}{2x_k}$

$x_{k+1}=\frac{1}{2}\left(x_k+\frac{a}{x_k}\right)$

Získáváme tak rekurentní rovnici, u které jako počáteční podmínku můžeme zvolit $x_0=a$.

Výpočet $\sqrt{9}$ (druhé odmocniny z devíti) bude podle výše uvedeného algoritmu probíhat následovně.

a  = 9
x0 = 9
x1 = 5
x2 = 3.4
x3 = 3.02352941176471
x4 = 3.00009155413138
x5 = 3.00000000139698
x6 = 3.00000000000000
x7 = 3.00000000000000

Je vidět, že po několika málo krocích se hodnota $x_k$ nemění a ustálí se (konverguje) na hodnotě 3, což odpovídá správnému výsledku.

Poznámky

Aproximace derivace

Pokud známe pouze funkci f(x) a neznáme její derivaci f'(x), můžeme se pokusit derivaci nahradit numerickou derivací. Případně je možné řešit úlohu metodou sečen, která znalost derivace nevyžaduje.

Vektory

Je-li funkce f(x) skalární funkcí vektorového argumentu („z vektoru vypočte skalár“), je nutné hledat x_k+1 proti směru gradientu. Předpis pro iteraci lze potom napsat následovně.

$x_{k+1}=x_k-<!--\; -\; -->\mathrm{\Delta <!--\; \Delta \; -->}{x}_k$

$\mathrm{\Delta <!--\; \Delta \; -->}x={\left[\begin{array}{c}\frac{f\left(x\right)}{\frac{\mathrm{\partial <!--\; \partial \; -->}f}{\mathrm{\partial <!--\; \partial \; -->}{x}_1}},\dots <!--\; \dots \; -->,\frac{f\left(x\right)}{\frac{\mathrm{\partial <!--\; \partial \; -->}f}{\mathrm{\partial <!--\; \partial \; -->}{x}_n}}\end{array}\right]}^T$

Pokud je funkce f(x) vektorovou funkcí vektorového argumentu („z vektoru vypočte vektor“), lze předpis pro iteraci napsat následovně.

$x_{k+1}=x_k-<!--\; -\; -->J\left(x_k{)}^{-<!--\; -\; -->1}f\right(x_k)$

Matice J je takzvaná Jacobiho matice obsahující parciální derivace.

Související články

• Fraktál Newton – fraktál generovaný Newtonovou metodou
• Linearizace – jedním ze způsobů linearizace je nahrazení části křivky její tečnou

Externí odkazy

• Obrázky, zvuky či videa k tématu metoda tečen na Wikimedia Commons
• Newtonova metoda: http://math.fce.vutbr.cz/vyuka/matematika/numericke_metody/node10.html
• Řešení rovnic, Newtonova metoda: https://web.archive.org/web/20071031001614/http://vydra.troja.mff.cuni.cz/bobo/fyzika/num3_metodanewtonova.cz

Autoritní data	LCCN: sh92005466 NLI: 987007536791105171