Kapitola 5
Derivace složené funkce, Taylorův vzorec

Stejně jako u funkce jedné proměnné potřebujeme u funkcí více proměnných určit parciální derivace složené funkce. To je obsahem prvního odstavce, kde také ukážeme použití odvozených vzorců. Druhý odstavec této kapitoly je věnován Taylorovu vzorci pro funkci více proměnných. Podrobnější srovnání s funkcí jedné proměnné provedeme v každém odstavci zvlášť.

5.1 Parciální derivace složených funkcí

Vzorce pro parciální derivace složených funkcí jsou jedním z nejdůležitějších nástrojů řešení rovnic matematické fyziky. Tyto rovnice jsou tzv. parciální diferenciální rovnice – to jsou rovnice, které obsahují parciální derivace neznámé funkce a jejichž řešení jsou funkce dvou či více proměnných. Odvozené vzorce umožňují transformovat tyto rovnice na jednodušší tvar, z něhož buď již umíme najít řešení, nebo alespoň můžeme vyvodit řadu důležitých vlastností řešení rovnice.

Na úvod připomeňme, jak se derivuje složená funkce jedné proměnné. Nechť funkce u = g(x) má derivaci v bodě x0. Označme u0 = g(x0). Má-li funkce y = f(u) derivaci v bodě u0, pak složená funkce y = F(x) = f(g(x)) má derivaci v bodě x0 a platí: y′(x0) = f′(u0)g′(x0).

Nyní odvodíme podobné vztahy pro parciální derivace složené funkce dvou proměnných. Bude nás především zajímat případ, kdy vnější funkce f není explicitně zadána (obvykle je to hledané řešení parciální diferenciální rovnice).

Věta 5.1. Nechť funkce u = u(x,y), v  = v(x,y) mají parciální derivace prvního řádu v bodě [x0,y0], označme u0 = u(x0,y0), v0 = v(x0,y0). Je-li funkce z = f(u,v) diferencovatelná v bodě [u0,v0], pak složená funkce z = F(x,y) = f(u(x,y),v(x,y)) má parciální derivace 1. řádu v bodě [x0,y0] a platí:

 ∂F- ∂f- ∂u- ∂f- ∂v- ∂x (x0,y0) = ∂u (u0,v0)∂x(x0,y0)+ ∂v(u0,v0)∂x (x0,y0) ∂F ∂f ∂u ∂f ∂v ---(x0,y0) = ---(u0,v0)---(x0,y0)+ ---(u0,v0)--(x0,y0). ∂y ∂u ∂y ∂v ∂y
(5.1)

Zkráceně píšeme

zx = zuux + zvvx, zy = zuuy + zvvy
(5.2)

nebo také

∂z-= ∂z-∂u-+ ∂z∂v-, ∂z-= ∂z-∂u-+ ∂z-∂v. ∂x ∂u ∂x ∂v∂x ∂y ∂u ∂y ∂v ∂y
(5.3)

 

 

Důkaz. Dokážeme pouze první vzorec v (5.1 ), druhý se dokáže zcela analogicky. Vyjdeme přímo z definice parciální derivace.

 ∂F- F-(x0 +-t,y0)−-F-(x0,y0) ∂x (x0,y0) = lit→m0 t = f (u(x + t,y ),v(x + t,y )) − f(u(x ,y ),v(x ,y )) = lim ------0-----0-----0-----0---------0--0-----0-0--. t→0 t
(5.4)

 

Označíme-li u(t) = u(x0 + t,y0), v(t) = v(x0 + t,y0), z diferencovatelnosti funkce f plyne existence funkce τ splňující (4.3) takové, že

 f(u(t),v(t))− f (u0,v0) = = fu(u0,v0)(u(t)− u0) + fv(u0,v0)(v(t)− v0)+ τ (u(t)− u0, v(t)− v0).
Dosazením tohoto vztahu do (5.4 ) dostáváme
∂F- ∂x(x0,y0) = limt→01- t [fu(u0,v0)(u(t)− u0) + fv(u0,v0)(v(t)− v0)+
+ τ(u(t) − u0,v(t) − v0)] = fu(u0,v0)limt→0u(x0 + t,y0)− u(x0,y0) ----------t-----------+
+ fv(u0,v0)limt→0v(x + t,y )− v(x ,y ) ---0-----0------0--0-- t + limt→0τ(u(t)− u ,v(t)− v ) ---------0--------0- t =
= fu(u0,v0)ux(x0,y0) + fv(u0,v0)vx(x0,y0) + limt→0τ(u(t)−-u0,v(t)−-v0) t.
K dokončení důkazu nyní stačí ukázat, že poslední limita je nulová:
limt→0τ(u(t) − u0,v(t) − v0) ---------t---------- = limt→0 τ(u(t)− u0,v(t)− v0) ∘-----------2------------2 (u(t)− u0) + (v(t)− v0)⋅
⋅ ∘ ---------------------------- (u(t) − u0)2 (v(t) − v0)2 --------- + -------- t t =
= ∘ -2----------2------- ux(x0,y0) + vx(x0,y0) ⋅ limt→0---τ(u(t)-−-u0,v(t)-−-v0)--- ∘(u(t)−--u0)2-+-(v(t)-−-v0)2- = 0.
V posledním výpočtu jsme využili faktu, že limt→0u(t) = u0,limt→0v(t) = v0, neboť funkce u(t) = u(x0 + t,y0),v(t) = v(x0 + t,y0) jsou spojité v bodě t = 0 – to plyne z existencí parciálních derivací funkcí u,v v bodě t = 0 a pro funkci jedné proměnné plyne z existence derivace spojitost.

 

Příklad 5.1. i) Je dána funkce z = eu sinv, kde u = xy a v = x + y. Vypočtěte zx a zy.

Řešení. Protože vnitřní i vnější složky mají spojité parciální derivace v celém ℝ2, má složená funkce parciální derivace v každém bodě tohoto prostoru. Dosazením do (5.2 ) dostáváme

 u u zx = zuux + zvvx = (e sin v)y + (e cosv), zy = zuuy + zvvy = (eu sin v)x + (eucosv).
Zbývá dosadit za u a v, u = xy a v = x + y a dostaneme
 xy xy zx = e (ysin(x+ y) + cos(x + y)), zy = e (x sin(x + y)+ cos(x + y)).

ii) Pomocí transformace do nových nezávisle proměnných u = x + y, v = x y najděte všechny diferencovatelné funkce f : ℝ2 → ℝ splňující rovnost

fx(x,y) + fy(x,y) = 0.
(5.5)

Řešení. Označme z = f(x,y). Pak zx = zuux + zvvx = zu + zv, zy = zuuy + zvvy = zu zv. Dosazením dostáváme zu + zv + zu zv = 2zu = 0, tedy zu = 0. To znamená, že funkce z = z(u,v) nezávisí na proměnné u, a tedy z(u,v) = g(v), kde g je libovolná diferencovatelná funkce jedné proměnné. Dosazením za v vidíme, že všechny diferencovatelné funkce dvou proměnných, které splňují (5.5 ), jsou tvaru f(x,y) = g(xy), kde g je libovolná diferencovatelná funkce jedné proměnné.

iii) Proveďte totéž jako v předchozím příkladě zavedením polárních souřadnic ϕ = arctg yx, r = ∘ ------- x2 + y2 do rovnice

yf (x,y) − xf (x,y) = 0. x y
(5.6)

Řešení. Vypočtěme nejprve parciální derivace funkcí r a ϕ:

 x y rx = ∘--------, ry = ∘--------, x2 + y2 x2 + y2 ϕ = − ---y---, ϕ = ---x---. x x2 + y2 y x2 + y2
Označíme-li opět z = f(x,y) a dosadíme-li do vzorečků pro derivace složené funkce prvního řádu, dostáváme
zx = zr----x---- ∘ x2 + y2 zϕ---y--- x2 + y2,
zy = zr y ∘--2----2 x + y + zϕ x -2----2 x + y,
což po dosazení do (5.6 ) a úpravě dává rovnici zϕ = 0, a tedy z(r,ϕ) = h(r). Všechny funkce dvou proměnných splňující rovnici (5.6) jsou tedy tvaru f(x,y) = h(∘ ------- x2 + y2), kde h je libovolná diferencovatelná funkce jedné proměnné.

 

Na předchozích příkladech vidíme, že zavedením nových nezávisle proměnných můžeme dosáhnout značného zjednodušení dané parciální diferenciální rovnice, což se velmi často využívá především při řešení diferenciálních rovnic popisujících různé fyzikální děje. Protože tyto rovnice jsou většinou 2. řádu (obsahují parciální derivace 2. řádu neznámé funkce), zvláště důležité jsou vzorce pro parciální derivace 2. řádu složených funkcí.

Dříve než si tyto vzorce pro parciální derivace 2. řádu uvedeme, připomeňme opět pro srovnání vzorec pro derivace 2. řádu složené funkce jedné proměnné. Derivováním rovnosti y′ = f′(u(x))g′(x) dostáváme

 ′′ ′ ′ ′ ′′ ′2 ′ ′′ y = (f (u)g (x))= f (u)g (x)+ f(u)g (x).

 

Věta 5.2. Nechť funkce u = u(x,y), v  = v(x,y) mají parciální derivace 2. řádu v bodě [x0,y0], označme u0 = u(x0,y0), v0 = v(x0,y0). Má-li funkce z = f(u,v) spojité parciální derivace 2. řádu v bodě [u0,v0], pak složená funkce z = F(x,y) = f(u(x,y),v(x,y)) má parciální derivace 2. řádu v bodě [x0,y0] a platí:

zxx =zuuu2x + 2zuvuxvx + zvvv2x + zuuxx + zvvxx zxy =zuuuxuy + zuvvyux + zvuuyvx + zvvvyvx + zuuxy + zvvxy zyy =zuuu2y + 2zuvuyvy + zvvv2y + zuuyy + zvvyy.
(5.7)

Funkce z a její parciální derivace mají argument (u0,v0), funkce u,v a jejich parciální derivace mají argument (x0,y0).

Důkaz. Dokážeme pouze rovnost pro zxx, důkaz zbývajících dvou vzorců je zcela analogický. Platí

zxx = ∂-- ∂x(zx) = -∂- ∂x(zuux + zvvx) = -∂- ∂x(zuux) + ∂-- ∂x(zvvx) =
= ∂ --- ∂x(zu)ux + zuuxx + ∂ --- ∂x(zv)vx + zvvxx =
= (zuuux + zuvvx)ux + zuuxx + (zvuux + zvvvx)vx+
+ zvvxx = zuuux2 + zuvvxux + zvvvx2 + zvuuxvx + zuuxx + zvvxx =
= zuuux2 + 2zuvuxvx + zvvvx2 + zuuxx + zvvxx.
K výpočtu ∂ ∂x-zu∂ ∂xzv jsme využili skutečnosti, že zu = zu(u(x,y),v(x,y))zv = zv(u(x,y),v(x,y)) jsou opět složené funkce proměnných x,y, a proto můžeme k výpočtu jejich derivací využít vztahů (5.1 ), ve kterých místo z dosadíme zu, resp. zv.

 

Poznámka 5.1. K zapamatování vzorců (5.7) můžeme použít formální umocnění, o kterém jsme se již zmínili u výpočtu diferenciálů vyšších řádů (Poznámka 4.5). Například pro výpočet zxx formálně umocníme pravou stranu rovnosti zx = zuux + zvvx. Dostaneme zu2ux2 + 2zuzvuxvx + zv2vx2, a nahradíme-li druhé mocniny, resp. součin prvních derivací funkce z odpovídajícími druhými derivacemi, obdržíme zuuux2+2zuvuxvx+zvvvx2, což jsou právě první tři členy v (5.7 ).

 

Příklad 5.2. i) Pomocí transformace do nových nezávisle proměnných u = x + ay, v = x ay najděte obecné řešení tzv. vlnové rovnice

 2 a zxx − zyy = 0
(tato rovnice popisuje např. chvění struny na hudebním nástroji, z(x,y) udává velikost výchylky struny ve vzdálenosti x od jednoho z bodů upevnění struny v čase t = y).

Řešení. Využitím vzorečků pro parciální derivace 1. a 2. řádu dostáváme

zx = zuux + zvvx = zu + zv, zy = zuuy + zvvy = azu − azv, 2 2 2 zxx = zuu + 2zuv + zvv, zyy = a zuu − 2a zuv + a zvv.
Dosazením a úpravou obdržíme rovnici zuv = 0, kterou řešíme takto: Označíme-li zu(u,v) = w(u,v), pak řešením rovnice wv = 0 je libovolná funkce nezávisející na v, tedy w = w(u). Řešení rovnice zu = w(u) je tvaru z(u,v) = w(u)du + g(v) (podobně jako při hledání kmenové funkce je „integrační konstantou“ funkce g(v) proměnné v, viz odst. 4.3). Označíme-li f(u) = w(u)du, dostáváme řešení rovnice ve tvaru z(u,v) = f(u) + g(v) a po dosazení za u a v
z(x,y) = f(x+ ay)+ g(x − ay),
kde f,g jsou libovolné funkce jedné proměnné mající derivaci 2. řádu.

Je-li ještě zadána počáteční poloha a rychlost chvějící se struny, tj. je daná dvojice funkcí ϕ,ψ jedné proměnné popisující počáteční stav struny, pak dvojice počátečních podmínek

z(x,0) = ϕ(x), zy(x,0) = ψ(x),
určuje jednoznačně funkci z(x,y) popisující chvění struny, viz např. [T-S].

ii) Pomocí transformace nezávisle proměnných u = xy, v = xy najděte všechny funkce dvou proměnných splňující rovnici

x2zxx + y2zyy − 2xyzxy + xzx + yzy = 0.
Řešení. Podobně jako v předchozím příkladu
 1 x zx = zuux + zvvx = zuy + zv-, zy = zuuy + zvvy = zux− zv--2, y y z = y2z + 2z + z -1, z = xyz − x-z + z − -1z , xx uu uv vvy2 xy uu y3 vv u y2 v x2 x2 x zyy = x2zuu − 2-2zuv +-4zvv + 2zv-3. y y y
Dosadíme do rovnice a po úpravě dostáváme
zuu(x2y2 + x2y2 2x2y2) + zuv(2x2 2x2) + zvv(x2- y4 + x2- y4 + 2x2- y4)+
+ zu(2xy + xy + xy) + zv(x -- y + 2x -- y x -- y) = 4x2 -2- yzvv + x -- yzv = 0,
odtud zvv + 1- 2vzv = 0. Řešením této rovnice je zv(u,v) = f√(u)- v, kde f je libovolná (diferencovatelná) funkce jedné proměnné, a odtud z(u,v) = 2f(u)√ -- v + g(u), kde g je libovolná funkce jedné proměnné se spojitou druhou derivací, což po dosazení za u,v dává
 ∘ -- z(x,y) = 2f(xy) x+ g(xy). y
Přesvědčte se zkouškou, že tato funkce je opravdu řešením dané rovnice.

iii) Transformujte tzv. Laplaceovu rovnici1ℝ2

zxx + zyy = 0
do polárních souřadnic x = r cosϕ, y = r sinϕ, za předpokladu, že funkce z má spojité parciální derivace 2. řádu.

Řešení. Podle (5.2 ) platí

zr = zx xr + zy yr = zx cosϕ + zy sinϕ,
zϕ = zx xϕ + zy yϕ = zx r sinϕ + zy r cosϕ.
Pro derivace 2. řádu dostáváme
zrr = zxxxr2 + 2zxy xryr + zyy yr2 + zx xrr + zy yrr =
= zxx cos2ϕ + zxy sin2ϕ + zyy sin2ϕ,
zϕϕ = zxxxϕ2 + 2zxyxϕyϕ + zyyyϕ2 + zx zϕϕ + zy yϕϕ =
= zxxr2 sin2ϕ zxyr2 sin2ϕ + zyyr2 cos2ϕ zxr cosϕ zyr sinϕ.
Vynásobíme-li vzorec pro zrr výrazem r2 a sečteme se vzorcem pro zϕϕ, dostáváme
r2zrr + zϕϕ = r2zxx(cos2ϕ + sin2ϕ) + r2zyy(cos2ϕ + sin2ϕ)+
+ zxy(sin2ϕ sin2ϕ) r[zx cosϕ + zy sinϕ] =
= r2(zxx + zyy) rzr.
Laplaceova rovnice v polárních souřadnicích má tedy tvar
 2 r zrr + zϕϕ + rzr = 0.

 

Poznámka 5.2. Ve všech řešených příkladech, které jsme zde uvedli, byla transformace do nových proměnných dána již v zadání. V rovnicích matematické fyziky se vyšetřují rovnice typu

a(x,y)z + 2b(x,y)z + c(x,y)z + f(x,y,z,z ,z ) = 0. xx xy yy x y
Chceme-li najít řešení této rovnice, je třeba ji vhodnou transformací do nových nezávisle proměnných zjednodušit – převést na tzv. kanonický tvar. Tuto „vhodnou“ transformaci najdeme prostřednictvím řešení obyčejné diferenciální rovnice
a(x,y)y′2 − 2b(x,y)y′ + c(x,y) = 0.
Podrobnější informace o tomto postupu lze nalézt například v [T-S].

 

Doposud jsme uvažovali pouze funkce dvou proměnných, ale situace pro funkce více proměnných je zcela analogická, včetně důkazu následujícího tvrzení.

Věta 5.3. Nechť je dána funkce f : ℝm → ℝ a m-tice funkcí gi : ℝn → ℝ, které mají spojité parciální derivace 2. řádu. Označme uk = gk(x1,,xn),k = 1,,m. Pak pro složenou funkci F(x1,,xn) = f(g1(x1,xn),,gm(x1,,xn)) platí

 ∂ ∑m ∂ ∂ ---F (x1 ...,xn) = ---f(u1,...,um) ---gk(x1,...,xn) ∂xi k=1 ∂uk ∂xi
(5.8)

 

 m -∂2-F (x ,...,x ) = ∑ -∂2--f(u) ∂--g(x)-∂-g (x)+ ∂xixj 1 n k,l=1 ∂ukul ∂xi k ∂xj l n 2 +∑ -∂-f(u)-∂---gk(x), k=1 ∂uk ∂xixj
(5.9)

kde i,j = 1,2,,n a ve vzorci (5.9 ) je u = (u1,un), x = (x1,,xn).

 

Poznámka 5.3. i) Jsou-li funkce gi ve Větě 5.3 lineární, pak všechny členy v druhé sumě v (5.9) jsou nulové (neboť druhá derivace lineární funkce je nulová). Metoda formálního vynásobení derivací prvního řádu a následná náhrada součinů prvních derivací odpovídajícími druhými derivacemi pak dává přímo vztahy pro druhou derivaci. Takto je tomu např. v Příkladu 5.2 i).

ii) Uvedli jsme si zde pouze vzorce pro parciální derivace složené funkce 1. a 2. řádu, které jsou potřeba v rovnicích matematické fyziky. Metodou stejnou jako v důkazu Věty 5.2 lze odvodit vztahy pro třetí a vyšší derivace, nebudeme je zde však již uvádět, neboť jsou formálně poměrně složité.

 

Příklad 5.3. i) Transformujte Laplaceovu rovnici v ℝ3

uxx + uyy + uzz = 0
do sférických souřadnic x = r cosϕsinϑ, y = r sinϕsinϑ, z = r cosϑ.

Řešení. Mohli bychom postupovat podobně jako při řešení Příkladu 5.2 iii), zde však pro ilustraci různých možných metod postupujeme odlišně. Vyjádříme nejprve r,ϕ,ϑ pomocí x,y,z. Jednoduchými úpravami dostáváme

 ∘ ------- r = ∘x2-+-y2-+-z2, ϕ = arctg-y, ϑ = arctg-x2 +-y2. x z
Nyní vypočteme všechny potřebné parciální derivace funkcí r,ϕ,ϑ. Platí
rx =  x ∘--2---2---2- x + y + z = x r-, ry = y r-, rz = z r,
rxx = 1 r x2 r3, ryy = 1 r y2 r3, rzz = 1 r z2 r3,
ϕx = ---y--- x2 + y2, ϕy = --x---- x2 + y2, ϕz = 0,
ϕxx = --−-2xy--- (x2 + y2)2, ϕyy = --−-2xy-- (x2 +y2)2,ϕzz = 0,
ϑx = ----1---- 1+ x2+y22- z----x----- z∘x2-+-y2- = ---------xz--------- ∘x2-+-y2(x2 + y2 + z2) = ----xz---- r2∘x2-+-y2-,
ϑy = --∘-yz----- r2 x2 + y2, ϑz = ---12--2- 1+ x-+z2y- ∘ -2---2- −---x-+-y- z2 = ∘ -2---2- --x-+-y-- r2,
ϑxx =  z -2∘--2---2- r x + y 2 x2z 4∘--2---2- r x +y  x2z -2--2---2-32 r (x + y ),
ϑyy = --∘-z------ r2 x2 + y2 2--∘-y2z---- r4 x2 + y2 ----y2z---- r2(x2 + y2)32,
ϑzz = 2z∘x2-+-y2- -----4---- r .
Podle vzorců pro derivace složené funkce
ux = urrx + uϕϕx + uϑϑx = urx- r + uϕ---y--- x2 + y2 + uϑ----xz--- r∘x2--+y2-,
uy = ury r- uϕ x x2 +-y2 + uϑ yz -∘--2---2 r x + y,
uz = urz r uϑ∘------- -x2-+-y2- r,

 

uxx = urrx2 r2 + uϕϕ y2 (x2 +-y2)2 + uϑϑ x2z2 r4(x2 +-y2) + 2urϕxyr(x2 + y2)+
+ 2urϑ 2 --∘-x-z---- r2 x2 + y2 2uϕϑ---2xyz23∕2 r(x + y ) + ur( 2) 1 − x3 r r
uϕ 2xy (x2-+-y2)2- + uϑ( ) z x2z x2z -2∘--2----2 − 2-4∘-2---2-− -2-2----232- r x + y r x + y r (x + y),
uyy = urry2 r2 + uϕϕ---x2---- (x2 + y2)2 + uϑϑ---y2z2--- r4(x2 + y2) 2urϕ---xy---- r(x2 + y2)+
+ 2urϑ y2z r2∘x2-+-y2- 2uϕϑ xyz r(x2-+y2)3∕2 + ur(1 y2) r − r3
uϕ---2xy--- (x2 + y2)2 + uϑ( ) --∘--z----- --∘-y2z----- ---y2z---- r2 x2 +y2 − 2r4 x2 + y2 − (x2 + y2)32,
uzz = urrz2 r2 + uϑϑx2 +-y2 r2 2urϑ ∘ ------- z--x2-+y2- r2+
+ ur(1 z2) r − r3 + uϑ2z∘x2-+-y2- ----r4----.
Odtud
 2 2 uxx + uyy + uzz = urr(x2 + y2 +z2)+ uϕϕ-x-+-y-+ ( r2 ) (x2 + y2)2 uϑϑ- --x2z2--- --y2z2--- 2 2 ---urϕ--- + r4 (x2 + y2)2 + (x2 + y2)2 + (x + y ) + 2r(x2 + y2)(xy− xy)+ ( 2 2 ∘ ------) +2 urϑ2- ∘-x-z---+ ∘-y-z---− z x2 + y2 + 2----uϕϑ---3(xyz − xyz)+ r x2 + y2 x2 + y2 r2(x2 + y2)2 (3 x2 +y2 + z2) uϕ +ur r − ----r3----- + (x2-+-y2)2(− 2xy+ 2xy)+ ( ) +u ϑ -∘--2z----− 2z∘x2--+-y2 + 2z∘x2-+-y2-−-∘--z----- = r2 x2 + y2 r4 r4 r2 x2 + y2 2 1 1 cotgϑ = urr + r ur + r2sin2-ϑuϕϕ + r2uϑϑ +-r2 uϑ.

ii) Určete řešení Laplaceovy rovnice v ℝ3, které je sféricky symetrické, tj. závisí pouze na vzdálenosti od počátku.

Řešení. Nechť funkce u závisí pouze na proměnné r = ∘ ----------- x2 + y2 + z2 a nikoliv na proměnných ϕ a ϑ, tj. u = u(r) (tento předpoklad je „rozumný“ vzhledem k fyzikálnímu významu Laplaceovy rovnice). Pak všechny parciální derivace podle ϕ, ϑ jsou rovny nule a dostáváme rovnici

u + 2u = 0. rr r r
Položíme-li ur = v, dostáváme dále rovnici vr + 2 rv = 0 a po úpravě r2 vr + 2rv = 0, což je ekvivalentní rovnici ∂- ∂r(r2v) = 0. Řešením této rovnice je např. v(r) = 1- r2, a tedy u = 1 r, tj.
u(x,y,z) = ∘----1------ x2 + y2 + z2
je jedním z řešení Laplaceovy rovnice (srov. Příklad 3.3 ii) ).

 

 

5.2 Taylorova věta

Nejprve připomeňme, co to je Taylorův polynom a Taylorova věta2 pro funkci jedné proměnné. Nechť f : ℝ → ℝ, x0,x ∈ ℝ a h = x x0. Taylorův polynom (mnohočlen) stupně n ∈ ℕ funkce f se středem v bodě x0 je polynom

 (k) Tn(x; x0) = a0 + a1(x− x0) + ⋅⋅⋅+ an(x − x0)n, ak = f-(x0), k!
k = 0,,n. Koeficienty ak určíme z požadavku, aby polynom Tn měl v bodě x0 stejnou funkční hodnotu a hodnotu prvních n derivací jako funkce f.

Taylorův polynom používáme k přibližnému výpočtu funkčních hodnot funkce f v okolí bodu x0. Taylorova věta udává velikost chyby, které se dopustíme, aproximujeme-li funkci Taylorovým polynomem.

Obdobně je tomu u funkce více proměnných. Taylorův polynom funkce f : ℝn → ℝ je polynom více proměnných, který má s funkcí f v daném bodě x = [x1,,xn] ∈ ℝn stejnou funkční hodnotu a stejnou hodnotu všech parciálních derivací až do řádu n, kde n je stupeň polynomu. Pro funkce dvou proměnných dostáváme toto tvrzení.

Věta 5.4. (Taylorova) Nechť funkce f : ℝ2 → ℝ má v bodě [x0,y0] a nějakém jeho okolí spojité parciální derivace až do řádu n + 1 včetně. Pak pro každý bod [x,y] z tohoto okolí platí

f (x, y) = Tn(x,y) + Rn(x, y),
(5.10)

kde

Tn(x,y) = f(x0,y0) + ∂f- ∂x(x0,y0)h + ∂f- ∂yy(x0,y0)k+
+ 1- 2!( 2 2 2 ) ∂-f-(x0,y0)h2 + 2-∂-f-(x0,y0)hk + ∂-f(x0,y0)k2 ∂x2 ∂x ∂y ∂y2 + +
+ -1 n! j=0n( ) n j---∂nf--- ∂xn −j∂yj(x0,y0)hnjkj,
 n+1 ( ) ---1----∑ n + 1 --∂n+1f---- n+1− j j Rn(x, y) = (n + 1)! j ∂xn+1−j∂yj (x0 + ϑh,y0 + ϑk)h k j=0
a kde h = x x0, k = y y0, ϑ ∈ (0,1).

 

Poznámka 5.4. Vzorec (5.10) se nazývá Taylorův vzorec, polynom Tn Taylorův polynom a Rn zbytek v Taylorově vzorci.

Taylorův vzorec lze zapsat pomocí diferenciálů takto:

f(x,y) = f(x0,y0) + df(x0,y0)(h,k) + 1- 2d2f(x0,y0)(h,k) + +
+ 1- n!dnf(x0,y0)(h,k) + ---1---- (n + 1)!dn+1f(x0 + ϑh,y0 + ϑk)(h,k).

 

Důkaz Věty 5.4. Zaveďme pomocnou funkci jedné proměnné F(t) = f(x0 + th,y0 + tk). Platí F(1) = F(x0 + h,y0 + k) = F(x,y), F(0) = f(x0,y0). Pomocí Taylorova vzorce pro funkci jedné proměnné dostáváme

 1 1 1 F (1) = F (0)+ F ′(0) +--F′′(0)+ ⋅⋅⋅+ --F(n)(0)+ --------F(n+1)(ϑ), 2! n! (n + 1)!
kde ϑ ∈ (0,1). Pro výpočet derivací funkce F využijeme vztahů pro parciální derivace složených funkcí. Dostáváme
 d ∂ ∂ F′(0) = --f(x0 + th, y0 + tk)∣t=0 =---f(x0,y0)h+ --f(x0,y0)k, dt ∂x ∂y
F′′(0) =  2 -d- dt2F(t)∣t=0 =  2 d-- dt2f(x0 + th,y0 + tk)∣t=0 =
= fxx(x0,y0)h2 + 2fxy(x0,y0)hk + fyy(x0,y0)k2
a analogicky obdržíme
 m ( ) (m) ∑ m --∂mf--- m−j j F (0) = j ∂xm −jyj(x0,y0)h k . j=0
Stejně postupujeme i při výpočtu zbytku Rn.

 

Příklad 5.4. i) Určete Taylorův polynom 2. stupně se středem v bodě [x0,y0] = [1,1] pro funkci f(x,y) = x y.

Řešení. Vypočteme nejprve všechny potřebné parciální derivace

f = 1-, f = − x-, f = 0, f = − 1-, f = 2x-. x y y y2 xx xy y2 yy y3
Podle Věty 5.4
 T2(x, y) = f (1,1) + fx(1,1)(x − 1)+ fy(1,1)(y − 1)+ 1- 2 2 + 2 [fxx(1,1)(x− 1) + 2fxy(1,1)(x − 1)(y − 1)+ fyy(1,1)(y − 1) ] = = 1+ (x − 1)− (y − 1)− (x − 1)(y − 1)− (y − 1)2 = 2 = − y − xy + 2x+ 2y − 1.

ii) Pomocí Taylorova polynomu 2. stupně vypočtěte přibližně:

a) ∘(2,-98)2 +-(4,05)2; b) 1,042,02.

 

Výsledek porovnejte s hodnotou získanou pomocí diferenciálu z Příkladu 4.2 ii).

Řešení. a) Přibližnou hodnotu vypočteme pomocí Taylorova polynomu 2. stupně funkce z = f(x,y) = ∘ -2---2- x + y v bodě [x0,y0] = [3,4] a diferencemi h = 0,02, k = 0,05. Parciální derivace funkce z jsou

 2 zx = ∘---x----, zy = ∘---y----, zxx =-----y-----, x2 + y2 x2 + y2 (x2 + y2)3∕2 xy x2 zxy = −--2----2-3∕2, zyy = --2---2-3∕2, (x + y ) (x + y )
Taylorův polynom je roven
 T2(x, y) = f (3,4) + fx(3,4)(x − 3)+ fy(3,4)(y − 4)+ 1 + --[fxx(3,4)(x− 3)2 + 2fxy(3,4)(x − 3)(y − 4)+ fyy(3,4)(y − 4)2] = 2 = 5 + 1[3(x − 3)+4(y − 4)]+ -1-[16(x− 3)2− 24(x − 3)(y− 4)+9(y − 4)2]. 5 250
Odtud
 ∘(2,98)2-+-(4,05)2=.5 + 1(− 0,06 + 0,2)+ 5 -1-- +250 (16⋅0,0004− 24⋅0,001+ 9 ⋅0,0025) = 5,0281332.
V příkladu 4.1 ii) jsme pomocí diferenciálu dostali výsledek ∘ ------2-------2- (2,98) + (4,05)5,028.
b) V Taylorově vzorci pro funkci z = xy položme [x0,y0] = [1,2], h = 0,04,k  = 0,02. Nejprve vypočtěme všechny potřebné parciální derivace. Platí zx = yxy1, zx(1,2) = 2, zy = xy lnx, zy(1,2) = 0, zxx = y(y1)xy2, zxx(1,2) = 2, zxy = xy1+yxy1 lnx = xy1(1+y lnx),zxy(1,2) = 1, zyy = xy lnxlnx = xy ln2x,zyy(1,2) = 0. Pak
T2(x, y) = 1 + 2(x− 1) + (x− 1)2 + (x − 1)(y − 1).
Odtud
 . 1,042,02= 1+ 2 ⋅0,04 + 0,0016+ 0,0008 = 1,0824.
V Příkladu 4.1 ii) jsme pomocí diferenciálu obdrželi přibližný výsledek 1,042,02 1,08.

iii) Mnohočlen P(x,y) = x3 + 3y3 + xy2 + 2x2 + xy + x 2y napište jako polynom v proměnných u = x 1,v = y + 2.

Řešení. Nechť T3(x,y) Taylorův polynom 3. stupně funkce P se středem x0 = 1,y0 = 2. Pak ve zbytku R3(x,y) vystupují 4. derivace funkce P, které jsou však všechny nulové, neboť P je polynom 3. stupně. Tedy T3(x,y) = P(x,y) a stačí určit koeficienty v T3(x,y). Postupně dostáváme P(1 ,2) = 20 Px = 3x2 + y2 + 4x + y + 1,Px(1,2) = 10, Py = 9y2 + 2xy + x 2,Py(1,2) = 31, Pxx = 6x + 4, Pxx(1,2) = 10, Pxy = 2y + 1,Pxy(1,2) = 3, Pyy = 18y + 2x, Pyy(1,2) = 34, Pxxx = 6,Pxxy = 0, Pxyy = 2,Pyyy = 18. Odtud

T3(x, y) = − 14+ 10(x − 1) + 31(y + 2)+ 5(x − 1)2 − 3(x − 1)(y + 2)− 2 3 2 3 − 17(y + 2) + (x− 1) + (x− 1)(y + 2) + 3(y + 2).
Jestliže ve výsledku provedeme umocnění, po úpravě samozřejmě dostáváme polynom P. Tuto kontrolu výsledku necháváme čtenáři jako cvičení.

 

Zformulujme na závěr kapitoly ještě Taylorův vzorec pro obecný případ funkcí n proměnných. Důkaz tohoto tvrzení neuvádíme, neboť je v podstatě stejný jako pro dvě proměnné.

Věta 5.5. Nechť funkce f : ℝn → ℝ má v bodě x = [x1,,xn] a nějakém jeho okolí spojité parciální derivace až do řádu m + 1. Pak pro h = [h1,,hn] platí

 ∗ ∗ ∗ 1 2 ∗ -1-m ∗ f(x + h) = f (x )+ df(x )(h) + 2d f(x )(h)+ ⋅⋅⋅+ m!d f(x )(h) +Rm(x),
kde
 1 Rm(x) = (m-+-1)!dm+1f (x∗ + ϑh)(h), ϑ ∈ (0,1)
je zbytek v Taylorově vzorci a
 k dkf(x ∗)(h) = ∑ ----k!--------∂-f---(x∗)hj1...hjn j1+⋅⋅⋅+jn=kj1!j2!...jn!∂xj11 ...xjnn 1 n
je k-tý diferenciál funkce f v bodě x.

 

Cvičení PIC

5.1.  Využitím uvedené substituce najděte všechny funkce splňující danou rovnost:
a) yzx xzy = 0, u = x, v =  ------- ∘ x2 + y2
b) xzx + yzy = 0, u = x, v = y x
c) ux + uy + uz = 0, ξ = x + y 2z, η = x 2y + z, χ = z

5.2.  Diferenciální rovnice transformujte do nových proměnných u,v. V případech, kdy po transformaci vyjde jednoduchý výsledek, se pokuste najít jejich řešení:
a) zxx yzyy 1 2zy = 0, u = x 2√ -- y,v = x + 2√ -- y
b) y2zxx + x2zyy 2xyzxy xzx yzy = 0,u = ∘ ------- x2 + y2,v = xy
c) x2zxx (x2 + y2)zxy + y2zyy = 0, u = x + y,v = 1 x + 1 y
d) zxx 2zxy + zyy = 0, u = x + y,v = -1-- x−y
e) xyzxx (x2 + y2)zxy + xyzyy + yzx + xzy = 0, u = 1 2(x2 + y2), v = xy
f) xzxx yzyy = 0, u = √ -- x + √-- y,v = √ -- x √-- y
g) xzxx + yzxy + zx = 0,u = x + y,v = xy+y-
h) x2zxx 2xyzxy + y2zyy + xzx + yzy = 0, u = xy, v = y
i) x2zxx y2zyy = 0, u = xy, v = y x

5.3.  Ukažte, že daná transformace do nových proměnných nemění tvar rovnice zxx + zyy = 0, x = ϕ(u,v),y = ψ(u,v), kde ϕ,ψ jsou funkce dvou proměnných splňující identity ϕu = ψv, ϕv = ψu.

5.4.  Určete Taylorův polynom 2. stupně se středem [x0,y0] následujících funkcí:

a) ∘ ----------- 1− x2 − y2, [x0,y0] = [12,12] e) arcsin√--x2--2 x +y, [x0,y0] = [0,1]

 

b) arctg 1+x+y- 1−x+y, [x0,y0] = [0,0] f) ln∘ -2----2 x + y, [x0,y0] = [1,1]

 

c) cosx- cosy, [x0,y0] = [0,0] g) xyz, [x0,y0,z0] = [1,1,1]

 

d) arctg yx, [x0,y0] = [1,1] h) sinxsiny, [x0,y0] = [0,0]

 

 

5.5.  Pomocí Taylorova polynomu 2. stupně vypočtěte přibližně funkční hodnoty:
a) arctg 1,04 0,98 b) sin29∘tg 46∘

 

Zkušenost není to, co člověk potká, ale co člověk udělá s tím, co ho potkalo. (A. Huxley)

1Pierre Simon Laplace (1749–1827), francouzský matematik, fyzik a astronom

2Brook Taylor (1685–1731), anglický matematik