MB152 Diferenciálny a integrálny počet, skupina 10, 10. cvičenie, 9.12.2020 Viacnásobný a viacrozmerný integrál Fubiniho veta je prostriedok, ktorý nám umožňuje previesť viacrozmerý integrál na viacnásobný integrál. Jej tvrdenie sa dá ilustrovať nasledujúcim obrázkom: Predpokladajme, že funkcia f(x, y) je spojitá na množine M = {(x, y) ∈ R2 | a ≤ x ≤ b, ϕ(x) ≤ y ≤ ψ(x)}, kde ϕ a ψ sú spojité funkcie na intervale [a, b]. Definujme funkciu F : [a, b] → R nasledujúcim predpisom: F(x) = ψ(x) ϕ(x) f(x, y) dy. Pre každé pevne zvolené x ∈ [a, b] je f(x, y) spojitá funkcia (jednej) premennej y definovaná na intervale [ϕ(x), ψ(x)], takže predošlá definícia má zmysel. Fubiniho veta hovorí o platnosti nasledujúcej rovnosti: M f(x, y) dx dy 1 = b a F(x) dx 2 = b a ψ(x) ϕ(x) f(x, y) dy dx 3 = b a ψ(x) ϕ(x) f(x, y) dy dx. Táto veta teda tvrdí, že objem telesa, ktoré je ohraničené rovinou xy a grafom funkcie f, vypočítame súčtom plôch rezov tohto telesa všetkými nekonečne veľa rovinami x = x0 ∈ [a, b]. Samozrejme, tvrdenie platí aj v prípade zámeny úloh premenných x a y. Analogické tvrdenie platí pre funkcie ľubovoľného počtu premenných. V takom prípade musíme byť schopní usporiadať jednotlivé premenné tak, aby množina, cez ktorú integrujeme bola vyjadriteľná v tvare xi ∈ [ai, bi], pre každú premennú xi, kde hranice tohto intervalu môžu závisieť len od premmenných, ktoré sú v spomínanom usporiadaní predchodcami premennej xi. Špeciálne, po správnom použití Fubiniho vety má posledný integrál konštantné medze. Opak svedčí o výskyte chyby. Uveďme ešte jednu poznámku. V prípade viacrozmerného integrálu chápeme diferenciály na jeho konci ako jeden symbol, a ich poradie pre naše účely nie je podstatné. Avšak pri práci s viacnásobným integrálom, špeciálne po použití Fubiniho vety, je každý diferenciál samostatný symbol a ich vzájomné poradie je podstatné. Viažu sa totiž na integračné medze zapísané pri integrálnom znaku a pripomínajú nám poradie integrovania jednotlivých premenných, teda kedy máme ktorú premennú integrovať. 1 Tvrdenie Fubiniho vety. 2 Definícia funkcie F. 3 Zátvorky sa zvyknú vynechávať. 1 Príklad 1. Vypočítajte dvojný integrál M y dx dy, kde M = {(x, y) ∈ R2 | y ≥ 0, x2 + y2 ≤ 1, (x − 1)2 + y2 ≥ 1}. Najprv nájdeme vyjadrenie množiny M, ktoré spĺňa predpoklady Fubiniho vety. Množina M je časťou prieniku kruhu so stredom v bode (0, 0) s polomerom 1 s doplnkom kruhu so stredom v bode (1, 0) s polomerom 1, ktorá leží nad osou x. Ak ju rozdelíme, na dve časti, podľa toho, či je x ≤ 0, alebo x > 0 dostaneme dve množiny, ktoré sú elementárne vzhľadom k ose x: A = M ∩ {(x, y) ∈ R2 | x ≤ 0} = {(x, y) ∈ R2 | −1 ≤ x ≤ 0, 0 ≤ y ≤ √ 1 − x2}, B = M ∩ {(x, y) ∈ R2 | x > 0} = {(x, y) ∈ R2 | 0 < x ≤ 1 2 , 1 − (x − 1)2 ≤ y ≤ √ 1 − x2}. Môžeme počítať: M y dx dy = A y dx dy + B y dx dy = 0 −1 √ 1−x2 0 y dy dx + 1 2 0 √ 1−x2 √ 1−(x−1)2 y dy dx = 0 −1 y2 2 √ 1−x2 0 dx + 1 2 0 y2 2 √ 1−x2 √ 1−(x−1)2 dx = 1 2 0 −1 1 − x2 dx + 1 2 1 2 0 1 − x2 − (1 − (x − 1)2 ) dx = 1 2 x − x3 3 0 −1 + 1 2 1 2 0 −2x + 1 dx = 1 3 + 1 2 −x2 + x 1 2 0 = 1 3 + 1 8 = 11 24 . v Príklad 2. Zameňte poradie integrácie nasledujúceho dvojnásobného integrálu: 1 0 1+ √ y 1− √ y f(x, y) dx dy. 2 Podľa Fubiniho vety platí 1 0 1+ √ y 1− √ y f(x, y) dx dy = M f(x, y) dx dy, kde M = {(x, y) ∈ R2 | 0 ≤ y ≤ 1, 1 − √ y ≤ x ≤ 1 + √ y}, teda množina M je elementárna vzhľadom k ose y. Našou úlohou je, ak je to možné, nájsť vyjadrenie množiny M, ktoré z nej urobí množinu elementárnu vzhľadom k ose x, a potom opäť použiť Fubiniho vetu, ktorá prevedie dvojný integrál M f(x, y) dx dy na dvojnásobný s tým rozdielom, že tentokrát budeme integrovať najprv premennú y a až potom premennú x. Môžeme si pomôcť obrázkom. 0 ≤y ≤ 1 1 − √ y ≤x ≤ 1 + √ y ⇔ 0 ≤y ≤ 1 1 − √ y ≤x x ≤ 1 + √ y ⇔ 0 ≤x ≤ 2 (x − 1)2 ≤y y ≤ 1 ⇔ 0 ≤x ≤ 2 (x − 1)2 ≤y ≤ 1 Preto platí 1 0 1+ √ y 1− √ y f(x, y) dx dy = 2 0 1 (x−1)2 f(x, y) dy dx. Všimnite si opačné poradie diferenciálov. v Príklad 3. Vypočítajte dvojnásobný integrál √π 2 0 √π 2 y y2 sin x2 dx dy. Hneď po prvom kroku √π 2 0 √π 2 y y2 sin x2 dx dy = √π 2 0 y2 √π 2 y sin x2 dx dy zistíme, že tadiaľto cesta nevedie, keďže funkcia sin x2 dx je vyššia transcendentná. Musíme zameniť poradie integrácie. 0 ≤y ≤ π 2 y ≤x ≤ π 2 ⇔ 0 ≤x ≤ π 2 0 ≤y ≤ x 3 Počítajme: √π 2 0 √π 2 y y2 sin x2 dx dy = √π 2 0 x 0 y2 sin x2 dy dx = √π 2 0 y3 sin x2 3 x 0 dx = √π 2 0 x3 sin x2 3 dx t = x2 dt = 2x dx 0 0 π 2 π 2 = 1 6 π 2 0 t sin t dt u = sin t u = − cos t v = t v = 1 = 1 6 [−t cos t] π 2 0 + 1 6 π 2 0 cos t dt = 1 6 [sin t] π 2 0 = 1 6 . v Príklad 4. Vypočítajte plochu útvaru M, ktorý je ohraničený krivkami x = 0, y = 1/x, y = 4x a y = 8. Obsah útvaru alebo presnejšie miera množiny M je daná vzťahom m(M) = M 1 dx dy. Preto nám stačí zistiť, či sa dá množina M poskladať z elementárnych množín vzhľadom k niektorej zo súradnicových osí a následne použiť Fubiniho vetu. Pomôžeme si obrázkom. Krivky y = 1 x a y = 4x sa pretínajú v bode (1 2 , 2). Množinu M môžeme vyjadriť ako zjednotenie M = A ∪ B, kde A = {(x, y) ∈ R2 | 0 ≤ y ≤ 2, 0 ≤ x ≤ y 4 }, B = {(x, y) ∈ R2 | 2 ≤ y ≤ 8, 0 ≤ x ≤ 1 y }. Obidve z uvedených množín sú elementárne vzhľadom k ose y a ich prienik, úsečka, je množinou nulovej miery, takže sa môžeme pustiť do integrovania. M 1 dx dy = A 1 dx dy + B 1 dx dy = 2 0 y 4 0 1 dx dy + 8 2 1 y 0 1 dx dy 4 = 2 0 [x] y 4 0 dy + 8 2 [x] 1 y 0 dy = 2 0 y 4 dy + 8 2 1 y dy = y2 8 2 0 + [log y]8 2 = 1 2 + log 4. v Príklad 5. Vypočítajte trojný integrál M (x + y)z dx dy dz, kde M je osmina gule x2 + y2 + z2 ≤ 1 ležiaca v prvom oktante. Množinu M môžeme určiť prostredníctvom nasledujúcich nerovností: M = {(x, y, z) ∈ R3 | 0 ≤ x ≤ 1, 0 ≤ y ≤ √ 1 − x2, 0 ≤ z ≤ 1 − x2 − y2}. Množina M je elementárna vzhľadom k rovine xy, čo nám opäť umožní použitie Fubiniho vety: M (x + y)z dx dy dz = 1 0 √ 1−x2 0 √ 1−x2−y2 0 (x + y)z dz dy dx Pozor na poradie diferenciálov. = 1 0 √ 1−x2 0 (x + y)z2 2 √ 1−x2−y2 0 dy dx = 1 2 1 0 √ 1−x2 0 (x + y)(1 − x2 − y2 ) dy dx = 1 2 1 0 √ 1−x2 0 x − x3 + y − x2 y − xy2 − y3 dy dx = 1 2 1 0 (x − x3 )y + (1 − x2 ) y2 2 − x y3 3 − y4 4 √ 1−x2 0 dx = 1 2 1 0 x(1 − x2 ) √ 1 − x2 + (1 − x2 ) 1 − x2 2 − x (1 − x2 ) √ 1 − x2 3 − (1 − x2 )2 4 dx = 1 2 1 0 2x(1 − x2 ) √ 1 − x2 3 + (1 − x2 )2 4 dx = 1 3 1 0 x(1 − x2 ) √ 1 − x2 dx t = 1 − x2 dt = −2x dx 0 1 1 0 + 1 8 1 0 1 − 2x2 + x4 dx = 1 6 1 0 t 3 2 dt + 1 8 x − 2x3 3 + x5 5 1 0 = 1 6 2 5 t 5 2 1 0 + 1 8 1 − 2 3 + 1 5 = 2 30 + 1 8 · 8 15 = 2 15 . v 5 Príklad 6. Vypočítajte objem kolmého ihlanu so štvorcovou podstavou. Označme symbolom v výšku ihlanu a symbolom 2a dĺžku strany jeho podstavy. Umiestnime tento ihlan do priestoru tak, aby vrcholy jeho podstavy boli body (a, a, 0), (−a, a, 0), (−a, −a, 0) a (a, −a, 0) a jeho vrchol mal súradnice (0, 0, v). Ihlan môžeme rozdeliť na štyri rovnaké časti tak, že ho rozrežeme rovinami , v ktorých leží vrchol ihlanu a jedna z uhlopriečok jeho podstavy. Vypočítame objem jednej takejto štvrtiny. Jedna z týchto štvrtín môže byť popísaná nerovnosťami 0 ≤ x ≤ a, −x ≤ y ≤ x a z-ová súradnica je zdola ohraničená 0 a zhora stenou ihlanu. Preto potrebujeme túto stenu vyjadriť ako graf funkcie z(x, y) = f(x, y). Vzhľadom k tomu, že sa jedná o rovinu, má funkcia z tento tvar z = bx + cy + d. Koeficienty b, c a d dopočítame zo sústavy lineárnych rovníc, ktoré vzniknú dosadením súradníc bodov (a, a, 0), (a, −a, 0), (0, 0, v) ležiacich na tejto stene. 0 = ba + ca + d 0 = ba − ca + d v = d ⇒ d = v b = − v a c = 0 Označme symbolom M skúmanú štvrtinu ihlanu. Zistili sme, že M je určená nasledujúcimi nerovnosťami: M = (x, y, z) ∈ R3 | 0 ≤ x ≤ a, −x ≤ y ≤ x, 0 ≤ z ≤ − v a x + v Miera tejto množiny je opäť daná integrálom funkcie 1 cez množinu M: m(M) = M 1 dx dy dz = a 0 x −x − v a x+v 0 1 dz dy dx = a 0 x −x [z] − v a x+v 0 dy dx = a 0 x −x − v a x + v dy dx = a 0 − v a x + v [y]x −x dx = a 0 2x − v a x + v dx = a 0 −2x2 v a + 2xv dx = −2x3 v 3a + x2 v a 0 = −2a3 v 3a + a2 v = a2 v 3 . Objem ihlanu je preto 4a2v 3 . Ak namiesto polovice dĺžky strany podstavy ihlanu a, použijeme dĺžku b = 2a, teda a = b 2 dostaneme známy vzorec 1 3 b2 v. v Transformácia integrálu Náš jediný prostriedok, ako vypočítať viacrozmerný integrál je jeho redukcia prostredníctvom Fubiniho vety na integrál viacnásobný. Táto veta požaduje špeciálny tvar množiny M, cez ktorú integrujeme. Jedným zo spôsobov ako vyhovieť tejto podmienke je prechod k novému 6 súradnicovému systému, v ktorom má množina M jednoduchší a hlavne požadovaný tvar. Tomuto postupu sa hovorí transformácia súradníc, respektíve integrálu. Zdôraznime, že zmenu súradnicového systému uskutočňujeme v prvom rade kvôli oboru integrácie a až v druhom rade kvôli zjednodušeniu integrovanej funkcie. Pri zmene súradníc nesmieme zabudnúť na absolútnu hodnotu jakobiánu použitej transformácie. Príklad 7. Pomocou vhodnej transformácie vypočítajte dvojný integrál M √ xy dx dy, kde M je množina ohraničená krivkami y2 = x, y2 = 2x, xy = 1 a xy = 2. Prepíšme uvedené nerovnosti do nasledujúceho tvaru: 1 ≤ y2 x ≤ 2, 1 ≤ xy ≤ 2. To nám ponúka možnosť zaviesť nové súradnice u a v nasledujúcim spôsobom: u = y2 x , v = xy. V takom prípade by množina, cez ktorú integrujeme, bola určená nerovnosťami 1 ≤ u ≤ 2 a 1 ≤ v ≤ 2 spĺňajúcimi prepdoklady Fubiniho vety. Vyjadrime staré premenné v tých nových: uv = y3 ⇒ y = 3 √ uv, v = xy ⇒ x = v y = v 3 √ uv = 3 v2 u . Spočítajme Jakobián transformácie: det xu xv yu yv = det −1 3 3 v2 u4 2 3 3 1 uv 1 3 3 v u2 1 3 3 u v2 = − 1 3u = 0, pre 1 ≤ u ≤ 2, 1 ≤ v ≤ 2, takže navrhnutú transformáciu skutočne môžeme použiť. M √ xy dx dy = [1,2]×[1,2] √ v · − 1 3u du dv = 2 1 2 1 √ v · − 1 3u du dv = 2 1 2 1 √ v 3u du dv = 1 3 2 1 √ v[log u]2 1 dv = 1 3 2 1 √ v log 2 dv = log 2 3 · 2 3 v 3 2 2 1 = 2 log 2 9 (2 √ 2 − 1). v Príklad 8. Vypočítajte dvojný integrál M e−x2−y2 dx dy, kde M = {(x, y) ∈ R2 | x2 + y2 ≤ 1, x ≥ 0}. 7 Vzhľadom k tomu, že množina cez ktorú integrujeme je polkruh nachádzajúci sa napravo od osi y so stredom v počiatku, má zmysel transformovať integrál do polárnych súradníc r a ϕ, ktoré sú dané rovnosťami x = r sin ϕ a y = r cos ϕ. V nich je daná množina určená nerovnosťami −π 2 ≤ ϕ ≤ π 2 a 0 ≤ r ≤ 1. Pripomeňme, že jakobián tejto transformácie je r. Môžeme počítať: M e−x2−y2 dx dy = [−π 2 , π 2 ]×[0,1] e−r2 |r| dϕ dr = π 2 −π 2 1 0 re−r2 dr dϕ t = r2 dt = 2r dr 0 0 1 1 = 1 2 π 2 −π 2 1 0 e−t dt dϕ = 1 2 π 2 −π 2 −e−t 1 0 dϕ = 1 2 π 2 −π 2 1 − e−1 dϕ = 1 − e−1 2 [ϕ] π 2 −π 2 = π(1 − e−1 ) 2 . v Príklad 9. Vypočítajte obsah oblasti M nachádzajúcej sa v štvrťrovine −x ≤ y ≤ x určenej nerovnosťou (x2 + y2 )3 ≤ (x2 − y2 )2 . Štvrťrovina −x ≤ y ≤ x je v polárnych súradniciach určená nerovnosťami −π 4 ≤ ϕ ≤ π 4 . Nerovnosť (x2 + y2 )3 ≤ (x2 − y2 )2 má v tejto štvrťrovine v polárnych súradniciach tiež jednoduchší tvar: (x2 + y2 )3 ≤ (x2 − y2 )2 ⇔ (r2 )3 ≤ (r2 cos2 ϕ − r2 sin2 ϕ)2 ⇔ r6 ≤ r4 (cos2 ϕ − sin2 ϕ)2 ⇔ r2 ≤ (cos 2ϕ)2 ⇔ r ≤ cos 2ϕ. Uvedená množina teda má v polárnych súradniciach tvar vyhovujúci prepdokladom Fubiniho vety, takže jej plochu môžeme vypočítať ako integrál cez túto množinu z funkcie 1: M 1 dx dy = {(r,ϕ)|− π 4 ≤ϕ≤π 4 , 0≤r≤cos 2ϕ} |r| dr dϕ = π 4 −π 4 cos 2ϕ 0 r dr dϕ = π 4 −π 4 r2 2 cos 2ϕ 0 dϕ = π 4 −π 4 (cos 2ϕ)2 2 dϕ = π 4 −π 4 cos 4ϕ + 1 4 = sin 4ϕ 16 + ϕ 4 π 4 −π 4 = π 16 + π 16 = π 8 . Takže plocha obrazca, ktorý ohraničuje krivka (x2 + y2 )3 = (x2 − y2 )2 v celej rovine R2 je π 2 . v Príklad 10. Vypočítajte integrál M x2 + y2 dx dy dz, kde M je množina, ktorá vznikne rotáciou obdĺžnika s vrcholmi (1, 0, 0), (2, 0, 0), (1, 0, 2) a (2, 0, 2) okolo osi z. 8 Obr. 1: Príklad 9 Vo valcových súradniciach4 má množina, cez ktorú integrujeme, veľmi jednoduchý tvar: 0 ≤ ϕ ≤ 2π, 1 ≤ r ≤ 2 a 0 ≤ z ≤ 2. Pripomeňme, že jakobián transformácie do valcových súradníc je r. Môžeme počítať: M x2 + y2 dx dy dz = [0,2π]×[1,2]×[0,2] r2 · |r| dϕ dr dz = 2π 0 2 0 2 1 r3 dr dz dϕ = 2π 0 2 0 r4 4 2 1 dz dϕ = 15 4 2π 0 2 0 dz dϕ = 15 4 2π 0 [z]2 0 dϕ = 30 4 2π 0 dϕ = 30 4 [ϕ]2π 0 = 30 4 · 2π = 15π. Vypočítali ste moment zotrvačnosti telesa M vzhľadom k jeho ose. v Príklad 11. Určte objem telesa M v R3 , ktoré je ohraničené časťou kúžeľa x2 + y2 = (z − 2)2 a paraboloidom x2 + y2 = 4 − z. Upravme uvedené rovnosti tak, aby sme boli schopní si ich lepšie predstaviť. Paraboloid je grafom funkcie z = 4−x2 −y2 , čiže je otočený hore bruchom a jeho vrchol je posunutý do bodu [0, 0, 4]. Kúžeľ je daný rovnicou z = 2 ± x2 + y2, jeho vrchol je teda v bode [0, 0, 2]. Nájdeme objem telesa M, ktoré je ohraničené paraboloidom a len tou „kladnou časťou kúžeľa, ktorej 4 x = r cos ϕ, y = r sin ϕ, z = z. Okrem materiálov z prednášky sú tieto súradnice popísané aj s obrázkom v [1] na strane 74 alebo v [3] na strane 140. 9 hrot smeruje dole, teda plochou, ktorá je grafom funkcie z = 2 + x2 + y2. Teleso M je zrejme rotačné, konkrétne, vzniklo rotáciou istej plochy v rovine, povedzme zx, okolo osi z. Pre jeho popis bude teda vhodné zvoliť valcové súradnice. Ohraničenie pre súradnice ϕ a z je zrejmé: 0 ≤ ϕ ≤ 2π a 2 + x2 + y2 ≤ z ≤ 4 − x2 − y2 ⇔ 2 + r ≤ z ≤ 4 − r2 . Ohraničenie pre premennú r dostaneme tak, že nájdeme polomer valca, v ktorom sa celé uvažované teleso nachádza. K tomu nám pomôže prienik plôch x2 +y2 = (z −2)2 a x2 +y2 = 4−z. Z druhej rovnice môžeme dosadiť do prvej, čo nás privádza k rovnici 4 − z = (z − 2)2 ⇔ z2 − 3z = 0 ⇔ z ∈ {0, 3}. To nás privádza k dvom kružniciam z = 3, x2 + y2 = 1 a z = 0, x2 + y2 = 4. Nás zaujíma len tá prvá z nich. Jej polomer je 1, a preto posledné nerovnosti, ktoré určujú teleso M sú 0 ≤ r ≤ 1. Môžeme integrovať: m(M) = M 1 dx dy dz = (r,ϕ,z) 0≤ϕ≤2π, 0≤r≤1, 2+r≤z≤4−r2 |r| dr dϕ dz = 2π 0 1 0 4−r2 2+r r dz dr dϕ = 2π 0 1 0 r[z]4−r2 2+r dr dϕ = 2π 0 1 0 r(4 − r2 − 2 − r) dr dϕ = 2π 0 1 0 −r3 − r2 + 2r dr dϕ = 2π 0 − r4 4 − r3 3 + r2 1 0 dϕ = 5 12 2π 0 dϕ = 5 12 · 2π = 5 6 π. v Príklad 12. Vypočítajte integrál M yz dx dy dz, kde M = {(x, y, z) ∈ R3 | x2 + y2 + z2 ≤ 1, z ≤ − x2 + y2}. Vzhľadom k tomu, že množina, cez ktorú integrujeme, je časťou gule, má zmysel transformovať integrál do sférických súradníc5 . Nerovnosť x2 +y2 +z2 ≤ 1 určuje guľu so stredom v počiatku 5 x = r cos ϕ sin ϑ, y = r sin ϕ sin ϑ, z = r cos ϑ. Okrem materiálov z prednášky sú tieto súradnice popísané aj s obrázkom v [1] na strane 90 alebo v [3] na strane 141. 10 a polomerom 1. My z nej však berieme len tie body, ktoré ležia pod kúžeľom z = − x2 + y2. Takáto množina je určená nasledujúcimi nerovnosťami: 0 ≤ ϕ ≤ 2π, 0 ≤ r ≤ 1 a 3 4 π ≤ ϑ ≤ π. Pripomeňme, že ϑ je uhol, ktorý zviera spojnica bodu s počiatkom s kladnou časťou osi z. Jakobián transformácie do sférických súradníc je −r2 sin ϑ. Môžeme počítať: M yz dx dy dz = [0,2π]×[0,1]×[ 3 4 π,π] r sin ϕ sin ϑ · r cos ϑ −r2 sin ϑ dϕ dr dϑ = 2π 0 1 0 π 3 4 π r3 sin ϕ sin2 ϑ cos ϑ dϑ dr dϕ. Ďalej by sme kľudne mohli pokračovať rovnako ako v predošlých príkladoch. Toto je však vhodná situácia na ukážku jedného užitočného pravidla. Predpokladajme, že počítame viacnásobný integrál s konštantnými medzami a navyše premenné v integrande vieme separovať, teda máme dočinenia s nasledujúcim integrálom: x2 x1 y2 y1 z2 z1 f(x)g(y)h(z) dz dy dx. Potom platí x2 x1 y2 y1 z2 z1 f(x)g(y)h(z) dz dy dx = x2 x1 f(x) dx · y2 y1 g(y) dy · z2 z1 h(z) dz. Takéto pravidlo funguje pre ľubovoľný počet dimenzií. V našom prípade sa dostaneme k súčinu nasledujúcich troch integrálov: = 2π 0 sin ϕ dϕ · 1 0 r3 dr · π 3 4 π sin2 ϑ cos ϑ dϑ t = sin ϑ dt = cos ϑ dϑ 3 4 π √ 2 2 π 0 = [− cos ϕ]2π 0 =0 · r4 4 1 0 · − √ 2 2 0 t2 dt = 0 11 Všimnime si, že množina, cez ktorú sme integrovali bola symetrická vzhľadom k rovine zx a že integrovaná funkcia f(x, y, z) bola nepárna v premennej y, teda spĺňala rovnosť f(x, −y, z) = −f(x, y, z). To znamená, že počítaný integrál musel byť nulový. Evidentne občasné zamyslenie môže ušetriť kopec roboty. v Príklad 13. Vypočítajte objem guľového odseku M, ktorý odrezáva rovina z = 1 z gule x2 + y2 + z2 ≤ 2. V prvom rade vypočítame prienik guľovej plochy x2 + y2 + z2 = 2 s rovinou z = 1. Tým je kružnica x2 + y2 = 1 vo výške z = 1, ktorej polomer je 1. Vďaka tomu vieme, že body uvažovaného odseku spĺňajú nerovnosti 0 ≤ ϕ ≤ 2π a 0 ≤ ϑ ≤ π 4 . Stačí vyšetriť súradnicu r. Keďže sa jedná o guľový odsek, je zrejme zhora ohraničená polomerom gule, teda √ 2. Dolné ohraničenie závisí od odklonu od osi z, teda od súradnice ϑ, o čom sa ľahko presvedčíme použitím akéhokoľvek rezu rovinou, v ktorej leží os z. Pre daný uhol ϑ je dolné ohraničenie súradnice r preponou v pravouhlom trojuholníku, v ktorom má priľahlá odvesna k uhlu ϑ dĺžku 1. Z vlastností goniometrických funkcií dostávame nerovnosť 1 cos ϑ ≤ r. Môžeme integrovať: M 1 dx dy dz =    (r,ϕ,ϑ) 0≤ϕ≤2π, 0≤ϑ≤π 4 , 1 cos ϑ ≤r≤ √ 2    −r2 sin ϑ dr dϕ dϑ = 2π 0 π 4 0 √ 2 1 cos ϑ r2 sin ϑ dr dϑ dϕ = 2π 0 π 4 0 sin ϑ r3 3 √ 2 1 cos ϑ dϑ dϕ = 2π 0 π 4 0 sin ϑ 2 √ 2 3 − 1 3 cos3 ϑ dϑ dϕ = 2π 0 dϕ · π 4 0 sin ϑ 2 √ 2 3 − 1 3 cos3 ϑ dϑ = 2π 2 √ 2 3 π 4 0 sin ϑ dϑ − 1 3 π 4 0 sin ϑ cos3 ϑ dϑ t = cos ϑ dt = − sin ϑ dϑ 0 1 π 4 √ 2 2 = 2π 2 √ 2 3 [− cos ϑ] π 4 0 − 1 3 1 √ 2 2 1 t3 dt = 4π √ 2 3 − √ 2 2 + 1 + 2π 3 1 2t2 1 √ 2 2 = 4π √ 2 3 − √ 2 2 + 1 + 2π 3 1 2 − 1 = π 3 (4 √ 2 − 5). Tento príklad bol vyriešený v Matematika drsne a svižne na strane 467 iným spôsobom. v 12 Príklad 14. Vypočítajte moment zotrvačnosti plného torusu M jednotkovej hustoty pri rotácii okolo jeho osi. V prvom rade vhodne umiestnime torus do priestoru. Predpokladajme, že vznikol rotáciou kruhu ležiaceho v rovine xz so stredom v bode [a, 0, 0] a polomerom b, kde b < a, okolo osi z. V takom prípade je os z jeho osou, a preto je moment zotrvačnosti daný nasledujúcim integrálom: M x2 + y2 dx dy dz. Nie je ťažké sa presvedčiť o tom, že aj vo valcových aj v sférických súradniciach má torus tvar, ktorý nám po transformácii integrálu umožní použiť Fubiniho vetu. V obidvoch prípadoch však v integrandoch vzniknú nepríjemné odmocniny. Skúsme preto použiť iný, vhodnejší súradnicový systém. Nové súradnice označíme symbolmi ϕ, r a ϑ. Význam symbolu ϕ je rovnaký ako v prípade valcových a sférických súradníc, teda je to odchýlka projekcie spojnice bodu (x0, y0, z0) s počiatkom do roviny xy s kladnou časťou osi x proti smeru hodinových ručičiek. Uvažujme polrovinu P určenú bodom (x0, y0, z0) a osou z. Táto polrovina pretína kružnicu z = 0, x2 +y2 = a2 v práve jednom bode, povedzme p. Symboly r a ϑ reprezentujú polárne súradnice v polrovine P so stredom v bode p. Táto geometrická predstava nám umožní odvodiť vzťah medzi starými kartézskymi súradnicami (x, y, z) a novými súradnicami (ϕ, r, ϑ) a v konečnom dôsledku aj popísať torus M v týchto nových súradniciach: x = (a + r cos ϑ) cos ϕ y = (a + r cos ϑ) sin ϕ z = r sin ϑ M = g({(ϕ, r, ϑ) | 0 ≤ ϕ ≤ 2π, 0 ≤ r ≤ b, 0 ≤ ϑ ≤ 2π}), kde g označuje transformáciu do nových súradníc. Pred tým, než sa pustíme do integrovania, musíme ešte spočítať determinant jakobiánu použitej transformácie: det   xϕ xr xϑ yϕ yr yϑ zϕ zr zϑ   = det   −(a + r cos ϑ) sin ϕ cos ϑ cos ϕ −r sin ϑ cos ϕ (a + r cos ϑ) cos ϕ cos ϑ sin ϕ −r sin ϑ sin ϕ 0 sin ϑ r cos ϑ   = − sin ϑ −(a + r cos ϑ) sin ϕ −r sin ϑ cos ϕ (a + r cos ϑ) cos ϕ −r sin ϑ sin ϕ + r cos ϑ −(a + r cos ϑ) sin ϕ cos ϑ cos ϕ (a + r cos ϑ) cos ϕ cos ϑ sin ϕ = − sin ϑ r(a + r cos ϑ) sin ϑ sin2 ϕ + r(a + r cos ϑ) sin ϑ cos2 ϕ + r cos ϑ −(a + r cos ϑ) cos ϑ sin2 ϕ − (a + r cos ϑ) cos ϑ cos2 ϕ = −r(a + r cos ϑ) sin2 ϑ − r(a + r cos ϑ) cos2 ϑ = −r(a + r cos ϑ). Môžeme integrovať: M x2 + y2 dx dy dz = [0,2π]×[0,b]×[0,2π] (a + r cos ϑ)2 |−r(a + r cos ϑ)| dϕ dr dϑ 13 = 2π 0 b 0 2π 0 r(a + r cos ϑ)3 dϑ dr dϕ = 2π 0 b 0 2π 0 a3 r + 3a2 r2 cos ϑ + 3ar3 cos2 ϑ + r4 cos3 ϑ dϑ dr dϕ = 2π 0 b 0 a3 rϑ + 3a2 r2 sin ϑ + 3ar3 sin 2ϑ 4 + ϑ 2 + r4 sin ϑ + sin3 ϑ 3 2π 0 dr dϕ = 2π 0 b 0 a3 rϑ + 3ar3 ϑ 2 2π 0 dr dϕ = 2π 0 b 0 2πa3 r + 3πar3 dr dϕ = 2π 0 πa3 r2 + 3πar4 4 b 0 dϕ = 2π 0 πa3 b2 + 3πab4 4 dϕ = 2π πa3 b2 + 3πab4 4 = 2π2 b2 a a2 + 3b2 4 . v Literatúra [1] PLCH, Roman, ŠARMANOVÁ, Petra a SOJKA, Petr. Integrální počet funkcí více proměnných, https://is.muni.cz/do/rect/el/estud/prif/js12/integral/publikace/ Integraly_funkce_promenne_v2.pdf [2] HILSCHER, Roman Šimon, HASIL, Petr, VESELÝ, Michal a ZEMÁNEK, Petr. Přednášky z matematické analýzy na FI, http://www.math.muni.cz/~hasil/Data/CZ/Teach/MU/ MB152/prednasky_MB152.pdf [3] KALAS, Josef a KUBEN, Jaromír. Integrální počet funkcí více proměnných 14