aboutsummaryrefslogtreecommitdiffstats
path: root/buch/papers/kreismembran/teil2.tex
blob: ec27bd34c1f10f29f57c9cc8b2ba19afd8ca1157 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
%
% (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil
%

\section{Die Hankel Transformation \label{kreismembran:section:teil2}}
\rhead{Die Hankel Transformation}

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 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.

\subsubsection{Definition der Hankel-Transformation \label{subsub:hankel_tansformation}}
Wir führen die Definition der Hankel-Transformation \cite{lokenath_debnath_integral_2015} aus der zweidimensionalen Fourier-Transformation und ihrer Umkehrung ein, die durch:
\begin{align}
	\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}
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:
\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}
\end{align}
wo $\phi_{0}=(\frac{\pi}{2}-\phi)$.

Unter Verwendung der Integraldarstellung
\begin{equation*}
	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 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),
	\label{equation:F_mit_bessel_step_2}
\end{align}
wo $\tilde{f}_n(\kappa)$ ist die \textit{Hankel-Transformation} von $f(r)$ und ist formell definiert durch:
\begin{align}
	\mathscr{H}_n\{f(r)\}=\tilde{f}_n(\kappa)=\int_{0}^{\infty}rJ_n(\kappa r) f(r) \; dr.
	\label{equation:hankel}
\end{align}

\subsubsection{Inverse Hankel-Transformation \label{subsub:inverse_hankel_tansformation}}
Wie bei der Entwicklung der Hankel-Transformation können auch für die Umkehrformel Analogien zur Fourier-Transformation hergestellt werden. Vergleicht man die beiden Transformationen, so stellt man fest, dass sie sehr ähnlich sind, wenn man den Term $J_n(\kappa r)$ der Hankel-Transformation durch $e^{-i( \bm{\kappa}\cdot \mathbf{r})}$ der Fourier-Transformation ersetzt. Diese beide Funktionen sind orthogonal, und bei orthogonalen Matrizen genügt bekanntlich die Transponierung, um sie zu invertieren. Da das Skalarprodukt der Bessel-Funktionen jedoch nicht dasselbe ist wie das der Exponentialfunktionen, muss man durch $\kappa\; d\kappa$ statt nur durch $d\kappa$ integrieren, um die Umkehrfunktion zu erhalten.

Von \eqref{equation:hankel} also ist, die inverse \textit{Hankel-Transformation} so definiert:
\begin{align}
	\mathscr{H}^{-1}_n\{\tilde{f}_n(\kappa)\}=f(r)=\int_{0}^{\infty}\kappa J_n(\kappa r) \tilde{f}_n(\kappa) \; d\kappa.
	\label{equation:inv_hankel}
\end{align}

Anstelle von $\tilde{f}_n(\kappa)$, wird häufig einfach $\tilde{f}(\kappa)$ für die Hankel-Transformation verwendet, indem die Ordnung angegeben wird.
Die Integrale \eqref{equation:hankel} und \eqref{equation:inv_hankel} existieren für bestimmte grosse Klassen von Funktionen, die normalerweise in physikalischen Anwendungen vorkommen.

Alternativ dazu kann die berühmte Hankel-Integralformel 

\begin{align*}
	f(r) = \int_{0}^{\infty}\kappa J_n(\kappa r) \; d\kappa \int_{0}^{\infty} p J_n(\kappa p)f(p) \; dp,
	\label{equation:hankel_integral_formula}
\end{align*}
verwendet werden, um die Hankel-Transformation \eqref{equation:hankel} und ihre Umkehrung \eqref{equation:inv_hankel} zu definieren.

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. 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 gilt:
	
	\begin{equation*}
		\mathscr{H}_n\{f(ar)\}=\frac{1}{a^{2}}\tilde{f}_n \left(\frac{\kappa}{a}\right), \quad a>0.
	\end{equation*}
\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 gilt:

\begin{equation*}
	\int_{0}^{\infty}rf(r)g(r) \; dr = \int_{0}^{\infty}\kappa\tilde{f}(\kappa)\tilde{g}(\kappa) \; d\kappa.
\end{equation*}
\end{satz}

\begin{satz}{Hankel-Transformationen von Ableitungen:}
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$.
\end{satz}

\begin{satz}
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$.
\end{satz}