20.11.2013 Views

Leonhard EULER

Leonhard EULER

Leonhard EULER

SHOW MORE
SHOW LESS

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>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!