Ausblendeeigenschaft

Die Delta-Distribution (auch δ-Funktion; Dirac-Funktion, -Impuls, -Puls, -Stoß (nach Paul Dirac); Stoßfunktion; sowie Einheitsimpulsfunktion genannt) wird in der Naturwissenschaft durch ein kleines Delta δ dargestellt und symbolisiert eine spezielle irreguläre Distribution, die in der Mathematik und Physik von grundlegender Bedeutung ist.

Inhaltsverzeichnis

Definition

Sei \phi \in C_c^\infty(\R^n) eine glatte Funktion und sei

φε(x) = ε nφ(ε − 1x)

eine Funktionenfolge in ε. Diese trägt den Namen Dirac-Folge. So ist die Delta-Distribution für alle Testfunktionen \psi \in C^\infty(\R^n) \cong \mathcal{E}(\R^n) definiert als

\begin{align}
\delta\colon\,\mathcal{E}&\to\mathbb{R}\\
\delta(\psi) &:= \int_{\R^n} \delta(x)\psi(x)\mathrm{d}x := \lim_{\epsilon \to \infty} \int_{\R^n} \phi_\epsilon (x) \psi(x) \mathrm{d} x = \psi(0)
\end{align}.

Der erste Integralausdruck ist nur symbolisch zu verstehen. Es gibt keine Funktion δ, welche der obigen Definition genügt. Denn man kann im zweiten Integralausdruck den Grenzwert nicht in das Integral hereinziehen.

Definition über Dirac-Maß

Das durch ein positives Radon-Maß μ erzeugte Funktional \langle\mu,f\rangle=\int f(x)\,\mathrm{d}\mu (für f\in\mathcal{D}) ist eine Distribution.

Die Delta-Distribution wird von folgenden Radon-Maß - man spricht hier speziell vom Diracmaß - erzeugt:

\delta(A)=\begin{cases}
1\  & \text{falls }0\in A\\
0\  & \text{sonst}\end{cases}\ ,\quad A\subset\mathbb{R}

Ein Maß lässt sich physikalisch interpretieren, z.B. als Massendichte oder Ladungsdichte des Raums. Dann entspricht die Delta-Distribution einem Massenpunkt der Masse 1 oder einer Punktladung der Ladung 1 im Ursprung.

\langle\delta,f\rangle=\int f(x)\,\mathrm{d}\delta=f(0)

Befinden sich an den Stellen x_i\in\mathbb{R} Punktladungen qi, wobei die Summe über alle Ladungen endlich bleibt, dann wird für A\subset\mathbb{R} ein Maß auf der σ-Algebra aller Teilmengen von \mathbb{R} definiert, die der Ladungsdichte entspricht (iA = Summe über alle i mit x_{i}\in A):

\rho(A):=\sum_{i_{A}}q_{i}

Für dieses Maß ist dann die zugehörige Distribution:

\langle\rho,f\rangle=\int f(x)\,\mathrm{d}\rho=\sum_{i_{A}}f(x_{i})q_{i}


Alternative Schreibweisen

Der Wert, den die Delta-Distribution nach Anwendung auf eine Testfunktion f\in\mathcal{D} liefert, schreibt man auch als

\delta(f) = \langle\delta,f\rangle = f(0)

bzw. formal auch als

\delta(f) = \int_{-\infty}^{\infty}\delta(x)\, f(x)\,\mathrm{d}x=f(0)

Diese Schreibweise ist eigentlich nicht richtig, weil die Delta-Distribution eine irreguläre (= singuläre) Distribution ist, d. h. sie lässt sich nicht durch eine lokal integrierbare Funktion in obiger Weise darstellen. Bei Verwendung der Integral-Schreibweise ist zu beachten, dass es sich nicht um ein Riemann-Integral oder Lebesgue-Integral, sondern um die Auswertung des Funktionals δ an der Stelle f, also δ(f) = f(0), handelt.

Ableitungen

Die Delta-Distribution kann wie jede Distribution beliebig oft distributiv differenziert werden:

\langle\delta',f\rangle=-\langle\delta,f'\rangle=-f'(0)

Und die n-te distributive Ableitung:

\langle\delta^{(n)},f\rangle=(-1)^{n}\langle\delta,f^{(n)}\rangle=(-1)^{n}f^{(n)}(0)

Die Heaviside-Funktion Θ(x) ist nicht stetig differenzierbar, aber die distributive Ableitung existiert, diese ist nämlich die Delta-Distribution:

\langle\Theta',f\rangle=-\langle\Theta,f'\rangle=-\int_{-\infty}^{\infty}\Theta(x)\, f'(x)\,\mathrm{d}x=-\int_{0}^{\infty}f'(x)\,\mathrm{d}x=-\underbrace{f(\infty)}_{=0}+f(0)=\langle\delta,f\rangle

Da f\in\mathcal{D} muss die Funktion im Unendlichen verschwinden.

Irregularität

Die Irregularität der Delta-Distribution lässt sich mit einem Widerspruchsbeweis zeigen:

Angenommen δ wäre regulär, dann gäbe es eine lokal integrierbare Funktion \delta(x)\in L^{1}_{lok}, also eine Funktion, die über jeden kompakten Intervall [a,b] bzgl. des Lebesgue-Maßes integrierbar ist

\int_{a}^{b}|\delta(x)|\mathrm{d}x<\infty

so dass für alle Testfunktionen f(x) gilt:

\langle\delta,f\rangle = \int_{-\infty}^{\infty}\delta(x)\, f(x)\,\mathrm{d}x=f(0)

Insbesondere muss dies für obige Testfunktion φb(x) gelten. Die Wirkung der Delta-Distribution auf diese ist:

\langle\delta,\phi_{b}\rangle = \phi_{b}(0)=\exp(-1)=\mathrm{const}_{\,b}

Mit der angenommenen regulären Distribution

\langle\delta,\phi_{b}\rangle=\int_{-\infty}^{\infty}\delta(x)\,\phi_{b}(x)\,\mathrm{d}x=\int_{-b}^{b}\delta(x)\,\phi_{b}(x)\,\mathrm{d}x

lässt sich folgende Abschätzung durchführen:

\phi_{b}(0)=|\langle\delta,\phi_{b}\rangle|=|\int_{-b}^{b}\delta(x)\,\phi_{b}(x)\,\mathrm{d}x|\leq\underbrace{\|\phi_{b}(x)\|_{\infty}}_{\phi_{b}(0)}\,\int_{-b}^{b}|\delta(x)|\,\mathrm{d}x\underset{(b<b_{c})}{<}\phi_{b}(0)

Weil \delta(x)\in L^{1}_{lok} wird das Integral \int_{-b}^{b}|\delta(x)|\,\mathrm{d}x für b < bc (wobei bc ein von der Funktion δ(x) abhängiger kritischer Wert ist) kleiner 1 (und konvergiert gegen 0 für b gegen 0). Man erhält φb(0) < φb(0), also einen Widerspruch; somit ist die Delta-Distribution nicht durch eine lokal integrierbare Funktion darstellbar. Der Widerspruch ergibt sich, weil die Menge {0} für das Lebesgue-Maß vernachlässigbar ist, nicht aber für das Dirac-Maß.

Praktische Definition

Insbesondere Physiker verwenden gerne die folgende praktische Definition, die allerdings mathematisch nicht ganz sauber ist, da ein einzelner Punkt das Lebesgue-Maß null hat und damit beide Bedingungen nicht unter einen Hut zu bringen sind:

\delta (x)=\begin{cases} 0 &amp;amp;amp; x\ne 0\\\infty &amp;amp;amp; x=0\end{cases}   und   \int_{-\infty}^\infty \delta(x) \, \mathrm dx = 1.

Man kann sich leicht überzeugen, dass es keine gewöhnliche Funktion gibt, die die obigen Bedingungen erfüllt. Deshalb ist eine exakte Definition nur im Rahmen der Theorie der Distributionen möglich (man spricht auch von Verallgemeinerten Funktionen). Die Delta-Distribution kann nicht wie eine gewöhnliche Funktion behandelt werden (z.B. kann man sie nicht quadrieren).

Anschaulich stellt man sich die Delta-Distribution meistens als eine beliebig hohe und beliebig schmale Funktion vor, deren Fläche den Wert 1 besitzt. Man lässt nun die Funktion immer schmaler und dafür immer höher werden - die Fläche darunter muss konstant 1 bleiben. Am Ende dieses Gedankenexperiments erhält man einen Graphen, den man wegen der unendlichen Amplitude nicht mehr zeichnen kann. Es existieren auch mehrdimensionale Dirac-Distributionen, diese werden anschaulich zu mehrdimensionalen „Keulen“ mit dem Volumen 1.

Für verschwindende Breite und zunehmende Höhe, bei Erhaltung des Flächeninhalts, gewinnt man aus dem Rechteckimpuls - und auch aus geeignet geglätteten, beliebig-oft differenzierbar gemachten Rechteckimpulsen - die Delta-Distribution einschließlich ihrer sämtlichen Ableitungen. Man muss sich klarmachen, dass viele Approximationen der Delta-Distributionen durch glatte Funktionen möglich sind und alle zum selben Resultat führen.

Verbesserte Heuristik

Standardnormalverteilung  \delta_{a}(x)=\frac {1}{\sqrt{2\pi a}} \cdot e^{-\frac {x^2}{2a}}.
Für a\to 0 wird die Funktion immer höher und schmaler, der Flächeninhalt, 1, bleibt jedoch unverändert. (Unterbrechung mit Esc-Taste.)

Statt der unsauberen praktischen Definition ist es sinnvoller, δ durch Limes-Bildungen mit einer Folge lokal integrierbarer Funktionen δε zu definieren (ε > 0), die folgende Eigenschaften besitzen:

  • \delta_{\epsilon}(-x)=\delta_{\epsilon}\,(x)
  • \int_{-\infty}^{\infty}\delta_{\epsilon}(x)\,\mathrm{d}x=1
  • \lim_{\epsilon\to 0}\delta_{\epsilon}(x)=\begin{cases} 0 &amp;amp;amp; x\ne 0\\\infty &amp;amp;amp; x=0\end{cases}

Mit den Funktionen δε(x) sind nun reguläre Distributionen δε gegeben, also Distributionen, die sich als Integral mit Kern δε(x) schreiben lassen:

\langle\delta_{\epsilon},f\rangle=\int_{-\infty}^{\infty}\delta_{\epsilon}(x)\, f(x)\,\mathrm{d}x

Nur im Limes \epsilon\to 0 erhält man das ungewöhnliche Verhalten der Delta-Distribution, wobei zu beachten ist, dass die Limes-Bildung nicht unter dem Integral, sondern davor erfolgt.

\lim_{\epsilon\to 0}\langle\delta_{\epsilon},f\rangle=f(0)=\langle\delta,f\rangle

Die Ableitungen der regulären Distributionen können mittels partieller Integration berechnet werden (hier exemplarisch für erste Ableitung, analog für höhere)

\begin{align}
\langle\delta_{\epsilon}^{\prime},f\rangle &amp;amp;amp;= \int_{-\infty}^{\infty}\delta_{\epsilon}^{\prime}(x)\, f(x)\,\mathrm{d}x \\
 &amp;amp;amp;= \underbrace{\left[\delta_{\epsilon}(x)\, f(x)\right]_{-\infty}^{\infty}}_{=0}-\int_{-\infty}^{\infty}\delta_{\epsilon}(x)\, f^{\prime}(x)\,\mathrm{d}x \\
 &amp;amp;amp;= -\int_{-\infty}^{\infty}\delta_{\epsilon}(x)\, f^{\prime}(x)\,\mathrm{d}x \\
 &amp;amp;amp;= -\langle\delta_{\epsilon},f^{\prime}\rangle
\end{align}

und ergeben im Limes \epsilon\to 0 das Verhalten der distributiven Ableitung:

\lim_{\epsilon\to0}\langle\delta_{\epsilon}^{\prime},f\rangle=-f^{\prime}(0)=\langle\delta^{\prime},f\rangle

Beispiele für Approximationen

Im Folgenden werden verschiedene Approximationen δε(x) angegeben, zunächst stetig differenzierbare:

\delta_{\epsilon}(x)=\frac{1}{\sqrt{2\pi\varepsilon}}\,\exp\left(-\frac{x^{2}}{2\varepsilon}\right)
Die angegebenen Funktionen sind Glockenfunktionen (Normalverteilungen) mit einem sehr schmalen und sehr hohen Maximum bei x=0, die Breite ist ~ \sqrt{\epsilon}\to 0 und die Höhe ~ 1/\sqrt{\epsilon}\to\infty. Der Flächeninhalt unter der Funktion ist aber immer 1, für alle ε.
\delta_{\epsilon}(x)=\frac{1}{\pi}\frac{\epsilon}{x^{2}+\epsilon^{2}}
\delta_{\epsilon}(x)=\frac{1}{\pi x}\sin\left(\frac{x}{\epsilon}\right)
\delta_{\epsilon}(x)=\frac{1}{\sqrt{i\pi\epsilon}}\,\exp\left(\frac{ix^{2}}{\epsilon}\right)
die man sich vorstellen kann, als eine Linie, die auf einen Zylinder gewickelt ist, und deren Wicklungen durch das x2 immer enger werden; die Grundfläche (in x-y-Ausrichtung) des Zylinders wird aus dem Imaginär- und Realteil der Funktion gebildet, die Funktion entwickelt sich dann in z-Richtung.

Es sind aber auch Approximationen möglich, die nur stückweise stetig differenzierbar sind:

\delta_{\epsilon}(x)=\frac{\textrm{rect}(x/\epsilon)}{\epsilon}=\begin{cases}
\frac{1}{\epsilon} &amp;amp;amp; \ |x|\leq\frac{\epsilon}{2}\\
0 &amp;amp;amp; \ \text{sonst}\end{cases}
\delta_{\epsilon}(x)=\begin{cases}
\frac{\epsilon+x}{\epsilon^{2}} &amp;amp;amp; \ -\epsilon\leq x\leq0\\
\frac{\epsilon-x}{\epsilon^{2}} &amp;amp;amp; \ 0&amp;amp;lt;x\leq\epsilon\\
0 &amp;amp;amp; \ \text{sonst}\end{cases}
\delta_{\epsilon}(x)=\frac{1}{2\epsilon}\exp\left(-\frac{|x|}{\epsilon}\right)

Eigenschaften

  • Definierende Eigenschaft der Delta-Distribution: Faltungseigenschaft, auch Ausblendeigenschaft, Siebeigenschaft genannt
\langle\delta,f\rangle=\int_{-\infty}^{\infty}\delta(x)\, f(x)\,\mathrm{d}x=f(0)
bzw. mit den Eigenschaften Translation und Skalierung (siehe unten) folgt:
\int_{- \infty}^\infty f(x)\,\delta (x-a)\,\mathrm{d}x=\int_{- \infty}^\infty f(x)\,\delta (a-x)\,\mathrm{d}x=f(a),
speziell für den Fall der konstanten Funktion 1:
\int_{- \infty}^\infty \delta (x-a)\,\mathrm{d}x=1
  • Linearität:
\langle\delta,f+g\rangle=\langle\delta,f\rangle+\langle\delta,g\rangle=f(0)+g(0)
  • Translation:
\langle\delta(\cdot-a),f\rangle=\langle\delta,f(\cdot+a)\rangle=f(a)
für \delta(\cdot-a) ist auch die Bezeichnung δa gebräuchlich.
  • Skalierung:
\langle\delta(a\cdot),f\rangle=\frac{1}{|a|}\langle\delta,f(\tfrac{\cdot}{a})\rangle=\frac{1}{|a|}f(0)
 \delta(\alpha x) = \frac{1}{|\alpha|} \delta(x)
d. h. die Delta-Distribution ist homogen vom Grad −1.
  • Dimension
Eine direkte Folgerung aus der Skalierungseigenschaft ist die Dimension bzw. Maßeinheit der Delta-Distribution. Sie entspricht genau der reziproken Dimension ihres Arguments. Hat x beispielsweise die Dimension einer Länge, so hat δ(x) die Dimension (1/Länge).
  • Hintereinanderausführung:
\int_{-\infty}^{\infty}\phi(x)\,\delta(g(x))\,\mathrm{d}x=\sum_{i=1}^{n}\int_{-\infty}^{\infty}\phi(x)\frac{\delta(x-x_{i})}{|g'(x_{i})|}\,\mathrm{d}x=\sum_{i=1}^{n}\frac{\phi(x_{i})}{|g'(x_{i})|},
\delta(g(x))=\sum_{i=1}^{n}\frac{\delta(x-x_{i})}{|g^{\prime}(x_{i})|}
wobei xi die einfachen Nullstellen von g(x) sind (sofern g(x) nur endlich viele und nur einfache Nullstellen hat).
  • damit folgt als ein Spezialfall die Rechenregel
 \delta(x^2 - \alpha^2) = \frac{1}{2|\alpha|} [\delta(x - \alpha)+ \delta(x + \alpha) ],
 \mathcal{F} \{ \delta(x) \}(k) = \int_{-\infty}^\infty e^{-ikx}\,\delta(x)\,\mathrm{d}x = 1
Anschaulich bedeutet das, dass in der Delta-Distribution alle Frequenzen enthalten sind, und zwar mit gleicher Stärke.
Die Fourierdarstellung \delta(x)=\mathcal{F}^{-1}\{1\} bzw. \delta(x-a)=\mathcal{F}^{-1}\{e^{-ika}\} ist eine in der Physik wichtige Darstellung der Delta-Distribution:
\delta(x-a) = \frac{1}{2\pi}\int_{-\infty}^\infty e^{ik(x-a)}\,\mathrm{d}k = \int_{-\infty}^\infty e^{i 2\pi\tilde{\nu}(x-a)}\,\mathrm{d}\tilde{\nu}
\mathcal{L}\left\{ \delta(x)\right\} (s)=\int_{0}^{\infty}\mathrm{e}^{-sx}\,\delta(x)\,\mathrm{d}x=1

Praktische Anwendung

Praktische Bedeutung hat der Dirac-Stoß bei der Ermittlung der Impulsantwort in der Akustik (in anderen Sparten der Physik spricht man auch von einer δ-Größe, wenn man meint, dass die betreffende Größe einer schmalst-möglichen Verteilung genügt). So hat jeder Raum ein eigenes Schallverhalten. Mit einem Dirac-Impuls (angenähert durch ein Klatschen mit den Händen) kann dieses Verhalten (durch Messen des „Echos“, also der Systemantwort) ermittelt werden.

Typische, technisch realisierbare Dirac-Werte:

Eine wichtige Anwendung der Delta-Distribution ist die Lösung inhomogener linearer Differentialgleichungen mit der Methode der Greenschen Funktion.

Mehrdimensionale Delta-Distribution

Im Mehrdimensionalen ist der Raum der Testfunktionen \mathcal{D} gleich C_{0}^{\infty}(\mathbb{R}^{n}), der Raum der unendlich oft stetig partiell differenzierbaren Funktionen f\colon\,\mathbb{R}^{n}\to\mathbb{R} mit kompaktem Träger, versehen mit einem sehr starken Konvergenzbegriff, nämlich der gleichmäßigen Konvergenz nicht nur der Funktion, sondern aller partiellen Ableitungen.

Die Delta-Distribution hat auf die Testfunktion f\colon\,\mathbb{R}^{n}\to\mathbb{R} die folgende Wirkung:

\delta\colon\,\mathcal{D}\to\mathbb{R}\,,\ f\mapsto f(\vec{0})

In der Integralschreibweise unter Verwendung von Translation und Skalierung:

\int f(\vec{x})\,\delta(\vec{x}-\vec{a})\,\mathrm{d}^{n}x=\int f(\vec{x})\,\delta(\vec{a}-\vec{x})\,\mathrm{d}^{n}x=f(\vec{a})

Die "mehrdimensionale" Delta-Distribution lässt sich als Produkt von "eindimensionalen" Delta-Distributionen schreiben:

\delta(\vec{x}-\vec{a})=\delta(x_{1}-a_{1})\,\delta(x_{2}-a_{2})\, ... \,\delta(x_{n}-a_{n})

Speziell im Dreidimensionalen gibt es eine Darstellung der Delta-Distribution, die häufig in der Elektrodynamik eingesetzt wird um Punktladungen darzustellen:

\delta(\vec{x}-\vec{a})=-\frac{1}{4\pi}\Delta\frac{1}{\|\vec{x}-\vec{a}\|_{2}}

Literatur

  • W. Walter: Einführung in die Theorie der Distributionen, BI-Wissenschaftsverlag.
  • G. Friedlander, M. Joshi: Introduction to the Theory of Distributions, Cambridge University Press.

Wikimedia Foundation.

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”