Besselfunktionen erster Art
Inhaltsverzeichnis
Programmbeschreibung
Dieses Applet ermöglicht die Berechnung und graphische Darstellung der Besselfunktionen erster Art und $n$–ter Ordnung entsprechend der Reihendarstellung:
- $${\rm J}_n (x) = \sum\limits_{k=0}^{\infty}\frac{(-1)^k \cdot (x/2)^{n \hspace{0.05cm} + \hspace{0.05cm} 2 \hspace{0.02cm}\cdot \hspace{0.05cm}k}}{k! \cdot (n+k)!} \hspace{0.05cm}.$$
- Graphisch dargestellt können die Funktionen ${\rm J}_n (x)$ für die Ordnung $n=0$ bis $n=9$ in verschiedenen Farben werden.
- Die linke Ausgabe liefert die Funktionswerte ${\rm J}_0 (x = x_1)$, ... , ${\rm J}_9 (x = x_1)$ für einen per Slider einstellbaren Wert $x_1$ im Bereich $0 \le x_1 \le 15$ mit Schrittweite $0.5$.
- Die rechte Ausgabe liefert die Funktionswerte ${\rm J}_0 (x = x_2)$, ... , ${\rm J}_9 (x = x_2)$ für einen per Slider einstellbaren Wert $x_2$ (gleicher Wertebereich und Schrittweite wie links).
Theoretischer Hintergrund
Allgemeines zu den Besselfunktionen
Besselfunktionen (oder auch Zylinderfunktionen) sind Lösungen der Besselschen Differentialgleichung der Form
- $$x^2 \cdot \frac{ {\rm d}^2}{{\rm d}x^2}\ {\rm J}_n (x) \ + \ x \cdot \frac{ {\rm d}}{{\rm d}x}\ {\rm J}_n (x) \ + \ (x^2 - n^2) \cdot {\rm J}_n (x)= 0. $$
Hierbei handelt es sich um eine gewöhnliche lineare Differentialgleichung zweiter Ordnung. Der Parameter $n$ ist meistens ganzzahlig, so auch in diesem Programm. Diese bereits 1844 von Friedrich Wilhelm Bessel eingeführten mathematischen Funktionen können auch in geschlossener Form als Integrale dargestellt werden:
- $${\rm J}_n (x) = \frac{1}{2\pi}\cdot \int_{-\pi}^{+\pi} {{\rm e}^{\hspace{0.05cm}{\rm j}\hspace{0.05cm}\cdot \hspace{0.05cm}[\hspace{0.05cm}x \hspace{0.05cm}\cdot \hspace{0.05cm}\sin(\alpha) -\hspace{0.05cm} n \hspace{0.05cm}\cdot \hspace{0.05cm}\alpha \hspace{0.05cm}]}}\hspace{0.1cm}{\rm d}\alpha \hspace{0.05cm}.$$
Die Funktionen ${\rm J}_n (x)$ gehören zur Klasse der Besselfunktionen erster Art (englisch: Bessel Functions of the First Kind). Den Parameter $n$ nennt man die Ordnung.
Anmerkung: Es gibt eine Vielzahl von Abwandlungen der Besselfunktionen, unter anderem die mit ${\rm Y}_n (x)$ benannten Besselfunktionen zweiter Art. Für ganzzahliges $n$ lässt sich ${\rm Y}_n (x)$ durch ${\rm J}_n (x)$–Funktionen ausdrücken. In diesem Applet werden jedoch nur die Besselfunktionen erster Art ⇒ ${\rm J}_n (x)$ betrachtet.
Eigenschaften der Besselfunktionen
Die Anwendungen der Besselfunktionen in den Natur– und Ingenieurswissenschaften sind vielfältig. Wir beschränken uns hier auf einige Gebiete, die in unserem Lerntutorial $\rm LNTwww$ angesprochen werden.
Anwendungen der Besselfunktionen
For example:
Electromagnetic waves in a cylindrical waveguide Pressure amplitudes of inviscid rotational flows Heat conduction in a cylindrical object Modes of vibration of a thin circular (or annular) acoustic membrane (such as a drum or other membranophone) Diffusion problems on a lattice Solutions to the radial Schrödinger equation (in spherical and cylindrical coordinates) for a free particle Solving for patterns of acoustical radiation Frequency-dependent friction in circular pipelines Dynamics of floating bodies Angular resolution
Bessel functions also appear in other problems, such as signal processing (e.g., see FM synthesis, Kaiser window, or Bessel filter).
Sie spielen eine wichtige Rolle in der Physik. Man trifft u.a. bei der Untersuchung von Eigenschwingungen von zylindrischen Resonatoren, der Analyse des Frequenzspektrums von frequenzmodulierten Signalen und dem Sättigungsverhalten von Klystrons auf die Besselfunktionen. Man zählt die Besselfunktionen wegen ihrer vielfältigen Anwendungen in der mathematischen Physik zu den speziellen Funktionen. Die Funktionswerte entnimmt man in der Praxis aus Bessel-Funktionstabellen. ${\rm Q}(x)$ und $1/2\cdot {\rm erfc}(x)$, die für die Fehlerwahrscheinlichkeitsberechnung von großer Bedeutung sind.
- Sowohl die Abszisse als auch der Funktionswert kann entweder linear oder logarithmisch dargestellt werden.
- Für beide Funktionen wird jeweils eine obere Schranke (englisch: Upper Bound ) und eine untere Schranke (englisch: Lower Bound) angegeben.
Diese bereits 1844 von Friedrich Wilhelm Bessel eingeführten mathematischen Funktionen sind wie folgt definiert (erste Gleichung) und können gemäß der zweiten Gleichung durch eine Reihe angenähert werden:
Die nebenstehende Grafik zeigt die jeweils ersten drei Summanden $(k = 0, 1, 2)$ der Reihen ${\rm J}_0(η)$, ... , ${\rm J}_3(η).$ Der rot umrandete Term – gültig für $n = 3$ und $k = 2$ – lautet beispielsweise in ausgeschriebener Form:
- $$\frac{(-1)^2 \cdot (\eta/2)^{3 \hspace{0.05cm} + \hspace{0.05cm} 2 \hspace{0.02cm}\cdot \hspace{0.05cm}2}}{2\hspace{0.05cm}! \cdot (3+2)\hspace{0.05cm}!} = \frac{1}{240}\cdot (\frac{\eta}{2})^7 \hspace{0.05cm}.$$
- Die Besselfunktionen $J_n(η)$ findet man aber auch in Formelsammlungen oder mit Hilfe des von uns bereitgestellten Berechnungsmodul Besselfunktion erster Art und $n$–ter Ordnung.
- Sind die Funktionswerte für $n = 0$ und $n = 1$ bekannt, so können daraus die Besselfunktionen für $n ≥ 2$ iterativ ermittelt werden:
- $${\rm J}_n (\eta) ={2 \cdot (n-1)}/{\eta} \cdot {\rm J}_{n-1} (\eta) - {\rm J}_{n-2} (\eta) \hspace{0.05cm}.$$
Bei der Untersuchung digitaler Übertragungssysteme muss oft die Wahrscheinlichkeit bestimmt werden, dass eine (mittelwertfreie) gaußverteilte Zufallsgröße $x$ mit der Varianz $σ^2$ einen vorgegebenen Wert $x_0$ überschreitet. Für diese Wahrscheinlichkeit gilt:
- $${\rm Pr}(x > x_0)={\rm Q}(\frac{x_0}{\sigma}) = 1/2 \cdot {\rm erfc}(\frac{x_0}{\sqrt{2} \cdot \sigma}).$$
Die Funktion ${\rm Q}(x )$
Die Funktion ${\rm Q}(x)$ bezeichnet man als das Komplementäre Gaußsche Fehlerintegral. Es gilt folgende Berechnungsvorschrift:
- $${\rm Q}(x ) = \frac{1}{\sqrt{2\pi}}\int_{x}^{ +\infty}\hspace{-0.2cm}{\rm e}^{-u^{2}/\hspace{0.05cm} 2}\,{\rm d} u .$$
- Dieses Integral ist nicht analytisch lösbar und muss – wenn man dieses Applet nicht zur Verfügung hat – aus Tabellen entnommen werden.
- Speziell für größere $x$–Werte von (also für kleine Fehlerwahrscheinlichkeiten) liefern die nachfolgend angegebenen Schranken eine brauchbare Abschätzung für das Komplementäre Gaußsche Fehlerintegral, die auch ohne Tabellen berechnet werden können.
- Eine obere Schranke (englisch: Upper Bound ) des Komplementären Gaußschen Fehlerintegrals lautet:
- $${\rm Q}_{\rm UB}(x )=\text{Upper Bound }\big [{\rm Q}(x ) \big ] = \frac{ 1}{\sqrt{2\pi}\cdot x}\cdot {\rm e}^{- x^{2}/\hspace{0.05cm}2} > {\rm Q}(x).$$
- Entsprechend gilt für die untere Schranke (englisch: Lower Bound ):
- $${\rm Q}_{\rm LB}(x )=\text{Lower Bound }\big [{\rm Q}(x ) \big ] =\frac{1-1/x^2}{\sqrt{2\pi}\cdot x}\cdot {\rm e}^{-x^ 2/\hspace{0.05cm}2} ={\rm Q}_{\rm UB}(x ) \cdot (1-1/x^2)< {\rm Q}(x).$$
In vielen Programmbibliotheken findet man allerdings die Funktion ${\rm Q}(x )$ nicht.
Die Funktion $1/2 \cdot {\rm erfc}(x )$
In fast allen Programmbibliotheken findet man dagegen die Komplementäre Gaußsche Fehlerintegral (englisch: Complementary Gaussian Error Function)
- $${\rm erfc}(x) = \frac{2}{\sqrt{\pi}}\int_{x}^{ +\infty}\hspace{-0.2cm}{\rm e}^{-u^{2}}\,{\rm d} u .$$
die mit ${\rm Q}(x)$ wie folgt zusammenhängt: ${\rm Q}(x)=1/2\cdot {\rm erfc}(x/{\sqrt{2}}).$ Da bei fast allen Anwendungen diese Funktion mit dem Faktor $1/2$ verwendet wird, wurde in diesem Applet genau diese Funktion realisiert:
- $$1/2 \cdot{\rm erfc}(x) = \frac{1}{\sqrt{\pi}}\int_{x}^{ +\infty}\hspace{-0.2cm}{\rm e}^{-u^{2}}\,{\rm d} u .$$
- Auch für diese Funktion kann wieder eine obere und eine untere Schranke angegeben werden:
- $$\text{Upper Bound }\big [1/2 \cdot{\rm erfc}(x) \big ] = \frac{ 1}{\sqrt{\pi}\cdot 2x}\cdot {\rm e}^{- x^{2}} ,$$
- $$\text{Lower Bound }\big [1/2 \cdot{\rm erfc}(x) \big ] = \frac{ {1-1/(2x^2)}}{\sqrt{\pi}\cdot 2x}\cdot {\rm e}^{- x^{2}} .$$
Wann bietet welche Funktion Vorteile?
$\text{Beispiel 1:}$ Wir betrachten die binäre Basisbandübertragung. Hier lautet die Bitfehlerwahrscheinlichkeit $p_{\rm B} = {\rm Q}({s_0}/{\sigma_d})$, wobei das Nutzsignal die Werte $\pm s_0$ annehmen kann und der Rauscheffektivwert $\sigma_d$ ist.
Es wird vorausgesetzt, dass Tabellen zur Verfügung stehen, in denen das Argument der Gaußschen Fehlerfunktionen im Abstand $0.1$ aufgelistet sind. Mit $s_0/\sigma_d = 4$ erhält man für die Bitfehlerwahrscheinlichkeit gemäß der Q–Funktion:
- $$p_{\rm B} = {\rm Q} (4) \approx 0.317 \cdot 10^{-4}\hspace{0.05cm}.$$
Nach der zweiten Gleichung ergibt sich:
- $$p_{\rm B} = {1}/{2} \cdot {\rm erfc} ( {4}/{\sqrt{2} })= {1}/{2} \cdot {\rm erfc} ( 2.828)\approx {1}/{2} \cdot {\rm erfc} ( 2.8)= 0.375 \cdot 10^{-4}\hspace{0.05cm}.$$
- Richtiger ist der erste Wert. Bei der zweiten Berechnungsart muss man runden oder – noch besser – interpolieren, was aufgrund der starken Nichtlinearität dieser Funktion sehr schwierig ist.
- Bei den gegebenen Zahlenwerten ist demnach Q–Funktion besser geeignet. Außerhalb von Übungsbeispielen wird $s_0/\sigma_d$ in der Regel einen „krummen” Wert besitzen. In diesem Fall bietet ${\rm Q}(x)$ natürlich keinen Vorteil gegenüber $1/2 \cdot{\rm erfc}(x)$.
$\text{Beispiel 2:}$ Mit der Energie pro Bit $(E_{\rm B})$ und der Rauschleistungsdichte $(N_0)$ gilt für die Bitfehlerwahrscheinlichkeit von Binary Phase Shift Keying (BPSK):
- $$p_{\rm B} = {\rm Q} \left ( \sqrt{ {2 E_{\rm B} }/{N_0} }\right ) = {1}/{2} \cdot {\rm erfc} \left ( \sqrt{ {E_{\rm B} }/{N_0} }\right ) \hspace{0.05cm}.$$
Für die Zahlenwerte $E_{\rm B} = 16 \ \rm mWs$ und $N_0 = 16 \ \rm mW/Hz$ erhält man:
- $$p_{\rm B} = {\rm Q} \left (4 \cdot \sqrt{ 2} \right ) = {1}/{2} \cdot {\rm erfc} \left ( 4\right ) \hspace{0.05cm}.$$
- Der erste Weg führt zum Ergebnis $p_{\rm B} = {\rm Q} (5.657) \approx {\rm Q} (5.7) = 0.6 \cdot 10^{-8}\hspace{0.05cm}$, während $1/2 \cdot{\rm erfc}(x)$ hier den richtigeren Wert $p_{\rm B} \approx 0.771 \cdot 10^{-8}$ liefert.
- Wie im ersten Beispiel erkennt man aber auch hier: Die Funktionen ${\rm Q}(x)$ und $1/2 \cdot{\rm erfc}(x)$ sind grundsätzlich gleich gut geeignet. Vor– oder Nachteile der einen oder anderen Funktion ergeben sich nur bei konkreten Zahlenwerten.
Zur Handhabung des Applets
(A) Vorauswahl für blauen Parametersatz
(B) Parametereingabe $I$ und $p$ per Slider
(C) Vorauswahl für roten Parametersatz
(D) Parametereingabe $\lambda$ per Slider
(E) Graphische Darstellung der Verteilungen
(F) Momentenausgabe für blauen Parametersatz
(G) Momentenausgabe für roten Parametersatz
(H) Variation der grafischen Darstellung
$\hspace{1.5cm}$„$+$” (Vergrößern),
$\hspace{1.5cm}$ „$-$” (Verkleinern)
$\hspace{1.5cm}$ „$\rm o$” (Zurücksetzen)
$\hspace{1.5cm}$ „$\leftarrow$” (Verschieben nach links), usw.
( I ) Ausgabe von ${\rm Pr} (z = \mu)$ und ${\rm Pr} (z \le \mu)$
(J) Bereich für die Versuchsdurchführung
Andere Möglichkeiten zur Variation der grafischen Darstellung:
- Gedrückte Shifttaste und Scrollen: Zoomen im Koordinatensystem,
- Gedrückte Shifttaste und linke Maustaste: Verschieben des Koordinatensystems.
Über die Autoren
Dieses interaktive Berechnungstool wurde am Lehrstuhl für Nachrichtentechnik der Technischen Universität München konzipiert und realisiert.
- Die erste Version wurde 2006 von Markus Elsberger und Slim Lamine im Rahmen von Abschlussarbeiten mit „FlashMX–Actionscript” erstellt (Betreuer: Günter Söder).
- 2018 wurde das Programm von Xiaohan Liu (Bachelorarbeit, Betreuer: Tasnád Kernetzky ) auf „HTML5” umgesetzt.
