Naučili jsme se pracovat s derivacemi, tedy s rychlostí změny. Známe-li funkci a zderivujeme ji, dostaneme rychlost změny. Pokud potom původní funkci “ztratíme” a zůstane nám jenom derivace, je otázka, jestli dokážeme původní funkci z této derivace najít. Odpověď je zní, že v jistém smyslu ano. Spojení “v jistém smyslu” naznačuje, že souvislost nebude tak snadná jako je souvislost u navzájem inverzních funkcí. Derivováním totiž můžeme ztratit aditivní konstanty, které v derivaci dávají nulu a zpětně není možné rekonstruovat, derivováním jaké konstanty jsme tuto nulu dostali. A protože problém uchopíme poněkud obecněji, uvedeme si dokonce hned tři různé “protijedy” na derivování.

Jeden představíme jako opak derivace (neurčitý integrál), druhý jako změnu funkce vypočtenou ze zadané rychlosti změny (Newtonův určitý integrál) a třetí jako náhradu součtu pro případ, kdy potřebujeme sčítat nekonečně mnoho příspěvků, z nichž každý má v podstatě nulovou hodnotu (Riemannův určitý integrál).

Intervalem $I$ budeme rozumět otevřený interval.

Motivace: Jak z rychlosti změny vypočítat změnu?

Derivace umožní z veličiny v prvním sloupci získat veličinu v pravém sloupci. Pohledem na tyto příklady věříme, že bude fungovat i něco, co naopak z rychlosti zrekonstruuje původní veličinu, která se touto rychlostí mění.

Závislá proměnná	Derivace podle času
veličina $x$	rychlost růstu veličiny $x$
výška stromu	rychlost růstu do výšky
objem kmene stromu (smrk)	rychlost růstu ve smyslu přírůstu dřevní hmoty
dráha	rychlost
rychlost	zrychlení
všeobecná cenová hladina (cca náklady na živobytí)	inflace

Neurčitý integrál

Představíme nástroj, který nám umožní odpovědět na následující otázky.

• Je znám směr křivky v každém bodě (tj. směr tečny, derivace). Jaká je rovnice křivky?
• Je známa rychlost, s jakou se mění veličina $f$. Jaká je rovnice udávající závislost veličiny $f$ na čase?

Definice (neurčitý integrál).

Řekneme, že funkce $F$ je primitivní funkcí k funkci $f$ na intervalu $I$, jestliže platí $$F^{\prime }(x)=f(x)$$ na intervalu $I$. Množina všech primitivních funkcí k funkci $f$ se nazývá neurčitý integrál funkce $f$ a značí $$\int f(x)\mathrm{d}x.$$

Otázkou existence primitivní funkce se budeme zabývat na další přednášce. Otázku (ne-)jednoznačnosti řeší následující věta.

Věta (jednoznačnost primitivní funkce).

Primitivní funkce je dána jednoznačně, až na aditivní konstantu.

• Je-li $F$ primitivní funkcí k funkci $f$ na intervalu $I$, platí totéž i pro funkci $G(x)=F(x)+c$, kde $c\in \mathbb{R}$.
• Jsou-li $F$ a $G$ primitivní funkce k téže funkci $f$ na intervalu $I$, existuje $c\in \mathbb{R}$ takové, že $$F(x)=G(x)+c$$ na $I$.

Příklad. Funkce $x^2$ má primitivní funkce například $\frac{1}{3}{x}^3$, nebo $\frac{1}{3}{x}^3+7$, nebo $\frac{1}{3}{x}^3+\pi$, protože derivace všech těchto tří funkcí je $x^2$. Platí $$\int x^2\mathrm{d}x=\frac{1}{3}{x}^3+c,c\in \mathbb{R}.$$

Věta (linearita neurčitého integrálu).

Neurčitý integrál zachovává součet a násobení konstantou. Tedy pro libovolné funkce $f$, $g$ a libovolnou konstantu $c$ platí $$\begin{array}{rl}\int f+g\mathrm{d}x& =\int f\mathrm{d}x+\int g\mathrm{d}x,\\ \int cf\mathrm{d}x& =c\int f\mathrm{d}x.\end{array}$$

Příklad.

$$\int 2x^4-e^{4x}+\frac{1}{x}\mathrm{d}x=\frac{2}{5}{x}^5-\frac{1}{4}{e}^{4x}+\ln |x|+C$$

Integrování si také můžete procvičit v následujících cvičeních.

Testová otázka Testová otázka Testová otázka Testová otázka Testová otázka Testová otázka

Funkční předpis z rychlosti změny a výchozího stavu

Uvažujme těleso, jehož teplota klesá známou rychlostí. Derivace teploty podle času je $\frac{\mathrm{d}T}{\mathrm{d}t}=-0.1e^{-0.01t}{}^{\circ }\mathrm{C}/\min .$ Cílem je najít teplotu jako funkci času. Dodatečná informace je, že počáteční teplota je ${28}^{\circ }\mathrm{C}$.

Použijeme skutečnost, že integrál konstantního násobku je konstantní násobek integrálu a vzorec $$\int e^{ax}\mathrm{d}x=\frac{1}{a}{e}^{ax}+c.$$ Teplota jako funkce času je dána integrálem $$T=\int -0.1e^{-0.01t}\mathrm{d}t=\frac{-0.1}{-0.01}{e}^{-0.01t}+C=10e^{-0.01t}+C.$$ Hodnota $C$ souvisí s počáteční teplotou. Protože počáteční teplota je ${28}^{\circ }\mathrm{C}$, dosadíme do vztahu pro $T$ hodnoty $T={28}^{\circ }\mathrm{C}$ a $t=0$ a ze vzniklé rovnice určíme $C$. Dostáváme takto podmínku $$28=10e^0+C,$$ která implikuje $C={18}^{\circ }\mathrm{C}$. Funkce udávající závislost teploty místnosti na čase je $$T=(18+10e^{-0.01t}){}^{\circ }\mathrm{C}.$$

Poznámka (vlhkost dřeva elektrickou metodou). Podobný výpočet se využívá u měření elektrického odporu dřeva pro stanovení vlhkosti. Protože elektrický odpor dřeva je velký, není vhodné pro určení elektrického odporu použít Ohmův zákon a změřený proud a napětí. Jedna z možností je měření času nutného k nabití nebo vybití kondenzátoru přes odpor. V případě nabíjení proud exponenciálně klesá (zdůvodníme později v přednášce věnované diferenciálním rovnicím) a proto (díky elektrickým vlastnostem kondenzátoru) exponenciálně klesá i rychlost, s jakou roste napětí na kondenzátoru. Toto napětí je nutné pro výpočet odporu. Pokud známe rychlost, s jakou se napětí mění, určíme napětí integrováním a znalostí napětí na začátku nabíjení.

Poznámka (veličina vypočtená z rychlosti své změny).

Pokud se veličina $f(t)$ mění v čase rychlostí $r(t)$, platí $$f(t)=\int r(t)\mathrm{d}t,$$ přičemž pravá strana je dána jednoznačně až na aditivní konstantu. To koresponduje s pozorováním, že rychlost změn k jednoznačné identifikaci časového průběhu měnící se veličiny nestačí. Je potřeba mít zadán ještě výchozí stav.

Příklad. V úvodu přednášky je popsáno, že křivka, která je přirozená pro nosné lano zavěšeného mostu, splňuje rovnici $$y^{\prime }=\frac{\mu g}{T}x.$$ Pouze za této podmínky bude lano namáháno ve směru své nejvyšší pevnosti, tj. v podélném směru, ve směru své osy. Integrací získáme $$y=\int \frac{\mu g}{T}x\mathrm{d}x=\frac{\mu g}{2T}{x}^2+C.$$ Lano tedy bude v každém bodě namáháno ve směru své osy pokud má tvar paraboly. Prohnutí paraboly (koeficient u $x^2$) je dáno hmotností mostu a tahem napínajícím lano.

Určitý integrál (Newtonův)

Představíme si mírnou modifikaci neurčitého integrálu. Rychlost změny nebudeme používat k hledání předpisu funkce, ale budeme hledat změnu funkce na zadaném intervalu.

Definice (Newtonův určitý integrál).

Buď $f$ funkce a $F$ její primitivní funkce na intervalu $I$. Buď $[a,b]\subset I$ podinterval v $I$. Určitým integrálem funkce $f$ na intervalu $[a,b]$ rozumíme veličinu označenou a definovanou vztahem $${\int }_a^{b}f(x)\mathrm{d}x:=F(b)-F(a).$$

Označení. Výraz $F(b)-F(a)$, tj. změnu funkce $F(x)$ na intervalu $[a,b]$, označujeme také $[F(x){]}_a^b$. Tento zápis se často používá jako mezivýpočet při výpočtu určitého integrálu. $${\int }_0^1{x}^2\mathrm{d}x={\left[\frac{1}{3}{x}^3\right]}_0^1=\frac{1}{3}(1{)}^3-\frac{1}{3}(0{)}^3=\frac{1}{3}$$

Věta (linearita určitého integrálu).

Určitý integrál zachovává součet a násobení konstantou. Tedy pro libovolné funkce $f$, $g$ a libovolnou konstantu $c$ platí $$\begin{array}{rl}{\int }_a^{b}f+g\mathrm{d}x& ={\int }_a^{b}f\mathrm{d}x+{\int }_a^{b}g\mathrm{d}x,\\ {\int }_a^{b}cf\mathrm{d}x& =c{\int }_a^{b}f\mathrm{d}x.\end{array}$$

Snadným důsledkem definice určitého integrálu je následující věta.

Věta (záměna mezí a rovnost mezí v určitém integrálu).

Platí $$\begin{array}{rl}{\int }_a^{a}f(x)\mathrm{d}x& =0,\\ {\int }_a^{b}f(x)\mathrm{d}x& =-{\int }_b^{a}f(x)\mathrm{d}x.\end{array}$$

Změna funkce z rychlosti změny (časová změna teploty)

Uvažujme těleso, jehož teplota klesá známou rychlostí. Derivace teploty podle času je $\frac{\mathrm{d}T}{\mathrm{d}t}=-0.1e^{-0.01t}{}^{\circ }\mathrm{C}/\min .$ Chceme určit pokles teploty za první hodinu a pokles teploty za druhou hodinu.

Neurčitý integrál $$\int -0.1e^{-0.01t}\mathrm{d}t=10e^{-0.01t}+C$$ jsme vypočítali v podkapitole s neurčitým integrálem. Potřebovali jsme ještě znát počáteční hodnotu teploty a našli jsme teplotu jako funkci času.

Nyní zapojíme určitý integrál. Nepotřebujeme informaci o počáteční teplotě, ale zato jsme schopni určit jenom změnu teploty za daný časový interval. Za první hodinu bude změna teploty $${\int }_0^{60}-0.1e^{-0.01t}\mathrm{d}t={\left[10e^{-0.01t}\right]}_0^{60}=10e^{-0.01\cdot 60}-10e^{-0.01\cdot 0}\approx -{4.5}^{\circ }\mathrm{C}.$$ Za druhou hodinu bude změna teploty $${\int }_{60}^{120}-0.1e^{-0.01t}\mathrm{d}t={\left[10e^{-0.01t}\right]}_{60}^{120}=10e^{-0.01\cdot 120}-10e^{-0.01\cdot 60}\approx -{2.5}^{\circ }\mathrm{C}.$$

Online výpočet.

Integrování určitým integrálem si také můžete procvičit v následujících cvičeních.

Testová otázka Testová otázka Testová otázka

Poznámka (změna veličiny vypočtená pomocí rychlosti).

Pokud se veličina $f(t)$ mění v časovém intervalu od $t=a$ do $t=b$ rychlostí $r(t)$, je změna veličiny $f$ za tento časový okamžik rovna $$\mathrm{\Delta }f=f(b)-f(a)={\int }_a^{b}r(t)\mathrm{d}t.$$

Slovní úlohy kde se hledaná veličina mění nekonstantní rychlostí jsou v následujících dvou cvičeních (první je s volbou s nabízených odpovědí).

Testová otázka Testová otázka Testová otázka

Volitelný text. Klikněte pro rozbalení.

Tato ukázka demonstruje převod rychlosti změny teploty v prostoru na změnu teploty mezi dvěma místy.

Změna funkce z rychlosti změny (prostorová změna teploty)

Mějme ustálené proudění tepelnou izolací mezi dvěma soustřednými válcovými plochami. Délka izolace je $L$, vnitřní a vnější poloměr jsou $r$ a $R$. Teploty uvnitř a vně jsou $T_1$ a $T_2$. Izolací prostupuje teplo rychlostí $Q$, tj. každou myšlenou válcovou plochou o poloměru $x$ projde za jednotku času teplo $Q$. Cílem je najít vztah dávající uvedené veličiny do souvislosti. Odvodíme vztah, který jsme použili v přednášce o lokálních extrémech a slíbili dokázat později.

Z Fourierova zákona plyne, že teplo, které projde jednotkovou plochou za jednotku času, je úměrné záporně vzatému gradientu prostorové změně teploty, tj. $$\frac{Q}{2\pi xL}=-k\frac{\mathrm{d}T}{\mathrm{d}x},$$ kde $k$ je tepelná vodivost. Tento zákon definuje rychlost, s jakou se mění teplota v prostoru a my jej použijeme ke stanovení vztahu mezi tokem tepla a teplotním rozdílem mezi vnější a vnitřní stranou izolace. Snadnou úpravou dostáváme $$\frac{\mathrm{d}T}{\mathrm{d}x}=-\frac{Q}{2k\pi xL}$$ a odsud $$\mathrm{\Delta }T={\int }_r^R-\frac{Q}{2k\pi xL}\mathrm{d}x=-\frac{Q}{2k\pi L}{\int }_r^R\frac{1}{x}\mathrm{d}x.$$ Protože platí $${\int }_r^R\frac{1}{x}\mathrm{d}x=[\ln x{]}_r^R=\ln R-\ln r=\ln \frac{R}{r},$$ dostáváme po dosazení a vynásobení minus jedničkou vztah $$-\mathrm{\Delta }T=T_1-T_2=\frac{Q}{2k\pi L}\ln \frac{R}{r},$$ který jsme použili v přednášce o lokálních extrémech a slíbili dokázat později.

Stejný princip funguje pro libovolné ustálené proudění radiálním směrem při konstantní materiálové charakteristice. Stejný přístup je možné použít pro proudění podzemní vody popsané Darcyho zákonem (namísto Fourierova zákona) pro zvodeň namísto izolace (zvodeň je prostor kde se nachází a teče podzemní voda, tj. něco jako podzemní rybník nebo řeka zasypaná pískem nasáklým vodou) a piezometrickou výšku $h$ namísto teploty (piezometrická výška udává, jak vysoko vystoupá voda ve zkušebním vrtu). Pokud máme zvodeň s napjatou hladinou (voda je pod tlakem sevřena mezi dvěma nepropustnými vrstvami), je vodivost konstantní. Rovnice popisující tuto situaci má tvar $$h-h_0=\frac{Q}{2\pi T}\ln \frac{r}{r_0}$$ a nazývá se Thiemova rovnice.

Pokud sledujeme prostup tepla izolací, jejíž teplotní vodivost se mění s teplotou, není veličina $k$ konstantní a proto výše uvedený postup není možné realizovat a odvozený vzorec pro takový případ neplatí. Stejná situace nastává u podzemní vody a proudění s volnou hladinou (není horní nepropustná vrstva zvodně). Takové úlohy vedou na jinou problematiku, kterou se naučíme řešit v kapitole s diferenciálními rovnicemi.

Poznámka (změna veličiny vypočtená pomocí gradientu).

Pokud se veličina $f$ mění podél přímky v závislosti na veličině $x$ na intervalu od $x=a$ do $x=b$ rychlostí $r(x)$ (tj. $r(x)=\frac{\mathrm{d}f(x)}{\mathrm{d}x}$), je změna veličiny $f$ na intervalu $[a,b]$ rovna $$\mathrm{\Delta }f=f(b)-f(a)={\int }_a^{b}r(x)\mathrm{d}x.$$

Další motivace

Ze středoškolské fyziky dobře známe vzorce pro dráhu, práci a tlakovou sílu. Ovšem jenom v extrémně pěkných případech.

• Dráha rovnoměrného pohybu je určena vzorcem $$\begin{array}{}\text{(1)}& s=vt.\end{array}$$ Tento vzorec není použitelný pro pohyb proměnnou rychlostí. Z kapitoly o neurčitém integrálu víme, že obecný vzorec je $$\begin{array}{}\text{(2)}& s=\int v\mathrm{d}t.\end{array}$$ Pokud je $v$ konstantní, vzorec (1) je důsledkem vzorce (2).
• Hydrostatická tlaková síla $F$ působící ve vodě v hloubce $h$ na plochu o velikosti $S$ se určí podle vztahu $$F=Sh\rho g,$$ kde $\rho$ je hustota vody a $g$ tíhové zrychlení. Tento vzorec však není možné použít, pokud různé části plochy jsou v různých hloubkách. Například není možné pomocí tohoto vzorce určit celkovou sílu na svislou stěnu reprezentující hráz přehrady.
• Práce vykonaná konstantní silou $F$ po dráze $s$ je $$\begin{array}{}\text{(3)}& W=Fs.\end{array}$$ Co když se ale síla nebo dráha mění? Pokud nás zajímá práce nutná k navinutí visícího řetězu na rumpál, síla se během namotávání plynule zmenšuje, protože visící kus řetězu se při namotávání zkracuje. Pokud nás zajímá práce nutná k vyčerpání vodní nádrže, musíme každý litr vody, který je na dně, “tahat” po delší dráze než každý litr vody, který je na hladině a proto se mění dráha. Vzorec (3) selhává v obou případech. Jednou kvůli nekonstantní síle, podruhé kvůli dráze.

• Obsah obrazce mezi konstantní funkcí $f$ a osou $x$ nad intervalem $[a,b]$ se vypočte snadno, protože se jedná o obdélník se stranami $f$ a $\mathrm{\Delta }x=b-a$. Proto $$S=f\cdot \mathrm{\Delta }x.$$ Tento přístup však není možné použít, pokud se funkce $f$ na intervalu $[a,b]$ mění. Formálně je tato úloha stejná jako ostatní úlohy výše, má však snadnou geometrickou interpretaci. Právě tuto interpretaci využijeme v následujícím k definici druhého typu určitého integrálu (Riemannova).

Určitý integrál (Riemannův)

Úloha 1. Snadným důsledkem vzorce pro obsah obdélníka je obsah obrazce mezi grafem konstantní funkce a osou $x$. $$S=f\mathrm{\Delta }x$$

Úloha 2. Obsah pod funkcí složené ze dvou konstantních funkcí napojených na sebe se vypočte jako součet obsahů dvou obdélníků. $$S=f_1\mathrm{\Delta }{x}_1+f_2\mathrm{\Delta }{x}_2$$ Toto se dá snadno zobecnit na libovolný počet intervalů a pro libovolnou po částech konstantní funkci.

Prostředky matematické analýzy je možné “zjemňovat dělení do nekonečna”, přesněji, můžeme použít limitní přechod podobný limitnímu přechodu, který v definici derivace převedl podíl (průměrnou rychlost) na derivaci (okamžitou rychlost). Díky tomu není nutné se omezovat na po částech konstantní funkce, ale postup bude fungovat i pro velmi obecné funkce. Výsledným produktem je Riemannův integrál.

Riemannův integrál je velmi názorný, ale poměrně obtížně se počítá, pokud postupujeme přímo podle definice. Pokud však je funkce v určitém smyslu pěkná (má primitivní funkci na intervalu, který uvnitř obsahuje interval $[a,b]$) jsou Riemannův a Newtonův integrál stejné. Proto mezi nimi nerozlišujeme, používáme jeden pojem určitý integrál a počítáme jej pomocí definice Newtonova integrálu. Obsah obrazce pod křivkou $f(x)$ je roven $$S={\int }_a^{b}f(x)\mathrm{d}x.$$

V teorii Riemannova integrálu má vzorec $${\int }_a^{b}f(x)\mathrm{d}x={[F(x)]}_a^b=F(b)-F(a)$$ postavení věty nazývané Newtonova–Leibnizova věta a je to věta udávající, jak vypočteme určitý integrál pomocí neurčitého. Zajímavé je, že v některých případech je vhodné postupovat naopak a určit neurčitý integrál pomocí integrálu určitého, což si ukážeme v následující přednášce.

Nasčítání příspěvků k celkové dráze

1. Těleso pohybující se po dobu $\mathrm{\Delta }t$ konstantní rychlostí $v$ po přímce urazí dráhu $$s=v\mathrm{\Delta }t.$$
2. Těleso pohybující se po dobu $\mathrm{\Delta }{t}_1$ konstantní rychlostí $v_1$ po přímce a poté po dobu $\mathrm{\Delta }{t}_2$ rychlostí $v_2$ urazí celkovou dráhu $$s=v_1\mathrm{\Delta }{t}_1+v_2\mathrm{\Delta }{t}_2.$$ Toto je možné zobecnit na libovolný pohyb skládající se z konečného počtu úseků, kdy se těleso pohybuje konstantní rychlostí. $$s=v_1\mathrm{\Delta }{t}_1+v_2\mathrm{\Delta }{t}_2+\cdots +v_k\mathrm{\Delta }{t}_k=\sum _{i=1}^k{v}_i\mathrm{\Delta }{t}_i$$ Příspěvek za každou část pohybu, kdy je rychlost konstantní, je $$\mathrm{\Delta }s=v\mathrm{\Delta }t,$$ kde $v$ a $\mathrm{\Delta }t$ jsou příslušná rychlost a doba pohybu, po kterou je rychlost konstantní.
3. Pokud se rychlost mění spojitě a $a$ a $b$ jsou počáteční a koncový okamžik pohybu, platí $$s={\int }_a^{b}v(t)\mathrm{d}t.$$

Volitelný text. Klikněte pro rozbalení.

Zde je ukázáno, jak vypočítat tlak na plochu, která je rozložena přes celý interval hloubky. Na každém místě totiž působí jiný tlak a příspěvky přes všechny hodnoty hloubky je nutné posčítat integrálem.

Nasčítání příspěvků k celkové síle na přehradu

Mojžíšův most (Holandsko, pevnost Fort de Roovere) je v celosvětovém měřítku unikátním mostem. Je postavený ze dřeva a zanořený do vodního příkopu okolo pevnosti tak, aby splýval s krajinou. Představme si zjednodušeně, že vodní masu drží svislá dřevěná stěna a budeme se snažit najít celkovou sílu působící na tuto stěnu tlakem vodní masy. (Ve skutečnosti most leží na dně a dno se zvedá směrem ke stěnám mostu. Google umí najít stavební plán mostu.) Délku mostu označíme $L$, výšku stěny (přesněji vzdálenost ode dna po hladinu vody) označíme $H$.

1. Tlaková síla na rovinnou plochu o obsahu $S$ vyvolaná tlakem $p$ je rovna $$F=pS.$$ Tlak v hloubce $h$ je dán vzorcem $$p=h\rho g,$$ kde $\rho$ je hustota vody a $g$ tíhové zrychlení.
2. Myšlenkově rozdělíme celou stěnu na části. Tlaková síla na celou stěnu je rovna součtu tlakových sil, které působí na jednotlivé části. Má smysl volit části tak, aby na nich byl tlak konstantní. Myšlenkově tedy stěnu rozřežeme na vodorovné pásky.
3. Na myšlený vodorovný pás, který má výšku $\mathrm{\Delta }x$ a je v hloubce $x$, působí tlak $p=x\rho g$. Obsah pásu je podle vzorce pro obsah obdélníka $\mathrm{\Delta }S=L\mathrm{\Delta }x$. Celková síla působící na tento pás je $$\mathrm{\Delta }F=p\mathrm{\Delta }S=L\rho gx\mathrm{\Delta }x.$$
4. Celkovou sílu na celou stěnu najdeme sečtením všech příspěvků. Formálně $$F=\sum L\rho gx\mathrm{\Delta }x.$$ Protože těchto příspěvků je nekonečně mnoho, sečteme je integrálem $$F={\int }_0^{H}L\rho gx\mathrm{d}x.$$
5. Po výpočtu dostáváme $$F={\int }_0^{H}L\rho gx\mathrm{d}x=L\rho g{\int }_0^{H}x\mathrm{d}x=L\rho g{\left[\frac{1}{2}{x}^2\right]}_0^H=L\rho g[\frac{1}{2}{H}^2-\frac{1}{2}{0}^2]=\frac{1}{2}L{H}^2\rho g.$$ Tento vztah je stejný, jako kdyby na celou plochu o velikosti $LH$ působila tlaková síla vyvolaná tlakem $\frac{1}{2}H\rho g$, tj. tlakem v poloviční hloubce.

Volitelný text. Klikněte pro rozbalení.

Ukázka demonstruje, jak stanovit tok trubkou. Pokud by v celém průřezu byla stejná rychlost, počítal by se tok násobením. Protože v reálném proudění rychlost směrem ke stěnám trubky klesá, je nutné celkový tok určit jako příspěvky z jednotlivých míst trubky a všechno posčítat integrálem.

Nasčítání příspěvků k celkovému toku potrubím

V předchozím příkladě jsme “krájeli” přehradu na vodorovné pásy, protože ve vodorovném směru byl konstantní parametr, který jsme potřebovali mít konstantní pro výpočet k celkovému příspěvku. V následujícím případě je obdobný parametr konstantní na kružnicích a proto budeme dělit a sčítat příspěvky pomocí mezikruží.

V potrubí o poloměru $R$ teče viskozní tekutina tak, že uprostřed má maximální rychlost a u stěn nulovou. Rychlost ve vzdálenosti $r$ od osy potrubí je dána vztahem $$v(r)=v_{max}\frac{R^2-r^2}{R^2}.$$ Střední rychlost dostaneme jako celkový tok potrubím dělený obsahem $S=\pi R^2$. Tok v případě konstantní rychlosti je součinem rychlosti a obsahu průřezu. V případě rychlosti nekonstantní rozdělíme průřez na části $\mathrm{\Delta }S$, na každé části určíme příspěvek k celkovému toku a poté vše sečteme, tj. $$Q=\sum _{\text{průřez}}v(r)\mathrm{\Delta }S$$ Protože rychlost závisí na $r$, je stejná na kružnicích a proto se jeví vhodné rozdělit průřez na mezikruží a sečíst přes poloměr těchto mezikruží. Budeme tedy integrovat přes proměnnou $r$. Vyjádření $Q$ přepíšeme do tvaru $$Q=\sum _{\text{průřez}}v(r)\frac{\mathrm{\Delta }S}{\mathrm{\Delta }r}\mathrm{\Delta }r$$ a limitním přechodem se změní suma na integrál a podíl změn na derivaci, tj. $$Q={\int }_0^{R}v(r)\frac{\mathrm{d}S}{\mathrm{d}r}\mathrm{d}r={\int }_0^R{v}_{max}\frac{R^2-r^2}{R^2}\frac{\mathrm{d}S}{\mathrm{d}r}\mathrm{d}r.$$ Odsud pomocí $S=\pi r^2$, $\frac{\mathrm{d}S}{\mathrm{d}r}=2\pi r$ dostáváme po vytknutí konstant, roznásobení závorek a integraci $$Q=\frac{2\pi v_{max}}{R^2}{\int }_0^R(R^2-r^2)r\mathrm{d}r=\frac{2\pi v_{max}}{R^2}{[R^2\frac{r^2}{2}-\frac{r^4}{4}]}_0^R=\frac{2\pi v_{max}}{R^2}\frac{R^4}{4}=\frac{v_{max}}{2}\pi R^2.$$ Tok je tedy formálně stejný, jako by voda tekla v celém průřezu rychlostí poloviční ve srovnání s maximální rychlostí ve středu trubice. Proto je $\frac{v_{max}}{2}$ nazývána střední profilová rychlost průřezu.

(Volně podle Dana Říhová a Jana Marková, Poznámky k přednáškám z Hydrauliky, přednáška č. 3.)

Volitelný text. Klikněte pro rozbalení.

Volitelná ukázka aplikace (výpočet momentu setrvačnosti pomocí integrálu)

Nasčítání příspěvků k celkovému momentu setrvačnosti tyče

• Kinetická energie tělesa o hmotnosti $m$ pohybujícího se posuvným pohybem rychlostí $v$ je dána vztahem $E=\frac{1}{2}m{v}^2$. Kinetická energie tělesa o momentu setrvačnosti $J$ pohybujícího se otáčivým pohybem úhlovou rychlostí $\omega$ je dána vztahem $E=\frac{1}{2}J{\omega }^2$. Odsud vidíme, že energie potřebná k vyvolání rotačního pohybu je úměrná momentu setrvačnosti. Moment setrvačnosti je tedy jakousi mírou odolnosti tělesa vůči silám, které se jej snaží uvést do rotačního pohybu. Zjednodušeně, větší moment setrvačnosti znamená, že těleso je stabilnější.
• Moment setrvačnosti hmotného bodu o hmotnosti $m$ vzhledem k ose otáčení vzdálené $r$ od tohoto bodu je $J=m{r}^2$. Pro soustavu hmotných bodů stačí příspěvky sečíst. Pro případ tělesa se spojitě rozloženou hmotností bychom museli “sečíst nekonečně mnoho nekonečně malých příspěvků” a proto sčítáme integrálem.

Budeme studovat rotaci tyče o hmotnosti $m$ a délce $L$ okolo osy kolmé k tyči. Nechť je tyč položena podél osy $x$ a rotuje okolo osy $y$. Kousek tyče o délce $\mathrm{\Delta }x$ má hmotnost $\frac{\mathrm{\Delta }x}{L}m$ a pokud je jeho vzdálenost od osy $y$ rovna $x$, příspěvek k celkovému momentu setrvačnosti je $$\mathrm{\Delta }J=\frac{\mathrm{\Delta }x}{L}m{x}^2=\frac{m}{L}{x}^2\mathrm{\Delta }x.$$ Celkový moment setrvačnosti je dán integrálem, ale závisí na poloze tyče vzhledem k ose otáčení.

1. Pro tyč umístěnou levým koncem v počátku dostáváme moment vzhledem k ose procházející koncem tyče ve tvaru $$J={\int }_0^L\frac{m}{L}{x}^2\mathrm{d}x=\frac{m}{L}{\left[\frac{1}{3}{x}^3\right]}_0^L=\frac{m}{L}\frac{1}{3}{L}^3=\frac{1}{3}m{L}^2.$$
2. Pro tyč umístěnou středem v počátku dostáváme moment vzhledem k ose procházející středem ve tvaru $$J={\int }_{-\frac{L}{2}}^{\frac{L}{2}}\frac{m}{L}{x}^2\mathrm{d}x=\frac{m}{L}{\left[\frac{1}{3}{x}^3\right]}_{-\frac{L}{2}}^{\frac{L}{2}}=\frac{m}{L}[\frac{1}{3}\frac{L^3}{8}-\frac{1}{3}(-1{)}^3\frac{L^3}{8}]=\frac{1}{12}m{L}^2.$$

Závěr.

• Na roztočení tyče okolo konce je potřeba více energie, než na roztočení okolo středu. Čtyřikrát více. (Z praxe víme, že s dlouhým žebřem se manipuluje nejlépe, pokud jej držíme uprostřed.)
• Tyč o konstantní délkové hustotě $\tau$ (dané použitým průřezem a materiálem) má hmotnost $m=\tau L$ a moment setrvačnosti vzhledem ke středu $$J=\frac{1}{12}\tau L^3.$$ Vidíme, že moment setrvačnosti roste dramaticky při zvětšování délky, s třetí mocninou. Proto provazochodci nosí na laně dlouhou tyč a proto při extrémních výkonech, jako je přechod Grand Canyon, bývá použita extrémně dlouhá tyč (pro Grand Canyon 9.1 metrů a 20 kilogramů, viz Nik Wallenda).

Shrnutí, hlavní myšlenky

• Někdy máme zadánu rychlost, s jakou se mění veličina a potřebujeme znát funkční předpis pro tuto veličinu, tj. hodnotu v libovolném čase. To je úloha inverzní k derivaci a řeší ji neurčitý integrál.
• Při zadané rychlosti změny není možné bez zadání výchozího stavu určit hodnotu veličiny, která se mění. Je možné vypočítat jenom změnu této veličiny za určitý časový úsek (Newtonův určitý integrál) anebo je řešení dáno až na počáteční stav vyjádřený integrační konstantou v neurčitém integrálu.
• Někdy potřebujeme veličinu, která nás zajímá, najít posčítáním nekonečně mnoha příspěvků. Toto je v situaci, kdy se “za běhu” mění parametry úlohy, například se během pohybu mění rychlost pohybu. V tomto případě používáme Riemannův určitý integrál, který je definovaný jinak než Newtonův, ale v prakticky zajímavých úlohách se počítá stejně.
• Další aplikací procesu opačného k derivování je úloha, kdy jsou vlastnosti křivky popsány pomocí derivace a hledáme rovnici pro tuto křivku. Příkladem jsou úlohy ve stavitelství a studiu materiálu (ohybová čára nosníku).