Plošný a křivkový integrál Jan Malý, Luboš Pick a Miroslav Zelený Obsah 1. Parametrické plochy 2. Křivkový integrál 3. Elementy teorie pole 4. Plošný integrál kodimenze 1 5. Věta o divergenci 6. Stokesova věta 7. Důkaz věty o divergenci a Stokesovy věty 8. Praktické hledání parametrizace a určování orientace
1 4 7 9 12 14 15 20
1. Parametrické plochy 1.1. Definice (Plocha). Plocha, přesněji k-rozměrná plocha v Rn , je množina, kterou lze lokálně popsat jako množinu všech řešení regulární soustavy n−k nelineárních rovnic. Tedy řekneme, že neprázdná množina G ⊂ Rn je k-rozměrná plocha, jestliže ke každému bodu x ∈ G existuje okolí U tak a C 1 zobrazení g : U → Rn−k tak, že g 0 má hodnost n−k na G a G ∩ U = {x : g(x) = 0} ∩ U. 1.2. Definice (Parametrická plocha). Plošný integrál bývá integrál přes plochu nebo parametrickou plochu. Zatímco plochy jsou množiny, parametrické plochy se definují jako zobrazení. Budeme definovat k-rozměrnou parametrickou plochu v Rn , k ≥ 1, jako spojitě diferencovatelné zobrazení Φ otevřené množiny G ⊂ Rk do Rn . Množinu G nazveme referenčním oborem plochy Φ. Parametrické plochy jsou definované jako zobrazení, ale intuitivně je často vnímáme jako množiny, tj. parametrickou plochu Φ vnímáme jako množinu Φ(G). Při takové intuitivní představě je třeba zachovávat opatrnost, například v definici parametrické plochy jsme nepožadovali prostotu, tj. plocha se může “křížit sama se sebou” nebo někde dokonce třeba “zdvojit”. Nějčastěji se však prostota objeví v dodatečných předpokladech. Je-li Φ : G → Rn k-rozměrná parametrická plocha, často se stává, že Φ(G) je k-rozměrná plocha, ale nemusí to být a ne každá plocha se dá “parametrizovat”. 1.3. Definice (Gramův determinant). Nechť Φ : G → Rn je k-rozměrná parametrická plocha v Rn . Matice Φ0 (t) se nazývá Jacobiho matice zobrazení Φ v bodě t. Máme k ∂Φ ∂Φ 0 T 0 Φ (t) Φ (t) = (t) · (t) . ∂ti ∂tj i,j=1 Determinant z této matice se nazývá Gramův determinant, jeho odmocnina se používá jako “jakobián” pro plošné integrály druhého druhu a značí |JΦ(t)|. Tedy q |JΦ(t)| := Φ0 (t)T Φ0 (t). Symbol | . . . | zde je použit k zdůraznění faktu, že jde o nezápornou veličinu, na rozdíl od “obyčejného” objemového jakobiánu. Také si lze správně myslet, že samotnému výrazu JΦ lze také přiřadit smysl, tím se však zde nebudeme zabývat. Totiž, jakobián parametrické plochy je složitý algebraický objekt, jehož souřadnice jsou minory k × k Jacobiho matice. 1.4. Definice (Regulární parametrická plocha). Říkáme, že k-rozměrná parametrická plocha Φ : G → Rn je regulární, jestliže Φ0 (t) má v každém bodě t ∈ G hodnost k. To je právě když |JΦ(t)| = 6 0. 1.5. Problém parametrizace. Při počítání obsahu jednotkového kruhu B = {x ∈ R2 : |x| < 1} používáme s výhodou polární souřadnice Φ : G → R2 , kde Φ(r, α) = (r cos α, r sin α). Vycházíme-li z přirozených požadavků, že zobrazení Φ má být prosté a G má být otevřená množina, vychází nám optimální volba G = (0, 1) × (−π, π). Přitom se musíme vzdát třetího přirozeného požadavku, totiž aby zobrazení Φ bylo na B. Místo toho se smíříme s tím, že množina B \ Φ(G) je lebesgueovsky nulová. 1
Se stejným problémem se s ještě větší naléhavostí setkáváme u parametrizací k-rozměrných ploch v Rn . Zatímco při výpočtu obsahu kruhu nás nikdo nenutí používat polární souřadnice, pro kulovou sféru S v R3 lze dokázat, že žádná prostá regulární parametrická plocha nemá S jako obor hodnot, viz. cvičení 1.3. Proto je vhodné zavést třídu množin, které lze považovat za zanedbatelné pro k-rozměrnou integraci. 1.6. Definice (Nulové množiny). Nechť 0 ≤ k ≤ n. (V aplikacích je zpravidla k celé číslo, ale v definici připouštíme i necelé hodnoty k.) Řekneme, že množina N ⊂ Rn je k-nulová, jestliže pro každé ε > 0 existují koule B(xj , rj ), j ∈ N, tak, že [ X N⊂ B(xj , rj ) a rjk < ε. j
j
Jako příklady k-nulových množin slouží např. plochy nižší dimenze a jejich spočetná sjednocení, nebo obrazy Φ(A), kde A je Lebesgueovsky k-nulová a Φ je k-rozměrná parametrická plocha. Je-li N k-nulová, pak všechny její k-rozměrné projekce jsou Lebesgueovsky k-nulové. Pro k = n pojmy k-nulovosti a lebesgueovské nulovosti splývají. Viz. cvičení 1.4–1.8 1.7. Definice (Parametrizace). Nechť M ⊂ Rn je libovolná množina (typicky k-rozměrná plocha) a Φ : G → Rn je k-rozměrná parametrická plocha. Řekneme, že Φ je (k-rozměrná) lokální parametrizace M , jestliže Φ je prosté, otevřené vzhledem k M , regulární a Φ(G) ⊂ M . (Zobrazení Φ : G → Rn je otevřené vzhledem k M , pokud obraz každé otevřené množiny je relativně otevřenou podmnožinou M . Pak Φ(G) je otevřená a inverzní zobrazení je spojité vzhledem k Φ(G).) Řekneme-li, že Φ je globální parametrizace M , znamená to, že navíc Φ(G) = M . Užitečný kompromis mezi lokální a globální parametrizací je zobecněná parametrizace, to je taková lokální parametrizace M , že M \ Φ(G) je k-nulová množina. Pokud existuje (lokální, globální, zobecněná) parametrizace množiny M , říkáme též, že M připouští (lokální, globální, zobecněnou) parametrizaci. 1.8. Věta. Nechť G ⊂ Rn je k-rozměrná plocha, k ≥ 1. Potom (a) ke každému bodu z ∈ G existuje okolí U tak, že U ∩ G přípouští globální parametrizaci (která je ovšem lokální parametrizace G). (b) G připouští k-rozměrnou zobecněnou parametrizaci. Důkaz. (a) Zvolme z ∈ G, najděme okolí W bodu z a spojitě diferencovatelnou funkci g : W → Rn−k tak, že {x ∈ W : g(x) = 0} = W ∩ G a g 0 (z) má hodnost n − k. Potom existuje n−k sloupců matice g 0 (z), které jsou lineárně nezávislé. BÚNO se jedná o sloupce s pořadovými čísly k + 1, k + 2, . . . , n. Podle věty o implicitních funkcích existují okolí G bodu (z1 , . . . , zk ), okolí V bodu (zk+1 , . . . , zn ) a spojitě diferencovatelná funkce ξ : G → V tak, že C × V ⊂ W , G ∩ (G × V ) = {(t, ξ(t)) : t ∈ G}. Potom t 7→ (t, ξ(t)) : t ∈ G je hledaná lokální parametrizace. (b) Podle části (a), ke každému bodu z ∈ G najdeme lokální parametrizaci Φz : Gz → G a s ∈ Gz tak, že Φz (s) = z. Zúžíme každé Φz na kouli B(s, r) takovou, že B(s, r) ⊂ Gz , takto zúžené zobrazení označíme Ψz a jeho definiční obor H z . Jelikož G je lokálně kompaktní, najdeme spočetný systém lokálních parametrizací Φj : Gj → G a jejich zúžení Ψj : Hj → G tak, že ∂Hj ⊂ Gj jsou k-nulové množiny (neboť to jsou sféry) a sjednocení množin Ψj (Gj ) pokrývá G. Pro každé j položme [ Aj = H j \ Ψ−1 j (Ψi (H i )). i<j
Potom každý bod z ∈ G je v obraze Ψj (Aj ) pro nejvýše jedno j, přitom pokud není v žádném z těchto obrazů, pak je v obrazu některé z množin ∂Hi a tento obraz je k-nulový (viz. cvičení 1.6). Nyní každou z parametrizací Ψj bAj složíme s vhodným posunutím κj prostoru Rk tak, aby výsledné definiční obory S 0 0 A0j = κ−1 j (Aj ) byly po dvou disjunktní. Položíme G = j Aj a Φ(t) = Φj (κj (t)), t ∈ Aj . Potom Φ : G → G je hledaná zobecněná parametrizace. 1.9. Příklad. Sférické souřadnice Φ : G → R3 , Φ = (x, y, z), x = cos γ cos α y = cos γ sin α , z = sin γ
(α, γ) ∈ G = (−π, π) × (−π/2, π/2) 2
tvoří zobecněnou parametrizaci sféry S = {[x, y, z] ∈ R3 : x2 + y 2 + z 2 = 1}. 2-nulová množina {[x, y, z] ∈ S : y = 0, x ≤ 0} je nepokryta. Porovnejte se cvičením 1.3. 1.10. Poznámka. Vždy si rozmyslete, zda symbol x je bod nebo pouze první souřadnice bodu. Značení x = (x1 , x2 , x3 ) se zpravidla používá v teorii a značení [x, y, z] v příkladech. První způsob je pohodlnější a pro n-rozměrnou analýzu jediný možný. Pokud si však nezvykneme i na druhé značení, můžeme zahodit většinu sbírek úloh a rezignovat na porozumění s fyziky, inženýry a dalšími uživateli matematiky. 1.11. Definice (Plošný integrál prvého druhu). Nechť Φ : G → Rn je k-rozměrná parametrická plocha v Rn a u je funkce na Φ(G). Definujeme Z Z Z (1) u dS = u(x) dS(x) := u(Φ(t)) |JΦ(t)| dt. Φ
Φ
G
Definici (a podobným definicím v dalším) rozumíme tak, že integrál vlevo má smysl, když má smysl integrál vpravo. 1.12. Lemma. Nechť Φ : G → Rn je prostá regulární k-rozměrná parametrická plocha a E ⊂ G je měřitelná množina. Nechť Φ(E) je k-nulová množina. Potom E je Lebesgueovsky nulová. Důkaz. Stačí ukázat, že ke každému bodu s ∈ G existuje okolí V tak, že V ∩ E je Lebesgueovsky nulová. Zvolme tedy bod s ∈ G a označme z = Φ(s). Najdeme k lineárně nezávislých řádků matice Φ0 (s). BÚNO jde o řádky odpovídající indexům 1, . . . , k. Buď π projekce x 7→ (x1 , . . . , xk ). Potom najdeme okolí V bodu s tak, že π ◦ Φ je difeomorfismus V na okolí bodu π(z). Jelikož množina Φ(E ∩ V ) je k-nulová, tím spíš je π(Φ(E ∩ V )) také k-nulová, tedy i Lebesgueovsky. Z věty o substituci dostaneme Z |J(π ◦ Φ)(t)| dt = λk (π(Φ(E ∩ V )) = 0. E∩V
Jelikož J(π ◦ Φ) > 0 na V , znamená to, že E ∩ V je Lebesgueovsky nulová.
n
1.13. Věta (nezávislost plošného integrálu na parametrizaci). Nechť M ⊂ R a u : M → R je funkce. Nechť Φ : G → Rn , Ψ : H → Rn jsou zobecněné parametrizace M . Potom Z Z u dS = u dS, Φ
Ψ
pokud aspoň jedna strana má smysl. Důkaz. Položme M1 = Φ(G), M2 = Ψ(H), M 0 = M1 ∩ M2 , G0 = Φ−1 (M 0 ), H 0 = Ψ−1 (M 0 ). Potom M1 , M2 , M 0 jsou relativně otevřené podmnožiny M a G0 je otevřená podmnožina G. Také víme, že M \ M 0 je k-nulová množina. Ukážeme, že zobrazení Ψ−1 ◦ Φ je difeomorfismus G0 na H 0 . Zvolme s0 ∈ G0 a t0 ∈ H 0 tak, že Φ(s0 ) = Ψ(t0 ). Jelikož Ψ0 (t0 ) má hodnost k, existuje k-tice indexů tak, že příslušné řádky matice Ψ0 (t0 ) jsou lineárně nezávislé. BÚNO jde o indexy 1, . . . , k. Buď π projekce x 7→ (x1 , . . . , xk ). Na G0 × H 0 uvažujme implicitní funkci Ξ(s, t) = π(Φ(s)) − π(Ψ(t)). Snadno ověříme splnění předpokladů, tedy na součinu okolí G00 bodu s a H 00 bodu t0 je nulová úrovňová množina Ξ grafem spojitě diferencovatelné funkce t = ξ( s). Jelikož Φ je spojité a Ψ je otevřené vzhledem k M , je Ψ−1 ◦ Φ spojitá funkce, jejíž graf na okolí (s0 , t0 ) též vyhovuje uvažované implicitní rovnici. Z jednoznačnosti plyne, že ξ = Ψ−1 ◦ Φ na okolí s0 . Tím dostaneme, že složené zobrazení Ψ−1 ◦ Φ je spojitě diferencovatelné na G0 a na celém G0 ho též nyní budeme značit ξ. Máme Φ = Ψ ◦ ξ a Φ0 (s) = Ψ0 (ξ(s))ξ 0 (s). Z věty o násobení determinantů snadno odvodíme |JΦ(s)| = |JΨ(ξ(s))| |Jξ(s)|. Tedy ξ je regulární. Z věty o substituci dostaneme Z Z Z u(Φ(s)) |JΦ(s)| ds = u(Ψ(ξ(s)) |JΨ(ξ(s))| |Jξ(s)| ds = u(t)|JΨ(t)| dt. G0
G0 0
0
H0 0
0
Víme, že Φ(G \ G ) ⊂ M \ M a Ψ(H \ H ) ⊂ M \ M jsou k-nulové množiny. Použijeme lemma 1.12 a dostaneme Z Z Z Z Z u dS = u(Φ(s)) |JΦ(s)| ds = u(Φ(s)) |JΦ(s)| ds = u(t)|JΨ(t)| dt = u(t)|JΨ(t)| dt Φ G0 H0 H ZG = u dS. Ψ
3
1.14. Definice (Integrál prvého druhu přes množinu). Řekneme, že M ⊂ Rn je přípustná k-rozměrná množina, pokud M připouští k-rozměrnou zobecněnou parametrizaci. Přes takové množiny budeme definovat k-rozměrnou integraci. Každá k-rozměrná plocha je přípustná pokud k ≥ 1, ale jsou i jiné důležité příklady. Například hranice n-rozměrné krychle, N ≥ 2, je přípustná (n−1)-rozměrná množina, ale není to (n−1)-rozměrná plocha. Nechť M ⊂ Rn je přípustná k-rozměrná množina a u : M → R je funkce. Potom definujeme (krozměrný) plošný integrál u přes M předpisem Z Z u dS = u dS, M
Φ
kde Φ je zobecněná k-rozměrná parametrizace M . Pokud integrál vpravo nemá smysl, zůstává integrál vlevo nedefinovaný. Z předchozí věty plyne, že taková definice je korektní. Pokud k = 0, z definice vychází jako přípustné množiny jen jednobodové množiny, avšak v tomto případě není zapotřebí definovat integrál pomocí parametrizace; odpovídající pojem je prostě suma přes množinu. 1.1. Cvičení. Nechť G ⊂ Rn je k-rozměrná plocha a Φ : G → G je prosté a regulární. Potom Φ je otevřené vzhledem ke G, tedy lokální parametrizace. Návod. Zvolme z = Φ(s) ∈ Φ(G) a okolí G0 bodu s. Na okolí W bodu z najdeme spojitě diferencovatelné zobrazení g do Rn−k tak, že g = 0 na G a g 0 (z) má hodnost n−k. Potom existuje (n − k)-tice indexů tak, že příslušné sloupce matice g 0 (z) jsou lineárně nezávislé. BÚNO jede o indexy k+1, . . . , n. Podle věty 1.8 a jejího důkazu pak existuje okolí H bodu (z1 , . . . , zk ), okolí V bodu (zk+1 , . . . , zn ) a spojitě diferencovatelná funkce ξ : H → V tak, že H × V ⊂ W a G je grafem ξ na G × V . Potom (Φk+1 , . . . , Φn ) = ξ ◦ (Φ1 , . . . , Φk ) na okolí s. Tedy ∇Φj jsou lineárními kombinacemi ∇Φ1 , . . . , ∇Φk pro j > k. Kdyby ∇Φ1 (s), . . . , ∇Φk (s) byly lineárně závislé, nedosáhli bychom hodnosti k pro Φ0 (s) a to by byl spor s regularitou Φ. Zobrazení (Φ1 , . . . , Φk ) splňuje předpoklady věty o lokálním difeomorfismu v bodě s a na okolí H 0 ⊂ H bodu (z1 , . . . , zk ) tedy existuje inverzní spojité zobrazení Ξ k restrikci (Φ1 , . . . , Φk ) na vhodné okolí G00 ⊂ G0 bodu s. Tedy (H 0 × V ) ∩ G ⊂ Φ(G0 ). 1.2. Cvičení. Nechť G = (−1, π) a ( (1, t), Φ(t) = (cos t, sin t),
(−1 < t ≤ 0), 0 < t < π.
Potom Φ je prosté a regulární na G, ale není otevřené vzhledem k Φ(G). Inverzní zobrazení není spojité v bodě (1, 0); na každém okolí tohoto bodu nabývá hodnot blízkých π, ale funkční hodnota inverzního zobrazení je 0. 1.3. Cvičení. Nechť S = {x ∈ Rn : |x| = 1}, n ≥ 2. Dokažte, že S je (n−1)-rozměrná plocha v Rn , která nepřipouští globální parametrizaci. Návod. Nechť G ⊂ Rn−1 je otevřená množina a Φ : G → S je globální parametrizace S. Potom podle cvičení 1.1 je Φ homeomorfismus a Φ−1 také. Tudíž G je spojitý obraz kompaktní množiny S, takže G je neprázdná kompaktní otevřená podmnožina Rn−1 a to je spor. 1.4. Cvičení. Spočetné sjednocení k-nulových množin je k-nulové. 1.5. Cvičení. Nechť 0 ≤ m < k ≤ n jsou celá čísla a N je m-rozměrná plocha v Rn . Ukažte, že N je k-nulová. 1.6. Cvičení. Nechť A ⊂ G ⊂ Rk , A je Lebesgueovsky k-nulová, G je otevřená a a Φ : G → Rn je k-rozměrná parametrická plocha. Pak Φ(A) je k-nulová v Rn . 1.7. Cvičení. Nechť N ⊂ U ⊂ Rn , N je k-nulová množina a U je otevřená. Nechť g : U → Rk je spojitě diferencovatelné zobrazení. Potom g(N ) je Lebesgueovsky k-nulová. 1.8. Cvičení. Nechť N ⊂ Rn . Potom N je N -nulová, právě když je Lebesgueovsky nulová.
2. Křivkový integrál Pojmy “plocha” a “křivka” se v matematice používají v mnoha různých významech. V tomto kurzu je zavedeme tak, jak se hodí pro účely integrace. 2.1. Definice (Křivka). Křivka (přesněji C 1 -křivka) v Rn je spojitě diferencovatelné zobrazení ψ intervalu [a, b] ⊂ R do Rn . Derivace ψ v krajních bodech a, b chápeme jako jednostranné. Interval [a, b] nazveme referenčním intervalem křivky ψ. 4
2.2. Definice (Zobecněná křivka). 1-rozměrnou parametrickou plochu budeme nazývat zobecněnou křivkou. Každé křivce ψ : [a, b] → Rn odpovídá zobecněná křivka ψ ◦ = ψb(a, b), tedy “ořízneme” hodnoty v krajních bodech referenčního intervalu. Ztráta informace je jen zdánlivá, “chybějící” krajní body křivky můžeme znovu zrekonstruovat jako limity v krajních bodech referenčního intervalu. Zobecněná křivka zobecňuje pojem křivky ve dvou směrech: • v krajních bodech nepožadujeme existenci jednostranných limit, • referenční obor nemusí být souvislý. Zdůrazněme, že (zobecněné) křivky pokládáme za zvláštní případ parametrických ploch. Tedy zformulujeme-li tvrzení (definici, poznámku,. . .) pro parametrické plochy, máme tím na mysli i aplikaci na křivky. Pojem křivka používáme jen mluvíme-li o specifikách jednorozměrného případu. 2.3. Definice (Křivkový integrál prvého druhu). Nechť Φ : G → Rn je zobecněná křivka. Potom Jacobiho matice Φ0 (t) má n řádků a jen jeden sloupec, je to tedy vlastně “jen svislý vektor”. Potom JΦ = Φ0 a pro křivkový integrál prvého druhu funkce u platí vzorec Z Z Z u ds = u(x) ds(x) = u(Φ(t)) |Φ0 (t)| dt. Φ
Φ
G
Všimněte si, že pro křivkovou integraci se zpravidla píše “diferenciál” ds místo dS. Integrál Z Z ds = |Φ0 (t)| dt Φ
G
má geometrický význam délky (zobecněné) křivky. (Bez ohledu na to, zda křivka je prostá či ne, může se i “protínat” či dokonce “probíhat některé úseky vícekrát”. V takovém případě se ovšem i délka příslušného úseku objeví ve výsledku vícekrát a délka křivky se může lišit od “délky” množiny Φ(G).) 2.4. Definice (Křivkový integrál druhého druhu). Nechť Φ = (Φ1 , . . . , Φn ) : G → Rn je zobecněná křivka a f = (f1 , . . . , fn ) : Φ(G) → Rn je vektorové pole. Definujeme Z Z (2) f · ds = f (Φ(t)) · Φ0 (t) dt. Φ
G
Také pro index i a skalární funkci u : Φ(G) → Rn píšeme Z Z (3) u dxi := u(Φ(t)) · Φ0i (t) dt. Φ
G
Zřejmě můžeme přepsat Z
Z f · ds =
Φ
f1 dx1 + · · · + fn dxn . Φ
Křivkový integrál (2) má velký význam ve fyzice, křivkovým integrálem druhého druhu se integrují veličiny, u nichž není zajímavý úhrn celkové velikosti, ale úhrn tečné složky. Například práce je křivkový integrál druhého druhu síly po dráze. 2.5. Pole. Pojmy skalární pole, vektorové pole se používají jako synonyma pro skalární, resp. vektorovou funkci. Jejich používání v některých situacích je dáno zvyklostmi. 2.6. Definice (Tečné pole). Nechť Φ : G → Rn je prostá regulární zobecněná křivka. Je-li x = Φ(t), t ∈ G, definujeme (4)
τ Φ (x) =
Φ0 (t) . |Φ(t)|
Funkce τ Φ : Φ(G) → Rn se nazývá pole jednotkových tečných vektorů (zkráceně tečné pole) ke křivce Φ. Zahrnuje v sobě informaci o tečném prostoru v každém bodě a směru probíhání křivky. 2.7. Definice (Orientovaná jednorozměrná plocha). Terminologie “plocha”דparametrická plocha” nám umožňuje rozlišovat mezi pojetím plochy jako množiny a zobrazení. Analogický model “křivka”דparametrická křivka” se nehodí, protože v současné matematice je již prakticky jasná shoda na tom, že křivka je zobrazení. Proto pro “křivku-množinu” budeme raději používat termín “jednorozměrná plocha”. Nechť Γ ⊂ Rn je jednorozměrná plocha a τ je spojité pole jednotkových vektorů na Γ. Řekneme, že (Γ, τ ) je orientovaná jednorozměrná plocha, jestliže pro každý bod x ∈ Γ existuje okolí U bodu x a globální jednorozměrná parametrizace Φ množiny U ∩ Γ tak, že τ = τΦ
na 5
U ∩ Γ.
Je-li (Γ, τ ) orientovaná jednorozměrná plocha, pak (Γ, −τ ) je také orientovaná jednorozměrná plocha. Souvislé jednorozměrné plochy se dají orientovat jen takto dvěma způsoby. Viz. cvičení 2.1-2.2. 2.8. Definice (Kladná parametrizace). Nechť (Γ, τ ) je orientovaná jednorozměrná plocha a Φ : G → Γ je lokální parametrizace Γ. Řekneme, že Φ je kladná, jestliže pro každé x ∈ Φ(G) je τ Φ (x) = τ (x). 2.9. Věta (o kladné parametrizaci). Nechť Γ ⊂ Rn je orientovaná jednorozměrná plocha. Pak Γ má kladnou zobecněnou parametrizaci. Jestliže f : Γ → Rn je vektorové pole a Φ : G → Rn , Ψ : H → Rn jsou kladné zobecněné parametrizace Γ, pak Z Z f · ds = f · ds, Φ
Ψ
pokud aspoň jedna strana má smysl. Důkaz. Ze samotné definice orientované jednorozměrné plochy je zřejmé, že ke každému bodu existuje okolí, které lze kladně parametrizovat. Nyní můžeme postupovat stejně, jako v důkazu věty 1.8 (b) a dostaneme existenci kladné zobecněné parametrizace. Jsou-li Φ : G → Γ a Ψ : H → Γ kladné zobecněné parametrizace, potom Z Z Z Z Φ0 (t) 0 |Φ (t)| dt = f (Φ(t)) · τ Γ (t)|Φ0 (t)| dt f · ds = f (Φ(t)) · Φ0 (t) dt = f (Φ(t)) · 0 |Φ (t)| G Φ G G Z (5) = f · τ Γ ds Φ
a podobně Z
Z f · ds =
f · τ Γ ds
Ψ
Ψ
Podle věty 1.13 je Z
Z f · ds =
Φ
Z
Z
f · τ Γ ds = Φ
f · τ Γ ds = Ψ
f · ds Ψ
2.10. Definice (Křivkový integrál druhého druhu přes orientovanou jednorozměrnou plochu). Nechť (Γ, τ ) je orientovaná jednorozměrná plocha v Rn a f : Γ → Rn je vektorové pole. Definujeme Z Z f · ds = f · ds, Γ
Φ
kde Φ : G → Γ je kladná zobecněná parametrizace Γ. 2.11. Věta (Vztah mezi křivkovým integrálem prvého a druhého druhu). Nechť (Γ, τ ) je orientovaná jednorozměrná plocha v Rn a f : Γ → Rn je vektorové pole. Potom Z Z f · ds = f · τ ds, Γ
Γ
pokud aspoň jeden z integrálů má smysl. Důkaz. Vzorec dostaneme z definic a výpočtu (5).
2.1. Cvičení. Dokažte, že každou jednorozměrnou plochu Γ je možné orientovat. Návod. Důkaz je docela těžký. Naznačíme zde hlavní myšlenky. Nejprve jako v důkazu věty 1.8 vytvoříme systém lokálních parametrizací tak, aby parametrizované části pokrývaly Γ a referenční obory byly intervaly. Z tohoto pokrytí vybereme spočetné a lokálně konečné, máme konečnou či nekonečnou posloupnost (Φj )j , Φj : Gj → Γ. Vytvoříme diskrétní množinu bodů A tak, že V každé průniku Φi (Gi ) ∩ Φj (Gj ) zvolíme tři body, pokud je tento průnik neprázdný. Na množině A uvažujeme graf: řekneme, že a sousedí s b, pokud existuje lokální parametrizace Φ : G → Γ a α, β ∈ G tak, že (α, β) ⊂ G, Φ(α) = a, Φ(β) = b a žádný bod z intervalu (α, β) se nezobrazí do A. Lze dokázat (dá to trochu práce), že každý bod množiny A má právě dva sousedy. Nyní přichází ke slovu kombinatorika, množina A se rozloží na cykly a řetězce, přičemž řetězce jsou z obou stran neomezené. Na každém cyklu či a řetězci se zvolí směr probíhání a s jeho pomocí se vytvoří orientace na Γ. 2.2. Cvičení. Dokažte, že každou souvislou jednorozměrnou plochu Γ je možné orientovat jen dvěma způsoby. 6
Návod. Nechť τ a = τ˜ jsou orientace Γ. Zvolíme bod z ∈ Γ a uvažujeme lokální parametrizace Φ : G → Γ, Ψ : H → Γ a body s0 ∈ G, t0 ∈ H tak, že z = Φ(s0 ) = Ψ(t0 ), přičemž τ Φ (z) = τ (z) a τ Ψ (z) = τ˜ (z). Podobně jako v důkazu věty 1.8 najdeme okolí G0 bodu s0 a regulární funkci ξ : G0 → H tak, že ξ(s0 ) = t0 a Φ = Ψ ◦ ξ na G0 . Potom Φ0 (s0 ) = Ψ0 (t0 )ξ 0 (s0 ) a tudíž ( τ (z), pokud ξ 0 (s0 ) > 0, τ˜ (z) = −τ (z), pokud ξ 0 (s0 ) < 0. Množina všech bodů, kde τ = τ˜ , resp. τ = −˜ τ je otevřená i uzavřená v Γ, tedy ze souvislosti dostaneme tvrzení.
3. Elementy teorie pole 3.1. Definice (Divergence, gradient, rotace). Nechť U ⊂ Rn je otevřená množina, u : U → R je spojitě diferencovatelná funkce a f = (f1 , . . . , fd ) : U → Rn je spojitě diferencovatelné vektorové pole (vektorové pole znamená zobrazení s hodnotami v Rn ). Nechť (e1 , . . . , ed ) je kanonická báze v Rn . Definujeme ∇u = grad u :=
div f :=
d X ∂u ei , ∂x i i=1
d X ∂fi , ∂x i i=1
(gradient u),
(divergence f )
∂f1 ∂f2 − (rotace f , d = 2), curl f := ∂x1 ∂x2 ∂f ∂f ∂f ∂f2 ∂f3 ∂f1 3 1 2 curl f := − e1 + − e2 + − e3 ∂x2 ∂x3 ∂x3 ∂x1 ∂x1 ∂x2
(rotace f , d = 3).
3.2. Věta o potenciálu. Nechť W ⊂ Rn je otevřená množina. Nechť ψ : [a, b] → Rn je křivka, ψ([a, b]) ⊂ W , A = ψ(a), B = ψ(b). Nechť u : W → R je spojitě diferencovatelná funkce. Potom Z u(B) − u(A) = ∇u · ds, ψ
pokud integrál vpravo konverguje. Důkaz. Máme Z u(B) − u(A) = u(ψ(b)) − u(ψ(a)) =
b 0
Z
(u ◦ ψ) (t) dt = a
b 0
Z
∇u(ψ(t)) · ψ (t) dt = a
∇u · ds. ψ
3.3. Definice (Hvězdovitá množina). Řekneme, že množina U ⊂ Rn je hvězdovitá, jestliže existuje A ∈ U tak, že pro každý bod x ∈ U je celá úsečka {A + t(x − A) : t ∈ [0, 1]} podmnožinou U . Bod A se nazývá střed hvězdovitosti množiny W . Každá konvexní množina C je hvězdovitá, v tom případě je každý bod x ∈ C středem hvězdovitosti množiny C. 3.4. Věta (Hlavní věta teorie pole). Nechť W ⊂ Rn je otevřená množina a f = (f1 , . . . , fn ) : W → Rn je spojité vektorové pole. Uvažujme následující podmínky: (i) (Existence potenciálu.) Existuje spojitě diferencovatelná funkce u : W → R tak, že f = ∇u. (ii) (Nezávislost integrálu na dráze.) Pro každé dva body A, B ∈ W existuje číslo c = c(A, B) tak, že a každou křivku ψ : [a, b] → W s počátečním bodem A = ψ(a) a koncovým bodem B = ψ(b) je Z f · ds = c. ψ
(iii) (Nulová rotace.) Pro každou dvojici i, j indexů z {1, . . . , d} je ∂fj ∂fi = . ∂xj ∂xi Potom platí následující vztahy: (a) (i) ⇐⇒ (ii), (b) Je-li f spojitě diferencovatelná, pak (i) =⇒ (iii). (c) Je-li f spojitě diferencovatelná a W hvězdovitá, pak pak (iii) =⇒ (ii). 7
Důkaz. Implikace (i) =⇒ (ii) plyne z věty 3.2 a implikace (ii) =⇒ (iii) je důsledkem záměnnosti smíšených derivací. Nechť je splněno (ii). Pokud existuje potenciál u v každé komponentě souvislosti množiny W , pak existuje na W . Můžeme tedy předpokládat, že W je souvislá. Zvolme bod A ∈ W a položme Z u(x) = f · ds, ψ
kde ψ : [a, b] → W je křivka, která “vede” z A do x, neboli ψ(a) = A, ψ(b) = x. Ze souvislosti plyne, že taková křivka existuje a podle (ii) integrál na volbě ψ nezávisí. Zvolme x ∈ W . Najděme r > 0 tak, že B(x, r) ⊂ W a zvolme h ∈ B(0, r). Položme γ(t) = x+th, t ∈ [0, 1]. Najděme ještě křivku ψ : [−1, 0] → W tak, že ψ(−1) = A, ψ(0) = x a položme ( ψ(t), −1 ≤ t ≤ 0, ˜ ψ(t) = γ(t), 0 < t ≤ 1. Potom ψ˜ vede z A do x + h, tedy Z u(x + h) − u(x) =
Z
Z
f · ds − ˜ ψ
f · ds = ψ
Z f · ds =
γ
1
f (x + th) · h dt. 0
Odtud pro h → 0 máme u(x + h) − u(x) − f (x) · h = h
Z
1
0
(f (x + th) − f (x)) · h → 0, |h|
tedy f (x) je derivace u v x. Konečně, předpokládejme, že W je hvězdovitá se středem hvězdovitosti A a f = (f1 , . . . , fn ) je spojitě diferencovatelná. BÚNO máme A = 0. Položme Z u(x) = f · ds, ψx
kde tentokrát ψ x (t) = tx,
t ∈ [0, 1].
Potom Z
1
f (tx) · x dt
u(x) = 0
a derivováním za integračním znamením zjistíme Z 1 n X ∂u(x) ∂fj (tx) (6) = fi (tx) + t · xj dt. ∂xi ∂xi 0 j=1 Máme též Z (7)
fi (x) = 0
1
d tfi (tx) dt = dt
Z 1 n X ∂fi (tx) fi (tx) + t · xj dt. ∂xj 0 j=1
Porovnáním (6), (7) a (iii) dostaneme, že ∂u(x) = fi (x), ∂xi
i = 1, . . . , n.
3.5. Poznámka. Nulovost rotace je rovnost curl f = 0 v dimenzi 2 a rovnost curl f = 0 v dimenzi 3. 3.6. Příklady. Množina P = {x ∈ R2 : x 6= 0} není hvězdovitá. Množina D = {x ∈ R2 : x2 6= 0 nebo x1 > 0} je hvězdovitá ale není konvexní. Tyto příklady mají značný význam pro teeorii funkcí komplexní proměnné. Množina D je přirozený definiční obor hlavní větve komplexního logaritmu. Na množině P můžeme najít příklad vektorového pole, které má nulovou rotaci, ale není potenciální, tj. není globálně gradientem žádné funkce. Nejznámějším příkladem je pole −x x 1 2 f (x) = , . |x|2 |x|2 8
Lokálně se potenciály k f nalézt dají, například arg(x1 + ix2 ) = arctg(x2 /x1 ) na {x ∈ R2 : x1 > 0}. Pro komplexní analýzu to znamená, že existuje uzavřená křivka ψ v P , konkrétně ψ(t) = eit , t ∈ [−π, π], a holomorfní funkce h na P , konkrétně h(z) = 1/z, tak, že Z h(z) dz = 2πi 6= 0, ψ
což je zdánlivý spor s Cauchyovou větou. 4. Plošný integrál kodimenze 1 4.1. Definice (Vektorový součin). Vektorový součin vektorů u1 , . . . , un−1 ∈ Rn je vektor u1 × · · · × un−1 :=
n X
det(ei , u1 , . . . , un−1 ) ei .
i=1
Vektorový součin je kolmý na své činitele. Při liché permutaci činitelů změní vektorový součin znaménko, při sudé zůstane zachován. V dimenzi tři je vektorovým součinem vektorů u = (u1 , u2 , u3 ) a v = (v1 , v2 , v3 ) vektor ! u2 , v2 u3 , v3 u1 , v1 u × v = det , det , det . u3 , v3 u1 , v1 u2 , v2 V dimenzi 2 má vektorový součin jen jednoho činitele. Vektorový součin vektoru [u1 , u2 ] je vektor [u2 , −u1 ]. 4.2. Věta. Nechť v, u1 , . . . , un−1 ∈ Rn a w = u1 × · · · × un−1 . Potom v · w = det(v, u1 , . . . , un−1 ). Důkaz. Pro bázové vektory v = ej P je vzorec snadným důsledkem definice vektorového součinu. V obecném případě použijeme rozvoj v = i vi ei . 4.3. Věta (Cauchyova-Binetova formule). Nechť u1 , . . . , un−1 ∈ Rn a w = u1 × · · · × un−1 . Nechť A je matice, která má jako sloupce vektory u1 , . . . un−1 . Potom (8)
det(AT A) = |w|2 .
Důkaz. Nechť u1 , . . . un−1 jsou lineárně nezávislé. Buď M matice, která má sloupce w, u1 , . . . un−1 . Potom M T M je blokově diagonální, totiž 2 |w| , 0T T M M= , 0, AT A kde 0 je nulová matice typu (n−1) × 1. Tedy (9)
|w|2 det(AT A) = det M T M.
Z věty 4.2 plyne det M = w · w = |w|2 . Podle vět o součinu determinantů a determinantu transponované matice je (10)
det M T M = |w|4 .
Porovnáním (9) a (10) dostaneme (8). Jestliže u1 , . . . un−1 jsou lineárně závislé, z elementární lineární algebry vyvodíme, že na obou stranách (8) je nula. 4.4. Definice (Vektorový jakobián a plošný integrál kodimenze jedna). Vektorový jakobián (n−1)rozměrné parametrické plochy Φ : G → Rn v bodě t ∈ G ⊂ Rn−1 definujeme předpisem JΦ(t) =
∂Φ ∂Φ (t) × · · · × (t). ∂t1 ∂tn−1 9
Podle věty 4.3 je |JΦ(t)| = |JΦ(t)|, tedy jakobián pro kalkulus plošného integrálu prvého druhu lze v kodimenzi 1 počítat alternativním způsobem Z Z Z f dS = f (x) dS(x) := f (Φ(t)) |JΦ(t)| dt. Φ
Φ
G
Integrál Z
Z |JΦ(t)| dt
dS = Φ
G
má geometrický význam obsahu (area) plochy. Podobně jako u křivky, o obsahu parametrické plochy Φ můžeme mluvit i tehdy, když Φ není prostá, pak se ale může lišit od “obsahu” množiny Φ(G). 4.5. Definice (Plošný integrál druhého druhu). Nechť Φ : G → Rn je (n−1)-rozměrná parametrická plocha v Rn . Nechť f = (f1 , . . . , fn ) : Φ(G) → Rn je vektorové pole. Potom definujeme Z Z (11) f · dS := f (Φ(t)) · JΦ(t) dt. Φ
G
Integrály typu (11) se hojně vyskytují ve fyzice, mají např. význam toku plochou. 4.6. Zápis pomocí diferenciálů. Jestliže n = 2, je Φ zobecněná křivka a integrál uvedený výše lze přepsat ve tvaru Z Z f · dS = f1 dx2 − f2 dx1 . Φ
Φ
Zde R velikost symbolu S v diferenciálu hraje významnou roli. Musíme striktně rozlišovat mezi integrálem f · dS a integrálem Φ Z Z f · ds = f1 dx1 + f2 dx2 . Φ
Φ
V dimenzi tři, pro dvourozměrnou parametrickou plochu Φ = (Φ1 , Φ2 , Φ3 ) : G → R3 , skalární pole u a dvojici indexů (i, j) ∈ {1, 2, 3}2 definujeme Z Z u dxi dxj = u(Φ(t)) det ∇Φi , ∇Φj dt Φ
G
Všimněme si, že takový integrál závisí znaménkem na pořadí diferenciálů a pro i = j je nulový! Pak lze psát Z Z f · dS = f1 dx2 dx3 − f2 dx1 dx3 + f3 dx1 dx2 . Φ
Φ
Podobně lze zapisovat různé integrály ve vyšších dimenzích, a nejen pro plochy dimenze či kodimenze jedna. Je-li Φ : G → Rn k-rozměrná parametrická plocha, u : Φ(G) → R skalární funkce a α = (α1 , . . . , αk ) je uspořádaná k-tice indexů, pak definujeme Z Z ∂Φ k αi dt. u(x) dxα1 . . . dxαk = u(Φ(t)) det ∂tj i,j=1 Φ G 4.7. Definice (Normálové pole). Nechť Φ : G → Rn je prostá regulární parametrická plocha kodimenze 1. Je-li x = Φ(t), t ∈ G, definujeme (12)
ν Φ (x) =
JΦ(t) . |JΦ(t)|
Funkce ν Φ : Φ(G) → Rn se nazývá pole jednotkových normálových vektorů (zkráceně normálové pole) k ploše Φ. Zahrnuje v sobě informaci o tečném prostoru v každém bodě a orientaci (rub-líc). 4.8. Definice (Orientovaná plocha kodimenze 1). Nechť Γ ⊂ Rn a ν je spojité pole jednotkových vektorů na Γ. Řekneme, že (Γ, ν) je orientovaná plocha kodimenze 1, jestliže pro každý bod x ∈ Γ existuje okolí U bodu x a globální jednorozměrná parametrizace Φ množiny U ∩ Γ tak, že ν = νΦ
na
U ∩ Γ.
Je-li (Γ, ν) orientovaná plocha kodimenze 1, pak (Γ, −ν) je také orientovaná plocha kodimenze 1. Souvislé plochy kodimenze 1 se dají orientovat jen takto dvěma způsoby (pokud vůbec, tzv. Möbiova páska se orientovat nedá, viz. cvičení 4.2). 10
4.9. Poznámka. Nechť Γ je plocha kodimenze 1 a z ∈ Γ. Potom existuje okolí U bodu z a spojitě diferencovatelná funkce g : U → R tak, že g 0 6= 0 na U a U ∩ Γ = {x ∈ U : g(x) = 0}. Nechť existuje globální parametrizace Φ : G → Γ plochy U ∩ Γ. Potom zderivováním identity g ◦ Φ = 0 dostaneme, že vektor ∇g(Φ(t)) je kolmý na všechny parciální derivace Φ v bodě t, t ∈ G. Tedy ∇g(Φ(t)) je reálný násobek JΦ(t). Položíme-li ν(x) =
∇g(x) , |∇g(x)|
x ∈ Γ ∩ U,
pak (Γ ∩ U, ν) je orientovaná plocha kodimenze 1. 4.10. Definice (Kladná parametrizace). Nechť (Γ, ν) je orientovaná plocha kodimenze 1 a Φ : G → Γ je lokální parametrizace Γ. Řekneme, že Φ je kladná, jestliže pro každé x ∈ Φ(G) je ν Φ (x) = ν(x). Jinými slovy, Φ je kladná, jestliže pro každé t ∈ G je JΦ(t) kladný násobek ν(Φ(t)). 4.11. Věta (o zobecněné parametrizaci). Nechť Γ ⊂ Rn je orientovaná plocha kodimenze 1. Pak Γ má kladnou zobecněnou parametrizaci. Jestliže f : Γ → Rn je vektorové pole a Φ : G → Rn , Ψ : H → Rn jsou kladné zobecněné parametrizace Γ, pak Z Z f · dS = f · dS, Φ
Ψ
pokud aspoň jedna strana má smysl. Důkaz. Dokazuje se analogicky jako věta 2.9. Namísto (5) počítame Z Z Z Z JΦ(t) |JΦ(t)| dt = f (Φ(t)) · ν Γ (t)|Φ0 (t)| dt f · dS = f (Φ(t)) · JΦ(t) dt = f (Φ(t)) · |JΦ(t)| G Φ G G Z (13) = f · ν Γ dS. Φ
4.12. Definice (Plošný integrál druhého druhu přes orientovanou plochu kodimenze 1). Nechť (Γ, ν) je orientovaná plocha kodimenze 1 v Rn a f : Γ → Rn je vektorové pole. Definujeme Z Z f · dS = f · dS, Γ
Φ
kde Φ : G → Γ je kladná zobecněná parametrizace Γ. 4.13. Věta (Vztah mezi plošným integrálem prvého a druhého druhu). Nechť (Γ, ν) je orientovaná plocha kodimenze 1 v Rn a f : Γ → Rn je vektorové pole. Potom Z Z f · dS = f · ν dS, Γ
Γ
pokud aspoň jeden z integrálů má smysl. Důkaz. Vzorec dostaneme z definic a výpočtu (13).
4.14. Příklad. Uvažujme sférické souřadnice jako v příkladu 1.9. Zobecněná parametrizace je kladná pro (S, ν) když ν(x, y, z) = [x, y, z], [x, y, z] ∈ S. Vektorový jakobián v bodě (α, γ) je ∂x ∂x − cos γ sin α − sin γ cos α cos γ cos α ∂γ ∂α ∂y ∂y × ∂γ = cos γ cos α × − sin γ sin α = cos γ cos γ sin α . ∂α ∂z ∂z 0 cos γ sin γ ∂α ∂γ Přesvědčili jsme se, že vektorový jakobián je kladný násobek normály, test orientace prošel. Kdybychom zaměnili pořadí proměnných na (γ, α) parametrizace by byla záporná. 4.1. Cvičení. Dokažte, že každou souvislou jednorozměrnou plochu Γ je možné orientovat nejvýše dvěma způsoby. Návod. Postupujte analogicky jako ve cvičení 2.2.
11
4.2. Cvičení. Uvažujme zobrazení Φ = (x, y, z) : G → R3 , x = 2 + t cos(α/2) cos α, y = 2 + t cos(α/2) sin α,
[t, α] ∈ G,
z = t sin(α/2), kde G = ((−1, 1) × R). Möbiova páska je množina M = Φ(G). Ukažte, M je dvourozměrná plocha v R3 a že se nedá orientovat. Totiž, kdyby ν byla orientace M , ze spojitosti bychom dostali ν(Φ(t)) =
JΦ(t) , |JΦ(t)|
t ∈ G,
nebo ν(Φ(t)) = −
JΦ(t) , |JΦ(t)|
t ∈ G.
Ale JΦ((0, 0)) = −JΦ(0, 2π), to je spor.
5. Věta o divergenci 5.1. Definice (Vnější normála). Buď n > 1. Nechť Ω ⊂ Rn je omezená otevřená množina. Řekneme, že v bodě z ∈ ∂Ω existuje vnější normála k Ω, jestliže existuje okolí U bodu z a rozhraničující funkce h : U → R tak, že ∇h(z) 6= 0, x ∈ U =⇒ [x ∈ Ω ⇐⇒ h(x) < 0] . V tom případě definujeme vnější normálu předpisem ν Ω (z) =
∇h(z) . |∇h(z)|
Vnější normála nezávisí na volbě rozhraničující funkce. Je to jednotkový vektor, směřuje vždy “ven z Ω” a je “kolmá na hranici Ω”. Řekneme, že z ∈ ∂Ω je regulární bod hranice Ω, pokud na okolí U bodu z najdeme spojitě diferencovatelnou rozhraničující funkci. Nechť Γ je množina všech regulárních bodů hranice Ω. Potom (Γ, ν Ω ) je orientovaná (n−1)-rozměrná plocha, viz. poznámka 4.9. 5.2. Poznámka. Ne každý bod existence vnější normály je regulárním bodem hranice, v definici vnější normály se nepředpokládá hladkost a dokonce ani spojitost rozhraničující funkce. 5.3. Věta (O hraniční parametrizaci). Nechť Ω ⊂ Rn je otevřená množina a z ∈ ∂Ω. Potom z je regulární bod hranice Ω, právě když existuje otevřená množina G ⊂ Rn a prosté regulární zobrazení Ψ : G → Rn tak, že platí (14)
0 ∈ G,
(15)
pro všechna y ∈ G je Ψ(y) ∈ Ω ⇐⇒ y1 < 0,
(16)
JΨ > 0 na G,
(17)
y ∈ G, y1 = 0 =⇒
Ψ(0) = z,
∂Ψ ∂Ψ (y) × · · · × (y) je kladný násobek ν Ω (Ψ(y)). ∂y2 ∂yn
Důkaz. Existuje-li hraniční parametrizace, pak najdeme rozhraničující funkci ve tvaru h(Ψ(y)) = y1 , ověření ponecháme čtenáři (ostatně věta se zpravidla aplikuje v obráceném směru). Nechť existuje rozhraničující funkce h : U → R. Potom existuje index i tak, že ∂h(z) 6= 0. ∂xi BÚNO předpokládejme, že i = 1 a ∂h(z) > 0. ∂x2 Z věty o implicitních funkcích dostaneme existenci okolí V bodu (z2 , . . . , zn ), okolí W bodu z1 a spojitě diferencovatelné funkce ξ : V → W tak, že h i x ∈ W × V =⇒ h(x) = 0 ⇐⇒ x1 = ξ(x2 , . . . , xn ) . 12
Případným zmenšením okolí V a W docílíme též h i x ∈ W × V =⇒ h(x) < 0 ⇐⇒ x1 < ξ(x2 , . . . , xn ) . Položme U 0 = W × V a na U 0 definujme funkci g předpisem g(x) = (x1 − ξ(x2 , . . . , xn ), x2 − z2 , . . . , xn − zn ). 0
Nyní položme G = g(U ) a definujme Ψ : G → Rn jako inverzní zobrazení g −1 . Vlastnosti (14), (15) jsou zřejmě splněny. Nechť y ∈ G a x = Ψ(y). Potom (JΨ(y))−1 = Jg(x) = 1, tedy máme (16). Konečně, pro y ∈ G a x = Ψ(y) máme Ψ(y) = (y1 + ξ(z2 + y2 , . . . , zn + yn ), z2 + y2 , . . . , zn + yn ), takže 1, Ψ0 (y) = 0, 0,
∂ξ(x2 ,...,xn ) , ∂x2
..., ..., ...,
1, 0,
∂ξ(x2 ,...,xn ) ∂xn
0 1
a ∂Ψ ∂Ψ ∂ξ(x2 , . . . , xn ) ∂ξ(x2 , . . . , xn ) (y) × · · · × (y) = 1, − ,...,− . ∂y2 ∂yn ∂x2 ∂xn Z věty o implicitních funkcích máme ∂h(x)
∂ξ(x2 , . . . , xn ) ∂xi , = − ∂h(x) ∂xi
i = 2, . . . , n,
∂x1
takže ∂h(x) ∂Ψ ∂Ψ (y) × · · · × (y) = ∇h(x). ∂x1 ∂y2 ∂yn Tím jsme ověřili (17).
5.4. Poznámka. Vlastnost popsaná ve větě 5.3 se často používá k definici regulárního bodu hranice. Naši definici jsme zvolili z toho důvodu, že bývá blíže zadání úlohy a snáze se ověřuje. Definice pomocí existence hraniční parametrizace bývá pohodlnější pro autory učebnic, protože se z ní snáze dokazuje věta o divergenci. Podmínka regularity se dá oslabit na lipschitzovskou regularitu, při níž se místo spojité diferencovatelnosti požaduje jen, aby hraniční parametrizace byla bilipschitzovské zobrazení, (14), (15) zůstávají, (16) má platit jen skoro všude a (17) (n−1)-skoro všude. Pak v regulárním bodě nemusí existovat vnější normála, ale na jeho relativním okolí vzhledem k ∂Ω existuje skoro všude ve smyslu (n−1)-nulovosti. Věta o divergenci s lipschitzovskou regularitou se využívá v kurzu parciálních diferenciálních rovnic. Její důkaz není myšlenkově jiný ani těžší, jsme-li zvyklí na práci s lipschitzovskými zobrazeními. 5.5. Věta o divergenci. Nechť Ω ⊂ Rn je omezená otevřená množina. Nechť všechny body ∂Ω až na (n−1)-nulovou množinu jsou regulární. Potom ∂Ω je přípustná (n−1)-rozměrná množina. Nechť vektorové pole f : Ω → Rn je spojitě diferencovatelné na Ω a spojité na Ω. Pak Z Z (18) f · ν Ω dS = div f (x) dx, ∂Ω
Ω
pokud integrály na obou stranách konvergují. Důkaz. Důkaz je složitější, budeme mu věnovat samostatnou kapitolu 7
5.6. Poznámka. Věta o divergenci se také nazývá Gaussova, Gaussova-Greenova nebo Ostrogradského. 5.7. Test orientace pro Greenovu větu. Následující varianta je důsledek věty o divergenci. Jedná se o to, že v dimenzi 2 je množina Ω ohraničena zobecněnou křivkou, takže integrál přes regulární část hranice můžeme vnímat i jako křivkový integrál. V tom případě je přirozenější integrovat s tečným polem než s normálovým, ale tomu se musí uzpůsobit diferenciální operátor na druhé straně rovnosti. Z pole ν Ω vyrobíme pole τ Ω otočením o pravý úhel doleva, tedy (ν1 , ν2 ) 7→ (τ1 , τ2 ), τ1 = −ν2 ,
τ2 = ν1 . 13
Test orientace ohraničující zobecněné křivky ψ v tomto případě v bodě z = ψ(t) je det(∇h(z), ψ 0 (t)) > 0. kde h je rozhraničující funkce v z. 5.8. Greenova věta. Nechť Ω ⊂ R2 je omezená otevřená množina. Nechť každý bod hranice Ω až na 1nulovou množinu je regulární. Potom ∂Ω přípustná 1-rozměrná množina. Nechť vektorové pole f : Ω → R2 je spojitě diferencovatelné na Ω a spojité na Ω. Pak Z Z f · τ Ω ds = curl f (x) dx, ∂Ω
Ω
pokud integrály na obou stranách konvergují. 5.9. Příklad (Koule a sféra). Buď Ω = {x ∈ R3 : |x| < 1}, Γ = ∂Ω = {x ∈ R3 : |x| < 1}. Rozhraničující funkce v bodě x ∈ Γ je h(x) = |x|2 − 1 = x21 + x22 + x23 − 1, tedy ∇h(x) = 2x a ν(x) = x. Každý bod Γ je regulární bod ∂Ω. 5.10. Příklad (Čtverec). Buď Ω = (0, 1)2 čtverec v R2 . Hranice Ω je sjednocením čtyř úseček (“stran”), jejichž relativní vnitřky vzhledem k ∂Ω tvoří regulární část hranice. Buď např. x = [1, x2 ] ∈ {1} × (0, 1). Rozhraničující funkce na (0, 2)×(0, 1) je h(x) = x1 −1. Normála v bodě x = [1, x2 ] je ν Ω (x) = e1 = [1, 0]. Ve vrcholech čtverce Ω neexistuje vnější normála, ale ty tvoří 1-nulovou množinu. Jako zobecněnou parametrizaci ∂Ω můžeme použít zobecněnou křivku Φ, která je definována na (0, 1)∪(1, 2)∪(2, 3)∪(3, 4) předpisem [t, 0], t ∈ (0, 1), [0, t − 1], t ∈ (1, 2), Φ(t) = [3 − t, 0], t ∈ (2, 3), [0, 4 − t], t ∈ (3, 4). 5.11. Příklad (Krychle). Hranici krychle Ω = (0, 1)3 tvoří sjednocení šesti čtverců (“stěn”), jejichž relativní vnitřky vzhledem k ∂Ω tvoří regulární část hranice. Do regulární části nepatří vrcholy a hrany, ale ty tvoří 2-nulovou množinu. Nechť například x ∈ (0, 1) × {0} × (0, 1) Rozhraničující funkce je h(x) = −x2 , x ∈ (0, 1) × (−1, 1) × (0, 1), tedy ν Ω (x) = −e2 = [0, −1, 0] pro takové x. 5.1. Cvičení. Alternativně se dá zavést ν(z) jako takový jednotkový vektor, že lim sup (19)
x→z,x∈Ω
lim inf
x→z,x∈Ω /
ν(z) · (x − z) ≤ 0, |x − z| ν(z) · (x − z) ≥ 0, |x − z|
Ukažte, že tato definice je ekvivalentní. Návod. Nechť platí (19). Najdeme spojitou funkci η na okolí bodu z tak, že η(z) = 0 a pro x 6= z je ν(z) · (x − z) < η(x), |x − z| ν(z) · (x − z) − < η(x), |x − z|
x ∈ Ω, x∈ / Ω,
pak ( ν(z) · (x − z) − η(x)|x − z|, x ∈ Ω, h(x) = ν(z) · (x − z) + η(x)|x − z|, x ∈ /Ω je rozhraničující funkce. Všimněme si, že jsme nikde nepředpokládali spojitost h na okolí z.
6. Stokesova věta 6.1. Ohraničení plochy. Nechť (G, ν) ⊂ R3 je 2-rozměrná orientovaná plocha, Ω je její omezená podplocha (tj. relativně otevřená podmnožina) a Γ = Ω \ Ω. Předpokládejme, že Γ ⊂ G. Potom Γ je relativní hranice Ω v G. Řekneme, že z je regulární bod Γ, jestliže existuje okolí U bodu z a spojitě diferencovatelná rozhraničující funkce h : U → R tak, že ν(z) × ∇h(z) 6= 0 a x ∈ U ∩ G =⇒ [x ∈ Ω ⇐⇒ h(x) < 0] . V takovém bodě definujeme τ Ω (z) =
ν(z) × ∇h(z) |ν(z) × ∇h(z)| 14
Tento vektor nezávisí na volbě rozhraničující funkce. 6.2. Poznámka. Mohli bychom uvažovat i případ, kdy Ω “přesáhne” ven z G. Taková obecnost by se hodila např. kdyby G byla kuželová plocha (ovšem bez vrcholu, s vrcholem nesplňuje požadavky kladené na plochu.) Pak bychom museli hlídat, aby přesah byl 1-nulový Pro jednoduchost však zde takový případ vyloučíme. 6.3. Stokesova věta. Nechť (G, ν) ⊂ R3 je 2-rozměrná orientovaná plocha, Ω je její omezená podplocha a Γ = Ω \ Ω. Předpokládejme, že Γ ⊂ G. Nechť každý bod Γ až na 1-nulovou množinu je regulární. Potom Γ je přípustná 1-rozměrná množina. Mějme okolí W množiny Ω a vektorové pole f : W ∪ Ω → Rn , které je je spojitě diferencovatelné na W a spojité na Ω. Potom Z Z f · τ Ω ds = curl f · ν dS, Γ
Ω
pokud integrály na obou stranách konvergují. Důkaz. Důkaz bude uveden v kapitole 7.
6.4. Poznámka. V širším smyslu, Stokesova věta je obecná věta o integrování diferenciálních forem na varietách. Jako taková zahrnuje větu o potenciálu, větu o divergenci, Greenovu větu i větu 6.3. Další aplikace už jsou v prostorech vyšší dimenze než tři a jejich nepovrchní studium vyžaduje pokročilou multilineární algebru. 7. Důkaz věty o divergenci a Stokesovy věty 7.1. Lemma. Nechť g2 , . . . , gn jsou dvakrát spojitě diferencovatelné funkce na otevřené množině W ⊂ Rn . Potom (20)
div(∇g2 × · · · × ∇gn ) = 0.
Důkaz. Uvažujme následující tvrzení: (20) platí, jsou-li gk , . . . , gn lineární polynomy. Toto tvrzení zřejmě platí pro k = 2, neboť pak ∇gi jsou konstanty, vektorový součin je konstantní vektor a jeho divergence je tudíž nulová. Indukcí dotáhneme platnost tvrzení až na k = n. Předpokládejme platnost pro k a chceme jej dokázat pro k+1, tedy g2 , . . . , gk nemusí být lineární polynomy. Máme n X ∂gk ei × · · · × ∇gn div(∇g2 × · · · × ∇gk × · · · × ∇gn ) = div ∇g2 × · · · × ∂xi i=1 = div
n X ∂g
k
i=1
=
n X i=1
+
∂xi
∇g2 × · · · × ei × · · · × ∇gn
∂gk div ∇g2 × · · · × ei × · · · × ∇gn ∂xi n X ∂gk ∇ · ∇g2 × · · · × ei × · · · × ∇gn . ∂xi i=1
Jelikož ei je gradient lineárního polynomu x 7→ xi , podle indukčního předpokladu je n X ∂gk div ∇g2 × · · · × ei × · · · × ∇gn = 0. ∂xi i=1 Podle věty o záměnnosti smíšených derivací je n n X ∂g X ∂gk k · ∇g2 × · · · × ei × · · · × ∇gn = det ∇ , ∇g2 , . . . , ei , . . . , ∇gn ∇ ∂xi ∂xi i=1 i=1 n X n X ∂ 2 gk = det ej , ∇g2 , . . . , ei , . . . , ∇gn = 0. ∂xj ∂xi i=1 j=1
Tím je završen indukční krok a celý důkaz lemmatu.
7.2. Lemma (Rozklad jednotky). Buď m ∈ N. Existuje systém funkcí {ωα : α ∈ I}, kde I = Zn , s následujícími vlastnostmi: (a) ωα , α ∈ I, jsou nekonečně diferencovatelné, 15
(b) (c) (d) (e) (f)
ω Pα ≥ 0, α ∈ I, α∈I ωα = 1 2 pro každé α ∈ I, nosič spt ωα je krychle o hraně menší než m , n n pro každý bod x ∈ R existuje nejvýš 2 indexů α tak, že x ∈ spt ωα , √ |∇ωα | ≤ 2m n, α ∈ I.
Důkaz. Buď η neklesající nekonečně diferencovatelná funkce, taková, že η 0 ≥ 2, η = 0 na okolí (−∞, 0] a η = 1 na okolí [1, ∞). Definujme ωk (t) = η(m(k + t − 1)) − η(m(k + t)),
t ∈ R.
Potom systém funkcí (ωk )k∈Z má požadované vlastnosti pro n = 1. Ve vyšší dimenzi použijeme α ∈ I.
ωα (x) = ωα1 (x1 ) . . . ωαn (xn ),
7.3. Věta (O aproximaci). Nechť W ⊂ Rn je otevřená množina, f : W → Rn je spojitě diferencovatelná funkce, K ⊂ W je kompaktní a ε > 0. Potom existuje nekonečně diferencovatelná funkce g : W → Rn tak, že |g(x) − f (x)| < ε, |∇g(x) − ∇f (x)| < ε, x ∈ K. Důkaz. Pokud tuto větu znáte odjinud (např. pomocí zhlazování), tento důkaz je možno přeskočit. Zvolme ε0 > 0. Najdeme konstantu M omezující ∇f na okolí K. Ze stejnoměrné spojitosti a kompaktnosti najdeme m ∈ N tak, že pro každou krychli Q o hraně 2/m protínající K platí Q ⊂ W,
oscQ ∇f < ε0 .
Nechť (ωα )α∈Zn je rozklad jednotky jako v Lemmatu 7.2. V průběhu důkazu budeme využívat jeho vlastnosti (a)–(f), někdy mlčky. Buď A množina všech α ∈ Zn takových, že spt ωα ∩ K 6= ∅. Označme Qα = spt ωα . Ke každému α ∈ A najdeme z = zα ∈ Qα ∩ K a položíme pα (x) = f (z) + ∇f (z)(x − z). Buď X g= ωα pα . α∈A
Potom (vlastnost (c)) g−f =
X
ωα (pα − f )
α∈A
a ∇g − ∇f =
X
∇ωα (pα − f ) +
X
α
ωα (∇pα − ∇f ).
α
Máme sup |∇pα − ∇f | ≤ oscQα (∇pα − ∇f ) ≤ ε0 . Qα
a
√ sup |pα − f | ≤ oscQα |pα − f | ≤ ε0 diam Qα ≤ 2 nε0 /m, Qα
tedy (vlastnost (f)) |∇ωα (pα − f )| + |ωα (∇pα − ∇f )| ≤ (4n + 1)ε0 a
√ |ωα (pα − f )| ≤ 2 nε0 /m
na Qα
na Qα .
Sečtením přes α (vlastnost (e)) dostaneme požadované odhady pro vhodně malé ε0 .
7.4. Lemma. Nechť W ⊂ Rn je otevřená množina a H = {x ∈ Rn : x1 < 0}. Nechť f je spojité vektorové pole na H ∩ W s kompaktním nosičem ve W . Nechť f je spojitě diferencovatelné na H ∩ W . Potom Z Z f · e1 dS = div f (x) dx, ∂H∩W
H∩W
pokud integrál vpravo konverguje. 16
Důkaz. Můžeme f = (f1 , . . . , fn ) dodefinovat nulou na H \ W a předpokládat, že W = Rn . Předpokládejme nejprve, že f je spojitě diferencovatelná dokonce na okolí H. Nechť (0, x2 , . . . , xn ) ∈ ∂H. Potom lim
x1 →−∞
f1 (x1 , x2 , . . . , xn ) = 0,
tedy 0
Z (21)
f1 (0, x2 , . . . , xn ) = −∞
∂f1 (x1 , . . . , xn ) dx1 . ∂x1
Vyintegrováním podle proměnných (x2 , . . . , xn ) (Fubiniova věta) dostanene Z Z ∂f1 (x) f · e1 dS = dx. ∂H H ∂x1 Nechť nyní i ∈ {2, . . . , n} a (x1 , . . . , xi−1 , xi+1 , . . . , xn ) ∈ Rn−1 , přičemž x1 < 0. Potom lim fi (x1 , . . . , xi , . . . , xn ) = 0,
xi →±∞
tedy Z
∞
0= −∞
∂fi (x1 , . . . , xi , . . . , xn ) dxi . ∂xi
Vyintegrováním podle zbývajících proměnných xj , j 6= i dostaneme Z ∂fi (x) (22) 0= dx. H ∂xi Sečteme-li (21) a (22) přes i = 2, . . . , n, dostaneme Z Z f · e1 dS = div f dx. ∂H
H
V obecném případě uvažujme posloupnost funkcí
(fk )∞ k=1
danou předpisem
fk (x) = f (x − 2−k e1 ). Potom z předchozího dostaneme Z
Z
Z
div f (x) dx = {x1 <−2−k }
fk · e1 dS.
div fk dx dx = ∂H
H
Jelikož f má kompaktní nosič, fk → f stejnoměrně na ∂H a integrál z div f přes H konverguje podle předpokladů, limitní přechod pro k → ∞ dává požadovaný výsledek. 7.5. Lemma. Nechť Ω ⊂ Rn a z ∈ Ω nebo z je regulární bod ∂Ω. Potom existuje okolí U bodu z tak, že pro každé spojité vektorové pole f na Ω s nosičem v U , spojitě diferencovatelné v Ω, platí Z Z f · ν Ω dS = div f (x) dx, ∂Ω
Ω
pokud integrál vpravo konverguje. Důkaz. Jestliže z ∈ Ω, okolí bodu z zvolíme celé v Ω a tvrzení dostaneme z lemmatu 7.4. Předpokládejme tedy, že z je regulární bod ∂Ω. Podle věty 5.3 najdeme hraniční parametrizaci, totiž otevřená množina G ⊂ Rn a prosté regulární zobrazení Ψ : G → Rn tak, že platí (14)–(17). Položme U = Ψ(G). Podle věty o lokálním difeomorfismu je U otevřená a Ψ je difeomorfismus. Nechť f je spojité vektorové pole na Ω s nosičem v U , spojitě diferencovatelné v Ω. Položme g(y) = f1 (Ψ(y)) ∇Ψ2 (y) × · · · × ∇Ψn (y),
y ∈ G.
Předpokládejme, že Ψ je dvakrát spojitě diferencovatelné. Potom div g = ∇(f1 ◦ Ψ) · ∇Ψ2 × · · · × ∇Ψn + (f1 ◦ Ψ) div ∇Ψ2 × · · · × ∇Ψn . Podle lemmatu 7.1 je div ∇Ψ2 × · · · × ∇Ψn = 0. 17
Zbývající člen upravíme ∇(f1 ◦ Ψ) · ∇Ψ2 × · · · × ∇Ψn = det(∇(f1 ◦ Ψ), ∇Ψ2 , . . . , ∇Ψn ) =
n X ∂f1
∂xi
i=1
= Pišme ν Ω =
1 n (νΩ , . . . , νΩ ).
∂f
1
∂x1
◦ Ψ det(∇Ψi , ∇Ψ2 , . . . , ∇Ψn )
◦ Ψ det(∇Ψ1 , ∇Ψ2 , . . . , ∇Ψn ).
Položme ˜ = {t ∈ Rn−1 : (0, t) ∈ G ∩ ∂H}, G ˜ Φ(t) = Ψ(0, t), t∈G
Potom podle (17) je pro y = (0, t) vektor ν Ω kladný násobek JΦ(t), takže ν Ω (Ψ(y))|JΦ(t)| = JΦ(t) a ∂Ψ (y) i 1 = e1 · ∇Ψ2 (y) × · · · × ∇Ψn (y) . νΩ (Ψ(y))|JΦ(t)| = det ∂yj i,j=2,...,n Podle lemmatu 7.4 s pomocí záměny proměnných a nezávislosti plošného integrálu na parametrizaci dostaneme Z Z Z 1 f1 νΩ dS = g · e1 dS = div g(y) dy ∂Ω∩W H∩G Z∂H∩G ∂f1 = Ψ(y) det ∇Ψ1 (y), ∇Ψ2 , . . . , ∇Ψn (y) dy ∂x ZH∩G 1 ∂f1 = (x) dx. Ω∩W ∂x1 Obecně s myšlenkově stejným důkazem Z Z i (23) fi νΩ dS = ∂Ω∩W
Ω∩W
∂fi (x) dx, ∂xi
i = 1, . . . , n
Sečteme-li (23) přes i dostaneme požadované, ale zatím v případě, že Ψ je dvakrát spojitě diferencovatelné. Pro obecný případ musíme aproximovat Ψ hladkými zobrazeními podle věty 7.3 a udělat limitní přechod v integrálu pro rovnost Z Z ∂f1 (f1 ◦ Ψ) ∇Ψ2 × · · · × ∇Ψn dS = Ψ(y) det ∇Ψ1 (y), ∇Ψ2 , . . . , ∇Ψn (y) dy. ∂x 1 ∂H∩G H∩G 7.6. Lemma. Nechť Ω ⊂ Rn , W ⊂ Rn jsou otevřené množiny, Ω je omezená. Nechť všechny body ∂Ω ∩ W jsou regulární pro Ω. Nechť f je spojité vektorové pole na Ω s kompaktním nosičem ve W , spojitě diferencovatelné v Ω. Potom Z Z f · ν Ω dS = div f (x) dx, ∂Ω
Ω
pokud integrál vpravo konverguje. Důkaz. Buď K ⊂ W ∩ Ω kompaktní. Podle lemmatu 7.5, ke každému bodu z ∈ K najdeme okolí U tak, že tvrzení platí pro každou f s nosičem v U . Z kompaktnosti odvodíme, že existuje δ > 0 tak, √ že pro každý bod z ∈ spt f platí tvrzení pro každou f s nosičem v B(z, δ). Najdeme m ∈ N tak, že n < mδ. Uvažujme rozklad jednotky (ωα )α∈Zn podle lemmatu 7.2 a pozorujeme, že pro každé α platí tvrzení pro f s nosičem v spt ωα . Omezme se na α ∈ A = {α ∈ Zn : spt ωα ∩ K 6= ∅}. Množina A je konečná. Zvolme f s nosičem v K a položme α ∈ A.
fα = ωα f ,
Potom tvrzení platí pro každou fα a tudíž i pro f jakožto jejich konečný součet. 18
7.7. Lemma. Nechť N ⊂ Rn je (n−1)-nulová kompaktní množina nula a ε > 0. Potom existuje posloupnost (vm )m spojitě diferencovatelných funkcí s kompaktním nosičem na Rn tak, že 0 ≤ vm ≤ 1, vm = 1 na okolí N , vm → χN a Z (24) |∇vm (x)| dx → 0. Rn
Důkaz. Nejprve si sestrojíme pomocnou spojitě diferencovatelnou funkci η tak, že 0 ≤ η ≤ 1, η = 0 na B(0, 1) a η = 1 vně B(0, 2). Označme Z c= |∇η(y)| dy. Rn
Zvolme m pevně. Najdeme posloupnost (B(xj , rj ))j koulí pokrývajících N tak, že X rjn−1 < 2−m . j
Jelikož N je kompaktní, existuje p ∈ N tak, že koule B(xj , rj ), j = 1, . . . , p, pokrývají N . Položme x − x j , ηj (x) = η rj vm (x) = 1 − η1 (x) . . . ηp (x). Pro každé j dostaneme záměnou proměnných Z |∇ηj (x)| dx = crjn−1 . Rn
Potom z věty o derivování součinu snadno dostaneme |∇vm | ≤
p X
|∇ηj |,
j=1
tedy Z |∇vm (x)| dx ≤ Rn
p Z X j=1
|∇ηj (x)| dx ≤ c
Rn
p X
rjn−1 ≤ c2−m .
j=1
Zřejmě vm = 0Pna okolí N . Jestliže x ∈ / N , existuje δ > 0 tak, že B(x, δ) ∩ N = ∅. Jestliže při konstrukci funkce vm je j rjn−1 < 2−m , znamená to rjn−1 < 2−m pro každou jednotlivou kouli B(xj , rj ). To ale znamená, že pro dostatečně velké m platí rj < δ, j = 1, . . . , p, takže pro všechna j máme |xj −x| ≥ δ > rj a ηj (x) = 1. Tedy vm (x) = 0. Tím jsme dokázali, že vm → χN . Důkaz věty 5.5. Nechť N je množina všech bodů ∂Ω, které nejsou regulární a W = Rn \ N . Podle lemmatu 7.7 najdeme posloupnost (vm )m spojitě diferencovatelných funkcí s kompaktním nosičem tak, že 0 ≤ vm ≤ 1, vm = 1 na okolí N a vm → χN . Buď m ∈ N pevné. Najděme otevřené okolí V množiny N tak, že vm = 1 na V a položme K = Ω \ V. Buď f m = (1 − vm )f . Potom f m má kompaktní nosič v K ⊂ W . Podle lemmatu 7.6 je Z Z Z Z Z (1−vm )f · ν Ω dS = f m · ν Ω dS = div f m (x) dx = (1−vm (x)) div f (x) dx − ∇vm (x) · f (x) dx. ∂Ω
∂Ω
Ω
Ω
Ω
Nyní provedeme limitní přechod. Jelikož integrály na obou stranách (18) konvergují, z Leesgueovy věty dostaneme Z Z (1−vm )f · ν Ω dS → f · ν Ω dS, ∂Ω Z Z∂Ω (1 − vm (x)) div f (x) dx → div f (x) dx. Ω
Ω
Zbývá odhadnout integrál Z ∇vm (x) · f (x) dx. Ω
Jelikož však Ω je omezená množina a f je omezená funkce (ze spojitosti na kompaktní množině Ω), tvrzení plyne z (??). 19
Důkaz věty 6.3. Stokesova věta se dá přepsat jako 3 Z 3 X 3 Z X X ∂fj fj dxj = dxi dxj , ∂xi j=1 Γ j=1 i=1 Ω kde používáme zápis f = (f1 , f2 , f3 ). Zvolme bod z ∈ G a najděme okolí U bodu z a globální parametrizaci Φ : G → G plochy U ∩ G. Nechť vektorové pole f má nosič v U . Buď G0 = G ∩ Φ−1 (Ω). Zvolme pevně j ∈ {1, 2, 3} a uvažujme na G vektorové pole g(t) = fj (Φ(t))∇Φj (t). Položme γ = G ∩ Φ−1 (Γ). Potom γ = G ∩ ∂G0 a 1-skoro všechny body γ jsou regulární vzhledem k G0 (rozhraničující funkce h ◦ Φ). Potom Z Z Z fj dxj = g · ds = (fj ◦ Φ) ∇Φj · ds. Γ
γ
γ
(Pro úpravy integrálu vlevo si uvědomme, že nezávisí na parametrizaci, můžeme tedy vzít parametrizaci tvaru Φ ◦ ψ, kde ψ je parametrizace γ.) Na pravé straně je Z X Z X 3 3 ∂fj ∂fj dxi dxj = (Φ(t)) det(∇Φi (t), ∇Φj (t)) dt. G0 i=1 ∂xi Ω i=1 ∂xi Potřebujeme tedy ověřit Z Z (25) (fj ◦ Φ) ∇Φj · ds =
3 X ∂fj
G0 i=1
γ
∂xi
(Φ(t)) det(∇Φi (t), ∇Φj (t)) dt.
Je-li Φ dvakrát spojitě diferencovatelná, použijeme greenovskou variantu lemmatu 7.6 a máme Z Z Z (fj ◦ Φ) ∇Φj · ds = g · ds = curl g(t) dt. γ
G0
γ
Máme curl g = det(∇(fj ◦ Φ), ∇Φj ) + (fj ◦ Φ) curl ∇Φj . Úpravou dostaneme det(∇(fj ◦ Φ)(t), ∇Φj (t)) =
3 X ∂fj i=1
∂xi
(Φ(t)) det(∇Φi (t), ∇Φj (t))
a curl ∇Φj = 0 (srov. věta 3.4), tedy (25) platí. Pokud Φ je pouze jednou spojitě diferencovatelné, pomocí věty 7.3 najdeme posloupnost (Φ(k) )k hladkých funkcí, které konvergují spolu se svými prvními derivacemi k Φ na okolí Φ−1 (z) a limitním přechodem v integrálu dostaneme též požadované. Uvědomme si, že při ověřování (25) jsme nikde nepoužili speciální vlastnosti zobrazení Φ jako prostotu nebo regularitu. Zatím jsme ověřili Stokesovu větu v případě, že f má nosič v malém okolí bodu z. Obecný případ odtud plyne metodou rozkladu jednotky jako v důkazu lemmatu 7.6. 8. Praktické hledání parametrizace a určování orientace 8.1. Poznámka. Sehranost orientací pro Greenovu a Stokesovu větu se heuristicky kontroluje pomocí názorných pomůcek. Kladná parametrizace regulární části hranice otevřené množiny Ω ⊂ R2 je “křivka”, která obíhá G proti směru hodinových ručiček. Kladná parametrizace regulární části relativní hranice 2rozměrné plochy G ⊂ R2 orientované pomocí normály se pozná podle pravidla pravé ruky: směřuje-li palec ve směru normály příslušné G, pak zakřivené prsty ukazují směr obíhání “křivky”, která parametrizuje regulární část hranice. Zde používáme konvenci, že osa x směřuje doprava, osa y dozadu a osa z nahoru. Tyto pomůcky nemůžou nahradit výpočet, ale mohou nám naznačit, zda jsme při výpočtu neudělali numerickou chybu. 8.2. Poznámka. Je-li Φ : G → Rn plocha, o níž chceme rozhodnout, zda ohraničuje množinu nebo zda je kladnou lokální parametrizací plochy, pak platí, že pokud G je souvislá (např. interval), stačí provést test orientace v jednom bodě. V některých následujících cvičeních budeme ze cvičných důvodů ověřovat znaménko ve všech bodech. Samostatně zkontrolujte, zda nalezené parametrizace jsou prostá regulární zobrazení do dané množiny a nepokrytá čast je nulová. 20
8.3. Poznámka. Pokud nám test orientace dá, že nalezená parametrizace je záporná, nezoufejme. Kladnou parametrizaci lze vyrobit prohozením pořadí proměnných (u plochy) nebo záměnou proměnných s = −t (u křivky). Také můžeme příklad dopočítat a u výsledku otočit znaménko. p 8.4. Cvičení. Nechť 0 < r < R a M = {( x2 + y 2 − R)2 + z 2 = r2 }. Najděte zobecněnou parametrizaci a rozhodněte o znaménku, víte-li, že kladná jednotková normála v bodě [R + r, 0, 0] je [1, 0, 0]. Řešení. Použijeme-li válcové souřadnice Φ: x = ρ¯ cos α ¯, y = ρ¯ sin α ¯,
[¯ ρ, α ¯ , z¯] ∈ (0, ∞) × (−π, π) × R,
z = z¯, rovnice se nám převede na (¯ ρ − R)2 + z¯2 = r2 . Tuto plochu můžeme parametrizovat posunutými polárními souřadnicemi ψ: ρ¯ = R + r cos β, [α, β] ∈ (−π, π)2 ,
z¯ = r sin β, α ¯ = α, takže složením parametrizací Φ ◦ ψ dostáváme Φ: x = (R + r cos β) cos α, y = (R + r cos β) sin α,
[α, β] ∈ (−π, π)2 .
z = r sin β, Znaménko parametrizace určíme z pravidla, že vektorový jakobián kladné parametrizace v [α, β] je kladným násobkem jednotkové normály v Φ(α, β). Stačí tedy kontrolovat x-ovou souřadnici JΦ(α, β), a to je ∂(y, z) (R + r cos β) cos α, −r sin β sin α (α, β) = det . 0, r cos β ∂(α, β) Jelikož náš bod [R + r, 0, 0] je Φ(0, 0), počítáme ∂(y, z) R + r, 0 (0, 0) = det = r(R + r) > 0. 0, r ∂(α, β) Tedy nalezená parametrizace je kladná.
8.5. Cvičení. Nechť G = {[x, y, z] : x2 + y 2 + z 2 = 1 a ν(x, y, z) = [x, y, z], [x, y, z] ∈ G. Spočtěte R x dy dz. G Řešení. Použijeme-li sférické souřadnice x = r¯ cos γ¯ cos α ¯, y = r¯ cos γ¯ sin α ¯, z = r¯ sin γ¯ ,
π π [¯ r, α ¯ , γ¯ ] ∈ (0, ∞) × (−π, π) × (− , ), 2 2
dostaneme z rovnice g = 0 podmínku r¯ = 1. Tedy zvolíme zobecněnou parametrizaci Φ, x = cos γ cos α, π π y = cos γ sin α, [α, γ] ∈ G := (−π, π) × (− , ). 2 2 z = sin γ, Potom G \ Φ(G) = G ∩ {x < 0} ∩ {y = 0}, což je 2-nulová množina. Pro [x, y, z] = Φ(α, γ) máme ∇g(x, y, z) = [2x, 2y, 2z] = [2 cos γ cos α, 2 cos γ sin α, 2 sin γ]. Počítejme 2 − cos γ sin α − sin γ cos α cos γ cos α ∂Φ ∂Φ cos γ cos α × − sin γ sin α = cos2 γ sin α , × = ∂α ∂γ 0 cos γ cos γ sin γ 21
takže
∂Φ
∂Φ (α, γ) je kladný násobek ∇g(Φ(α, γ)). ∂α ∂γ Test orientace prošel, parametrizace je kladná. Z Z ∂(cos γ sin α, sin γ) x dy dz = cos γ cos α dα dγ ∂(α, γ) G G Z π/2 Z π 4 = cos3 γ cos2 α dα dγ = π . 3 −π/2 −π ×
8.6. Cvičení. Nechť G = {[x, y, z] : x2 + y 2 + z 2 = 1 a ν(x, y, z) = [x, y, z], [x, y, z] ∈ G. Nechť √ 3 h(x, y, z) = 2 − x. Nechť Ω = {[x, y, z] ∈ G : h(x) < 0}. Najděte zobecněnou křivku, která je kladnou parametrizací ∂Ω. Řešení. Plocha G je daná implicitně rovnicí g = x2 + y 2 + z 2 − 1 = 0, rozhraničující funkce je h. Hledaná křivka má parametrizovat plochu g = h = 0. Použijeme-li válcové souřadnice x=x ¯, y = ρ¯ cos α ¯,
[¯ ρ, α ¯, x ¯] ∈ (0, ∞) × (−π, π) × R,
z = ρ¯ sin α ¯, vyjádříme danou soustavu rovnic jako x ¯2 + ρ¯2 = 1, √ x ¯=
3 . 2
Odtud dostaneme zobecněnou parametrizaci Φ: √
x= y= z=
3 2 , 1 2 cos α, 1 2 sin α,
α ∈ (−π, π). √
Zobecněná křivka Φ pokrývá {g = h = 0} až na bod [ 23 , − 12 , 0]. Máme √ 1 1 ∇g = [2x, 2y, 2z] = [ 3, cos α, sin α], ∇h = [−1, 0, 0], Φ0 = [0, − sin α, cos α]. 2 2 Počítejme x −1 0 ν(x, y, z) × ∇h(x, y, z) = y × 0 = −z , 0 y z 0 tedy Φ (α) = ν(Φ(α)) × ∇h(Φ(α)) Test orientace prošel, takže nalezená parametrizace je kladná.
22