20.6 Integralsätze

Wir wollen Pendants zum Hauptsatz (18.2.2) für mehrfach-Integrale finden.




20.6.1 Proposition. Länge von Kurven.
Es sei \bgroup\color{demo}$ c:[a,b]\to E$\egroup eine Kurve. Unter der Länge \bgroup\color{demo}$ L_a^b(c)$\egroup der Kurve versteht man das Supremum der Länge aller interpolierenden Polygonzüge, d.h.

$\displaystyle L_a^b(c)$ $\displaystyle := \sup\{L(c, Z): Z$ ist Zerlegung $\displaystyle [a,b]\}$    
$\displaystyle L(c, Z)$ $\displaystyle := \sum_{i=1}^m \vert c(t_i)-c(t_{i-1})\vert$ für $\displaystyle Z=\{a=t_0<\dots<t_m=b\}$    

Falls \bgroup\color{demo}$ c$\egroup \bgroup\color{demo}$ C^1$\egroup ist, so ist

\bgroup\color{demo}$\displaystyle L_a^b(c) = \int_a^b \vert c'\vert.
$\egroup

Image ..//pic-6-05.gif

Beweis. Wegen dem Hauptsatz (18.2.2) und dem Mittelwertsatz (18.1.10) der Integralrechnung ist

$\displaystyle L(c, Z)$ $\displaystyle =\sum_{i=1}^m \vert c(t_i)-c(t_{i-1})\vert =\sum_{i=1}^m \vert\int_{t_{i-1}}^{t_i} c'(t)\vert dt$    
  $\displaystyle \leq \sum_{i=1}^m \int_{t_{i-1}}^{t_i} \vert c'(t)\vert dt = \int_{a}^{b} \vert c'(t)\vert dt,$    

also

\bgroup\color{demo}$\displaystyle L_a^b(c) = \sup\{L(c, Z): Z\}\leq \int_a^b \vert c'\vert.
$\egroup

Somit ist

$\displaystyle \vert\frac{c(t+h)-c(t)}{h}\vert \leq \vert\frac1h L_t^{t+h}(c)\ve...
..._a^{t+h}(c)-L_a^t{c}}{h} \leq \vert\frac1h\int_t^{t+h}\vert c'(s)\vert ds\vert$    

und Grenzübergang für \bgroup\color{demo}$ t\to 0$\egroup liefert

\bgroup\color{demo}$\displaystyle \frac{d}{dt}L_a^t(c)=\lim_{h\to 0} \frac{L_a^{t+h}(c)-L_a^t{c}}{h}
= \vert c'(t)\vert.
$\egroup

Nun folgt das Resultat aus dem Hauptsatz.     []


Geometrische Interpretation des Kurvenintegrals.

Sei nun \bgroup\color{demo}$ f=(f_1,\dots,f_n):B\to \mathbb{R}^n$\egroup ein Vektorfeld und \bgroup\color{demo}$ \tilde f=\sum_i f_i(x) dx^i$\egroup die zugehörige 1-Form \bgroup\color{demo}$ \tilde f:\mathbb{R}^n\supseteq B\to L(\mathbb{R}^n,\mathbb{R})\cong\mathbb{R}^n$\egroup. Nach der Kettenregel (17.2.3) hängt das Kurvenintegral nicht von der Parametrisierung (sondern nur deren Orientierung) ab, wir können \bgroup\color{demo}$ \int_B \tilde f := \int_c \tilde f$\egroup setzen, und brauchen uns bei \bgroup\color{demo}$ B$\egroup nicht die Parametrisierung sondern nur die Durchlaufungsrichtung merken. Weiters ist

$\displaystyle \int_B f$ $\displaystyle := \int_B \tilde f = \int_c \sum_i f_i(x) dx_i = \int_a^b \sum_i f_i(c(t)) c_i'(t) dt$    
  $\displaystyle = \int_a^b \langle f(c(t)), c'(t)\rangle dt \overset{c'(t)\ne 0}...
...t}}\rangle  \undersetbrace{=:\operatorname{vol}_c(t)}\to{\Vert c'(t)\Vert dt}$    
  $\displaystyle = \int_c \langle f\o c,\tau _c \rangle  \operatorname{vol}_c$    
  $\displaystyle =: \int_B \langle f,\tau _B \rangle  \operatorname{vol}_B$ falls $ c$ injektiv ist$\displaystyle .
$    

Dabei bezeichnet \bgroup\color{demo}$ \tau _B(c(t)):=\tau _c(t):=\frac{c'(t)}{\Vert c'(t)\Vert}$\egroup den Einheitstangentialvektor an \bgroup\color{demo}$ B$\egroup im Punkte \bgroup\color{demo}$ c(t)$\egroup und \bgroup\color{demo}$ \operatorname{vol}_B$\egroup ist nur ein Symbol (wie \bgroup\color{demo}$ dt$\egroup im 1-dimensionalen Riemannintegral) und heißt des Längenelement oder auch 1-dimensionale Volumselement. Beachte, daß

\bgroup\color{demo}$\displaystyle \int_I \operatorname{vol}_c = \int_I \Vert c'(t)\Vert dt
$\egroup

die Länge von \bgroup\color{demo}$ c$\egroup ist, und \bgroup\color{demo}$ \operatorname{vol}_c(t)$\egroup somit die Infinitesimale Version der Länge ist:

\bgroup\color{demo}$\displaystyle \operatorname{vol}_{c}(t) = \Vert c'(t_0)\Vert...
...\vert\to 0, t_0\in I} \frac1{\vert I\vert} \int_I \Vert c'(t)\Vert dt.
$\egroup

Beachte weiters, daß \bgroup\color{demo}$ \tau _B$\egroup das Vorzeichen ändert, wenn die Parametrisierung umgekehrt wird.


20.6.2 Der Gauß'sche Integralsatz in der Ebene.
Es sei \bgroup\color{demo}$ B=M(\varphi _1,\varphi _2)$\egroup mit \bgroup\color{demo}$ \varphi _i:[a,b]\to\mathbb{R}$\egroup eine Ordinaten-Menge mit \bgroup\color{demo}$ C^1$\egroup-Grenzen \bgroup\color{demo}$ \varphi _1$\egroup und \bgroup\color{demo}$ \varphi _2$\egroup. Es sei \bgroup\color{demo}$ (p,q)$\egroup ein \bgroup\color{demo}$ C^1$\egroup-Vektorfeld auf \bgroup\color{demo}$ B$\egroup und \bgroup\color{demo}$ (x,y)\mapsto p(x,y) dx + q(x,y) dy$\egroup die zugehörige 1-Form. Der Rand von \bgroup\color{demo}$ B$\egroup wird durch folgende 4 Kurven positiv orientiert parametrisiert:

$\displaystyle c_1(t)$ $\displaystyle := \binom{t}{\varphi _1(t)},$    
$\displaystyle c_2(t)$ $\displaystyle := \binom{b}{t},$    
$\displaystyle c_3(t)$ $\displaystyle := \binom{t}{\varphi _2(t)}$ verkehrt durchlaufen$\displaystyle ,$    
$\displaystyle c_4(t)$ $\displaystyle := \binom{a}{t}$ verkehrt durchlaufen$\displaystyle .
$    

Also ist

$\displaystyle \int_{\d B} p(x,y) dx$ $\displaystyle = \int_{c_1} + \undersetbrace{=0}\to{\int_{c_2}} - \int_{c_3} - \undersetbrace{=0}\to{\int_{c_4}}$    
  $\displaystyle = \int_a^b p(t,\varphi _1(t))\cdot 1\cdot dt - \int_a^b p(t,\varphi _2(t))\cdot 1\cdot dt$    
  $\displaystyle = \int_a^b [p(t,y)]_{y=\varphi _2(t)}^{y=\varphi _1(t)}  dt$    
  $\displaystyle = \int_a^b \int_{\varphi _2(t)}^{\varphi _1(t)} \d _2 p(t,y) dy dt$    
  $\displaystyle = -\int_a^b \int_{\varphi _1(t)}^{\varphi _2(t)} \d _2 p(t,y) dy dt$    

Analog erhalten wir falls $ B$ auch Ordinatenmenge bzgl. der anderen Achse ist


  $\displaystyle = -\int_B \d _2 p.\int_{\d B} q(x,y) dy$ $\displaystyle = \int_B \d _1 q$    

und somit


$\displaystyle \int_{\d B} p(x,y) dx + q(x,y) dy$ $\displaystyle = \int_B \left(\frac{\d q}{\d x}-\frac{\d p}{\d y}\right)(x,y) d(x,y).$    

Das zeigt den




Gauß'scher Integralsatz in der Ebene.
Es sei \bgroup\color{demo}$ B$\egroup eine Ordinatenmenge bezüglich beiden Achsen mit \bgroup\color{demo}$ C^1$\egroup Rändern. Weiters sei \bgroup\color{demo}$ \d B$\egroup der positiv orientierte Rand. Das Vektorfeld \bgroup\color{demo}$ (p,q):B\to\mathbb{R}^2$\egroup sei \bgroup\color{demo}$ C^1$\egroup. Dann ist

\bgroup\color{demo}$\displaystyle \int_{\d B} p  dx + q  dy = \int_B \d _1 q - \d _2 p.{\rm\quad[]}
$\egroup

Insbesonders folgt mit \bgroup\color{demo}$ p(x,y):=-y$\egroup und \bgroup\color{demo}$ q(x,y):=x$\egroup




20.6.3 Folgerung. Fläche via Kurvenintegral.
Es sei \bgroup\color{proclaim}$ B$\egroup eine Ordinatenmenge bezüglich beiden Achsen und stetiger Funktionen beschränkter Variation. Weiters sei \bgroup\color{proclaim}$ \d B$\egroup der positiv orientierte Rand. Dann ist

\bgroup\color{proclaim}$\displaystyle \vert B\vert = \frac12 \int_{\d B} (x dy-y dx).{\rm\quad[]}
$\egroup

Ist \bgroup\color{proclaim}$ B$\egroup gegeben durch alle Punkte mit Polarkoordinaten \bgroup\color{proclaim}$ (r,\varphi )$\egroup die \bgroup\color{proclaim}$ 0\leq r\leq r(\varphi )$\egroup erfüllen, so ist

\bgroup\color{proclaim}$\displaystyle \vert B\vert = \frac12 \int_{0}^{2\pi} r(\varphi )^2 d\varphi .
$\egroup

Beachte, daß also \bgroup\color{proclaim}$ x dy-y dx$\egroup eine somit recht nützliche aber dennoch nicht-exakte 1-Form ist.

Image ..//pic-6-06.gif
Beweis Der Rand \bgroup\color{proclaim}$ \d B$\egroup wird in der letztgenannten Situation durch die Kurve \bgroup\color{proclaim}$ c:t\mapsto (r(t) \cos(t),r(t) \sin(t))$\egroup mit Ableitung

\bgroup\color{proclaim}$\displaystyle c'(t)=(r'(t) \cos(t)-r(t) \sin(t),r'(t) \sin(t)+r(t) \cos(t))
$\egroup

beschrieben. Somit ist

$\displaystyle \int_{\d B} (x dy-y dx)$ $\displaystyle = \int_0^{2\pi} r(t) \cos(t) \Bigl(r'(t) \sin(t)+r(t) \cos(t)\Bigr)$    
  $\displaystyle \qquad -r(t) \sin(t) \Bigl(r'(t) \cos(t)-r(t) \sin(t)\Bigr) dt$    
  $\displaystyle = \int_0^{2\pi} r(t)^2 (\cos(t)^2+\sin(t)^2) dt =\int_0^{2\pi} r(\varphi )^2 d\varphi .{\rm\quad[]}$    

Für die Kreisscheibe \bgroup\color{proclaim}$ B$\egroup mit Radius \bgroup\color{proclaim}$ r$\egroup ergibt sich daraus nun sehr einfach

\bgroup\color{proclaim}$\displaystyle \vert B\vert=\frac12\int_0^{2\pi} r^2 d\varphi =r^2  \pi.
$\egroup


Geometrische Interpretation.
Für das Vektorfeld \bgroup\color{proclaim}$ f=(f_1,f_2)=(p,q):\mathbb{R}^2\supseteq B\to\mathbb{R}^2$\egroup definieren wir eine Funktion

\bgroup\color{proclaim}$\displaystyle \operatorname{rot}f := \frac{\d q}{\d x} -...
...{\d y} = \d _1 f_2 - \d _2 f_1 :
\mathbb{R}^2\supseteq B\to \mathbb{R},
$\egroup

die Rotation von $ f$. Der Gauß'sche Integralsatz ließt sich dann wie folgt:

\bgroup\color{proclaim}$\displaystyle \int_{\d B} \langle f,\tau _{\d B}\rangle \operatorname{vol}_{\d B} = \int_B \operatorname{rot}f,
$\egroup

wobei \bgroup\color{proclaim}$ \tau _{\d B}$\egroup der Einheitstangentialvektor an \bgroup\color{proclaim}$ \d B$\egroup und \bgroup\color{proclaim}$ \operatorname{vol}_{\d B}$\egroup das 1-dimensionale Volumselement von \bgroup\color{proclaim}$ \d B$\egroup ist. Die linke Seite ist die Arbeit von \bgroup\color{proclaim}$ f$\egroup längs \bgroup\color{proclaim}$ \d B$\egroup (das Mittel des Tangentialanteils von \bgroup\color{proclaim}$ f$\egroup) und somit ist

\bgroup\color{proclaim}$\displaystyle \operatorname{rot}f(x) = \lim_{\vert B\vert\to 0,x\in B} \frac1{\vert B\vert} \int_B \operatorname{rot}f$\egroup die Wirbeldichte\bgroup\color{proclaim}$\displaystyle .
$\egroup

Diese mißt wie sehr der Fluß des Vektorfelds um den Punkt \bgroup\color{proclaim}$ x$\egroup herumfließt.
Image ..//pic-6-07.gif

Es sei \bgroup\color{proclaim}$ g:=(g_1,g_2)=(q,-p)$\egroup, d.h. \bgroup\color{proclaim}$ f=g^{\perp}$\egroup, der positiv gedrehte Normalvektor. Dann ist

\bgroup\color{proclaim}$\displaystyle \operatorname{rot}f = \frac{\d q}{\d x} - \frac{\d p}{\d y} = \d _1 g_1 + \d _2 g_2 =: \operatorname{div}g,
$\egroup

und

\bgroup\color{proclaim}$\displaystyle \langle f,\tau _{\d B} \rangle = \langle g...
...= \langle g, - \tau _{\d B}^\perp\rangle = \langle g,\nu_{\d B}\rangle,
$\egroup

wobei \bgroup\color{proclaim}$ \nu_{\d B}$\egroup der nach außen weisende Einheits-Normalvektor ist. Der Gauß'sche Integralsatz hat nun also die Form

\bgroup\color{proclaim}$\displaystyle \int_{\d B} \langle g,\nu_{\d B}\rangle \operatorname{vol}_{\d B} = \int_B \operatorname{div}f.
$\egroup

Die linke Seite mißt nun die Durchflußmenge durch \bgroup\color{proclaim}$ \d B$\egroup und somit ist

\bgroup\color{proclaim}$\displaystyle \operatorname{div}f(x) = \lim_{B\to 0,x\in B} \frac1{\vert B\vert} \int_B \operatorname{div}f$\egroup die Quelldichte\bgroup\color{proclaim}$\displaystyle .
$\egroup

Image ..//pic-6-08.gif

Das Vektorfeld \bgroup\color{proclaim}$ f$\egroup heißt wirbelfrei \bgroup\color{proclaim}$ :\Leftrightarrow$\egroup \bgroup\color{proclaim}$ \operatorname{rot}f=0$\egroup. Es heißt quellenfrei \bgroup\color{proclaim}$ :\Leftrightarrow$\egroup \bgroup\color{proclaim}$ \operatorname{div}f=0$\egroup.


20.6.4 Definition. Fläche in $ \mathbb{R}^3$.
Unter einer Fläche \bgroup\color{proclaim}$ S$\egroup im \bgroup\color{proclaim}$ \mathbb{R}^3$\egroup verstehen wir eine Teilmenge \bgroup\color{proclaim}$ S\subseteq \mathbb{R}^3$\egroup zusammen mit einer \bgroup\color{proclaim}$ C^1$\egroup-Parametrisierung \bgroup\color{proclaim}$ \Phi =(x,y,z):K\to S\subseteq\mathbb{R}^3$\egroup mit \bgroup\color{proclaim}$ K\subseteq \mathbb{R}^2$\egroup \bgroup\color{proclaim}$ J$\egroup-meßbar und kompakt.


20.6.5 Der Stokes'sche Integralsatz im Raum.
Es sei \bgroup\color{proclaim}$ f=(p,q,r)$\egroup ein \bgroup\color{proclaim}$ C^1$\egroup-Vektorfeld auf \bgroup\color{proclaim}$ \Phi (K)$\egroup, mit \bgroup\color{proclaim}$ \Phi $\egroup \bgroup\color{proclaim}$ C^2$\egroup (!) und der Rand \bgroup\color{proclaim}$ \d K$\egroup sei durch eine \bgroup\color{proclaim}$ C^1$\egroup-Kurve \bgroup\color{proclaim}$ c=(u,v):[a,b]\to \mathbb{R}^2$\egroup parametrisiert. Mit \bgroup\color{proclaim}$ \d S$\egroup bezeichnen wir die durch \bgroup\color{proclaim}$ \Phi \o c$\egroup parametrisierte Kurve \bgroup\color{proclaim}$ \Phi (\d K)$\egroup. Dies ist nicht der topologische Rand von \bgroup\color{proclaim}$ S\subseteq \mathbb{R}^3$\egroup. Dann ist

$\displaystyle \int_{\d S} p dx$ $\displaystyle = \int_{\Phi \o c} p(x,y,z) dx = \int_a^b p(x(c(t)),y(c(t)),z(c(t))) (x\o c)'(t) dt$    
  $\displaystyle = \int_a^b p\cdot \left(\frac{\d x}{\d u}\cdot \frac{du}{dt} +\fr...
...\int_c p\cdot \left(\frac{\d x}{\d u}\cdot du +\frac{\d x}{\d v}\cdot dv\right)$    
  $\displaystyle \overset{\text{Gauß{}}}= \int_K \left(\frac{\d }{\d u}\left(p\cdo...
...}{\d v}\right) - \frac{\d }{\d v}\left(p\cdot \frac{\d x}{\d u}\right)\right) =$    

Für den Integranden erhalten wir, da $ \Phi $ $ C^2$ ist, nach dem Satz von Schwarz


  $\displaystyle \frac{\d }{\d u}\left(p\cdot \frac{\d x}{\d v}\right) - \frac{\d }{\d v}\left(p\cdot \frac{\d x}{\d u}\right) =$    
  $\displaystyle \quad = \frac{\d p}{\d u}\cdot \frac{\d x}{\d v} + p\cdot \frac{\...
... v} -\frac{\d p}{\d v}\cdot \frac{\d x}{\d u} - p\cdot \frac{\d ^2 x}{\d v\d u}$    
  $\displaystyle \quad = \left(\frac{\d p}{\d x}\cdot \frac{\d x}{\d u} +\frac{\d ...
...}{\d u} +\frac{\d p}{\d z}\cdot \frac{\d z}{\d u}\right)\cdot \frac{\d x}{\d v}$    
  $\displaystyle \qquad -\left(\frac{\d p}{\d x}\cdot \frac{\d x}{\d v} +\frac{\d ...
...}{\d v} +\frac{\d p}{\d z}\cdot \frac{\d z}{\d v}\right)\cdot \frac{\d x}{\d u}$    
  $\displaystyle \quad = \frac{\d p}{\d y}\Bigl( \undersetbrace{-\det\frac{\d (x,y...
...{\frac{\d z}{\d u}\frac{\d x}{\d v} - \frac{\d z}{\d v}\frac{\d x}{\d u}}\Bigr)$    

somit gilt weiter


  $\displaystyle = \int_K \left(\frac{\d p}{\d z}\cdot\det\frac{\d (z,x)}{\d (u,v)} - \frac{\d p}{\d y}\cdot\det\frac{\d (x,y)}{\d (u,v)} \right)  d(u,v) =$    

Wir setzen


  $\displaystyle \quad dz\wedge dx := \det\frac{\d (z,x)}{\d (u,v)}, \quad dx\wedge dy := \det\frac{\d (x,y)}{\d (u,v)},\quad\dots$    

und erhalten schließlich


  $\displaystyle = \int_K \frac{\d p}{\d z} dz\wedge dx - \frac{\d p}{\d y}  {dx\wedge dy}$    

und analog oder durch zyklisches Vertauschen


$\displaystyle \int_{\d S} q dy$ $\displaystyle = \int_K \frac{\d q}{\d x}  dx\wedge dy -\frac{\d q}{\d z} dy\wedge dz$    
$\displaystyle \int_{\d S} r dz$ $\displaystyle = \int_K \frac{\d r}{\d y}  dy\wedge dz -\frac{\d r}{\d x} dz\wedge dx$    

und insgesamt also


$\displaystyle \int_{\d S} f$ $\displaystyle = \int_K \left(\tfrac{\d q}{\d x}-\tfrac{\d p}{\d y}\right)  dx\...
...dy\wedge dz + \left(\tfrac{\d p}{\d z}-\tfrac{\d r}{\d x}\right)  dz\wedge dx.$    

Zum \bgroup\color{proclaim}$ C^1$\egroup-Vektorfeld \bgroup\color{proclaim}$ f=(p,q,r):\mathbb{R}^3\supseteq S\to \mathbb{R}^3$\egroup definieren wir nun ein neues Vektorfeld

\bgroup\color{proclaim}$\displaystyle \operatorname{rot}f = \operatorname{rot}(p...
... p}{\d z}-\frac{\d r}{\d x}, \frac{\d q}{\d x}-\frac{\d p}{\d y}\Bigr),
$\egroup

die Rotation des Vektorfelds. Als Memotechnik ist

$\displaystyle \operatorname{rot}f$ $\displaystyle = \left(\det\left(\begin{matrix}\frac{\d }{\d y} & q  \frac{\d ...
...matrix}\frac{\d }{\d x} & p  \frac{\d }{\d y} & q \end{matrix}\right) \right)$    
  $\displaystyle = \det\left(\begin{matrix}e_1 & \frac{\d }{\d x} & p  e_2 & \fr...
...{\d y} & q  e_3 & \frac{\d }{\d z} & r \end{matrix}\right) = \nabla \times f,$    

wobei


$\displaystyle \nabla$ $\displaystyle := \left(\begin{matrix}\frac{\d }{\d x}  \frac{\d }{\d y}  \frac{\d }{\d z} \end{matrix}\right)$   und$\displaystyle \quad a\times b := \det \left(\begin{matrix}e_1 & a_1 & b_1  e_2 & a_2 & b_2  e_3 & a_3 & b_3 \end{matrix}\right)$    

das Kreuzprodukt von Vektoren im \bgroup\color{proclaim}$ \mathbb{R}^3$\egroup ist. Dies ist durch folgende Eigenschaften geometrisch charakterisiert: \bgroup\color{proclaim}$ a\times b\perp a,b$\egroup, weiters ist \bgroup\color{proclaim}$ (a\times b,a,b)$\egroup positiv orientiert, und \bgroup\color{proclaim}$ \vert a\times b\vert$\egroup ist die Fläche des Parallelogramms, das von \bgroup\color{proclaim}$ (a,b)$\egroup aufgespannt wird: Es ist

$\displaystyle \langle v,a\times b\rangle$ $\displaystyle = \det( v,a,b) \Rightarrow a\times b\perp a,b$    
$\displaystyle \det(a\times b,a,b)$ $\displaystyle = \det\left(\begin{matrix}\det\left(\begin{matrix}a_2 & b_2  a_...
...trix}a_1 & b_1  a_2 & b_2 \end{matrix}\right) & a_3 & b_3 \end{matrix}\right)$    
  $\displaystyle = \det\left(\begin{matrix}a_2 & b_2  a_3 & b_3 \end{matrix}\rig...
...right)^2 + \det\left(\begin{matrix}a_1 & b_1  a_2 & b_2 \end{matrix}\right)^2$    
  $\displaystyle = \vert a\times b\vert^2 \geq 0$    

und andererseits


$\displaystyle \det(a\times b,a,b)$ $\displaystyle = \vert a\times b\vert\cdot$Fläche$\displaystyle (a,b)$    
$\displaystyle \Rightarrow \vert a\times b\vert$ $\displaystyle =$Fläche$\displaystyle (a,b)$    

Es sei \bgroup\color{proclaim}$ g=(a,b,c):\mathbb{R}^3\supseteq S\to\mathbb{R}^3$\egroup ein stetiges Vektorfeld (wie z.B. \bgroup\color{proclaim}$ g=\operatorname{rot}F$\egroup). Wir definieren ein Integral von \bgroup\color{proclaim}$ g$\egroup über \bgroup\color{proclaim}$ \Phi :K\to S$\egroup als

$\displaystyle \int_\Phi g$ $\displaystyle := \int_K a\o\Phi   dy\wedge dz + b\o\Phi  dz\wedge dx + c\o\Phi   dx\wedge dy$    
mit z.B. $\displaystyle \quad dy\wedge dz$ $\displaystyle := \det\left(\frac{\d (y,z)}{\d (u,v)}\right).$    

Dieses Integral ist Reparametrisierungsinvariant, denn wenn \bgroup\color{proclaim}$ h:(r,s)\mapsto(u,v)$\egroup eine Orientierungs-erhaltende (d.h. \bgroup\color{proclaim}$ \det h'(r,s)\geq 0$\egroup) Substitutionsfunktion \bgroup\color{proclaim}$ h:K'\to K$\egroup ist, dann ist wegen der Kettenregel

\bgroup\color{proclaim}$\displaystyle \det\left(\frac{\d (y,z)}{\d (r,s)}\right)...
...,z)}{\d (u,v)}\right) \cdot \det\left(\frac{\d (u,v)}{\d (r,s)}\right)
$\egroup

und somit wegen der Transformationsformel

\bgroup\color{proclaim}$\displaystyle \int_{\Phi } g = \int_{\Phi \o h} g
$\egroup

und wir können

\bgroup\color{proclaim}$\displaystyle \int_S g := \int_\Phi g
$\egroup

setzen, müssen uns dabei aber die ``Orientierung'' von \bgroup\color{proclaim}$ S$\egroup merken. Damit erhalten wir schließlich den




Satz von Stokes für Flächen im 3-dimensionalen.
Es sei \bgroup\color{proclaim}$ K\subseteq \mathbb{R}^2$\egroup eine Ordinatenmenge bezüglich beider Achsen mit \bgroup\color{proclaim}$ C^1$\egroup-Grenzen und Rand \bgroup\color{proclaim}$ \d K$\egroup der durch eine \bgroup\color{proclaim}$ C^1$\egroup-Kurve \bgroup\color{proclaim}$ c:[a,b]\to\d K$\egroup parametrisiert werden kann. Weiters sei \bgroup\color{proclaim}$ \Phi :\mathbb{R}^2\supseteq K\to \mathbb{R}^3$\egroup \bgroup\color{proclaim}$ C^2$\egroup und \bgroup\color{proclaim}$ f$\egroup ein \bgroup\color{proclaim}$ C^1$\egroup-Vektorfeld auf \bgroup\color{proclaim}$ S:=\Phi (K)$\egroup. Dann ist

\bgroup\color{proclaim}$\displaystyle \int_{\d S} f = \int_S \operatorname{rot}f.{\rm\quad[]}
$\egroup


Geometrische Interpretation.
Nach obigen können wir die linke Seite geometrisch als

\bgroup\color{proclaim}$\displaystyle \int_{\d S}\langle f,\tau _{\d S}\rangle \operatorname{vol}_{\d S}
$\egroup

interpretieren. Nun zur rechten Seite: Es ist die zugehörige 2-Form zum Vektorfeld \bgroup\color{proclaim}$ \operatorname{rot}f$\egroup gerade

\begin{multline*}
\left(\tfrac{\d q}{\d x}-\tfrac{\d p}{\d y}\right)  dx\wedge ...
...c{\d\Phi }{\d u}\times \tfrac{\d\Phi }{\d v}\Bigr\rangle d(u,v),
\end{multline*}

denn

\bgroup\color{proclaim}$\displaystyle \Phi '(u,v) = \frac{\d (x,y,z)}{\d (u,v)}
...
...y}{\d v} \\
\frac{\d z}{\d u} & \frac{\d z}{\d v} \end{matrix}\right).
$\egroup

Somit ist

\bgroup\color{proclaim}$\displaystyle \int_S \operatorname{rot}f = \int_K \langle \operatorname{rot}f,\nu_\Phi \rangle \operatorname{vol}_\Phi ,
$\egroup

wobei

\bgroup\color{proclaim}$\displaystyle \nu_\Phi (u,v) := \frac{\tfrac{\d\Phi }{\d...
... := \Vert\tfrac{\d\Phi }{\d u}\times \tfrac{\d\Phi }{\d v}\Vert d(u,v)
$\egroup

die Einheitsnormale an die Fläche und das Oberflächen- oder 2-dimensionale Volumselement bezeichnet. Dies beiden Objekte sind wohldefiniert, falls \bgroup\color{proclaim}$ \{\d _1\Phi (u,v),\d _2\Phi (u,v)\}$\egroup linear unabhängig ist, d.h. \bgroup\color{proclaim}$ \Phi '(u,v)$\egroup (maximalen) Rang 2 hat. Beachte, daß für \bgroup\color{proclaim}$ C^1$\egroup-Parametrisierungen \bgroup\color{proclaim}$ \Phi :K\to S$\egroup und stetiges \bgroup\color{proclaim}$ \rho :S\to\mathbb{R}$\egroup folgendes Integral

\bgroup\color{proclaim}$\displaystyle \int_\Phi \rho \o\Phi  \operatorname{vol}...
...Vert\frac{\d\Phi }{\d u}\times \frac{\d\Phi }{\d v}\Bigr\Vert d(u,v)}.
$\egroup

reparametrisierungsinvariant ist, denn wenn \bgroup\color{proclaim}$ g:(r,s)\mapsto(u,v)$\egroup eine (nicht notwendigerweise Orientierungs-erhaltende) Substitutionsfunktion ist, dann ist

$\displaystyle \frac{\d\Phi }{\d r}\times \frac{\d\Phi }{\d s}$ $\displaystyle = \Bigl(\frac{\d\Phi }{\d u} \frac{\d u}{\d r} + \frac{\d\Phi }{...
...Phi }{\d u} \frac{\d u}{\d s} + \frac{\d\Phi }{\d v} \frac{\d v}{\d s} \Bigr)$    
  $\displaystyle = \frac{\d u}{\d r}  \frac{\d u}{\d s}   \undersetbrace{=0}\to{...
...u}{\d r}  \frac{\d v}{\d s}   \frac{\d\Phi }{\d u}\times \frac{\d\Phi }{\d v}$    
  $\displaystyle \quad+ \frac{\d v}{\d r}  \frac{\d u}{\d s}   \undersetbrace{- ...
...\d s}   \undersetbrace{=0}\to{\frac{\d\Phi }{\d v}\times \frac{\d\Phi }{\d v}}$    
  $\displaystyle = \det\left(\frac{\d (u,v)}{\d (r,s)}\right)\cdot \left(\frac{\d\Phi }{\d u}\times \frac{\d\Phi }{\d v}\right).$    

Es macht also Sinn folgende Definition zu geben:

\bgroup\color{proclaim}$\displaystyle \int_S \rho \operatorname{vol}_S := \int_\Phi \rho \o\Phi  \operatorname{vol}_\Phi
$\egroup

und dies kann als Masse von \bgroup\color{proclaim}$ S$\egroup bzgl. der Dichtefunktion \bgroup\color{proclaim}$ \rho $\egroup aufgefaßt werden, oder im Spezialfall \bgroup\color{proclaim}$ \rho =1$\egroup als Oberfläche von \bgroup\color{proclaim}$ S$\egroup. Dabei ist wieder \bgroup\color{proclaim}$ \operatorname{vol}_S$\egroup nur ein Symbol (das Oberflächenelement von \bgroup\color{proclaim}$ S$\egroup heißt) und dazu dient, die definierende rechte Seite zu imitieren. In der Tat sind die Riemann-Summen von \bgroup\color{proclaim}$ \int_K \rho \o\Phi  \operatorname{vol}_\Phi $\egroup gerade

$\displaystyle R(\rho \o\Phi  \operatorname{vol}_\Phi )$ $\displaystyle = \sum_{I,J} \rho (\Phi (\xi _I,\eta _J))\cdot \undersetbrace{=\t...
...Phi \times \d _2\Phi )(\xi _I,\eta _J) \Bigr\Vert  \vert I\vert \vert J\vert}$    
  $\displaystyle \approx \sum_{I,J} \rho (\Phi (\xi _I,\eta _J))\cdot ($``Fläche'' des Bildes von $ I\times J$ unter $ \Phi $$\displaystyle )$    
  $\displaystyle =$   ``Fläche'' des Bildes von $ K$ unter $ \Phi $$\displaystyle .
$    

Im Fall, daß \bgroup\color{proclaim}$ \Phi $\egroup injektiv ist (und natürlich noch immer der Rang von \bgroup\color{proclaim}$ \Phi '$\egroup gleich 2 ist) können wir die Einheitsnormale an die Fläche \bgroup\color{proclaim}$ S$\egroup als

\bgroup\color{proclaim}$\displaystyle \nu_S(\Phi (u,v)):=\nu_{\Phi }(u,v):= \fra...
...hi }{\d v}}
{\Vert\frac{\d\Phi }{\d u}\times \frac{\d\Phi }{\d v}\Vert}
$\egroup

definieren. Diese steht normal auf \bgroup\color{proclaim}$ \d _1\Phi $\egroup und \bgroup\color{proclaim}$ \d _2\Phi $\egroup und somit an jede Richtungsableitung (d.h. Tangentialvektor an die Fläche \bgroup\color{proclaim}$ S$\egroup) \bgroup\color{proclaim}$ d_w\Phi $\egroup für \bgroup\color{proclaim}$ w\in\mathbb{R}^2$\egroup. Der Satz von Stokes hat nun folgende geometrische Form:

\bgroup\color{proclaim}$\displaystyle \int_{\d S} \langle f, \tau _{\d S}\rangle...
...\int_S \langle \operatorname{rot}f,\nu_S\rangle  \operatorname{vol}_S.
$\egroup

Der Anteil der Rotation \bgroup\color{proclaim}$ \operatorname{rot}f$\egroup in Richtung eines Vektors \bgroup\color{proclaim}$ \nu$\egroup mißt also wieder die Wirbeldichte von \bgroup\color{proclaim}$ f$\egroup in der Ebene \bgroup\color{proclaim}$ S$\egroup die normal auf \bgroup\color{proclaim}$ \nu$\egroup steht.


Beispiel.
Die Sphäre wird durch \bgroup\color{proclaim}$ \Phi :(\varphi ,\th )\mapsto (\cos(\varphi )\cos(\th ),\sin(\varphi )\cos(\th ),\sin(\th ))$\egroup mit \bgroup\color{proclaim}$ K:=[0,2\pi]\times [-\pi/2,\pi/2]$\egroup parametrisiert. Somit ist die Oberfläche der Sphäre durch

\bgroup\color{proclaim}$\displaystyle \int_K \left\Vert\frac{\d\Phi }{\d\varphi }\times \frac{\d\Phi }{\d\th }\right\Vert d(\varphi ,\th )
$\egroup

gegeben. Wir erhalten:

$\displaystyle \Phi '(\varphi ,\th )$ $\displaystyle = \left(\begin{matrix}-\cos(\th )\sin(\varphi ) & -\sin(\th )\cos...
...cos(\varphi ) & -\sin(\th )\sin(\varphi )  0 & \cos(\th ) \end{matrix}\right)$    
$\displaystyle \Rightarrow \frac{\d\Phi }{\d\varphi }\times \frac{\d\Phi }{\d\th }$ $\displaystyle = \left(\begin{matrix}\cos(\th )^2\cos(\varphi )  \cos(\th )^2\...
...\sin(\th )\cos(\th ) \end{matrix}\right) = \cos(\th )\cdot \Phi (\varphi ,\th )$    

Somit ist die Oberfläche

$\displaystyle \int_K \left\Vert\frac{\d\Phi }{\d\varphi }\times \frac{\d\Phi }{\d\th }\right\Vert d(\varphi ,\th )$ $\displaystyle = \int_K \vert\cos(\th )\Vert \oversetbrace 1\to {\Vert\Phi \Vert} d(\varphi ,\th )$    
  $\displaystyle = \int_0^{2\pi}\int_{-\pi/2}^{\pi/2}\cos(\th ) d\th d\varphi = 2\pi\cdot [\sin(\th )]_{\th =-\pi/2}^{\pi/2} = 4\pi.$    


20.6.6 Definition. Normalbereich.
Unter einem $ C^1$-Normalbereich bezüglich der \bgroup\color{proclaim}$ xy$\egroup-Ebene verstehen wir eine Ordinatenmenge

\bgroup\color{proclaim}$\displaystyle B=\{(x,y,z):(x,y)\in K,\varphi _1(x,y)\leq z\leq\varphi _2(x,y)\}
$\egroup

wobei \bgroup\color{proclaim}$ K\subseteq \mathbb{R}^2$\egroup kompakt ( \bgroup\color{proclaim}$ J$\egroup-meßbar \bgroup\color{proclaim}$ \Leftarrow$\egroup) \bgroup\color{proclaim}$ \d K$\egroup \bgroup\color{proclaim}$ C^1$\egroup-parametrisierbar und \bgroup\color{proclaim}$ \varphi _1\leq\varphi _2:K\to\mathbb{R}$\egroup stetig sind. Die Deckel sind dann durch \bgroup\color{proclaim}$ S_i:=\operatorname{graph}(\varphi _i)$\egroup für \bgroup\color{proclaim}$ i=1,2$\egroup gegeben und der Mantel \bgroup\color{proclaim}$ S_0:=\{(x,y,z):(x,y)\in\d K, \varphi _1(x,y)\leq z\leq\varphi _2(x,y)\}$\egroup. Wir benötigen eine \bgroup\color{proclaim}$ C^1$\egroup-Parametrisierung der Oberfläche \bgroup\color{proclaim}$ \d B=S_0\cup S_1\cup S_2$\egroup. Da in wichtigen Beispielen (Kugel) die \bgroup\color{proclaim}$ \varphi _i$\egroup aber nicht überall differenzierbar sind, fordern wir zusätzlich die Existenz von Substitutionsfunktionen \bgroup\color{proclaim}$ g_i:K_i\to K$\egroup, s.d. \bgroup\color{proclaim}$ \varphi _i\o g$\egroup \bgroup\color{proclaim}$ C^1$\egroup ist. Dann ist \bgroup\color{proclaim}$ \Phi _i:K_i\to S_i$\egroup eine Parametrisierung von \bgroup\color{proclaim}$ S_i$\egroup für \bgroup\color{proclaim}$ i\in\{0,1,2\}$\egroup, wobei \bgroup\color{proclaim}$ K_0:=\{(u,v):u\in[a,b],\varphi _1(c(u))\leq v\leq\varphi _2(c(u))\}$\egroup, weiters \bgroup\color{proclaim}$ c:[a,b]\to\d K$\egroup eine positiv orientierte \bgroup\color{proclaim}$ C^1$\egroup-Parametrisierung von \bgroup\color{proclaim}$ \d K$\egroup ist, und

\bgroup\color{proclaim}$\displaystyle \Phi _i(u,v) := \left\{\begin{array}{ll} (...
...=0  (g_i(u,v),\varphi _i(g_i(u,v))) &\text{sonst.}
\end{array}\right.
$\egroup

Wir benötigen noch, daß die Orientierung dieser Parametrisierungen paßt, d.h. der Normalvektor \bgroup\color{proclaim}$ \nu_\Phi $\egroup nach außen weißt. Dafür folgende Rechnung: Es ist

$\displaystyle \Phi _i'(u,v)$ $\displaystyle = \left\{\begin{array}{ll} \left(\begin{matrix}x'(u) & 0  y'(u)...
... (\varphi _i\o g_i)'(u,v) \end{matrix}\right) &\text{sonst,}\end{array}\right.$    
$\displaystyle \Rightarrow \nu_{\Phi _i}(u,v)$ $\displaystyle = \left\{\begin{array}{ll} \left(\begin{matrix}y'(u)  -x'(u) \\...
...}*  *  \det(g_i'(u,v)) \end{matrix}\right) &\text{sonst.}\end{array}\right.$    

Also verlangen wir zusätzlich o.B.d.A. , daß \bgroup\color{proclaim}$ \det(g_1'(u,v))\leq 0$\egroup und \bgroup\color{proclaim}$ \det(g_2'(u,v))\geq 0$\egroup für alle \bgroup\color{proclaim}$ (u,v)$\egroup gilt.


Beispiel.
Der Quader \bgroup\color{proclaim}$ [a_1,a_2]\times [b_1,b_2]\times [c_1,c_2]$\egroup: Dabei ist \bgroup\color{proclaim}$ K=[a_1,a_2]\times [b_1,b_2]$\egroup und \bgroup\color{proclaim}$ \varphi _i:=c_i$\egroup. Hier müssen wir nicht mehr umparametrisieren, d.h. \bgroup\color{proclaim}$ K_i:=K$\egroup und \bgroup\color{proclaim}$ g_i:=\operatorname{id}$\egroup.

Die Kugel \bgroup\color{proclaim}$ \{(x,y,z):x^2+y^2+z^2\leq R^2\}$\egroup mit Radius \bgroup\color{proclaim}$ R$\egroup: Dabei ist \bgroup\color{proclaim}$ K:=\{(x,y):x^2+y^2\leq \mathbb{R}^2\}$\egroup, \bgroup\color{proclaim}$ \varphi _2(x,y):=\sqrt{R^2-(x^2+y^2)}$\egroup, \bgroup\color{proclaim}$ \varphi _1(x,y):=-\sqrt{R^2-(x^2+y^2)}$\egroup. Die \bgroup\color{proclaim}$ \varphi _i$\egroup sind nun auf \bgroup\color{proclaim}$ \d K$\egroup nicht mehr differenzierbar. Wir verwenden Kugelkoordinaten zur Umparametrisierung, d.h. \bgroup\color{proclaim}$ K_1:=[0,2\pi]\times [-\pi/2,0]$\egroup, \bgroup\color{proclaim}$ K_2:=[0,2\pi]\times [0,\pi/2]$\egroup und \bgroup\color{proclaim}$ g_i(u,v):=R\cos(v)(\cos(u),\sin(u))$\egroup. Dann ist \bgroup\color{proclaim}$ (\varphi _i\o g_i)(u,v)=\pm\sqrt{R^2(1-\cos^2 v)}=R\sin(v)$\egroup \bgroup\color{proclaim}$ C^{\infty}$\egroup in \bgroup\color{proclaim}$ (u,v)$\egroup.


20.6.7 Der Gauß'sche Integralsatz im Raum.
Für stetiges \bgroup\color{proclaim}$ r:B\to\mathbb{R}$\egroup ist

\bgroup\color{proclaim}$\displaystyle \int_{\d B} r dx\wedge dy = \sum_{i=1}^3 \int_{S_i} r(x,y,z) dx\wedge dy
$\egroup

Weiters ist

$\displaystyle \int_S r(x,y,z) dx\wedge dy$ $\displaystyle = \int_K r(\Phi (u,v)) \det\frac{\d (x,y)}{\d (u,v)} d(u,v).$    

Da

\bgroup\color{proclaim}$\displaystyle \frac{\d (x,y)}{\d (u,v)} = \left\{\begin{...
...0) &\text{für }i=0 \\
(g_i'(u,v)) &\text{für }i=1,2 \end{array}\right.
$\egroup

ist das Integral 0 für \bgroup\color{proclaim}$ i=0$\egroup und für \bgroup\color{proclaim}$ i=1,2$\egroup ist

$\displaystyle \int_{S_i} r(x,y,z)  dx\wedge dy$ $\displaystyle = \int_{K_i} r(g_i(u,v),\varphi _i(g_i(u,v))) \det g_i'(u,v) d(u,v)$    
  $\displaystyle = \mp \int_K r(x,y,\varphi _i(x,y)) d(x,y)$    

Folglich ist


$\displaystyle \int_{\d B} r dx\wedge dy$ $\displaystyle = \int_K [r(x,y,z)]_{z=\varphi _1(x,y)}^{\varphi _2(x,y)} d(x,y)$    
  $\displaystyle \overset{\text{HS}}{=} \int_K\int_{\varphi _1(x,y)}^{\varphi _2(x,y)} \frac{\d }{\d z} r(x,y,z) dz d(x,y)$    
  $\displaystyle = \int_B \frac{\d }{\d z} r(x,y,z) d(x,y,z)$    

Falls $ B$ ein $ C^1$-Normalbereich auch bezüglich der anderen Koordinatenebenen ist, so folgt analog


   
$\displaystyle \int_B \left( \frac{\d }{\d x}p +\frac{\d }{\d y}q+\frac{\d }{\d z}r\right)$ $\displaystyle =\int_{\d B} p\;dy\wedge dz + q\;dz\wedge dx + r\;dx\wedge dy$    




Satz von Gauß im 3-dimensionalen.
Es sei \bgroup\color{proclaim}$ B$\egroup ein \bgroup\color{proclaim}$ C^1$\egroup-Normalbereich bezüglich aller 3 Koordinatenebenen und \bgroup\color{proclaim}$ f=(p,q,r)$\egroup ein \bgroup\color{proclaim}$ C^1$\egroup-Vektorfeld auf \bgroup\color{proclaim}$ B$\egroup. Dann ist

\bgroup\color{proclaim}$\displaystyle \int_B \left( \frac{\d }{\d x}p +\frac{\d ...
...nt_{\d B} p dy\wedge dz + q dz\wedge dx + r dx\wedge dy {\rm\quad[]}
$\egroup


Geometrische Bedeutung.

Für ein \bgroup\color{proclaim}$ C^1$\egroup-Vektorfeld \bgroup\color{proclaim}$ f=(p,q,r):\mathbb{R}^3\supseteq B\to\mathbb{R}^3$\egroup definieren wir eine Funktion \bgroup\color{proclaim}$ \operatorname{div}f:B\to \mathbb{R}$\egroup durch \bgroup\color{proclaim}$ \operatorname{div}f:=\d _1 p + \d _2 q + \d _3 r$\egroup. Dann ist

\bgroup\color{proclaim}$\displaystyle \int_B \operatorname{div}f = \int_{\d B} \langle f,\nu_{\d B}\rangle  \operatorname{vol}_{\d B}.
$\egroup

Die rechte Seite beschreibt wie im 2-dimensionalen die Quellenergiebigkeit oder auch Quellenstärke des Bereichs \bgroup\color{proclaim}$ B$\egroup, und somit nennt man \bgroup\color{proclaim}$ \operatorname{div}f$\egroup die Quellendichte des Vektorfelds \bgroup\color{proclaim}$ f$\egroup.

glossary

index

Andreas Kriegl 2002-07-01