Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
15<br />
1. MILNÍKY MATEMATIKY PŘED <strong>EULER</strong>EM<br />
1.6. Diferenciální rovnice, mechanika a další rozvoj matematiky<br />
Proměnné, funkce a derivace mohou být svázány do relace, diferenciální rovnice.<br />
Např. zrychlení volného pádu tělesa v atmosféře závisí na výšce a odporu vzduchu,<br />
d y<br />
který je funkcí rychlosti tělesa, R(<br />
d t ). Pád z malé výšky y0 > 0 a při bezvětří lze<br />
aproximovat jednodušším pohybem s konstantním gravitačním zrychlením –g a<br />
d y<br />
lineárním odporem vzduchu R(<br />
d t<br />
dy<br />
)/<br />
dt<br />
= const (kvadratická závislost R(<br />
by samozřejmě byla realističtější [27], jak dosvědčí každý cyklista)<br />
d<br />
dt<br />
s počátečními podmínkami y(0) = y 0 ,<br />
2<br />
2<br />
d y dy<br />
) =<br />
dt<br />
dt<br />
d y<br />
(c 1 + c 2 )<br />
d t<br />
d<br />
y + g + a y = 0 (*)<br />
dt<br />
d y (0) = 0.<br />
dt<br />
------<br />
Relace (*) je lineární diferenciální rovnice a její řešení se najde snadno. Přepíšeme ji do tvaru<br />
d<br />
dt (<br />
d<br />
dt<br />
y + ay) = –g. Z něj dostaneme<br />
d<br />
dt<br />
y + ay = –gt + c, kde c = const. Pro t = 0 z počátečních<br />
podmínek plyne<br />
d<br />
dt y(0) + a×y(0) = –g×0 + c, tj. c = ay0 . Řešení nové diferenciální rovnice<br />
d<br />
dt y + ay = –gt + ay0<br />
hledejme ve tvaru y(t) = C(t) e –at . Po dosazení dostaneme<br />
d<br />
dt y(t) + a y(t) =<br />
t<br />
C(t) – C(0) = ∫<br />
0<br />
Takže y(t) = {C(0) + y 0 (e at – 1) –<br />
odkud pro t = 0 plyne:<br />
dC<br />
e –at – aC e –at + aC e –at = –gt + ay<br />
dt<br />
0 , tj.<br />
t<br />
(–gt + ay 0 ) e at dt = y 0 e at – g ∫<br />
0<br />
g [ t e at +<br />
a<br />
y(0) = y 0 = C(0).<br />
------<br />
Rovnice (*) má pro a ≠ 0 řešení<br />
Po rozvinutí exponenciály v řadu<br />
dC<br />
dt<br />
t e at dt = y 0 (e at – 1) –<br />
1 (1 – e at )]} e –at = C(0) e –at –<br />
a<br />
g g<br />
y(t) = y 0 – t +<br />
−at<br />
a<br />
( 1− e )<br />
a<br />
2<br />
= (–gt + ay 0 ) e at a<br />
.<br />
g [ t e at +<br />
a<br />
g t + (y0 +<br />
a<br />
1 (1 – e at )] .<br />
a<br />
g<br />
2<br />
a<br />
)(1 – e –at ),<br />
2<br />
gt ⎛ 3 4 2 5 ⎞<br />
y(t) = y 0 – + ag⎜<br />
t at a t<br />
− + − ... ⎟<br />
2<br />
⎝ 3! 4! 5! ⎠<br />
je ihned zřejmé, že pro a = 0 (nulový odpor) dostáváme elementární volný pád ve<br />
vakuu.<br />
Obr. 1.11 ilustruje pád tělesa v gravitačním poli na povrchu Země (g = 9.8m/s 2 )<br />
z výšky 10m při různém a, tedy s různým odporem prostředí. Funkce y pro nejhustší<br />
prostředí s a = 3s –1 je prakticky lineární už pro t >1s (kámen ve vodě), u funkce y pro<br />
a = 0.2s –1 by se linearity dosáhlo až po mnohem delším čase (>30s, kámen ve vzduchu).<br />
F. KOUTNÝ: <strong>Leonhard</strong> <strong>EULER</strong>