Kurs:Mathematik für Anwender (Osnabrück 2019-2020)/Teil II/Vorlesung 60/kontrolle

Die Transformationsformel für Integrale

Wir kommen zur Transformationsformel für Integrale, wofür wir noch eine Bezeichnung einführen.

Definition Referenznummer erstellen

Es sei ${}G\subseteq V$ eine offene Teilmenge in einem endlichdimensionalen reellen Vektorraum und sei

$\varphi \colon G\longrightarrow V$

eine total differenzierbare Abbildung. Dann nennt man die Determinante

$\det \left(D\varphi \right)_{P}$

die Jacobi-Determinante (oder Fundamental-Determinante) von ${}\varphi$ in ${}P\in G$.

Wir betrachten die Jacobi-Determinante als eine auf ${}G$ definierte reellwertige Funktion ${}P\mapsto J(\varphi )(P)=J(P):=\det \left(D\varphi \right)_{P}$. Bei einer stetig partiell differenzierbaren Abbildung ist sie stetig, da dann die Einträge in der Jacobi-Matrix stetige Funktionen sind. Bei einer linearen Abbildung

$\varphi \colon \mathbb {R} ^{n}\longrightarrow \mathbb {R} ^{n}$

stimmt das totale Differential in jedem Punkt mit ${}\varphi$ selbst überein, und daher ist die Jacobi-Determinante konstant. In Satz 58.9 haben wir gesehen, dass die Determinante der linearen Abbildung das Verhältnis zwischen dem Volumen von Bild und Urbild festlegt. Eine wesentliche Verallgemeinerung von dieser Beziehung wird durch die beiden folgenden Aussagen gegeben.

Satz

Es seien ${}G$ und ${}H$ offene Mengen im ${}\mathbb {R} ^{n}$ und es sei

$\varphi \colon G\longrightarrow H$

ein ${}C^{1}$-Diffeomorphismus mit der Jacobi-Determinante ${}(J(\varphi ))(x)=\det \left(D\varphi \right)_{x}$ für ${}x\in G$. Es sei ${}T\subseteq H$ eine (in ${}\mathbb {R} ^{n}$) kompakte Teilmenge und es sei

$f\colon T\longrightarrow \mathbb {R}$

eine stetige Funktion.

Dann ist ${}\varphi ^{-1}(T)$ ebenfalls kompakt und es gilt

${}\int _{T}f\,d\lambda ^{n}=\int _{\varphi ^{-1}(T)}(f\circ \varphi )\vert {J(\varphi )}\vert \,d\lambda ^{n}\,.$

Beweis

Dieser Beweis wurde in der Vorlesung nicht vorgeführt.

$\Box$

Häufig startet man auch mit einer kompakten Teilmenge ${}S\subseteq G$ und setzt ${}T=\varphi (S)$. Es kann auch auf ${}S$ eine stetige Funktion ${}{\tilde {f}}$ definiert sein, dann muss man ${}f={\tilde {f}}\circ \varphi ^{-1}$ setzen. Da ein Diffeomorphismus vorausgesetzt wird, ist die Aussage dieses Satzes grundsätzlich symmetrisch.

Korollar Korollar 60.3 ändern

Es seien ${}G$ und ${}H$ offene Mengen im ${}\mathbb {R} ^{n}$ und es sei

$\varphi \colon G\longrightarrow H$

ein ${}C^{1}$-Diffeomorphismus mit der Jacobi-Determinante

${}(J(\varphi ))(x)=\det \left(D\varphi \right)_{x}\,$

für ${}x\in G$. Es sei ${}T\subseteq H$ eine (in ${}\mathbb {R} ^{n}$) kompakte Teilmenge.

Dann gilt

${}\lambda ^{n}(T)=\int _{\varphi ^{-1}(T)}\vert {J(\varphi )}\vert \,d\lambda ^{n}\,.$

Beweis

Dies folgt aus Satz 57.2, angewendet auf die konstante Funktion ${}f=1$.

$\Box$

Beispiele zur Transformationsformel

Beispiel Referenznummer erstellen

Wir betrachten das komplexe Quadrieren

${\mathbb {C} }\longrightarrow {\mathbb {C} },\,z\longmapsto z^{2}.$

In reellen Koordinaten ist dies die differenzierbare Abbildung

$\varphi \colon \mathbb {R} ^{2}\longrightarrow \mathbb {R} ^{2},\,(x,y)\longmapsto (x^{2}-y^{2},2xy).$

Diese Abbildung ist wegen ${}\varphi (x,y)=\varphi (-x,-y)$ nicht injektiv. Allerdings ist die Einschränkung auf die positive Halbebene ${}G={\left\{(x,y)\mid x>0\right\}}$ injektiv, und das Bild davon ist ${}H=\mathbb {R} ^{2}\setminus \mathbb {R} _{-}$ (also die Ebene ohne die negative reelle Achse). Die Jacobi-Matrix von ${}\varphi$ ist

${}\operatorname {Jak} (\varphi )_{(x,y)}={\begin{pmatrix}2x&-2y\\2y&2x\end{pmatrix}}\,$

mit der Jacobi-Determinante

${}(J(\varphi ))(x,y)=4x^{2}+4y^{2}\,.$

Wir möchten den Flächeninhalt des Bildes ${}T=\varphi (S)$ des Einheitsquadrates ${}S=[0,1]\times [0,1]$ unter dieser Abbildung berechnen (die eine Seite des Einheitsquadrates gehört nicht zu ${}G$, dieser Rand ist aber eine Nullmenge nach Lemma 58.7 und daher für den Flächeninhalt und die Integration unerheblich). Aufgrund von Korollar 60.3 ist dann

${}{\begin{aligned}\lambda ^{2}(T)&=\int _{S}4x^{2}+4y^{2}d\lambda ^{2}\\&=\int _{0}^{1}\int _{0}^{1}{\left(4x^{2}+4y^{2}\right)}dxdy\\&=\int _{0}^{1}{\left({\frac {4}{3}}x^{3}+4xy^{2}\right)}{|}_{0}^{1}dy\\&=\int _{0}^{1}{\left({\frac {4}{3}}+4y^{2}\right)}dy\\&={\left({\frac {4}{3}}y+{\frac {4}{3}}y^{3}\right)}{|}_{0}^{1}\\&={\frac {8}{3}}.\end{aligned}}$

Korollar Korollar 60.5 ändern

Es sei

$\varphi \colon \mathbb {R} ^{2}\longrightarrow \mathbb {R} ^{2},\,(r,\theta )\longmapsto (r\cos \theta ,r\sin \theta ),$

die Polarkoordinatenauswertung und es seien ${}G$ und ${}H$ offene Mengen, auf denen ${}\varphi$ einen Diffeomorphismus induziert. Es sei ${}T\subseteq H$ eine (in ${}\mathbb {R} ^{2}$) kompakte Teilmenge und

$f\colon T\longrightarrow \mathbb {R}$

eine stetige Funktion.

Dann ist

$\int _{T}f(x,y)\,d\lambda ^{2}(x,y)=\int _{\varphi ^{-1}(T)}f(r\cos \theta ,r\sin \theta )\cdot r\,d\lambda ^{2}(r,\theta )\,.$

Dies gilt auch dann, wenn außerhalb von Nullmengen ein Diffeomorphismus vorliegt. Insbesondere gilt bei stetigem ${}f$ die Formel

$\int _{\mathbb {R} ^{2}}f(x,y)\,d\lambda ^{2}(x,y)=\int _{0}^{\infty }\int _{0}^{2\pi }f(r\cos \theta ,r\sin \theta )\cdot r\,d\theta \,dr\,.$

Beweis

Dies folgt wegen

${}\det {\left(D\varphi \right)}=\det {\begin{pmatrix}\cos \theta &-r\sin \theta \\\sin \theta &r\cos \theta \end{pmatrix}}=r\cos ^{2}\theta +r\sin ^{2}\theta =r\,$

direkt aus Satz 57.2.

$\Box$

Lemma Referenznummer erstellen

Es ist

${}\int _{-\infty }^{+\infty }e^{-x^{2}}\,dx={\sqrt {\pi }}\,.$

Beweis

Durch eine einfache Substitution ist die Aussage äquivalent zu

${}\int _{-\infty }^{+\infty }{\frac {1}{\sqrt {2\pi }}}e^{-{\frac {t^{2}}{2}}}\,dt=1\,.$

Nennen wir dieses Integral ${}I$. Nach einer Variante des Satzes von Fubini für uneigentliche Integrale ist

${}I^{2}={\left(\int _{-\infty }^{+\infty }{\frac {1}{\sqrt {2\pi }}}e^{-{\frac {t^{2}}{2}}}\,dt\right)}\cdot {\left(\int _{-\infty }^{+\infty }{\frac {1}{\sqrt {2\pi }}}e^{-{\frac {t^{2}}{2}}}\,dt\right)}=\int _{\mathbb {R} ^{2}}{\frac {1}{2\pi }}e^{-{\frac {x^{2}+y^{2}}{2}}}\,d\lambda ^{2}\,.$

Durch Einführung von Polarkoordinaten ${}x=r\cos \theta$ und ${}y=r\sin \theta$ ist dieses Integral nach Korollar 60.5 und nach Fubini gleich

${}{\begin{aligned}\,&={\frac {1}{2\pi }}\int _{[0,2\pi ]\times \mathbb {R} _{\geq 0}}e^{-{\frac {r^{2}}{2}}}\cdot r\,d\lambda ^{2}(r,\theta )\\&={\frac {1}{2\pi }}{\left(\int _{[0,2\pi ]}1\,d\lambda ^{1}(\theta )\right)}{\left(\int _{\mathbb {R} _{\geq 0}}e^{-{\frac {r^{2}}{2}}}\cdot r\,d\lambda ^{1}(r)\right)}\\&=\int _{\mathbb {R} _{\geq 0}}e^{-{\frac {r^{2}}{2}}}\cdot r\,d\lambda ^{1}(r)\\&=\int _{0}^{\infty }e^{-{\frac {r^{2}}{2}}}\cdot r\,dr\\&=-e^{-{\frac {r^{2}}{2}}}|_{0}^{\infty }\\&=1.\end{aligned}}$

Damit ist auch ${}I=1$.

$\Box$

Korollar Referenznummer erstellen

Für die Zylinderkoordinatenauswertung

$\Psi \colon G=\mathbb {R} _{+}\times ]0,2\pi [\times \mathbb {R} \longrightarrow \mathbb {R} ^{3},\,(r,\theta ,z)\longmapsto \left(r\cos \theta ,\,r\sin \theta ,\,z\right),$

eine kompakte Teilmenge ${}T\subseteq \Psi (G)$ und eine stetige Funktion

$f\colon T\longrightarrow \mathbb {R}$

gilt die Beziehung

${}\int _{T}f\,d\lambda ^{3}=\int _{\Psi ^{-1}(T)}(f\circ \Psi )\cdot r\,d\lambda ^{3}\,.$

Dies kann man auch als

$\int \int \int _{T}f(x,y,z)dxdydz=\int \int \int _{\Psi ^{-1}(T)}{\tilde {f}}(r,\theta ,z)\cdot rdrd\theta dz\,$

schreiben, wobei ${}{\tilde {f}}=f\circ \Psi$ bezeichnet.

Beweis

Dies folgt aus Satz 57.2, da die Jacobi-Determinante der Zylinderkoordinatenauswertung gleich ${}r$ ist.

$\Box$

Korollar Referenznummer erstellen

Für die Kugelkoordinatenauswertung

$\Psi \colon G=\mathbb {R} _{+}\times ]0,\pi [\times ]0,2\pi [\longrightarrow \mathbb {R} ^{3},\,(r,\theta ,\varphi )\longmapsto \left(r\cos \varphi \sin \theta ,\,r\sin \varphi \sin \theta ,\,r\cos \theta \right),$

eine kompakte Teilmenge ${}T\subseteq \Psi (G)$ und eine stetige Funktion

$f\colon T\longrightarrow \mathbb {R}$

gilt die Beziehung

${}\int _{T}f\,d\lambda ^{3}=\int _{\Psi ^{-1}(T)}(f\circ \Psi )r^{2}\sin \theta \,d\lambda ^{3}\,.$

Dies kann man auch als

$\int \int \int _{T}f(x,y,z)dxdydz=\int \int \int _{\Psi ^{-1}(T)}{\tilde {f}}(r,\theta ,\varphi )\cdot r^{2}\sin \theta drd\theta d\varphi \,$

schreiben, wobei ${}{\tilde {f}}=f\circ \Psi$ bezeichnet.

Beweis

Nach Beispiel 53.7 ist die Jacobi-Determinante von ${}\Psi$ im Punkt ${}\left(r,\,\theta ,\,\varphi \right)$ gleich ${}r^{2}\sin \theta$, sodass die Aussage aus Satz 57.2 folgt.

$\Box$

Beispiel Referenznummer erstellen

Es soll eine Straße in der Ebene der Breite ${}2a$ asphaltiert werden. Dabei wird die Straße durch den Verlauf des Mittelstreifen vorgegeben, der durch die Kurve

$[0,s]\longrightarrow \mathbb {R} ^{2},\,t\longmapsto \psi (t)={\begin{pmatrix}f(t)\\g(t)\end{pmatrix}},$

bestimmt ist. Dabei sei ${}\psi$ zweimal stetig differenzierbar und bogenparametrisiert, d.h. es sei ${}f'(t)^{2}+g'(t)^{2}=1$, was bedeutet, dass die Mittelstreifenkurve mit normierter Geschwindigkeit durchlaufen wird. Die Breite ist dabei senkrecht zum Mittelstreifen zu messen. Die zu asphaltierende Trasse wird dann durch die Abbildung

$\varphi \colon [0,s]\times [-a,a]\longrightarrow \mathbb {R} ^{2},\,(t,r)\longmapsto {\begin{pmatrix}f(t)\\g(t)\end{pmatrix}}+r{\begin{pmatrix}-g'(t)\\f'(t)\end{pmatrix}},$

parametrisiert. Wir nehmen an, dass diese Parametrisierung injektiv ist, was erfüllt ist, wenn die Mittelstreifenabbildung ${}\psi$ injektiv ist und die Straße nicht zu breit werden soll.

Die Jacobi-Matrix der Parametrisierung ist

${}\left(D\varphi \right)_{(t,r)}={\begin{pmatrix}f^{\prime }(t)-rg^{\prime \prime }(t)&-g'(t)\\g^{\prime }(t)+rf^{\prime \prime }(t)&f'(t)\end{pmatrix}}\,.$

Die Determinante davon ist

${}{\begin{aligned}f'(t)f'(t)+g'(t)g'(t)-r{\left(g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)\right)}&=1-r{\left(g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)\right)}\\\end{aligned}}$

Daher ist die Asphaltfläche nach der Transformationsformel gleich

$\int _{[0,s]\times [-a,a]}\vert {1-r{\left(g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)\right)}}\vert \,d\lambda ^{2}.$

Wenn wir weiter annehmen, dass

${}\vert {g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)}\vert \leq {\frac {1}{a}}\,$

ist (was bedeutet, dass die Straßenbreite nicht allzu groß ist), so ist dieses Integral nach Korollar 59.6 geich

${}{\begin{aligned}\int _{[0,s]\times [-a,a]}1-r{\left(g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)\right)}\,d\lambda ^{2}&=2as-{\left(\int _{-a}^{a}r\,dr\right)}{\left(\int _{0}^{s}g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)\,dt\right)}\\&=2as-0\cdot {\left(\int _{0}^{s}g^{\prime \prime }(t)f'(t)-f^{\prime \prime }(t)g'(t)\,dt\right)}\\&=2as.\,\end{aligned}}$

Dies bedeutet, dass die Asphaltfläche gleich der Mittelstreifenlänge mal der Straßenbreite ist.

Volumentreue Abbildungen

In einer früheren Vorlesung haben wir über Isometrien gesprochen, also lineare Abbildungen zwischen euklidischen Vektorräumen, die das Skalarprodukt und insbesondere die Norm respektieren. Aufgrund von Aufgabe 60.7 ist die Determinante einer Isometrie auf dem ${}\mathbb {R} ^{n}$ gleich ${}1$ oder ${}-1$, so dass sich gemäß Korollar 60.3 das Volumen von beliebigen kompakten Teilmengen ${}T\subseteq \mathbb {R} ^{n}$ unter der Abbildung nicht ändert, d.h. es ist stets ${}\lambda ^{n}(\varphi (T))=\lambda ^{n}(T)$.

Definition Referenznummer erstellen

Es seien ${}G$ und ${}H$ offene Mengen im ${}\mathbb {R} ^{n}$ und es sei

$\varphi \colon G\longrightarrow H$

ein ${}C^{1}$-Diffeomorphismus. Man sagt, dass ${}\varphi$ volumentreu ist, wenn

${}\vert {(J(\varphi ))(x)}\vert =\vert {\det \left(D\varphi \right)_{x}}\vert =1\,$

für alle ${}x\in G$ ist.

Im ebenen Fall spricht man natürlich von flächentreu. Für einen volumentreuen Diffeomorphismus ist

${}\lambda ^{n}(\varphi (T))=\lambda ^{n}(T)\,$

nach Korollar 60.3.

Beispiel Referenznummer erstellen

Es sei ${}h\in \mathbb {R} [y]$ ein beliebiges Polynom in der einen Variablen ${}y$. Dann ist die Abbildung

$\varphi \colon \mathbb {R} ^{2}\longrightarrow \mathbb {R} ^{2},\,(x,y)\longmapsto (x+h(y),y)$

ein flächentreuer Diffeomorphismus. Die Jacobi-Matrix von ${}\varphi$ ist ja

${}\operatorname {Jak} (\varphi )_{(x,y)}={\begin{pmatrix}1&h'(y)\\0&1\end{pmatrix}}\,,$

sodass die Jacobi-Determinante konstant gleich ${}1$ ist. Wenn man die Rollen von ${}x$ und ${}y$ vertauscht und die Hintereinanderschaltung von solchen Abbildungen betrachtet, so erhält man flächentreue Abbildungen, denen man es nicht auf den ersten Blick ansieht. Beispielsweise ist zu ${}\varphi (x,y)=(x+y^{2},y)$ und ${}\psi (x,y)=(x,y+x^{3})$ die Hintereinanderschaltung

${}{\begin{aligned}(\psi \circ \varphi )(x,y)&=\psi (\varphi (x,y))\\&=\psi \left(x+y^{2},\,y\right)\\&=\left(x+y^{2},\,y+{\left(x+y^{2}\right)}^{3}\right)\\&=\left(x+y^{2},\,y+x^{3}+3x^{2}y^{2}+3xy^{4}+y^{6}\right).\end{aligned}}$