Taylorov rad

Modrá krivka označuje Taylorov polynóm. V animácii sa stupeň polynómu postupne zväčšuje, čim sa aproximácia v okolí nuly spresňuje.

Taylorov rad funkcie $f$ premennej $x$ v bode $a$ je potenčný rad (mocninový rad) so stredom $a$ tvaru

$\sum_{n=0}^{\infin} \frac{f^{(n)}(a)}{n!}(x-a)^{n}$

pričom

• $f^{(n)}(a)$ je $n$-tá derivácia funkcie $f$ v bode $a$
• $f$ má v okolí bodu $a$ derivácie všetkých rádov

Funkcia sa nazýva analytická (v bode a) ak jej Taylorov rad sa v niektorom okolí bodu a zhoduje s danou funkciou. Toto nemusí byť automaticky tak, t. j. existujú funkcie, ktorých Taylorove rady sa s nimi nezhodujú. Príkladom takej funkcie je

$0\not=x\mapsto e^{-\frac{1}{|x|}},\ 0\mapsto 0$

Jej Taylorov rad so stredom v bode 0 je nulový rad pričom daná funkcia má hodnotu nula jedine keď jej argument je nula.

Taylorov rozvoj (funkcie f premennej x v bode a) je Taylorov rad, pre ktorý platí, že jeho súčet (teda výsledná hodnota) v okolí bodu a sa rovná f(x). Maclaurinov rad je Taylorov rad so stredom v a=0.

Účel

Mnoho značne zložitých funkcií je ťažké si predstaviť, zobraziť ich, prípadne odhadnúť ich funkčné hodnoty. Taktiež elementárne funkcie, ako napríklad sínus, kosínus, nadobúdajú najmä iracionálne hodnoty, ktoré nemožno presne vyčísliť, niekedy ani odhadnúť. Formálne definovať tieto základné goniometrické funkcie a mnohé iné umožňuje práve Taylorov rad. Napríklad pre funkciu sínus platí odhad v okolí nuly

$\sin x \approx x - \frac{x^3}{3!} + \frac{x^5}{5!} - \frac{x^7}{7!}$

Táto aproximácia je veľmi silná. Pri výpočte hodnôt funkcie sínus v okolí nuly, možno počítať s veľmi malou chybou. Na prelome 17. a 18. storočia sa viacerí matematici pokúšali nahradiť funkciu nejakou jednoduchšou. Za najjednoduchšie sa všeobecne považujú polynomické funkcie, respektíve polynómy. Vybudovanie teórie, ktorá umožňovala aproximáciu funkcií práve polynómami, však vyžadovala poznatky z vyššej matematiky, hlavne diferenciálneho počtu. Teóriu nezávisle od seba budovali Brook Taylor a Colin Maclaurin. Táto teória umožňuje zapísať, za určitých predpokladov, funkciu ako súčet nekonečného mocninového radu, ktorý sa nazýva Taylorov rad.

Intuitívne odvodenie

Hlavná myšlienka konštrukcie Taylorovho radu spočíva v rovnosti derivácií dvoch funkcií. Obmedzme sa na polynóm stupňa $n$. Zovšeobecnenie pre polynóm nekonečného stupňa - Taylorov rad bude presnejšie popísané v samotnej definícii. Majme dve funkcie, definované v okolí nejakého bodu ich definičného oboru. Ak sa ich funkčné hodnoty rovnajú a ich derivácie všetkých rádov sú v tomto bode rovnaké, potom možno považovať funkcie za rovnaké. Pre zjednodušenie uvažujme bod $x=0$. Vezmime funkciu $f$ a všeobecný polynóm

$p(x)=a_0+a_1x+a_2x^2+a_3x^3+\cdots+a_nx^n$

Treba však nájsť koeficienty $\{a_k\}_{k=1}^{n}$, aby nastala rovnosť $f(x)=p(x)$, pre $x\in\mathcal{O}(0)$. Jednoducho možno odvodiť, že $a_0=f(0)$. Ďalší koeficient možno osamostatniť derivovaním a dosadením nuly. Koeficient $a_1$ teda vypočítame prvou deriváciou

$f'(x)=a_1+2a_2x+3a_3x^2+\cdots+na_nx^{n-1}$

Po dosadení jednoducho $a_1=f'(0)$. Všeobecne pre $k$-ty člen polynómu platí

$a_k=\frac{f^{(k)}(0)}{k!}$

Týmto spôsobom sa jednoducho nájde tvar hľadaného polynómu. V tomto odvodení, ktoré je veľmi hrubé, nie sú zahrnuté všetky predpoklady pre existenciu takého polynómu. Nie pre každú funkciu existuje Taylorov polynóm, respektíve Taylorov rad. Všetky predpoklady sú zhrnuté spolu so všeobecnou formálnou definíciou v nasledujúcom odseku.

Definícia

Nech $f$ je (n+1)-krát diferencovateľná funkcia v bode $a$, definovaná na okolí $\mathcal{O}(a)$ bodu $a$. Potom platí

$f(x)=f(a)+\frac{f'(a)}{1!}(x-a)+\frac{f''(a)}{2!}(x-a)^2+\cdots+\frac{f^{(n)}(a)}{n!}(x-a)^n+R_{n+1}(x)$

kde výraz $R_{n+1}(x)$ označuje zvyšok. Tento vzorec sa nazýva Taylorov vzorec a jeho prvá časť (teda časť bez zvyšku) sa nazýva (n-tý) Taylorov polynóm alebo Taylorov aproximačný polynóm funkcie $f$ so stredom v bode $a$.

Aby bol polynóm $p(x)$ konečného stupňa $n$ Taylorovým polynómom, musí platíť:

$\lim_{x\to a}\frac{f(x)-p(x)}{(x-a)^n}=0$

Implikácia platí aj opačne, teda ak je uvedená limita nulová, potom $p(x)$ je Taylorovým polynómom funkcie $f$ stupňa $n$.

Ak sa n blíži k nekonečnu, nazýva sa Taylorov aproximačný polynóm Taylorov rad.

Tvary zvyšku

Cauchyho tvar zvyšku

Cauchyho tvar vychádza z Lagrangeovej vety o strednej hodnote diferenciálneho počtu. Oveľa praktickejšie je však jeho zovšeobecnenie - Lagrangeov tvar, ktorý je ľahšie zapamätateľný. Nech číslo $\zeta$ je medzi $x$ a stredom $a$. Teda v prípade $x<a$ je $\zeta\in(x,a)$, v opačnom prípade $x>a$ je $\zeta\in(a,x)$. Potom zvyšok možno zapísať v tvare

$R_{n+1}(x)=\frac{f^{(n+1)}(\zeta)}{n!}(x-\zeta)^n(x-a)$

Lagrangeov tvar zvyšku

Lagrangeov tvar je zovšeobecnením Cauchyho tvaru a využíva Cauchyho vetu o strednej hodnote diferenciálneho počtu. Tento tvar je jednoduchší ako Cauchyho, pretože sa podobá na nasledujúci (n+1)-vý člen rozvoja.

$R_{n+1}(x)=\frac{f^{(n+1)}(\zeta)}{(n+1)!}(x-a)^{n+1}$

Integrálny tvar zvyšku

$R_{n+1}(x)=\int\limits_{a}^{x}\frac{f^{(n+1)}(z)}{n!}(x-z)^{n}\,\mathrm{d}z$

Taylorove rady elementárnych funkcií

V tabuľke sú odvodené Taylorove rady niektorých elementárnych funkcií. Na základe znalosti týchto, možno odvodiť rozvoje iných, zložitejších funkcií.

Funkcia	Taylorov rad	Konvergencia
Exponenciálna funkcia	$\mathrm{e}^{x}=\sum^{\infty}_{n=0}\frac{x^n}{n!}$	$x\in\mathbb{R}$
Goniometrická funkcia sínus	$\sin x=\sum^{\infty}_{n=0}\frac{(-1)^nx^{2n+1}}{(2n+1)!}$	$x\in\mathbb{R}$
Goniometrická funkcia cosínus	$\cos x=\sum^{\infty}_{n=0}\frac{(-1)^nx^{2n}}{(2n)!}$	$x\in\mathbb{R}$
Cyklometrická funkcia arkussínus	$\arcsin x=\sum^{\infty}_{n=0}\frac{(2n)!x^{2n+1}}{4^n(n!)^2(2n+1)}$	$x\in\langle-1;1\rangle$
Cyklometrická funkcia arkustangens	$\arctan x=\sum^{\infty}_{n=0}\frac{(-1)^nx^{2n+1}}{2n+1}$	$x\in\mathbb{R}$
Nekonečný geometrický rad	$\frac{1}{1-x}=\sum^{\infty}_{n=0}x^n$	$|x|<1$
Prirodzený logaritmus	$\ln(1-x)=-\sum^{\infty}_{n=1}\frac{x^n}{n}$ $\ln(1+x)=\sum^{\infty}_{n=1}\frac{(-1)^{n+1}x^n}{n}$	$|x|\leq1,x\ne+1$ $|x|\leq1,x\ne-1$
Zovšeobecnený binomický dvojčlen	$(1+x)^\alpha = \sum_{n=0}^\infty {\alpha \choose n} x^n$	$|x|<1$

Príklad

Nájdime Taylorov rad funkcie $f(x)=\arctan3x-x\ln(1+x^2)$ so stredom v bode $a=0$. Najprv uvážime Taylorove rady funkcií arkustangens a prirodzený logaritmus. Pre arkustangens platí

$\left.\arctan u\right|_{u=3x}=\left.\sum^{\infty}_{n=0}\frac{(-1)^nu^{2n+1}}{2n+1}\right|_{u=3x}=\sum^{\infty}_{n=0}\frac{(-1)^n(3x)^{2n+1}}{2n+1}=\sum^{\infty}_{n=0}\frac{(-1)^n3^{2n+1}}{2n+1}x^{2n+1}$

Podobne odvodíme vzťah pre logaritmus

$\left.x\ln(1-u)\right|_{u=-x^2}=-\left.x\sum^{\infty}_{n=1}\frac{u^n}{n}\right|_{u=-x^2}=-\sum^{\infty}_{n=1}\frac{(-1)^n}{n}x^{2n+1}$

Súčtom týchto radov a úpravou do vhodného tvaru dostaneme Taylorov rad funkcie $f$

$\begin{array}{ll}f(x)&=\sum\limits^{\infty}_{n=0}\frac{(-1)^n3^{2n+1}}{2n+1}x^{2n+1}+\sum\limits^{\infty}_{n=1}\frac{(-1)^n}{n}x^{2n+1}=\\&= 3x+\sum\limits^{\infty}_{n=1}\frac{(-1)^n3^{2n+1}}{2n+1}x^{2n+1}+\sum\limits^{\infty}_{n=1}\frac{(-1)^n}{n}x^{2n+1}=\\&= 3x+\sum\limits^{\infty}_{n=1}(-1)^n\left(\frac{3^{2n+1}}{2n+1}+\frac{1}{n}\right)x^{2n+1} \end{array}$

Použitie

V praxi je využiteľný práve Taylorov polynóm, čo je špeciálny prípad Taylorovho radu, pričom sa vypočíta niekoľko prvých členov. Chybu, ktorej sa pri určovaní funkčnej hodnoty dopustíme, možno odhadnúť pomocou zvyškov. V teórii a pri teoretických dôkazoch je nutné používať rozvoje funkcií cez definíciu Taylorovho radu ako nekonečného súčtu. Ak používame pri výpočte radu len niekoľko $k$ prvých členov, ostatné je zvykom zapísať v tvare $o(x^k)$. Tento zápis znamená napríklad pre $p(x)=x+x^2+x^3+o(x^3)$, že ide o členy stupňa vyššieho, ako $3$. Zároveň platí, že

$\lim_{x\to0}\frac{o(x^n)}{x^n}=0$

Súčty radov

Nájdime súčet nekonečného číselného radu

$\sum_{n=0}^{\infty}\frac{(-1)^n\pi^{2n}}{4^{n}(2n+1)!}$

Na výpočet použijeme rozvoj funkcie sínus, pričom zvolíme $\frac{\pi^{2n}}{4^n}=\left(\frac{\pi}{2}\right)^{2n}=x^{2n}$

$\left.\sum_{n=0}^{\infty}\frac{(-1)^nx^{2n}}{(2n+1)!}\right|_{x=\frac{\pi}{2}}=\left.\frac{1}{x}\sum_{n=0}^{\infty}\frac{(-1)^nx^{2n+1}}{(2n+1)!}\right|_{x=\frac{\pi}{2}}=\left.\frac{1}{x}\sin x\right|_{x=\frac{\pi}{2}}=\frac{2}{\pi}$

Výpočet niektorých limít

Niektoré limity nie je možné vypočítať bežnými prostriedkami, prípadne sú príliš zdĺhavé. Napríklad použitie l'Hospitalovho pravidla nezaručuje vždy jednoduchý postup. Možno jednoducho vypočítať mnohé známe limity

$\lim_{x\to0}\frac{\mathrm{e}^x-1}{x}=\lim_{x\to0}\frac{1+x+o(x)-1}{x}=\lim_{x\to0}\left(1+\frac{o(x)}{x}\right)=1$

$\lim_{x\to0}\frac{x-\sin x}{x^3}=\lim_{x\to0}\frac{x-x+\dfrac{x^3}{3!}+o(x^3)}{x^3}=\frac{1}{6}$

Maticová exponenciála

Pomocou Taylorovho radu možno definovať aj abstraktný pojem z algebry - maticová exponenciála. Táto exponenciála sa využíva napríklad pri riešení systémov diferenciálnych rovníc. Uvažujme systém obyčajných diferenciálnych rovníc, ktorý možno napísať v tvare

$\mathbf{A}\mathbf{x}=\frac{\mathrm{d}\mathbf{x}}{\mathrm{d}t}$

Podobne ako obyčajná diferenciálna rovnica prvého rádu, má aj táto maticová diferenciálna rovnica riešenie v tvare

$\mathbf{x}=\mathrm{e}^{\mathbf{A}t}\mathbf{x}_0$

kde $\mathbf{x}_0$ je vektor s počiatočnými podmienkami pre partikulárne riešenie. Definovať maticovú exponenciálu umožňuje teoretický poznatok o diagonalizácii matice a Taylorovom rade. Exponenciálnu funkciu $\mathrm{e}^{\mathbf{A}t}$, bez ohľadu na to, aký objekt predstavuje $\mathrm{A}t$ (v tomto prípade matica) môžeme zapísať

$\mathrm{e}^{\mathbf{A}t}=\mathbf{I}+\mathbf{A}t+\frac{1}{2!}\mathbf{A}^2t^2+\frac{1}{3!}\mathbf{A}^3t^3+\cdots$

Z diagonalizácie matice platí

$\begin{array}{ll}\mathrm{e}^{\mathbf{A}t}&=\mathbf{V}\mathbf{V}^{-1}+\mathbf{V}\Lambda\mathbf{V}^{-1}t+\dfrac{1}{2!}\mathbf{V}\Lambda^2\mathbf{V}^{-1}t^2+\dfrac{1}{3!}\mathbf{V}\Lambda^3\mathbf{V}^{-1}t^3+\cdots=\\&= \mathbf{V}\left(\mathbf{I}+\Lambda t+\dfrac{1}{2!}\Lambda^2t^2+\dfrac{1}{3!}\Lambda^3t^3+\cdots\right)\mathbf{V}^{-1}=\\&= \mathbf{V}\mathrm{e}^{\Lambda t}\mathbf{V}^{-1} \end{array}$

Matica $\mathbf{V}$ je matica vlastných vektorov a matica $\Lambda$ je matica vlastných čísel matice $\mathbf{A}$. Podmienka diagonalizovateľnosti matice je však nutná. V prípade, že matica nie je diagonalizovateľná, je možné použiť na riešenie Jordanov tvar matice a rozklad $\mathbf{A}=\mathbf{MJM}^{-1}$. Pre všeobecné riešenie pôvodnej diferenciálnej sústavy platí

$\mathbf{x}=\mathbf{V}\mathrm{e}^{\Lambda t}\mathbf{V}^{-1}\mathbf{x}_0$

$\mathrm{e}^{\Lambda t}=\begin{bmatrix}\mathrm{e}^{\lambda_1t}&0&0&\cdots&0\\0&\mathrm{e}^{\lambda_2t}&0&\cdots&0\\0&0&\mathrm{e}^{\lambda_3t}&\cdots&0\\\vdots&\vdots&\vdots&\ddots&\vdots\\0&0&0&\cdots&\mathrm{e}^{\lambda_nt}\end{bmatrix}$

Pozri aj

• Nekonečný číselný rad
• Funkcionálny rad
• Infinitezimálny počet
• Derivácia
• Lagrangeova veta o strednej hodnote
• Cauchyho veta o strednej hodnote
• Fourierov rad