197 Kapitola 4 Aplikace vícerozměrných integrálů Integrální počet funkcí více proměnných má rozsáhlé použití jak v matematice a fyzice, tak i v oborech, které na nich staví. V této kapitole uvedeme některé geometrické a fyzikální aplikace, přičemž se omezíme jen na základní z nich. Vesměs se jedná o výpočet podobných veličin, které se vyskytují u aplikací jednoduchého určitého integrálu, nyní je však budeme umět určit v podstatně obecnější situaci (objem a hmotnost nepravidelných nehomogenních těles apod.). 4.1. Geometrické aplikace Z geometrických aplikací se zaměříme na výpočet obsahů rovinných množin, obsahů ploch v prostoru a objemů prostorových množin. 4.1.1. Míra (obsah) rovinné množiny Nechť M je rovinná měřitelná množina. Z definic 1.31 a 1.45, popřípadě z tvrzení a) věty 1.50 vyplývá, že pro míru (obsah) m2(M) této množiny platí m2(M) = / / áxáy. (4.1) M Podle důsledku 1.41 víme, že omezená množina je měřitelná právě tehdy, když její hranice má míru nula. Tuto vlastnost mají všechny „rovinné obrazce", se kterými se setkáváme v geometrii a v dalších běžných aplikacích. 198 Aplikace vícerozměrných integrálů Příklad 4.1. Vypočtěte obsah rovinné množiny M omezené křivkami xy = a2, x2 = ay, y = 2a, x = 0, kde a > 0 je konstanta. Řešení. Zjistíme nejprve průsečík křivek xy = = a2, x2 = ay.Z první z těchto rovnic vidíme, že obě souřadnice průsečíku musí být různé od nuly. Vypočítáme-li z první rovnice neznámou y a dosadíme do druhé rovnice, obdržíme x2 = a3/x, takže x3 = a3. Odtud plyne, že x = a, a z rovnic křivek snadno určíme y = a. Jediným průsečíkem uvažovaných dvou křivek je tedy bod [a, a]. Množina M je znázorněna na obr. 4.1. Obr. 4.1 Označíme-li Mi OS y S a, M2: a ^ y ^ 2a, OSxS «Jay, ' 0 ^ x ^ a2/y, máme M = Mi U M2, m2(Mi n M2) = 0. Užitím tvrzení c) věty 1.50 dostáváme m2(M) = dxdy = dxdy + / / dxdy = M Mi ay dx o dy + M2 2a r ľa ly dx dy = = *Ja 9 l a -y3'2 2 o + a2 ľln y li a = - a2 + a2 ln 2 = a + hi2 Příklad 4.2. Vypočtěte obsah množiny A omezené kardioidou, mající v polárních souřadnicích rovnici q = a{\ + cos
0 je konstanta. Řešení. Množina je znázorněna na obr. 4.2 a). Vzhledem ke způsobu zadání kardioidy je vhodné použít transformaci do polárních souřadnic. Množina A se transformuje na množinu B -TT ^ (p ^ It, 0 S Q S «U +cos
0. Řešení. Abychom si usnadnili výpočet, umístíme střed elipsy do počátku souřadnicové soustavy a osy souměrnosti do souřadnicových os (na obsah to nemá vliv) — viz 2 2 obr. 4.3. Rovnice elipsy pak bude ^ + fi = = 1. Nyní použijeme transformaci do eliptických souřadnic (3.7), konkrétně x = ag cos
0, dostáváme jako speciální případ obsah itr2 kruhu o poloměru r. Á 4.1.2. Míra (objem) měřitelné množiny v trojrozměrném prostoru Nechť A je měřitelná množina v trojrozměrném prostoru. Z definice míry v IR3 a z definice trojného integrálu vyplývá, že pro objem 1113(A) této množiny platí 1113(A) = ííí dxdy dz. (4.2) A Příklad 4.4. Vypočtěte objem tělesa A omezeného plochami 2(x2+y2) — z2 = 0 a x2 + y2 — z2 = —9, leží-li A v poloprostoru z ^ 0. Řešení. První plochou je rotační kužel s osou z, druhou rotační dvojdílný hyperboloid s osou z. To lze snadno nahlédnout pomocí řezů plochy rovinami x = 0, j = 0 a z = c, kde c G IR je konstanta. Množina A je znázorněna na obr. 4.4 a). Na obr. 4.4 b) je řez této množiny rovinou x = 0. Zadané plochy se protínají v kružnici. Její rovnici dostaneme odečtením rovnic obou kvadrik. Vyjde x2 + y2 = 9. Tedy kolmým průmětem množiny A do roviny x y je kruh K se středem v počátku o poloměru 3 (obr. 4.4 c)). Použijeme transformaci do cylindrických souřadnic. Nejprve vyjádříme průmět K v polárních souřadnicích, tj.0^
dg = j 2tzq [z] ^ +9 >ÚQ =
2tt / (g^Q2 + 9-Q2V2)dQ = Jo
2tt / q\Jq2 + 9de -2tt / £2\/2de = Jo Jo
ř?2 + 9 = ř 2^dg) = d/ £>dg> = ^ d/ 0^9, 3^ 18
"18
= TT / Vřdř -2ttV2
2tt
T
18 9
18tta/2 =
= -^ (54V2 - 27) - 18tta/2 = 18tt(a/2 - l).
▲ množina
Častý je případ, kdy těleso A, jehož objem hledáme, je elementární vzhledem k rovině xy, tj.
A = {[x, y, z]eť: [x, y] g M, f(x, y) ^ z ^ g(x, y)},
kde M c ]R2 je uzavřená měřitelná množina a/ag jsou spojité funkce na M přičemž f(x, y) S g(x, y) pro každé [x, y] g M (srovnejte obr. 2.2). Pak j< možné vzorec (4.2) pro výpočet objemu upravit pomocí Fubiniovy věty na tvai
m3(A) = / / (g(x, y) - f{x, y)) dxdy,
(4.3)
M
v němž figuruje jen dvojný integrál.
Příklad 4.5. Vypočtěte objem tělesa A určeného nerovnostmi z ^ l-x2az^ ^ y1 + 2, přičemž —l^jc^l,—l^j^l.
Řešení. Plochy, které omezují množinu A zdola resp. shora, jsou parabolické válce, jejichž povrchové přímky jsou rovnoběžné se souřadnicovou osou y resp. x — viz obr. 4.5. Kolmým průmětem množiny A do roviny xy je čtverec M: (— -l,l)x(-l,l).
Použijeme vzorec (4.3), v němž zvolíme f{x,y) = \—x2dLg{x,y) = y2+2. Vzniklý integrál vypočítáme pomocí Fubiniovy věty. Vyjde:
1113(A) = / / ((/ + 2) - (1 - x1)) dxdy =
M
4.1 Geometrické aplikace
203
Obr. 4.5
M
(x2 + y2 + l)dxdy =
dx =
v:
(x2 + y2 + l)ďy dx
Ĺ x2y + y3 3 + y i
2x . 3 3 8jc" -+3. i -1 z ,0 3"
\{2x2+í)dx =
4.1.3. Míra měřitelné množiny v «-rozměrném prostoru
Stejně jako ve dvojrozměrném a trojrozměrném prostoru plyne z definice míry a z definice integrálu pro míru měřitelné množiny M v Mn vzorec
mn(M) =
dxid^2 • • -dxn.
(4.4)
M
Příklad 4.6. Vypočtěte míru Tn «-rozměrného jehlanu (simplexu)
Mn = {[xux2, ...,xn]: x\ ^ 0, x2 ^ 0, ...,xn ^ 0, x\ +x2-\-------\-xn^h) ,
kde h > 0 je konstanta.
Řešení. Proveďme nejprve dilataci x\ = hu\,x2 = hu2,...,xn = hun. Pro jakobián J této transformace platí J = hn — viz (3.19). Množina Mn přejde po
204
Aplikace vícerozměrných integrálů
zmíněné transformaci v množinu
M* = {[u\, u2, ...,«„]: u\ ^ 0, u2 ^ 0,... ,un ^ 0, u\ + u2 H-------\-un ^ l}.
Množinu M* jakožto elementární množinu můžeme vymezit pomocí nerovností:
0Sun ^1,
0 ^ «w_i ^ 1 —un,
0 ^ Ww_2 ^ 1 -«„ —Un-\,
M*:
0 ^ «2 ^ 1 — «, 0 ^ «i ^ 1 — W,
Dostáváme pak
«3,
W3 — U2.
dxidx2 • • • dxn = • • • h du\dii2 ... dun = h an,
Mn M*
kde
an = / • • • / du\dii2 • • • dun.
M*
Označíme-li při pevném un € (0, 1)
Mn-\ = {[u\, u2,..., un-\\: u\ ^ 0, u2 ^ 0,..., un-\ ^ 0,
«i + «2 H-------V un-\ S 1 - ««},
vidíme, že
du\du2 • • • dun~\ )dww.
Provedeme-li ve vnitřním integrálu novou dilataci wi = (1 — un)v\,u2 = (1 — — un) v2,..., un-\ = (1 — un) vn-\ s jakobianem Jn-\ = (1 — un)n~x, obdržíme
i 2^ 17
ŕ11
208
Aplikace vícerozměrných integrálů
ji 4
ř3/2
3/2
17
JI
= -(17VÍ7- 1). 6 v ;
b) Jde o část hyperbolického paraboloidu. Platí f'x(x, y) = y a fý(x,y) = x, takže
S2(G) = // ^\+y2+x2dxdy =
Q
\/\ + q2 sin2 2 cos2 (p • q dgd(p =
= / / qV^ + Q2dgd(p = d(p ■ / g-y/l + £2d£ =
1 + Q2 = t2
2q dg = 2ř dt q dg = t dt
V5
= 2jt-
/•V5
= Wo"71 ř2dř =
= ^(5V5-l).
Všimněte si, že pro procvičení jsme na zcela obdobný jednoduchý určitý integrál vzhledem k proměnné q použili pokaždé poněkud odlišnou substituci. Á
4.1.5. Míra (obsah) (n — l)-rozměrné plochy v «-rozměrném prostoru
Vzorec pro (« — \)-rozměrnou míru (obsah) v «-rozměrném prostoru grafu
G = {[xux2, ...,xn-u f(xux2, ...,*„_i)] g M.n: [xux2,. ..,*„_i] g Q) funkce / o n — 1 proměnných je analogický vzorci (4.8):
/„ p re—1
• • •/ i+^2f?Mi'x2,• • •'Xn-^
1/2
dx\dx2- • -dxn-\. (4.9)
Přitom předpokládáme, že Q c Mw_1 je omezená oblast, funkce / je spojitá na Q a má spojité a ohraničené parciální derivace na oblasti Q, jejíž hranice je po částech hladká (tj. je sjednocením konečně mnoha grafů diferencovatelných funkcí n — 2 proměnných x7l,..., x/„_2, 1 ^ ji ^ • • • ^ jn-i = n — í, definovaných na kompaktních množinách, přičemž různé grafy mají společné nejvýše
4.1 Geometrické aplikace
209
body svých „okrajů"). Přesný výklad pojmů plocha v W a míra plochy v W viz např. [26] nebo [27].
Příklad 4.9. Vypočtěte obsah Sn-\ povrchu «-rozměrné koule Kn v «-rozměrném prostoru, která má daný poloměr R > 0.
Řešení. Výpočet nejprve provedeme pro « ^ 4. Určíme obsah „horní" poloviny povrchu (hranice) koule Kn. Užitím vzorce (4.9) dostáváme:
kde f{x\,x2, ...,xn-\) = V/?2 - x\ - x\---------x2n_x
1/2
dxid^2 • • -dxw_i,
^re-1 = {[*1, X2, . . . , Xn-\]
i>n — \ . „2 , „2
Xi + *o H--------hl„_, <
■\ ~ A2
n-\
R2}.
>2 __ v2 __ v2
1 x2
Protože f'x.{x\, X2,..., xn-i) = —Xj/vR2 — x
n-\
1 + S f*){XX' *2' • • • ' *w-l} = P2 Y2 '2
A Ai Ai * * *
x^_\, mame
R'
7 = 1
-1 A2
X
n-1
Dosazením do vzorce z úvodu řešení dostáváme
Sn-l =
7?
K„-i
V/?2 - JC? - x}
áx\áx2 • • -áxn-\.
-1 A2
X
n-\
Použitím transformace do sférických souřadnic
X\ = Q 20S (p sin ů\ sin Ů2
X2 = Q sin (p sin ů\ sin Ů2
x3 = Q cos ů\ sin Ů2
X4 = Q COS Ů2
Xn—3 — Q
Xfi—2 — Q
Xn—\ — Q
• sin ůn-5 sin ůn-4 sin #w_3,
• sin ůn-5 sin íV4 sin #w_3,
• sin#w_5sin#w_4sin#w_3,
• sin ůn-5 sin íV4 sin #w_3,
cos #w_5 sin #w_4 sin #w_3,
cos#„_4sin#n_3,
cos #w_3
210
Aplikace vícerozměrných integrálů
přejde množina Kn_\ y množinu
OSq < R,
n—l •
0 S 1V3 S it. Pro absolutní hodnotu jakobianu této transformace platí
\J\ = Qn~2 siní^i sin2 ů2 ■■ ■ sinw"3 ůn_3-
Užitím vztahu x2 + x\ -\-----+ x2_x = g2 dostáváme
1
Sn-l = R
1
K
v^3^
n-1
= R
R ( ľ2n
0 UO
Qn 2 sin ů\ sin2 Ů2 • • • sin" 3 ůn-3 x
x dg dep dů\ dí^2 • • • d#w-2 =
Jt n-2
sin ů\ sin2 #2 • • • sin""3 ůn-3 x
o
2jt
v^3^
x d#i • • • diV
d dg> =
"Ä ow~2 = 7? / dg I d n— 1 „ „ i?«— 1
1 iv I n—3 I I n — 2 I iv I n —1 I In
-5n_i =2tt-----------2L—J • jtL—J =-----------2L—J -jt^J,
2 (w-2)ü (« -2)!!
kde L* J značí celou část čísla x. Obsah povrchu celé koule je tedy roven číslu
pn-l
J\ I H+l I In
Sw_i =-----------2L—J • jtW.
(«-2)1!
Výsledek je zřejmě správný i pro n = 2 a n = 3. Poznamenejme ještě závěrem, že pro limitu podílu objemu «-rozměrné koule o poloměru R (viz příklad 4.7) a obsahu jejího povrchu platí
lim —— = lim — = 0.
n^oo Sn-\ re^oo n A
Poznámka 4.10. Pozorný čtenář si jistě všiml, že provedený výpočet nebyl zcela korektní. Integrovaná funkce totiž nebyla na Kn_\ ohraničená. Náš výpočet však lze ospravedlnit, předpokládáme-li, že obsah \S^_X (n — 1)-rozměrných kulových vrchlíků s osou v ose xn, středem v počátku, o středovém úhlu 2ý, kde \[r G (0, tt/2), je spojitou funkcí proměnné ý. Skutečně, podobným postupem jako dříve dostaneme
ls: 1=27t/?w-1 rsinB-2řdř.2L^U¥J
2 n~X Jo (n-m
pre—1 í-ip
= —---- / sin^ídí^LlV^J.
(«-3)!!7o
Vzhledem k tomu, že J0 sinw~21 dt -> Yn (w-2)n Pro ^ ~^ ^12—, odtud snadno zjistíme, že
1 1 , t\ , I n I I n —1 I /v I n —1 I In
-5n_i = lim -S* 1 =-----------k^^JjtL—J =-----------2L—J7TL2J.
2 f^x/2-2 w_1 (n-2)V/n (ra-2)!!
Lze ukázat, že vzorec (4.9) platí i v případě, že integrál stojící na jeho pravé straně je nevlastní (viz kapitola 5) a konvergentní (parciální derivace f'x. tedy nemusí být nutně ohraničené na Q). Srovnejte poznámku 5.27.
4.2. Fyzikální aplikace
Z fyzikálních aplikací uvedeme výpočet hmotnosti, souřadnic těžiště, momentu setrvačnosti a elektrického náboje. Odvození uvedených vzorců patří do teoretické fyziky.