aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers/kreismembran
diff options
context:
space:
mode:
authorntobler <nitobler@gmail.com>2022-08-15 09:38:05 +0200
committerntobler <nitobler@gmail.com>2022-08-15 09:38:05 +0200
commitb4b8e0a1e1c2ba5a8756c1b72bf2991eeb176f71 (patch)
tree18043b3c9912c787d5e05d2bc7dbc7daa75d0ef5 /buch/papers/kreismembran
parentcorrections (diff)
parentMerge pull request #49 from HeadAndToes/master (diff)
downloadSeminarSpezielleFunktionen-b4b8e0a1e1c2ba5a8756c1b72bf2991eeb176f71.tar.gz
SeminarSpezielleFunktionen-b4b8e0a1e1c2ba5a8756c1b72bf2991eeb176f71.zip
Merge branch 'master' of https://github.com/AndreasFMueller/SeminarSpezielleFunktionen
Diffstat (limited to '')
-rw-r--r--buch/papers/kreismembran/teil0.tex22
-rw-r--r--buch/papers/kreismembran/teil1.tex22
-rw-r--r--buch/papers/kreismembran/teil2.tex22
-rw-r--r--buch/papers/kreismembran/teil3.tex28
-rw-r--r--buch/papers/kreismembran/teil4.tex43
5 files changed, 71 insertions, 66 deletions
diff --git a/buch/papers/kreismembran/teil0.tex b/buch/papers/kreismembran/teil0.tex
index 6f55358..c6dac06 100644
--- a/buch/papers/kreismembran/teil0.tex
+++ b/buch/papers/kreismembran/teil0.tex
@@ -5,18 +5,18 @@
%
\section{Einleitung\label{kreismembran:section:teil0}}
\rhead{Membran}
-Eine Membran oder selten ein Schwingblatt ist laut Duden \cite{kreismembran:Duden:Membran} ein ``dünnes Blättchen aus Metall, Papier o. Ä., das durch seine Schwingungsfähigkeit geeignet ist, Schallwellen zu übertragen ...''.
+Eine Membran oder selten ein Schwingblatt ist laut Duden \cite{kreismembran:Duden:Membran} ein ``dünnes Blättchen aus Metall, Papier o. Ä., das durch seine Schwingungsfähigkeit geeignet ist, Schallwellen zu übertragen \dots''.
Ein dünnes Blättchen aus Metall zeig jedoch nicht die selben dynamischen Eigenschaften wie ein gespanntes Stück Papier.
Beschreibt man das dynamische Verhalten, muss zwischen einer dünnen Platte und einer Membrane unterschieden werden \cite{kreismembran:membrane_vs_thin_plate}.
Eine dünne Platte zum Beispiel aus Metall, wirkt selbst entgegen ihrer Deformation sobald sie gekrümmt wird.
Eine Membran auf der anderen Seite besteht aus einem Material, welches sich ohne Kraftaufwand verbiegen lässt wie zum Beispiel Papier.
-Bevor Papier als schwingende Membran betrachtet werden kann, wird jedoch noch eine Spannung $ T $ benötigt welche das Material daran hindert aus der Ruhelage gebracht zu werden.
+Bevor Papier als schwingende Membran betrachtet werden kann, wird jedoch noch eine Spannung $ T $ benötigt, welche das Material daran hindert, aus der Ruhelage gebracht zu werden.
Ein geläufiges Beispiel einer Kreismembran ist eine runde Trommel.
Sie besteht herkömmlicherweise aus einem Leder (Fell), welches auf einen offenen Zylinder (Zargen) aufgespannt wird.
Das Leder alleine erzeugt nach einem Aufschlag keine hörbaren Schwingungen.
Sobald das Fell jedoch über den Zargen gespannt wird, kann das Fell auf verschiedensten Weisen weiter schwingen, was für den Klang der Trommel verantwortlich ist.
-Wie genau diese Schwingungen untersucht werden können wird in der folgenden Arbeit diskutiert.
+Wie genau diese Schwingungen untersucht werden können, wird in der folgenden Arbeit diskutiert.
\subsection{Annahmen} \label{kreimembran:annahmen}
@@ -36,8 +36,8 @@ Das untersuchte Modell erfüllt folgende Eigenschaften:
\end{enumerate}
-\subsection{Wellengleichung} Um die Wellengleichung einer Membran herzuleiten wird vorerst eine schwingende Saite betrachtet.
-Es lohnt sich das Verhalten einer Saite zu beschreiben, da eine Saite das selbe Verhalten wie eine Membran aufweist mit dem Unterschied einer fehlenden Dimension.
+\subsection{Wellengleichung} Um die Wellengleichung einer Membran herzuleiten, wird vorerst eine schwingende Saite betrachtet.
+Es lohnt sich, das Verhalten einer Saite zu beschreiben, da eine Saite dasselbe Verhalten wie eine Membran aufweist, mit dem Unterschied einer fehlenden Dimension.
Die Verbindung zwischen Membran und Saite ist intuitiv ersichtlich, stellt man sich einen Querschnitt einer Trommel vor.
\begin{figure}
@@ -48,9 +48,10 @@ Die Verbindung zwischen Membran und Saite ist intuitiv ersichtlich, stellt man s
\end{center}
\end{figure}
-Abbildung \ref{kreismembran:im:Saite} ist ein infinitesimales Stück einer Saite mit Länge $ dx $ skizziert.
-Wie für die Membran ist die Annahme iii) gültig, keine Bewegung in die Richtung $ \hat{x} $.
-Um dies zu erfüllen muss der Punkt $ P_1 $ gleich stark in Richtung $ -\hat{x} $ gezogen werden wie der Punkt $ P_2 $ in Richtung $ \hat{x} $ gezogen wird. Ist $ T_1 $ die Kraft welche mit Winkel $ \alpha $ auf Punkt $ P_1 $ wirkt sowie $ T_2 $ und $ \beta$ das analoge für Punkt $ P_2 $ ist, so können die Kräfte
+In Abbildung \ref{kreismembran:im:Saite} ist ein infinitesimales Stück einer Saite mit Länge $ dx $ skizziert.
+Wie für die Membran ist die Annahme iii) gültig, es entsteht keine Bewegung entlang der $ x $-Achse.
+Um dies zu erfüllen, muss der Punkt $ P_1 $ gleich stark entgegen der $ x $-Achse gezogen werden wie der Punkt $ P_2 $ in Richtung der $ x $-Achse gezogen wird.
+Ist $ T_1 $ die Kraft, welche mit Winkel $ \alpha $ auf Punkt $ P_1 $ wirkt sowie $ T_2 $ und $ \beta$ das analoge für Punkt $ P_2 $ ist, so können die Kräfte
\begin{equation}\label{kreismembran:eq:no_translation}
T_1 \cos \alpha = T_2 \cos \beta = T
\end{equation}
@@ -81,9 +82,10 @@ Durch die Division mit $ dx $ entsteht
\begin{equation*}
\frac{1}{dx} \left[\frac{\partial u}{\partial x} \bigg|_{x_0 + dx} - \frac{\partial u}{\partial x} \bigg|_{x_0}\right] = \frac{\rho}{T}\frac{\partial^2 u}{\partial t^2}.
\end{equation*}
-Auf der linken Seite der Gleichung wird die Differenz der Steigungen durch die Intervalllänge geteilt, in anderen Worten die zweite Ableitung von $ u(x,t) $ nach $ x $ berechnet.
+Auf der linken Seite der Gleichung wird die Differenz der Steigungen durch die Intervalllänge geteilt.
+Wenn $ dx $ als unendlich kleines Stück betrachtet wird, ergibt sich als Grenzwert die zweite Ableitung von $ u(x,t) $ nach $ x $.
Der Term $ \frac{\rho}{T} $ wird durch $ c^2 $ ersetzt, da der Bruch für eine gegebene Membran eine positive Konstante sein muss.
-Somit resultiert die in der Literatur gebräuchliche Form
+Damit resultiert die in der Literatur gebräuchliche Form
\begin{equation}
\label{kreismembran:Ausgang_DGL}
\frac{1}{c^2}\frac{\partial^2u}{\partial t^2} = \Delta u.
diff --git a/buch/papers/kreismembran/teil1.tex b/buch/papers/kreismembran/teil1.tex
index f0d478f..f6ba7d1 100644
--- a/buch/papers/kreismembran/teil1.tex
+++ b/buch/papers/kreismembran/teil1.tex
@@ -23,7 +23,7 @@ Da es sich um eine Kreisscheibe handelt, werden Polarkoordinaten verwendet, so d
\frac1r
\frac{\partial}{\partial r}
+
- \frac{1}{r 2}
+ \frac{1}{r^2}
\frac{\partial^2}{\partial\varphi^2}
\label{buch:pde:kreis:laplace}
\end{equation*}
@@ -39,16 +39,16 @@ Daher ist die Membranabweichung im Punkt $(r,\varphi)$ $\in$ $\overline{\rm \Ome
\end{align*}
Um die Vergleichbarkeit der beiden nachfolgend vorgestellten Lösungsverfahren in Abschnitt \ref{kreismembran:vergleich} zu vereinfachen, werden keine Randbedingungen vorgegeben.
-Um eine eindeutige Lösung bestimmen zu können, werden die folgenden Anfangsbedingungen festgelegt zur zeit $t = \text{0}$:
+Um eine eindeutige Lösung bestimmen zu können, werden die folgenden Anfangsbedingungen festgelegt zur Zeit $t = \text{0}$:
\begin{align*}
u(r,\varphi, 0) &= f(r,\varphi)\\
u_t(r,\varphi, 0) &= g(r,\varphi).
\end{align*}
\subsection{Lösung\label{sub:lösung1}}
-Nun wird das in Abschnitt \ref{sub:aufgabestellung} vorgestellte Problem mit Hilfe der varibalen Trennungsmethode gelöst.
+Nun wird das in Abschnitt \ref{sub:aufgabestellung} vorgestellte Problem mit Hilfe der Separationsmethode gelöst.
\subsubsection{Ansatz der Separation der Variablen\label{subsub:ansatz_separation}}
-Bezug muss an dieser Stelle von einer Separation der Variablen ausgegangen werden:
+Hierfür wird folgenden Ansatz gemacht:
\begin{equation*}
u(r,\varphi, t) = F(r)G(\varphi)T(t)
\end{equation*}
@@ -64,26 +64,26 @@ Dies bedeutet, dass die Konstante negativ sein muss, also schreibt man $-\kappa^
T''(t) + c^2\kappa^2T(t) &= 0\\
r^2\frac{F''(r)}{F(r)} + r \frac{F'(r)}{F(r)} +\kappa^2 r^2 &= - \frac{G''(\varphi)}{G(\varphi)}.
\end{align*}
-In der zweiten Gleichung hängt die linke Seite nur von $r$ ab, während die rechte Seite nur von $\varphi$ abhängt. Sie müssen also wiederum gleich einer reellen Zahl $\nu$ sein. Also das:
+In der zweiten Gleichung hängt die linke Seite nur von $r$ ab, während die rechte Seite nur von $\varphi$ abhängt. Sie müssen also wiederum gleich einer reellen Zahl $\nu$ sein. Also:
\begin{align*}
- r^2F''(r) + rF'(r) + (\kappa^2 r^2 - \nu)F(r) &= 0 \\
- G''(\varphi) &= \nu G(\varphi).
+ r^2F''(r) + rF'(r) + (\kappa^2 r^2 - \nu)F(r) = 0 \quad \text{und} \quad
+ G''(\varphi) = \nu G(\varphi).
\end{align*}
\subsubsection{Lösung für $G(\varphi)$\label{subsub:lösung_G}}
Da für die zweite Gleichung Lösungen von Schwingungen erwartet werden, für die $G''(\varphi)=-\omega^2 G(\varphi)$ gilt, schreibt man die gemeinsame Konstante als $\nu=-\omega^2$, was die Formeln später vereinfacht. Also:
\begin{equation*}
- G(\varphi) = C_n \cos(\varphi) + D_n \sin(\varphi)
+ G(\varphi) = C_n \cos(\nu\varphi) + D_n \sin(\nu\varphi)
\label{eq:cos_sin_überlagerung}
\end{equation*}
\subsubsection{Lösung für $F(r)$\label{subsub:lösung_F}}
-Die Gleichung für $F$ hat die Gestalt (verweis auf \ref{buch:differentialgleichungen:bessel-operator})
+Die Gleichung für $F$ hat die Gestalt (Verweis auf \label{buch:differentialgleichungen:bessel-operator}
\begin{align}
r^2F''(r) + rF'(r) + (\kappa^2 r^2 - n^2)F(r) = 0
\label{eq:2nd_degree_PDE}
\end{align}
-Wir bereits in Kapitel \ref{buch:differntialgleichungen:section:bessel} gezeigt, sind die Besselfunktionen
+Wie bereits in Kapitel \ref{buch:differntialgleichungen:section:bessel} gezeigt, sind die Bessel-Funktionen
\begin{equation*}
J_{\nu}(x) = r^\nu \displaystyle\sum_{m=0}^{\infty} \frac{(-1)^m x^{2m}}{2^{2m+\nu}m! \Gamma (\nu + m+1)}
\end{equation*}
@@ -104,7 +104,7 @@ Durch Überlagerung aller Ergebnisse erhält man die Lösung
\end{align}
Dabei sind $m$ und $n$ ganze Zahlen, wobei $m$ für die Anzahl der Knotenkreise und $n$
-für die Anzahl der Knotenlinien steht. Es gibt bestimmte Bereiche auf der Membran, in denen es keine Bewegung oder Vibration gibt. Wenn der nicht schwingende Bereich ein Kreis ist, nennt man ihn einen Knotenkreis, und wenn er eine Linie ist, nennt man ihn ebenfalls eine Knotenlinie; siehe Abbildung \ref{buch:pde:kreis:fig:pauke}. $Jn(\kappa_{mn}r)$ ist die Besselfunktion $n$-ter Ordnung, wobei $\kappa mn$ die Wellenzahl und $r$ der Radius ist. $a_{mn}$ und $b_{mn}$ sind die zu bestimmenden Konstanten.
+für die Anzahl der Knotenlinien steht. Es gibt bestimmte Bereiche auf der Membran, in denen es keine Bewegung oder Vibration gibt. Wenn der nicht schwingende Bereich ein Kreis ist, nennt man ihn einen Knotenkreis, und wenn er eine Linie ist, nennt man ihn ebenfalls eine Knotenlinie (siehe Abbildung \ref{buch:pde:kreis:fig:pauke}). $J_n(\kappa_{mn}r)$ ist die Besselfunktion $n$-ter Ordnung, wobei $\kappa mn$ die Wellenzahl und $r$ der Radius ist. $a_{mn}$ und $b_{mn}$ sind die zu bestimmenden Konstanten.
\begin{figure}
\centering
diff --git a/buch/papers/kreismembran/teil2.tex b/buch/papers/kreismembran/teil2.tex
index 4fb139c..ec27bd3 100644
--- a/buch/papers/kreismembran/teil2.tex
+++ b/buch/papers/kreismembran/teil2.tex
@@ -7,7 +7,7 @@
Hermann Hankel (1839--1873) war ein deutscher Mathematiker, der für seinen Beitrag zur mathematischen Analysis und insbesondere für die nach ihm benannte Transformation bekannt ist.
Diese Transformation tritt bei der Untersuchung von Funktionen auf, die nur von der Entfernung des Ursprungs abhängen.
-Er studierte auch Funktionen, jetzt Hankel- oder Bessel- Funktionen genannt, der dritten Art.
+Er untersuchte auch Funktionen, jetzt Hankel- oder Bessel- Funktionen genannt, der dritten Art.
Die Hankel-Transformation, die die Bessel-Funktion enthält, taucht natürlich bei achsensymmetrischen Problemen auf, die in zylindrischen Polarkoordinaten formuliert sind.
In diesem Abschnitt werden die Theorie der Transformation und einige Eigenschaften der Grundoperationen erläutert.
@@ -17,12 +17,12 @@ Wir führen die Definition der Hankel-Transformation \cite{lokenath_debnath_inte
\mathscr{F}\{f(x,y)\} & = F(k,l)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{-i( \bm{\kappa}\cdot \mathbf{r})}f(x,y) \; dx \; dy,\label{equation:fourier_transform}\\
\mathscr{F}^{-1}\{F(x,y)\} & = f(x,y)=\frac{1}{2\pi}\int_{-\infty}^{\infty}\int_{-\infty}^{\infty}e^{i(\bm{\kappa}\cdot \mathbf{r})}F(k,l) \; dx \; dy \label{equation:inv_fourier_transform}
\end{align}
-wo $\mathbf{r}=(x,y)$ und $\bm{\kappa}=(k,l)$. Polarkoordinaten sind für diese Art von Problem am besten geeignet, mit $(x,y)=r(\cos\theta,\sin\theta)$ und $(k,l)=\kappa(\cos\phi,\sin\phi)$ findet man $\bm{\kappa}\cdot\mathbf{r}=\kappa r(\cos(\theta-\phi))$ und danach:
+definiert ist, wobei $\mathbf{r}=(x,y)$ und $\bm{\kappa}=(k,l)$. Polarkoordinaten sind für diese Art von Problem am besten geeignet. Mit $(x,y)=r(\cos\theta,\sin\theta)$ und $(k,l)=\kappa(\cos\phi,\sin\phi)$ findet man $\bm{\kappa}\cdot\mathbf{r}=\kappa r(\cos(\theta-\phi))$ und danach:
\begin{align}
F(k,\phi)=\frac{1}{2\pi}\int_{0}^{\infty}r \; dr \int_{0}^{2\pi}e^{-ikr\cos(\theta-\phi)}f(r,\theta) \; d\phi.
\label{equation:F_ohne_variable_wechsel}
\end{align}
-Dann wird angenommen dass, $f(r,\theta)=e^{in\theta}f(r)$, was keine strenge Einschränkung ist, weil die \textit{Fourier-Theorie} besagt, dass sich jede Funktion durch Überlagerung solcher Terme darstellen lässt. Es wird auch eine Änderung der Variabeln vorgenommen $\theta-\phi=\alpha-\frac{\pi}{2}$, um \eqref{equation:F_ohne_variable_wechsel} zu reduzieren:
+Dann wird angenommen, dass $f(r,\theta)=e^{in\theta}f(r)$, was keine strenge Einschränkung ist, weil die \textit{Fourier-Theorie} besagt, dass sich jede Funktion durch Überlagerung solcher Terme darstellen lässt. Es wird auch eine Änderung der Variabeln vorgenommen $\theta-\phi=\alpha-\frac{\pi}{2}$, um \eqref{equation:F_ohne_variable_wechsel} zu reduzieren:
\begin{align}
F(k,\phi)=\frac{1}{2\pi}\int_{0}^{\infty}rf(r) \; dr \int_{\phi_{0}}^{2\pi+\phi_{0}}e^{in(\phi-\frac{\pi}{2})+i(n\alpha-kr\sin\alpha)} \; d\alpha,
\label{equation:F_ohne_bessel}
@@ -34,7 +34,7 @@ Unter Verwendung der Integraldarstellung
J_n(\kappa r)=\frac{1}{2\pi}\int_{\phi_{0}}^{2\pi + \phi_{0}}e^{i(n\alpha-\kappa r \sin \alpha)} \; d\alpha
\label{equation:bessel_n_ordnung}
\end{equation*}
- der Besselfunktion vom Ordnung $n$ \eqref{buch:fourier:eqn:bessel-integraldarstellung} wird \eqref{equation:F_ohne_bessel} zu:
+ der Bessel-Funktion vom Ordnung $n$ \eqref{buch:fourier:eqn:bessel-integraldarstellung} wird \eqref{equation:F_ohne_bessel} zu:
\begin{align}
F(k,\phi)&=e^{in(\phi-\frac{\pi}{2})}\int_{0}^{\infty}rJ_n(\kappa r) f(r) \; dr \nonumber \\
&=e^{in(\phi-\frac{\pi}{2})}\tilde{f}_n(\kappa),
@@ -69,10 +69,10 @@ verwendet werden, um die Hankel-Transformation \eqref{equation:hankel} und ihre
Insbesondere die Hankel-Transformation der nullten Ordnung ($n=0$) und der ersten Ordnung ($n=1$) sind häufig nützlich, um Lösungen für Probleme mit der Laplace Gleichung in einer achsensymmetrischen zylindrischen Geometrie zu finden.
\subsection{Operatoreigenschaften der Hankel-Transformation \label{sub:op_properties_hankel}}
-In diesem Kapitel werden die operativen Eigenschaften der Hankel-Transformation aufgeführt. Der Beweis für ihre Gültigkeit wird jedoch nicht analysiert.
+In diesem Kapitel werden die operativen Eigenschaften der Hankel-Transformation aufgeführt. Die Beweise für ihre Gültigkeit werden jedoch nicht analysiert, diese sind im Buch \textit{Integral Tansforms and Their Applications} \cite{lokenath_debnath_integral_2015} zu finden.
\begin{satz}{Skalierung:}
- Wenn $\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)$, dann:
+ Wenn $\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)$, dann gilt:
\begin{equation*}
\mathscr{H}_n\{f(ar)\}=\frac{1}{a^{2}}\tilde{f}_n \left(\frac{\kappa}{a}\right), \quad a>0.
@@ -80,7 +80,7 @@ In diesem Kapitel werden die operativen Eigenschaften der Hankel-Transformation
\end{satz}
\begin{satz}{Parsevalsche Relation:}
-Wenn $\tilde{f}(\kappa)=\mathscr{H}_n\{f(r)\}$ und $\tilde{g}(\kappa)=\mathscr{H}_n\{g(r)\}$, dann:
+Wenn $\tilde{f}(\kappa)=\mathscr{H}_n\{f(r)\}$ und $\tilde{g}(\kappa)=\mathscr{H}_n\{g(r)\}$, dann gilt:
\begin{equation*}
\int_{0}^{\infty}rf(r)g(r) \; dr = \int_{0}^{\infty}\kappa\tilde{f}(\kappa)\tilde{g}(\kappa) \; d\kappa.
@@ -88,20 +88,20 @@ Wenn $\tilde{f}(\kappa)=\mathscr{H}_n\{f(r)\}$ und $\tilde{g}(\kappa)=\mathscr{H
\end{satz}
\begin{satz}{Hankel-Transformationen von Ableitungen:}
-Wenn $\tilde{f}_n(\kappa)=\mathscr{H}_n\{f(r)\}$, dann:
+Wenn $\tilde{f}_n(\kappa)=\mathscr{H}_n\{f(r)\}$, dann gilt:
\begin{align*}
&\mathscr{H}_n\{f'(r)\}=\frac{\kappa}{2n}\left[(n-1)\tilde{f}_{n+1}(\kappa)-(n+1)\tilde{f}_{n-1}(\kappa)\right], \quad n\geq1, \\
&\mathscr{H}_1\{f'(r)\}=-\kappa \tilde{f}_0(\kappa),
\end{align*}
-vorausgesetzt dass $[rf(r)]$ verschwindet wenn $r\to0$ und $r\to\infty$.
+vorausgesetzt, dass $rf(r)$ verschwindet wenn $r\to0$ und $r\to\infty$.
\end{satz}
\begin{satz}
-Wenn $\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)$, dann:
+Wenn $\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)$, dann gilt:
\begin{equation*}
\mathscr{H}_n \left\{ \left( \nabla^2 - \frac{n^2}{r^2} f(r)\right)\right\}= \mathscr{H}_n\left\{\frac{1}{r}\frac{d}{dr}\left(r\frac{df}{dr}\right) - \frac{n^2}{r^2}f(r)\right\}=-\kappa^2\tilde{f}_{n}(\kappa),
\end{equation*}
-bereitgestellt dass $rf'(r)$ und $rf(r)$ verschwinden für $r\to0$ und $r\to\infty$.
+bereitgestellt, dass $rf'(r)$ und $rf(r)$ verschwinden für $r\to0$ und $r\to\infty$.
\end{satz}
diff --git a/buch/papers/kreismembran/teil3.tex b/buch/papers/kreismembran/teil3.tex
index 276f911..a9dcd95 100644
--- a/buch/papers/kreismembran/teil3.tex
+++ b/buch/papers/kreismembran/teil3.tex
@@ -6,25 +6,22 @@
\section{Lösungsmethode 2: Transformationsmethode
\label{kreismembran:section:teil3}}
\rhead{Lösungsmethode 2: Transformationsmethode}
-Die Hankel-Transformation wird dann zur Lösung der Differentialgleichung verwendet. Es müssen jedoch einige Änderungen an dem Problem vorgenommen werden, damit es mit den Annahmen übereinstimmt, die für die Verwendung der Hankel-Transformation erforderlich sind. Das heisst, dass die Funktion $u$ nur von der Entfernung zum Ausgangspunkt abhängt.
+Die Hankel-Transformation kann hier zur Lösung der Differentialgleichung verwendet werden. Es müssen jedoch einige Änderungen an dem Problem vorgenommen werden, damit es mit den Annahmen übereinstimmt, die für die Verwendung der Hankel-Transformation erforderlich sind. Das heisst, dass die Funktion $u$ nur von der Entfernung zum Ausgangspunkt abhängt.
\subsubsection{Transformation und Reduktion auf eine algebraische Gleichung\label{subsub:transf_reduktion}}
Führt man also das Konzept einer unendlichen und achsensymmetrischen Membran ein:
-\begin{equation*}
+\begin{align}
\frac{\partial^2u}{\partial t^2}
=
c^2 \left(\frac{\partial^2 u}{\partial r^2}
+
\frac{1}{r}
- \frac{\partial u}{\partial r} \right), \quad 0<r<\infty, \quad t>0
- \label{eq:PDE_inf_membane}
-\end{equation*}
-
-\begin{align}
- u(r,0)=f(r), \quad u_t(r,0) = g(r), \quad \text{für} \quad 0<r<\infty
+ \frac{\partial u}{\partial r} \right), \quad 0<r<\infty, \quad t>0 \label{eq:PDE_inf_membane} \\
+ u(r,0)=f(r), \quad u_t(r,0) = g(r), \quad \text{für} \quad 0<r<\infty.
\label{eq:PDE_inf_membane_RB}
\end{align}
+
Mit Anwendung der Hankel-Transformation nullter Ordnung in Abhängigkeit von $r$ auf die Gleichungen \eqref{eq:PDE_inf_membane} und \eqref{eq:PDE_inf_membane_RB}:
\begin{align}
@@ -45,7 +42,7 @@ Die allgemeine Lösung für diese Gleichung lautet, wie in Abschnitt \eqref{eq:c
\begin{equation*}
\tilde{u}(\kappa,t)=\tilde{f}(\kappa)\cos(c\kappa t) + \frac{1}{c\kappa}\tilde{g}(\kappa)\sin(c\kappa t).
\end{equation*}
-Wendet man an nun die inverse Hankel-Transformation an, so erhält man die formale Lösung
+Wendet man nun die inverse Hankel-Transformation an, so erhält man die formale Lösung
\begin{align}
u(r,t)=\int_{0}^{\infty}\kappa\tilde{f}(\kappa)\cos(c\kappa t) J_0(\kappa r) \; d\kappa +\frac{1}{c}\int_{0}^{\infty}\tilde{g}(\kappa)\sin(c\kappa t)J_0(\kappa r) \; d\kappa.
@@ -53,7 +50,7 @@ Wendet man an nun die inverse Hankel-Transformation an, so erhält man die forma
\end{align}
\subsubsection{Erfüllung der Anfangsbedingungen\label{subsub:erfüllung_AB}}
-Es wird in Folgenden davon ausgegangen, dass sich die Membran verformt und zum Zeitpunkt $t=0$ freigegeben wird
+Es wird im Folgenden davon ausgegangen, dass sich die Membran verformt und zum Zeitpunkt $t=0$ freigegeben wird
\begin{equation*}
u(r,0)=f(r)=Aa(r^2 + a^2)^{-\frac{1}{2}}, \quad u_t(r,0)=g(r)=0
@@ -62,10 +59,17 @@ so dass $\tilde{g}(\kappa)\equiv 0$ und
\begin{equation*}
\tilde{f}(\kappa)=Aa\int_{0}^{\infty}r(a^2 + r^2)^{-\frac{1}{2}} J_0 (\kappa r) \; dr=\frac{Aa}{\kappa}e^{-a\kappa}.
\end{equation*}
+
+Aus der Laplace-Transformation und unter Verwendung der Skalierungseigenschaft ergibt sich, dass
+
+\begin{align*}
+ \int_{0}^{\infty}e^{-px} J_0(\kappa x) \; dx = \frac{1}{\sqrt{\kappa^2 + p^2}}.
+\end{align*}
+
Die formale Lösung \eqref{eq:formale_lösung} lautet also
\begin{align*}
u(r,t)&=Aa\int_{0}^{\infty}e^{-a\kappa} J_0(\kappa r)\cos(c\kappa t) \; dk=AaRe\int_{0}^{\infty}e^{-\kappa(a+ict)} J_0(\kappa r) \; dk\\
- &=AaRe\left\{r^2+\left(a+ict\right)^2\right\}^{-\frac{1}{2}}
+ &=AaRe\left\{r^2+\left(a+ict\right)^2\right\}^{-\frac{1}{2}}.
\end{align*}
Nimmt man jedoch die allgemeine Lösung durch Überlagerung,
@@ -80,6 +84,6 @@ kann man die Lösungsmethoden 1 und 2 vergleichen.
\label{kreismembran:vergleich}}
Bei der Analyse der Gleichungen \eqref{eq:lösung_endliche_generelle} und \eqref{eq:lösung_unendliche_generelle} fällt sofort auf, dass die Gleichung \eqref{eq:lösung_unendliche_generelle} nicht mehr von $m$ und $n$ abhängt, sondern nur noch von $n$ \cite{nishanth_p_vibrations_2018}.
Das macht Sinn, denn $n$ beschreibt die Anzahl der Knotenlinien, welche unter der Annahme einer rotationssymmetrischen Lösung nicht vorhanden sein können. Tatsächlich werden $a_{m0}$, $b_{m0}$ und $\kappa_{m0}$ in $a_m$, $b_m$ bzw. $\kappa_m$ umbenannt. Die beiden Termen $\cos(n\varphi)$ und $\sin(n\varphi)$ verschwinden ebenfalls, da für $n=0$ der $\cos(n\varphi)$ gleich 1 und der $\sin(n \varphi)$ gleich 0 ist.
-Die Funktion hängt also nicht mehr von der Besselfunktionen $n$-ter Ordnung ab, sondern nur von der $0$-ter Ordnung.
+Die Funktion hängt also nicht mehr von der Besselfunktionen $n$-ter Ordnung ab, sondern nur von der nullter Ordnung.
diff --git a/buch/papers/kreismembran/teil4.tex b/buch/papers/kreismembran/teil4.tex
index 74bb87d..01a6029 100644
--- a/buch/papers/kreismembran/teil4.tex
+++ b/buch/papers/kreismembran/teil4.tex
@@ -12,7 +12,7 @@ Jedes Element $ U_{ij} $ steh für die Auslenkung der Membran $ u(x,y,t) $ an d
Zwischen benachbarten Elementen in der Matrix $ U $ liegt immer der Abstand $ dh $, eine Inkrementierung von $ i $ oder $ j $ entspricht somit einem Schritt in Richtung $ x $ oder $ y $ von Länge $ dh $ auf der Membran.
Die zeitliche Dimension wird in Form des Array $ U[] $ aus $ z \times U $ Matrizen dargestellt, wobei $ z $ der Anzahl Zeitschritten entspricht.
Das Element auf Zeile $ i $, Spalte $ j $ der $ w $-ten Matrix von $ U[] $ also $ U[w]_{ij} $ entspricht somit der Auslenkung $ u(i,j,w) $.
-Da die DGL von Zweiter Ordnung ist, reicht eine Zustandsvariabel pro Membran-Element nicht aus.
+Da die DGL von zweiter Ordnung ist, reicht eine Zustandsvariabel pro Membran-Element nicht aus.
Es wird neben der Auslenkung auch die Geschwindigkeit jedes Membran-Elementes benötigt um den Zustand eindeutig zu beschreiben.
Dazu existiert neben $ U[] $ ein analoger Array $ V[] $ welcher die Geschwindigkeiten aller Membran-Elementen repräsentiert.
$ V[w]_{ij} $ entspricht also $ \dot{u}(i,j,w) $.
@@ -67,39 +67,39 @@ Die Filtermaske kann dann auf jedes Element einzeln angewendet werden mit einer
\subsection{Simulation: Kreisförmige Membran}
Als Beispiel soll nun eine schwingende kreisförmige Membran simuliert werden.
-\paragraph{Initialisierung}
-Die Anzahl der simulierten Elementen soll $ m \times n $ was dementsprechend die Dimensionen von $ U $ und $ V $ vorgibt.
+\subsubsection{Initialisierung}
+Die Anzahl der simulierten Elemente soll $ m \times n $ sein, was die Dimensionen von $ U $ und $ V $ vorgibt.
Als Anfangsbedingung wird eine Membran gewählt, welche bei $ t=0 $ mit einer Gauss-Kurve ausgelenkt wird.
Die Membran soll sich zu Beginn nicht bewegen, also wird $ V[0] $ mit Nullen initialisiert.
Die Auslenkung kann kompakt erreicht werden, wenn $ U[0] $ als Null-Matrix mit einer $ 1 $ in der Mitte initialisiert wird.
Diese Matrix wird anschliessend mit einer Filtermaske in Form einer Gauss-Glocke gefaltet.
-Die Faltung mit einer Gauss-Glocke ist in Programmen wie Matlab eine Standartfunktion, da dies einm Tiefpassfilter in der Bildverarbeitung entspricht.
+Die Faltung mit einer Gauss-Glocke ist in Programmen wie Matlab eine Standartfunktion, da dies einem Tiefpassfilter in der Bildverarbeitung entspricht.
-\paragraph{Rand}
+\subsubsection{Rand}
Bislang ist die definierte Matrix rechteckig.
-Um eine kreisförmige Membran zu simulieren muss der Rand angepasst werden.
+Um eine kreisförmige Membran zu simulieren, muss der Rand angepasst werden.
Da in den meisten Programme keine Möglichkeit besteht, mit runden Matrizen zu rechnen, wird der Rand in der Berechnung des Folgezustandes implementiert.
-Der Rand bedeutet, das Membran-Elemente auf dem Rand sich nicht Bewegen können.
-Die Position sowie die Geschwindigkeit aller Elemente welche nicht auf der definierten Membran sind müssen zu beliebiger Zeit $0$ entsprechen.
+Der Rand bedeutet, dass Membran-Elemente auf dem Rand sich nicht Bewegen können.
+Die Position, sowie die Geschwindigkeit aller Elemente, welche nicht auf der definierten Membran sind, müssen zu beliebiger Zeit $0$ sein.
Hierzu wird eine Maske $M$ erstellt.
Diese Maske besteht aus einer binären Matrix von identischer Dimension wie $ U $ und $ V $.
-Ist in der Matrix $M$ eine $1$ abgebildet so ist an jener stelle ein Element der Membran, ist es eine $0$ so befindet sich dieses Element auf dem Rand oder ausserhalb der Membran.
+Ist in der Matrix $M$ eine $1$ abgebildet, so ist an jener Stelle ein Element der Membran, ist es eine $0$ so befindet sich dieses Element auf dem Rand oder ausserhalb der Membran.
In dieser Anwendung ist $M$ eine Matrix mit einem Kreis voller $1$ umgeben von $0$ bis an den Rand der Matrix.
-Die Maske wird angewendet indem das Resultat des nächsten Zustandes noch mit der Maske elementweise multipliziert wird.
+Die Maske wird angewendet, indem das Resultat des nächsten Zustandes noch mit der Maske elementweise multipliziert wird.
Der Folgezustand kann also mit den Gleichungen
\begin{align}
\label{kreismembran:eq:folge_U}
- U[w+1] &= (U[w] + dt \cdot V[w])*M\\
+ U[w+1] &= (U[w] + dt \cdot V[w])\odot M\\
\label{kreismembran:eq:folge_V}
- V[w+1] &= (V[w] + dt \cdot \Delta_h u \cdot c^2)*M
+ V[w+1] &= (V[w] + dt \cdot \Delta_h u \cdot c^2)\odot M
\end{align}
berechnet werden.
-\paragraph{Simulation}
+\subsubsection{Simulation}
Mit den gegebenen Gleichungen \ref{kreismembran:eq:folge_U} und \ref{kreismembran:eq:folge_V} das Verhalten der Membran mit einem Loop über das zu untersuchende Zeitintervall berechnet werden.
In der Abbildung \ref{kreismembran:im:simres_rund} sind Simulationsresultate zu sehen.
-Die Erste Figur zeigt die Ausgangslage gefolgt von den Auslenkungen nach jeweils $ 50 $ weiteren Iterationsschritten.
+Die erste Figur zeigt die Ausgangslage gefolgt von den Auslenkungen nach jeweils $ 50 $ weiteren Iterationsschritten.
Es ist zu erkennen, wie sich die Störung vom Zentrum an den Rand ausbreitet.
-Erreicht die Störung den Rand wird sie reflektiert und nähert sich dem Zentrum.
+Erreicht die Störung den Rand, wird sie reflektiert und nähert sich dem Zentrum.
\begin{figure}
\begin{center}
@@ -117,19 +117,19 @@ Erreicht die Störung den Rand wird sie reflektiert und nähert sich dem Zentrum
\end{figure}
\subsection{Simulation: Unendliche Membran}
-Um eine unendlich grosse Membran zu simulieren könnte der unpraktische weg gewählt werden die Matrix unendlich gross zu definieren, dies wird jedoch spätestens bei der numerischen Berechnung seine Probleme mit sich bringen.
-Etwas geeigneter ist es die Matrix so gross wie möglich zu definieren wie es die Kapazitäten erlauben.
+Um eine unendlich grosse Membran zu simulieren, könnte der unpraktische Weg gewählt werden, die Matrix unendlich gross zu definieren, dies wird jedoch spätestens bei der numerischen Berechnung seine Probleme mit sich bringen.
+Etwas geeigneter ist es, die Matrix so gross wie möglich zu definieren, wie es die Kapazitäten erlauben.
Wenn anschliessend nur das Verhalten im Zentrum, bei der Störung beobachtet wird, verhaltet sich die Membran wie eine unendliche.
Dies aber nur bis die Störung am Rand reflektiert wird und wieder das innere zu beobachtende Zentrum beeinflusst.
Soll erst gar keine Reflexion entstehen, muss ein Absorber modelliert werden welcher die Störung möglichst ohne Reflexion aufnimmt.
-\paragraph{Absorber}
+\subsubsection{Absorber}
Sehr knapp formuliert entstehen Reflexionen, wenn eine Welle von einem Material in ein anderes Material mit unterschiedlichen Eigenschaften eindringen möchte.
Je unterschiedlicher und abrupter der Übergang zwischen den Materialien umso ausgeprägter die Reflexion.
In diesem Fall sind die Eigenschaften vorgegeben.
Im Zentrum soll sich die Membran verhalten, wie von der DGL vorgegeben, am Rand jedoch muss sich jedes Membran-Element in der Ausgangslage befinden.
Der Spielraum welcher dem Absorber übrig bleibt ist die Art der Überganges.
-Bei der endlichen kreisförmigen Membran hat die Maske $M$ ein binärer Übergang von Membran zu Rand bezweckt.
+Bei der endlichen kreisförmigen Membran hat die Maske $M$ einen binären Übergang von Membran zu Rand bezweckt.
Anstelle dieses abrupten Wechsels wird nun eine Maske definiert, welche graduell von Membran $1$ zu Rand-Element $0$ wechselt.
Die Elemente werden auf Basis ihres Abstand $r$ zum Zentrum definiert.
Der Abstand entspricht
@@ -156,11 +156,10 @@ In der Abbildung \ref{kreismembran:im:masks} ist der Unterschied der beiden Mask
\label{kreismembran:im:masks}
\end{center}
\end{figure}
-\paragraph{Simulation}
+\subsubsection{Simulation}
Bis auf die Absorber-Maske kann nun identisch zur endlichen Membran simuliert werden.
Auch hier wurde eine Gauss-Glocke als Anfangsbedingung gewählt.
Die Simulationsresultate von Abbildung \ref{kreismembran:im:simres_unendlich}
-
\begin{figure}
\begin{center}
@@ -183,7 +182,7 @@ Dieses Verhalten spricht für den Absorber-Ansatz, es soll jedoch erwähnt sein,
Die DGL \ref{kreismembran:Ausgang_DGL} welche simuliert wird geht jedoch von der Annahme \ref{kreimembran:annahmen} iv) aus, dass die Membran keine Art von Dämpfung erfährt.
\section{Schlusswort}
-Auch wenn ein Physikalisches Verhalten bereits durch Annahmen und Annäherungen deutlich vereinfacht wird, bestehen auch dann noch eine Vielzahl von Lösungsansätzen.
+Auch wenn ein physikalisches Verhalten bereits durch Annahmen und Annäherungen deutlich vereinfacht wird, bestehen auch dann noch eine Vielzahl von Lösungsansätzen.
Lösungen einer unendlich grosse Membran scheinen fern der Realität zu sein, doch dies darf es im Sinne der Mathematik.
Und wer weis, für eine Ameise auf einem Trampolin ist eine unendliche Membran vielleicht eine ganz gute Annäherung.