Spektraldarstellung stationärer stochastischer Prozesse

Die Spektraldarstellung stationärer stochastischer Prozesse ist in gewissem Sinne ein Analogon zur Fourierreihenentwicklung einer Funktion. Jede beschränkte stetige Funktion kann als additive Überlagerung harmonischer Schwingungen dargestellt werden. Auch ein stationärer stochastischer Prozess kann dargestellt werden als additive Überlagerung harmonischer Schwingungen, allerdings mit zufälliger Amplitude. Die Spektraldarstellung eines stationären Prozesses bietet in der Regel tiefere Einblicke in die Struktur des Prozesses, insbesondere wenn es sich um eine Mischung verschiedener periodischer Anteile handelt.

Mathematische Beschreibung

Sei $\mathbb{Z}$ die Menge der ganzen Zahlen und $X_t;t\in \mathbb{Z}$ ein zeitdiskreter stationärer stochastischer Prozess mit Erwartungswert $\mathrm{E}{X}_t=0$ und Kovarianzfunktion $\mathrm{Cov}(X_t,X_{t+h})=\mathrm{\Gamma }(t,t+h)$, die wegen der Stationarität nur von der Differenz der Zeitpunkte abhängt, also nur die Funktion einer Variablen ist:$\mathrm{\Gamma }(t,t+h)=:\gamma (h)$.

Spektraldarstellung von X_t

Jeder stationäre Prozess $X_t$ mit $\mathrm{E}{X}_t=0$ hat die sogenannte Spektraldarstellung

$X_t=\underset{-\pi }{\overset{\pi }{\int }}{e}^{it\lambda }dZ(\lambda )$.

Dies ist ein stochastisches Integral, und zwar bzgl. eines Prozesses $Z(\lambda )$ mit unkorrelierten Zuwächsen, d. h. für ${\lambda }_1<{\lambda }_2<{\lambda }_3<{\lambda }_4$ sind die Zuwächse $Y({\lambda }_2)-Y({\lambda }_1)$ und $Y({\lambda }_4)-Y({\lambda }_3)$ unkorreliert.

Wenn $Z(\lambda )$ nur endlich viele Zuwächse hat, z. B. $n$ Zuwächse $A({\lambda }_1),\dots ,A({\lambda }_n)$ bei $-\pi <{\lambda }_1<\dots <{\lambda }_n<\pi$, dann kann obiges Integral als Summe geschrieben werden:

$X_t=\sum _{k=1}^{n}A({\lambda }_k)e^{it{\lambda }_k}$.

Jeder Summand ist eine harmonische Schwingung mit Frequenz ${\lambda }_k$ und der zufälligen Amplitude $A({\lambda }_k)$.

Spektraldarstellung der Kovarianzfunktion

Die Kovarianzfunktion $\gamma (h)$ ist eine symmetrische und positiv semidefinite Funktion und hat damit nach dem Satz von Bochner (in diskreter Variante als Satz von Herglotz bezeichnet) die Darstellung

$\gamma (h)=\underset{-\pi }{\overset{\pi }{\int }}{e}^{ih\lambda }dF(\lambda )$.

Dabei heißt $F$ Spektralverteilungsfunktion. Sie ist auf $[-\pi ,\pi ]$ monoton nicht fallend und es gilt $F(-\pi )=0;F(\pi )=\gamma (0)=\mathrm{Var}{X}_t$. Die Beziehung

$\mathrm{E}|Z({\lambda }_2)-Z({\lambda }_1){|}^2=F({\lambda }_2)-F({\lambda }_1);{\lambda }_2\ge {\lambda }_1$

stellt die Verbindung zwischen der Spektraldarstellung von $X_t$ und der Spektraldarstellung von $\gamma (h)$ dar.

Spektraldichte

Wenn $\sum _{h=-\mathrm{\infty }}^{\mathrm{\infty }}|\gamma (h)|<\mathrm{\infty }$, dann kann die Spektraldarstellung von $\gamma (h)$ als Riemannsches Integral geschrieben werden:

$\gamma (h)=\underset{-\pi }{\overset{\pi }{\int }}f(\lambda )e^{ih\lambda }d\lambda$.

Die Funktion $f$ heißt Spektraldichte von $X_t$. Anschaulich gesprochen gibt $f(\lambda )$ an, mit welcher Intensität die Frequenz $\lambda$ im Spektrum von $X_t$ vorkommt. Die Spektraldichte selbst hat die Darstellung

$f(\lambda )=\frac{1}{2\pi }\sum _{h=-\mathrm{\infty }}^{\mathrm{\infty }}{e}^{-ih\lambda }\gamma (h)$.

$f$ ist also die Fouriertransformierte von $\gamma$, bzw. $\gamma$ ist die inverse Fouriertransformierte von $f$. Für $h=0$ gilt speziell

$\gamma (0)=\mathrm{Var}{X}_t=\underset{-\pi }{\overset{\pi }{\int }}f(\lambda )d\lambda$.

Dies kann als Streuungszerlegung (signaltechnisch Leistungsverteilung) auf die verschiedenen Frequenzen $\lambda$ interpretiert werden.

Zeitstetiger Fall

Sei nun $X_t$ ein stationärer Prozess mit reellwertigem $t\in ℝ$. Dann modifizieren sich obige Formeln zu:

$X_t=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{i\lambda t}dZ(\lambda );\gamma (h)=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{i\lambda h}dF(\lambda )$.

Dabei ist $Z(\lambda )$ wiederum ein stochastischer Prozess mit unkorrelierten Zuwächsen. Falls $\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}|\gamma (h)|dh<\mathrm{\infty }$, dann hat die Spektralverteilungsfunktion $F$ eine Spektraldichte $f$, und es gilt:

$\gamma (h)=\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}{e}^{i\lambda h}f(\lambda )d\lambda ;f(\lambda )=\frac{1}{2\pi }\underset{-\mathrm{\infty }}{\overset{\mathrm{\infty }}{\int }}\gamma (h)e^{-ih\lambda }dh$.

Beispiele

• Ein stationärer Prozess $X_t$ mit der häufig benutzten Kovarianzfunktion $\gamma (h)={\sigma }^2{e}^{-\alpha |h|}$, wobei $\alpha$ eine positive Konstante ist, hat die Spektraldichte $f(\lambda )=\frac{{\sigma }^2\alpha }{\pi ({\lambda }^2+{\alpha }^2)}$.
• Weißes Rauschen hat die Kovarianzfunktion $\gamma (h)=\{\begin{array}{cc}{\sigma }^2& \text{wenn }h=0\\ 0& \text{sonst }\end{array}$ und die Spektraldichte $f(\lambda )=\frac{1}{2\pi }{\sigma }^2$.

Die Spektraldichte ist also konstant. Alle Frequenzen sind gleichstark im Spektrum vertreten (Analogie zum weißen Licht).

Anwendungen

Spektraldarstellungen benötigt man in der Zeitreihenanalyse, in der Signalverarbeitung (siehe z. B. auch Spektrale Leistungsdichte), bei der Konstruktion geeigneter Filter (beispielsweise Tiefpass, Hochpass oder Bandpass).

Besonders wichtig in den Anwendungen sind geeignete Methoden zur Spektraldichteschätzung.