Hyperexponentialverteilung

Die durchgezogene, blaue Linie zeigt die Wahrscheinlichkeitsdichte einer Hyperexponentialverteilung am Beispiel p₁=0.9, p₂=0.1, λ₁=1 und λ₂=20.

Die Hyperexponentialverteilung ist eine stetige Wahrscheinlichkeitsverteilung. Anschaulich gesprochen ist sie eine Überlagerung mehrerer Exponentialverteilungen.

Definition

Seien $Y_i$ (mit $i=1,\dots ,N$) unabhängige, exponentialverteilte Zufallsvariablen mit Raten ${\lambda }_i$ und seien $p_i$ Wahrscheinlichkeiten, deren Summe 1 ergibt. Dann heißt die Zufallsvariable $X$ hyperexponentialverteilt, wenn sie folgende Wahrscheinlichkeitsdichte besitzt:

$f_X(x)=\sum _{i=1}^N{p}_i{f}_{Y_i}(x)=\{\begin{array}{cc}\sum _{i=1}^N{p}_i{\lambda }_i{e}^{-{\lambda }_{i}x}& x\ge 0,\\ 0& x<0.\end{array}$

Einordnung und Bemerkungen

Bei einer Exponentialverteilung ist der Variationskoeffizient (Standardabweichung geteilt durch Erwartungswert) gleich 1. Die Bezeichnung „hyper“-exponential rührt daher, dass der Variationskoeffizient hier größer als 1 ist (sofern verschiedene ${\lambda }_i$ auftreten). Im Unterschied dazu ist er bei der Hypoexponentialverteilung kleiner als 1. Während die Exponentialverteilung das stetige Analogon zur geometrischen Verteilung ist, ist die Hyperexponentialverteilung kein Analogon zur hypergeometrischen Verteilung. Die Hyperexponentialverteilung ist ein Beispiel für eine Mischverteilung.

Als Anwendungsbeispiel kann die Auslastung eines Internetanschlusses dienen, über welchen entweder (mit Wahrscheinlichkeit $p$ und Rate ${\lambda }_1$) Internettelefonie oder (mit Wahrscheinlichkeit $q$ und Rate ${\lambda }_2$) Dateiübertragungen laufen, wobei $p+q=1$. Die Gesamtauslastung ist dann hyperexponentialverteilt.

Eine gegebene Wahrscheinlichkeitsverteilung, inklusive endlastiger Verteilungen, kann durch eine Hyperexponentialverteilung angenähert werden, indem rekursiv verschiedene Zeitskalen (${\lambda }_i$) mittels der sogenannten Prony-Methode angefittet werden.

Eigenschaften

Aus der Linearität des Integrals ergibt sich:

$\mathrm{E}[X]=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}xf(x)dx=\sum _{i=1}^N{p}_i\underset{0}{\overset{\mathrm{\infty }}{\int }}x{\lambda }_i{e}^{-{\lambda }_{i}x}dx=\sum _{i=1}^N\frac{p_i}{{\lambda }_i}$

und

$\mathrm{E}\left[X^2\right]=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{x}^{2}f(x)dx=\sum _{i=1}^N{p}_i\underset{0}{\overset{\mathrm{\infty }}{\int }}{x}^2{\lambda }_i{e}^{-{\lambda }_{i}x}dx=\sum _{i=1}^N\frac{2}{{\lambda }_i^2}{p}_i.$

Mit Hilfe des Verschiebungssatzes ergibt sich daraus die Varianz:

$\mathrm{Var}[X]=\mathrm{E}\left[X^2\right]-\mathrm{E}{\left[X\right]}^2=\sum _{i=1}^N\frac{2}{{\lambda }_i^2}{p}_i-{\left[\sum _{i=1}^N\frac{p_i}{{\lambda }_i}\right]}^2={\left[\sum _{i=1}^N\frac{p_i}{{\lambda }_i}\right]}^2+\sum _{i=1}^N\sum _{j=1}^N{p}_i{p}_j{(\frac{1}{{\lambda }_i}-\frac{1}{{\lambda }_j})}^2.$

Sofern nicht alle ${\lambda }_i$ gleich groß sind, ist die Standardabweichung größer als der Erwartungswert.

Die momenterzeugende Funktion ist

$\mathrm{E}\left[e^{tX}\right]=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{tx}f(x)dx=\sum _{i=1}^N{p}_i\underset{0}{\overset{\mathrm{\infty }}{\int }}{e}^{tx}{\lambda }_i{e}^{-{\lambda }_{i}x}dx=\sum _{i=1}^N\frac{{\lambda }_i}{{\lambda }_i-t}{p}_i.$

Siehe auch

• Phasenverteilung