Taylorov rad

z Wikipédie, slobodnej encyklopédie
Prejsť na: navigácia, hľadanie
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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

Cauchyho tvar zvyšku[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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í[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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[upraviť | upraviť zdroj]

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}

Referencie[upraviť | upraviť zdroj]


Pozri aj[upraviť | upraviť zdroj]