diff options
Diffstat (limited to '')
| -rw-r--r-- | buch/SeminarSpezielleFunktionen.pdf | bin | 0 -> 22768314 bytes | |||
| -rw-r--r-- | buch/papers/parzyl/img/D_plot.png | bin | 0 -> 712446 bytes | |||
| -rw-r--r-- | buch/papers/parzyl/img/plane.pdf | bin | 0 -> 2072 bytes | |||
| -rw-r--r-- | buch/papers/parzyl/img/v_plot.png | bin | 0 -> 637451 bytes | |||
| -rw-r--r-- | buch/papers/parzyl/main.tex | 4 | ||||
| -rw-r--r-- | buch/papers/parzyl/references.bib | 24 | ||||
| -rw-r--r-- | buch/papers/parzyl/teil0.tex | 110 | ||||
| -rw-r--r-- | buch/papers/parzyl/teil1.tex | 135 | ||||
| -rw-r--r-- | buch/papers/parzyl/teil2.tex | 44 | ||||
| -rw-r--r-- | buch/papers/parzyl/teil3.tex | 81 |
10 files changed, 317 insertions, 81 deletions
diff --git a/buch/SeminarSpezielleFunktionen.pdf b/buch/SeminarSpezielleFunktionen.pdf Binary files differnew file mode 100644 index 0000000..36b612f --- /dev/null +++ b/buch/SeminarSpezielleFunktionen.pdf diff --git a/buch/papers/parzyl/img/D_plot.png b/buch/papers/parzyl/img/D_plot.png Binary files differnew file mode 100644 index 0000000..f76e35b --- /dev/null +++ b/buch/papers/parzyl/img/D_plot.png diff --git a/buch/papers/parzyl/img/plane.pdf b/buch/papers/parzyl/img/plane.pdf Binary files differnew file mode 100644 index 0000000..c52c336 --- /dev/null +++ b/buch/papers/parzyl/img/plane.pdf diff --git a/buch/papers/parzyl/img/v_plot.png b/buch/papers/parzyl/img/v_plot.png Binary files differnew file mode 100644 index 0000000..b8c803e --- /dev/null +++ b/buch/papers/parzyl/img/v_plot.png diff --git a/buch/papers/parzyl/main.tex b/buch/papers/parzyl/main.tex index 528a2e2..fd2aea7 100644 --- a/buch/papers/parzyl/main.tex +++ b/buch/papers/parzyl/main.tex @@ -6,13 +6,13 @@ \chapter{Parabolische Zylinderfunktionen\label{chapter:parzyl}} \lhead{Parabolische Zylinderfunktionen} \begin{refsection} -\chapterauthor{Thierry Schwaller, Alain Keller} +\chapterauthor{Alain Keller und Thierry Schwaller} \input{papers/parzyl/teil0.tex} \input{papers/parzyl/teil1.tex} \input{papers/parzyl/teil2.tex} - +\input{papers/parzyl/teil3.tex} \printbibliography[heading=subbibliography] \end{refsection} diff --git a/buch/papers/parzyl/references.bib b/buch/papers/parzyl/references.bib index 494ff7c..40be69a 100644 --- a/buch/papers/parzyl/references.bib +++ b/buch/papers/parzyl/references.bib @@ -33,3 +33,27 @@ url = {https://doi.org/10.1016/j.acha.2017.11.004} } +@book{parzyl:whittaker, + place={Cambridge}, + edition={4}, + series={Cambridge Mathematical Library}, + title={A Course of Modern Analysis}, + DOI={10.1017/CBO9780511608759}, + publisher={Cambridge University Press}, + author={Whittaker, E. T. and Watson, G. N.}, + year={1996}, + collection={Cambridge Mathematical Library}} + +@book{parzyl:abramowitz-stegun, + added-at = {2008-06-25T06:25:58.000+0200}, + address = {New York}, + author = {Abramowitz, Milton and Stegun, Irene A.}, + edition = {ninth Dover printing, tenth GPO printing}, + interhash = {d4914a420f489f7c5129ed01ec3cf80c}, + intrahash = {23ec744709b3a776a1af0a3fd65cd09f}, + keywords = {Handbook}, + publisher = {Dover}, + timestamp = {2008-06-25T06:25:58.000+0200}, + title = {Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables}, + year = 1972 +}
\ No newline at end of file diff --git a/buch/papers/parzyl/teil0.tex b/buch/papers/parzyl/teil0.tex index 4b251db..4a6f8f4 100644 --- a/buch/papers/parzyl/teil0.tex +++ b/buch/papers/parzyl/teil0.tex @@ -4,42 +4,65 @@ % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % \section{Einleitung\label{parzyl:section:teil0}} -\rhead{Teil 0} -Die Laplace-Gleichung ist eine wichtige Gleichung in der Physik. -Mit ihr lässt sich zum Beispiel das elektrische Feld in einem ladungsfreien Raum bestimmen. -In diesem Kapitel wird die Lösung der Laplace-Gleichung im -parabolischen Zylinderkoordinatensystem genauer untersucht. -\subsection{Laplace Gleichung} -Die partielle Differentialgleichung -\begin{equation} - \Delta f = 0 -\end{equation} -ist als Laplace-Gleichung bekannt. -Sie ist eine spezielle Form der Poisson-Gleichung +\rhead{Einleitung} +%Die Laplace-Gleichung ist eine wichtige Gleichung in der Physik. +%Mit ihr lässt sich zum Beispiel das elektrische Feld in einem ladungsfreien Raum bestimmen. +%In diesem Kapitel wird die Lösung der Laplace-Gleichung im +%parabolischen Zylinderkoordinatensystem genauer untersucht. +Die Helmholtz-Gleichung ist eine wichtige Gleichung in der Physik. Mit ihr lässt sich zum Beispiel das Verhalten von elektromagnetischen Wellen beschreiben. +In diesem Kapitel wird die Lösung der Helmholtz-Gleichung im parabolischen Zylinderkoordinatensystem, die parabolischen Zylinderfunktionen, genauer untersucht. + +\subsection{Helmholtz-Gleichung} +Die partielle Differentialgleichung \begin{equation} - \Delta f = g + \nabla f = \lambda f \end{equation} -mit g als beliebige Funktion. -In der Physik hat die Laplace-Gleichung in verschieden Gebieten -verwendet, zum Beispiel im Elektromagnetismus. -Das Gaussche Gesetz in den Maxwellgleichungen +ist als Helmholtz-Gleichung bekannt und beschreibt das Eigenwert Problem für den Laplace-Operator. Sie ist eine der Gleichungen welche auftritt wenn die Wellengleichung \begin{equation} - \nabla \cdot E = \frac{\varrho}{\epsilon_0} -\label{parzyl:eq:max1} + \left ( \nabla^2 - \frac{1}{c^2}\frac{\partial^2}{\partial t^2} \right ) u(\textbf{r},t) + = + 0 \end{equation} -besagt das die Divergenz eines Elektrischen Feldes an einem -Punkt gleich der Ladung an diesem Punkt ist. -Das elektrische Feld ist hierbei der Gradient des elektrischen -Potentials +mit Hilfe von Separation \begin{equation} - \nabla \phi = E. -\end{equation} -Eingesetzt in \eqref{parzyl:eq:max1} resultiert + u(\textbf{r},t) = A(\textbf{r})T(t) +\end{equation} +in zwei Differentialgleichungen aufgeteilt wird. Die Helmholtz-Gleichung ist der Teil, welcher Zeit unabhängig ist \begin{equation} - \nabla \cdot \nabla \phi = \Delta \phi = \frac{\varrho}{\epsilon_0}, + \nabla^2 A(\textbf{r}) = \lambda A(\textbf{r}). \end{equation} -was eine Possion-Gleichung ist. -An Ladungsfreien Stellen, ist der rechte Teil der Gleichung $0$. + +%\subsection{Laplace Gleichung} +%Die partielle Differentialgleichung +%\begin{equation} +% \Delta f = 0 +%\end{equation} +%ist als Laplace-Gleichung bekannt. +%Sie ist eine spezielle Form der Poisson-Gleichung +%\begin{equation} +% \Delta f = g +%\end{equation} +%mit $g$ als beliebiger Funktion. +%In der Physik hat die Laplace-Gleichung in verschiedenen Gebieten +%verwendet, zum Beispiel im Elektromagnetismus. +%Das Gaussche Gesetz in den Maxwellgleichungen +%\begin{equation} +% \nabla \cdot E = \frac{\varrho}{\epsilon_0} +%\label{parzyl:eq:max1} +%\end{equation} +%besagt, dass die Divergenz eines elektrischen Feldes an einem +%Punkt gleich der Ladungsdichte an diesem Punkt ist. +%Das elektrische Feld ist hierbei der Gradient des elektrischen +%Potentials +%\begin{equation} +% \nabla \phi = E. +%\end{equation} +%Eingesetzt in \eqref{parzyl:eq:max1} resultiert +%\begin{equation} +% \nabla \cdot \nabla \phi = \Delta \phi = \frac{\varrho}{\epsilon_0}, +%\end{equation} +%was eine Poisson-Gleichung ist. +%An ladungsfreien Stellen ist der rechte Teil der Gleichung $0$. \subsection{Parabolische Zylinderkoordinaten \label{parzyl:subsection:finibus}} Im parabolischen Zylinderkoordinatensystem bilden parabolische Zylinder die Koordinatenflächen. @@ -51,7 +74,7 @@ Die Koordinate $(\sigma, \tau, z)$ sind in kartesischen Koordinaten ausgedrückt z & = z. \label{parzyl:coordRelationse} \end{align} -Wird $\tau$ oder $\sigma$ konstant gesetzt resultieren die Parabeln +Wird $\tau$ oder $\sigma$ konstant gesetzt, resultieren die Parabeln \begin{equation} y = \frac{1}{2} \left( \frac{x^2}{\sigma^2} - \sigma^2 \right) \end{equation} @@ -67,7 +90,6 @@ und konstantes $\sigma$ und die grünen ein konstantes $\tau$.} \label{parzyl:fig:cordinates} \end{figure} - Abbildung \ref{parzyl:fig:cordinates} zeigt das Parabolische Koordinatensystem. Das parabolische Zylinderkoordinatensystem entsteht wenn die Parabeln aus der Ebene gezogen werden. @@ -75,8 +97,6 @@ Ebene gezogen werden. Um in diesem Koordinatensystem integrieren und differenzieren zu können braucht es die Skalierungsfaktoren $h_{\tau}$, $h_{\sigma}$ und $h_{z}$. -\dots - Wird eine infinitessimal kleine Distanz $ds$ zwischen zwei Punkten betrachtet kann dies im kartesischen Koordinatensystem mit \begin{equation} @@ -85,7 +105,7 @@ kann dies im kartesischen Koordinatensystem mit \label{parzyl:eq:ds} \end{equation} ausgedrückt werden. -Das Skalierungsfaktoren werden so bestimmt, dass +Die Skalierungsfaktoren werden so bestimmt, dass \begin{equation} \left(ds\right)^2 = \left(h_{\sigma}d\sigma\right)^2 + \left(h_{\tau}d\tau\right)^2 + \left(h_z dz\right)^2 @@ -106,7 +126,7 @@ von \eqref{parzyl:coordRelationsa} - \eqref{parzyl:coordRelationse} als dz &= \frac{\partial \tilde{z} }{\partial \sigma} d\sigma + \frac{\partial \tilde{z} }{\partial \tau} d\tau + \frac{\partial \tilde{z} }{\partial \tilde{z}} d \tilde{z} - = d \tilde{z} \\ + = d \tilde{z} \end{align} substituiert. Wird diese Gleichung in der Form von \eqref{parzyl:eq:dspara} @@ -120,7 +140,7 @@ geschrieben, resultiert Daraus ergeben sich die Skalierungsfaktoren \begin{align} h_{\sigma} &= \sqrt{\sigma^2 + \tau^2}\\ - h_{\sigma} &= \sqrt{\sigma^2 + \tau^2}\\ + h_{\tau} &= \sqrt{\sigma^2 + \tau^2}\\ h_{z} &= 1. \end{align} \subsection{Differentialgleichung} @@ -216,26 +236,12 @@ und + \mu \right ) - i(\tau) + i(z) = 0 \end{equation} führt. -Wobei die Lösung von \eqref{parzyl:sep_dgl_3} -\begin{equation} - i(z) - = - A\cos{ - \left ( - \sqrt{\lambda + \mu}z - \right )} - + - B\sin{ - \left ( - \sqrt{\lambda + \mu}z - \right )} -\end{equation} -ist und \eqref{parzyl:sep_dgl_1} und \eqref{parzyl:sep_dgl_2} die sogenannten Weberschen Differentialgleichungen sind, welche die parabolischen Zylinder Funktionen als Lösung haben. + diff --git a/buch/papers/parzyl/teil1.tex b/buch/papers/parzyl/teil1.tex index f297189..673fa7f 100644 --- a/buch/papers/parzyl/teil1.tex +++ b/buch/papers/parzyl/teil1.tex @@ -5,24 +5,143 @@ % \section{Lösung \label{parzyl:section:teil1}} -\rhead{Problemstellung} -Die Differentialgleichungen \eqref{parzyl:sep_dgl_1} und \eqref{parzyl:sep_dgl_2} können mit einer Substitution -in die Whittaker Gleichung gelöst werden. +\rhead{Lösung} + +\eqref{parzyl:sep_dgl_3} beschriebt einen ungedämpften harmonischen Oszillator. +Die Lösung ist somit +\begin{equation} + i(z) + = + A\cos{ + \left ( + \sqrt{\lambda + \mu}z + \right )} + + + B\sin{ + \left ( + \sqrt{\lambda + \mu}z + \right )}. +\end{equation} +Die Differentialgleichungen \eqref{parzyl:sep_dgl_1} und \eqref{parzyl:sep_dgl_2} werden in \cite{parzyl:whittaker} +mit Hilfe der Whittaker Gleichung gelöst. \begin{definition} Die Funktion \begin{equation*} W_{k,m}(z) = e^{-z/2} z^{m+1/2} \, - {}_{1} F_{1}(\frac{1}{2} + m - k, 1 + 2m; z) + {}_{1} F_{1} + ( + {\textstyle \frac{1}{2}} + + m - k, 1 + 2m; z) \end{equation*} heisst Whittaker Funktion und ist eine Lösung - von + von der Whittaker Differentialgleichung \begin{equation} \frac{d^2W}{d z^2} + \left(-\frac{1}{4} + \frac{k}{z} + \frac{\frac{1}{4} - m^2}{z^2} \right) W = 0. + \label{parzyl:eq:whitDiffEq} \end{equation} \end{definition} +Es wird nun die Differentialgleichung bestimmt, welche +\begin{equation} + w = z^{-1/2} W_{k,-1/4} \left({\textstyle \frac{1}{2}} z^2\right) +\end{equation} +als Lösung hat. +Dafür wird $w$ in \eqref{parzyl:eq:whitDiffEq} eingesetzt woraus +\begin{equation} + \frac{d^2 w}{dz^2} - \left(\frac{1}{4} z^2 - 2k\right) w = 0 +\label{parzyl:eq:weberDiffEq} +\end{equation} +resultiert. DIese Differentialgleichung ist dieselbe wie +\eqref{parzyl:sep_dgl_2} und \eqref{parzyl:sep_dgl_2}, welche somit +$w$ als Lösung haben. +Da es sich um eine Differentialgleichung zweiter Ordnung handelt, hat sie nicht nur +eine sondern zwei Lösungen. +Die zweite Lösung der Whittaker-Gleichung ist $W_{k,-m} (z)$. +Somit hat \eqref{parzyl:eq:weberDiffEq} +\begin{align} + w_1(k, z) & = z^{-1/2} W_{k,-1/4} \left({\textstyle \frac{1}{2}} z^2\right)\\ + w_2(k, z) & = z^{-1/2} W_{k,1/4} \left({\textstyle \frac{1}{2}} z^2\right) +\end{align} +als Lösungen. +Mit der Hypergeometrischen Funktion ausgeschrieben ergeben sich die Lösungen +\begin{align} + \label{parzyl:eq:solution_dgl} + w_1(k,z) &= e^{-z^2/4} \, + {}_{1} F_{1} + ( + {\textstyle \frac{1}{4}} + - k, {\textstyle \frac{1}{2}} ; {\textstyle \frac{1}{2}}z^2) \\ + w_2(k,z) & = z e^{-z^2/4} \, + {}_{1} F_{1} + ({\textstyle \frac{3}{4}} + - k, {\textstyle \frac{3}{2}} ; {\textstyle \frac{1}{2}}z^2). +\end{align} +In der Literatur gibt es verschiedene Standartlösungen für $w(k,z)$ präsentiert. +Whittaker und Watson zeigen in \cite{parzyl:whittaker} eine Lösung +\begin{equation} + D_n(z) = \frac{ + \Gamma \left( {\textstyle \frac{1}{2}}\right) 2^{\frac{1}{2}n + \frac{1}{2}} z^{-\frac{1}{2}} + }{ + \Gamma \left( {\textstyle \frac{1}{2}} \right) - {\textstyle \frac{1}{2}} n) + } + M_{\frac{1}{2} n + \frac{1}{4}, - \frac{1}{4}} \left(\frac{1}{2}z^2\right) + + + \frac{ + \Gamma\left(-{\textstyle \frac{1}{2}}\right) 2^{\frac{1}{2}n + \frac{1}{4}} z^{-\frac{1}{2}} + }{ + \Gamma\left(- {\textstyle \frac{1}{2}} n\right) + } + M_{\frac{1}{2} n + \frac{1}{4}, \frac{1}{4}} \left(\frac{1}{2}z^2\right) +\end{equation} +welche die Differentialgleichung +\begin{equation} + \frac{d^2D_n(z)}{dz^2} + \left(n + \frac{1}{2} - \frac{1}{4} z^2\right)D_n(z) = 0 +\end{equation} +löst. -Lösung Folgt\dots - - +In \cite{parzyl:abramowitz-stegun} sind zwei Lösungen $U(a, z)$ und $V(a,z)$ +\begin{align} + U(a,z) &= + \cos\left[\pi \left({\textstyle \frac{1}{4}} + {\textstyle \frac{1}{2}} a\right)\right] Y_1 + - \sin\left[\pi \left({\textstyle \frac{1}{4}} + {\textstyle \frac{1}{2}} a\right)\right] Y_2 \\ + V(a,z) &= \frac{1}{\Gamma \left({\textstyle \frac{1}{2} - a}\right)} \left\{ + \sin\left[\pi \left({\textstyle \frac{1}{4}} + {\textstyle \frac{1}{2}} a\right)\right] Y_1 + + \cos\left[\pi \left({\textstyle \frac{1}{4}} + {\textstyle \frac{1}{2}} a\right)\right] Y_2 + \right\} +\end{align} +mit +\begin{align} + Y_1 &= \frac{1}{\sqrt{\pi}} + \frac{\Gamma\left({\textstyle \frac{1}{4} - + {\textstyle \frac{1}{2}}a}\right)} + {2^{\frac{1}{2} a + \frac{1}{4}}} w_1\\ + Y_2 &= \frac{1}{\sqrt{\pi}} + \frac{\Gamma\left({\textstyle \frac{3}{4} - + {\textstyle \frac{1}{2}}a}\right)} + {2^{\frac{1}{2} a - \frac{1}{4}}} w_2 +\end{align} +der Differentialgleichung +\begin{equation} + \frac{d^2 y}{d z^2} - \left(\frac{1}{4} z^2 + a\right) y = 0 +\end{equation} +beschrieben. Die Lösungen $U(a,z)$ und $V(a, z)$ können auch mit $D_n(z)$ +ausgedrückt werden +\begin{align} + U(a,z) &= D_{-a-1/2}(z) \\ + V(a,z) &= \frac{\Gamma \left({\textstyle \frac{1}{2}} + a\right)}{\pi} + \left[\sin\left(\pi a\right) D_{-a-1/2}(z) + D_{-a-1/2}(-x)\right]. +\end{align} +TODO Plot +\begin{figure} + \centering + \includegraphics[scale=0.3]{papers/parzyl/img/D_plot.png} + \caption{$D_a(z)$ mit unterschiedlichen Werten für $a$.} + \label{parzyl:fig:dnz} +\end{figure} +\begin{figure} + \centering + \includegraphics[scale=0.3]{papers/parzyl/img/v_plot.png} + \caption{$V(a,z)$ mit unterschiedlichen Werten für $a$.} + \label{parzyl:fig:Vnz} +\end{figure}
\ No newline at end of file diff --git a/buch/papers/parzyl/teil2.tex b/buch/papers/parzyl/teil2.tex index 3f890d0..4af6860 100644 --- a/buch/papers/parzyl/teil2.tex +++ b/buch/papers/parzyl/teil2.tex @@ -5,16 +5,19 @@ % \section{Anwendung in der Physik \label{parzyl:section:teil2}} -\rhead{Teil 2} +\rhead{Anwendung in der Physik} - -\subsection{Elektrisches Feld einer semi-infiniten Platte -\label{parzyl:subsection:bonorum}} -Die parabolischen Zylinderkoordinaten tauchen auf, wenn man das elektrische Feld einer semi-infiniten Platte finden will. -Das dies so ist kann im zwei Dimensionalen mit Hilfe von komplexen Funktionen gezeigt werden. Wobei die Platte dann nur eine Linie ist. +Die parabolischen Zylinderkoordinaten tauchen auf, wenn man das elektrische Feld einer semi-infiniten Platte, wie in Abbildung \ref{parzyl:fig:leiterplatte} gezeigt, finden will. +\begin{figure} + \centering + \includegraphics[width=0.9\textwidth]{papers/parzyl/img/plane.pdf} + \caption{Semi-infinite Leiterplatte} + \label{parzyl:fig:leiterplatte} +\end{figure} +Das dies so ist kann im zwei Dimensionalen mit Hilfe von komplexen Funktionen gezeigt werden. Die Platte ist dann nur eine Linie, was man in Abbildung TODO sieht. Jede komplexe Funktion $F(z)$ kann geschrieben werden als \begin{equation} - F(z) = U(x,y) + iV(x,y) \qquad z \in \mathbb{C}; x,y \in \mathbb{R}. + F(s) = U(x,y) + iV(x,y) \qquad s \in \mathbb{C}; x,y \in \mathbb{R}. \end{equation} Dabei muss gelten, falls die Funktion differenzierbar ist, dass \begin{equation} @@ -35,7 +38,7 @@ Aus dieser Bedingung folgt \frac{\partial^2 U(x,y)}{\partial y^2} = 0 - }_{\nabla^2U(x,y)=0} + }_{\displaystyle{\nabla^2U(x,y)=0}} \qquad \underbrace{ \frac{\partial^2 V(x,y)}{\partial x^2} @@ -43,26 +46,35 @@ Aus dieser Bedingung folgt \frac{\partial^2 V(x,y)}{\partial y^2} = 0 - }_{\nabla^2V(x,y) = 0}. + }_{\displaystyle{\nabla^2V(x,y) = 0}}. \end{equation} -Zusätzlich zeigen diese Bedingungen auch, dass die zwei Funktionen $U(x,y)$ und $V(x,y)$ orthogonal zueinander sind. +Zusätzlich kann auch gezeigt werden, dass die Funktion $F(z)$ eine winkeltreue Abbildung ist. Der Zusammenhang zum elektrischen Feld ist jetzt, dass das Potential an einem quellenfreien Punkt gegeben ist als \begin{equation} \nabla^2\phi(x,y) = 0. \end{equation} -Da dies bei komplexen differenzierbaren Funktionen gilt, wie Gleichung \ref{parzyl_e_feld_zweite_ab} zeigt, kann entweder $U(x,y)$ oder $V(x,y)$ von einer solchen Funktion als das Potential angesehen werden. Im weiteren wird für das Potential $U(x,y)$ verwendet. -Da die Funktion, welche nicht das Potential beschreibt, in weiteren angenommen als $V(x,y)$, orthogonal zum Potential ist, zeigt dies das Verhalten des elektrischen Feldes. -Um nun zu den parabolische Zylinderkoordinaten zu gelangen muss nur noch eine geeignete komplexe Funktion $F(z)$ gefunden werden, welche eine semi-infinite Platte beschreiben kann. Man könnte natürlich auch nach anderen Funktionen suchen, welche andere Bedingungen erfüllen und würde dann auf andere Koordinatensysteme stossen. Die gesuchte Funktion in diesem Fall ist +Dies ist eine Bedingung welche differenzierbare Funktionen, wie in Gleichung \ref{parzyl_e_feld_zweite_ab} gezeigt wird, bereits besitzen. +Nun kann zum Beispiel $U(x,y)$ als das Potential angeschaut werden +\begin{equation} + \phi(x,y) = U(x,y). +\end{equation} +Orthogonal zum Potential ist das elektrische Feld +\begin{equation} + E(x,y) = V(x,y). +\end{equation} +Um nun zu den parabolische Zylinderkoordinaten zu gelangen muss nur noch eine geeignete komplexe Funktion $F(s)$ gefunden werden, +welche eine semi-infinite Platte beschreiben kann. +Die gesuchte Funktion in diesem Fall ist \begin{equation} - F(z) + F(s) = - \sqrt{z} + \sqrt{s} = \sqrt{x + iy}. \end{equation} Dies kann umgeformt werden zu \begin{equation} - F(z) + F(s) = \underbrace{\sqrt{\frac{\sqrt{x^2+y^2} + x}{2}}}_{U(x,y)} + diff --git a/buch/papers/parzyl/teil3.tex b/buch/papers/parzyl/teil3.tex index 4e44bd6..972fd33 100644 --- a/buch/papers/parzyl/teil3.tex +++ b/buch/papers/parzyl/teil3.tex @@ -3,6 +3,81 @@ % % (c) 2020 Prof Dr Andreas Müller, Hochschule Rapperswil % -\section{Teil 3 -\label{parzyl:section:teil3}} -\rhead{Teil 3} +\section{Eigenschaften +\label{parzyl:section:Eigenschaften}} +\rhead{Eigenschaften} + +\subsection{Potenzreihenentwicklung + \label{parzyl:potenz}} +Die parabolischen Zylinderfunktionen, welche in Gleichung \ref{parzyl:eq:solution_dgl} gegeben sind, können auch als Potenzreihen geschrieben werden +\begin{align} + w_1(k,z) + &= + e^{-z^2/4} \, + {}_{1} F_{1} + ( + {\textstyle \frac{1}{4}} + - k, {\textstyle \frac{1}{2}} ; {\textstyle \frac{1}{2}}z^2) + = + e^{-\frac{z^2}{4}} + \sum^{\infty}_{n=0} + \frac{\left ( \frac{1}{4} - k \right )_{n}}{\left ( \frac{1}{2}\right )_{n}} + \frac{\left ( \frac{1}{2} z^2\right )^n}{n!} \\ + &= + e^{-\frac{z^2}{4}} + \left ( + 1 + + + \left ( \frac{1}{2} - 2k \right )\frac{z^2}{2!} + + + \left ( \frac{1}{2} - 2k \right )\left ( \frac{5}{2} - 2k \right )\frac{z^4}{4!} + + + \dots + \right ) +\end{align} +und +\begin{align} + w_2(k,z) + &= + ze^{-z^2/4} \, + {}_{1} F_{1} + ( + {\textstyle \frac{3}{4}} + - k, {\textstyle \frac{3}{2}} ; {\textstyle \frac{1}{2}}z^2) + = + ze^{-\frac{z^2}{4}} + \sum^{\infty}_{n=0} + \frac{\left ( \frac{3}{4} - k \right )_{n}}{\left ( \frac{3}{2}\right )_{n}} + \frac{\left ( \frac{1}{2} z^2\right )^n}{n!} \\ + &= + e^{-\frac{z^2}{4}} + \left ( + z + + + \left ( \frac{3}{2} - 2k \right )\frac{z^3}{3!} + + + \left ( \frac{3}{2} - 2k \right )\left ( \frac{7}{2} - 2k \right )\frac{z^5}{5!} + + + \dots + \right ). +\end{align} +Bei den Potenzreihen sieht man gut, dass die Ordnung des Polynoms im generellen ins unendliche geht. Es gibt allerdings die Möglichkeit für bestimmte k das die Terme in der Klammer gleich null werden und das Polynom somit eine endliche Ordnung $n$ hat. Dies geschieht bei $w_1(k,z)$ falls +\begin{equation} + k = \frac{1}{4} + n \qquad n \in \mathbb{N}_0 +\end{equation} +und bei $w_2(k,z)$ falls +\begin{equation} + k = \frac{3}{4} + n \qquad n \in \mathbb{N}_0. +\end{equation} + +\subsection{Ableitung} +Es kann gezeigt werden, dass die Ableitungen $\frac{\partial w_1(z,k)}{\partial z}$ und $\frac{\partial w_2(z,k)}{\partial z}$ einen Zusammenhang zwischen $w_1(z,k)$ und $w_2(z,k)$ zeigen. Die Ableitung von $w_1(z,k)$ nach $z$ kann über die Produktregel berechnet werden und ist gegeben als +\begin{equation} + \frac{\partial w_1(z,k)}{\partial z} = \left (\frac{1}{2} - 2k \right ) w_2(z, k -\frac{1}{2}) - \frac{1}{2} z w_1(z,k), +\end{equation} +und die Ableitung von $w_2(z,k)$ als +\begin{equation} + \frac{\partial w_2(z,k)}{\partial z} = w_1(z, k -\frac{1}{2}) - \frac{1}{2} z w_2(z,k). +\end{equation} +Über diese Eigenschaft können einfach weitere Ableitungen berechnet werden. + |
