|
Autor |
Ableitung der Inversion |
|
Phoensie
Aktiv  Dabei seit: 11.04.2020 Mitteilungen: 442
Wohnort: Muri AG, Schweiz
 | Themenstart: 2020-10-02
|
Aufgabe:
Sei $\sigma: \mathbb{R}^n\setminus \{a\} \to \mathbb{R}^n\setminus \{a\}$ die Inversion bezüglich der Sphäre $S_r(a)=\{x \in \mathbb{R}^n \mid \|x-a\| = r\}$, definiert durch
\[ \begin{align*}
\sigma(x)
= a + \frac{r^2}{\|x-a\|^2}(x-a).
\end{align*}
\]
Berechne die Ableitung $\mathrm{D}\sigma(x): \mathbb{R}^n \to \mathbb{R}^n$ und die Jacobimatrix $J\sigma(x)$. Zeige, dass $\sigma$ konform ist.
Ich habe folgendes gemacht: Seien $x,y \in \mathbb{R}^n$. Dann gilt:
\[ \begin{align*}
\mathrm{D}\sigma(x)y
&= \frac{\mathrm{d}}{\mathrm{d}t}\Big|_{t=0} \sigma(x + ty) \\
&\quad\color{red}{\downarrow\;\text{Definition der Richtungsableitung.}}\\
&= \lim_{t \to 0} \frac{\sigma(x+ty) - \sigma(x)}{t} \\
&\quad\color{red}{\downarrow\;\text{Definition von $\sigma$.}}\\
&= \lim_{t \to 0} \frac{\left(a + \frac{r^2}{\|x+ty-a\|^2}(x+ty-a)\right) - \left(a + \frac{r^2}{\|x-a\|^2}(x-a)\right)}{t} \\
&= r^2\lim_{t \to 0} \frac{\frac{1}{\|x+ty-a\|^2}(x+ty-a) - \frac{1}{\|x-a\|^2}(x-a)}{t} \\
&\quad\color{red}{\downarrow\;\text{Distributivgesetz (Bruch aufspalten).}}\\
&= r^2\lim_{t \to 0} \left(\frac{x}{t\|x+ty-a\|^2} + \frac{ty}{t\|x+ty-a\|^2} - \frac{a}{t\|x+ty-a\|^2} - \frac{x}{t\|x-a\|^2} + \frac{a}{t\|x-a\|^2}\right) \\
&= r^2\lim_{t \to 0} \left(\frac{x}{t\|x+ty-a\|^2} + \frac{y}{\|x+ty-a\|^2} - \frac{a}{t\|x+ty-a\|^2} - \frac{x}{t\|x-a\|^2} + \frac{a}{t\|x-a\|^2}\right) \\
&= r^2\lim_{t \to 0} \left(\frac{\|x-a\|^2 - \|x+ty-a\|^2}{t\|x+ty-a\|^2\|x-a\|^2}x + \frac{y}{\|x+ty-a\|^2} - \frac{\|x-a\|^2 - \|x+ty-a\|^2}{t\|x+ty-a\|^2\|x-a\|^2}a\right) \\
&\quad\color{red}{\downarrow\;\|x-a\|^2=\|x\|^2 - 2\langle x,a \rangle + \|a\|^2.}\\
&\quad\color{red}{\downarrow\;\|x+ty-a\|^2= \|x\|^2 + 2t\langle x,y \rangle + t^2\|y\|^2 - 2\langle x,a \rangle - 2t\langle y,a \rangle + \|a\|^2.}\\
&\quad\color{red}{\downarrow\; \|x-a\|^2 - \|x+ty-a\|^2 = 2t(\langle y,a \rangle - \langle x,y \rangle) - t^2\|y\|^2.}\\
&= r^2\lim_{t \to 0} \left(\frac{2t(\langle y,a \rangle - \langle x,y \rangle) - t^2\|y\|^2}{t\|x+ty-a\|^2\|x-a\|^2}x + \frac{y}{\|x+ty-a\|^2} - \frac{2t(\langle y,a \rangle - \langle x,y \rangle) - t^2\|y\|^2}{t\|x+ty-a\|^2\|x-a\|^2}a\right) \\
&\quad\color{red}{\downarrow\;\text{Kürze wo möglich.}} \\
&= r^2\lim_{t \to 0} \left(\frac{2(\langle y,a \rangle - \langle x,y \rangle) - t\|y\|^2}{\|x+ty-a\|^2\|x-a\|^2}x + \frac{y}{\|x+ty-a\|^2} - \frac{2(\langle y,a \rangle - \langle x,y \rangle) - t\|y\|^2}{\|x+ty-a\|^2\|x-a\|^2}a\right) \\
&\quad\color{red}{\downarrow\;\text{Berechne den Grenzwert.}} \\
&= r^2\left(\frac{2(\langle y,a \rangle - \langle x,y \rangle)}{\|x-a\|^4}x + \frac{y}{\|x-a\|^2} - \frac{2(\langle y,a \rangle - \langle x,y \rangle)}{\|x-a\|^4}a\right) \\
&\quad\color{red}{\downarrow\;\text{Distributivgesetz.}} \\
&= \frac{2r^2(\langle y,a \rangle - \langle x,y \rangle)}{\|x-a\|^4}(x-a) + \frac{r^2}{\|x-a\|^2}y \\
&\quad\color{red}{\downarrow\;\text{Das Skalarprodukt ist symmetrisch.}} \\
&= \frac{2r^2(\langle a,y \rangle - \langle x,y \rangle)}{\|x-a\|^4}(x-a) + \frac{r^2}{\|x-a\|^2}y \\
&\quad\color{red}{\downarrow\;\text{Das Skalarprodukt ist bilinear.}} \\
&= \frac{2r^2\langle a-x,y \rangle}{\|x-a\|^4}(x-a) + \frac{r^2}{\|x-a\|^2}y \\
\end{align*}
\]
Mein Problem jetzt ist, dass ich hier die "lineare Abbildung" nicht erkenne...
Habe ich was übersehen?🤔
|
Profil
|
Triceratops
Wenig Aktiv  Dabei seit: 28.04.2016 Mitteilungen: 6472
Wohnort: Berlin
 | Beitrag No.1, eingetragen 2020-10-02
|
Vorab: Ich habe deine Rechnung nicht überprüft.
Die lineare Abbildung ist $y \mapsto D \sigma (x) y = $(das was du berechnet hast).
Man kann sich die Rechnung einfacher machen, indem man o.B.d.A. $a=0$ annimmt. Ganz am Ende kann man alles "nach $a$ verschieben", um den allgemeinen Fall abzudecken (Kettenregel). Denn es gilt $\sigma_a = \tau_a \circ \sigma_0 \circ \tau_{-a}$ mit $\tau_a(x) := x+a$.
|
Profil
|
Phoensie
Aktiv  Dabei seit: 11.04.2020 Mitteilungen: 442
Wohnort: Muri AG, Schweiz
 | Beitrag No.2, vom Themenstarter, eingetragen 2020-10-05
|
Lieber Triceratops
Danke für deinen Tipp. Nun habe ich folgendes gemacht:
Definiere die Abbildungen $\sigma_z(x) = z + \frac{r^2}{\|x-z\|^2}(x-z)$ und $\tau_z(x)=x+z$ für $z \in \mathbb{R}^n$. Damit gilt
\[
\begin{align*}
\sigma(x)
= \tau_a\left(\frac{r^2}{\|\tau_{-a}(x)\|^2}\tau_{-a}(x)\right) = (\tau_a \circ \sigma_0 \circ \tau_{-a})(x)
\end{align*}
\]
und für die Ableitung folgt gemäss der Kettenregel:
\[
\begin{align*}
\mathrm{D}\sigma(x)
= \mathrm{D}\tau_a\big(\sigma_0 (\tau_{-a}(x))\big) \circ \mathrm{D}\sigma_0 \big(\tau_{-a}(x)\big) \circ \mathrm{D}\tau_{-a}(x).
\end{align*}
\]
Somit können wir zur Berechnung des Differentials von $\sigma$ auch schlicht die Differentiale der Abbildungen $\tau_a,\tau_{-a},\sigma_0$ berechnen. Mit leichter Rechnerei (die ich hier nicht aufführe) erhält man
\[
\begin{align*}
\mathrm{D}\tau_a\big(\sigma_0 (\tau_{-a}(x))\big)h = h
&\implies \mathrm{D}\tau_{a}(\sigma_0 (\tau_{-a}(x))) = \mathrm{id}. \\
\mathrm{D}\tau_{-a}(x)h
= h
&\implies \mathrm{D}\tau_{-a}(x) = \mathrm{id}.
\end{align*}
\]
Damit haben wir als Zwischenresultat
\[
\begin{align*}
\mathrm{D}\sigma(x)
&= \mathrm{D}\tau_a\big(\sigma_0 (\tau_{-a}(x))\big) \circ \mathrm{D}\sigma_0 \big(\tau_{-a}(x)\big) \circ \mathrm{D}\tau_{-a}(x) \\
&= \mathrm{id} \circ \mathrm{D}\sigma_0 \big(\tau_{-a}(x)\big) \circ \mathrm{id} \\
&= \mathrm{D}\sigma_0 \big(\tau_{-a}(x)\big) \\
\end{align*}
\]
Bis hierhin ist alles klar.👍😁 Jedoch habe ich bei der Berechnung von $\mathrm{D}\sigma_0(x)$ wieder so meine Mühe:
\[
\begin{align*}
\mathrm{D}\sigma_0(x)h
&= \frac{\mathrm{d}}{\mathrm{d}t}\Big|_{t=0} \sigma_0(x + th) \\
&\quad\color{red}{\downarrow\;\text{Definition der Richtungsableitung.}}\\
&= \lim_{t \to 0} \frac{\sigma_0(x+th) - \sigma_0(x)}{t} \\
&\quad\color{red}{\downarrow\;\text{Definition von $\sigma_0$.}}\\
&= \lim_{t \to 0} \frac{\frac{r^2}{\|x+th\|^2}(x+th) - \frac{r^2}{\|x\|^2}x}{t} \\
&\quad\color{red}{\downarrow\;\text{Nenner gleichnamig machen.}} \\
&= \lim_{t \to 0} \frac{\frac{r^2\|x\|^2}{\|x+th\|^2\|x\|^2}(x+th) - \frac{r^2\|x+th\|^2}{\|x+th\|^2\|x\|^2}x}{t} \\
&\quad\color{red}{\downarrow\;\text{Brüche zusammenziehen.}}\\
&= \lim_{t \to 0} \left(\frac{r^2\|x\|^2x + r^2\|x\|^2 th - r^2\|x+th\|^2x }{t\|x+th\|^2\|x\|^2}\right) \\
&\quad\color{red}{\downarrow\;\text{Faktorisieren.}}\\
&= \frac{r^2}{\|x\|^2} \lim_{t \to 0} \left(\frac{\left(\|x\|^2 - \|x+th\|^2\right) x + \|x\|^2 th }{t\|x+th\|^2}\right) \\
&\quad\color{red}{\downarrow\;\|x\|^2 - \|x+th\|^2 = \|x\|^2 - \|x\|^2 - 2t\langle x,h \rangle - t^2\|h\|^2 = - 2t\langle x,h \rangle - t^2\|h\|^2.}\\
&= \frac{r^2}{\|x\|^2} \lim_{t \to 0} \left(\frac{\left(- 2t\langle x,h \rangle - t^2\|h\|^2\right) x + \|x\|^2 th }{t\|x+th\|^2}\right) \\
&\quad\color{red}{\downarrow\;\text{Bruch aufspalten.}}\\
&= \frac{r^2}{\|x\|^2} \lim_{t \to 0} \left(-2\frac{\langle x,h \rangle x}{\|x+th\|^2} - \frac{t\|h\|^2x}{\|x+th\|^2} + \frac{\|x\|^2h}{\|x+th\|^2}\right) \\
&\quad\color{red}{\downarrow\;\text{Evaluiere den Grenzwert.}}\\
&= \frac{r^2}{\|x\|^2} \left(-2\frac{\langle x,h \rangle x}{\|x\|^2} + h\right)
\end{align*}
\]
Wo versteckt sich denn hier die lineare Abbildung $\mathrm{D}\sigma_0(x)$?🥵
|
Profil
|
Triceratops
Wenig Aktiv  Dabei seit: 28.04.2016 Mitteilungen: 6472
Wohnort: Berlin
 | Beitrag No.3, eingetragen 2020-10-05
|
\quoteon(2020-10-02 15:51 - Triceratops in Beitrag No. 1)
Die lineare Abbildung ist $y \mapsto D \sigma (x) y = $(das was du berechnet hast).
\quoteoff
|
Profil
|
Phoensie
Aktiv  Dabei seit: 11.04.2020 Mitteilungen: 442
Wohnort: Muri AG, Schweiz
 | Beitrag No.4, vom Themenstarter, eingetragen 2020-10-05
|
Ich muss eben anschliessend die Jacobimatrix dazu aufstellen, und das klappt irgendwie nicht.
Ich weiss, dass dort, wo in meiner letzten Zeile das $h$ steht, die Identitätsmatrix hingehört. Aber mit dem ersten Summanden kann ich nichts anfangen...
|
Profil
|
Wally
Senior  Dabei seit: 02.11.2004 Mitteilungen: 9721
Wohnort: Dortmund, Old Europe
 | Beitrag No.5, eingetragen 2020-10-05
|
\(\begingroup\)\(\newcommand{\D}{\displaystyle}\)
Hallo, du musst Matrizensprech lernen - eine wichtige Mathe-Fachsprache.
1. Wenn man Produkte Skalar*Vektor in Matrixschreibweise aufstellen will, gehört die 1x1-Matrix mit dem Skalar RECHTS vom Vektor (aufschreiben, und überzeuge dich, dass dann die Dimensionen passen).
2. Übersetzen von Skalarprodukt im Matrizensprache:
\( \langle x,h\rangle = \langle h,x \rangle = x^T h\)
Viele Grüße
Wally \(\endgroup\)
|
Profil
|
Phoensie hat die Antworten auf ihre/seine Frage gesehen. Phoensie hat selbst das Ok-Häkchen gesetzt. | Phoensie wird per Mail über neue Antworten informiert. |
|
All logos and trademarks in this site are property of their respective owner. The comments are property of their posters, all the rest © 2001-2023 by Matroids Matheplanet
This web site was originally made with PHP-Nuke, a former web portal system written in PHP that seems no longer to be maintained nor supported. PHP-Nuke is Free Software released under the GNU/GPL license.
Ich distanziere mich von rechtswidrigen oder anstößigen Inhalten, die sich trotz aufmerksamer Prüfung hinter hier verwendeten Links verbergen mögen. Lesen Sie die
Nutzungsbedingungen,
die Distanzierung,
die Datenschutzerklärung und das Impressum.
[Seitenanfang]
|