Quelltext der Seite Gleitender Mittelwert

Der '''gleitende Durchschnitt''' (auch '''gleitender Mittelwert''') ist eine Methode zur [[Glätten (Mathematik)|Glättung]] von Zeit- bzw. Datenreihen. Die Glättung erfolgt durch das Entfernen höherer Frequenzanteile. Im Ergebnis wird eine neue Datenpunktmenge erstellt, die aus den [[Mittelwert]]en gleich großer Untermengen der ursprünglichen Datenpunktmenge besteht. In der [[Signaltheorie]] wird der gleitende Durchschnitt als [[Tiefpass]]filter mit [[Filter mit endlicher Impulsantwort|endlicher Impulsantwort]] (FIR-Tiefpass) beschrieben. In der gleichgewichteten Form stellt der gleitende Durchschnitt das einfachste FIR-Tiefpassfilter dar.

Anwendung findet der gleitende Durchschnitt beispielsweise bei der [[Zeitreihenanalyse|Analyse von Zeitreihen]]. Die gleichgewichtete Variante entspricht der [[Faltung (Mathematik)|Faltung]] mit einer [[Rechteckfunktion]] und führt zu einer Reihe von Problemen, denen beispielsweise mit speziellen [[Fensterfunktion|Gewichtungen]] entgegengewirkt werden kann.

== Prinzipielles Vorgehen ==
Die Menge der gleitenden Durchschnitte wird iterativ („gleitend“) über einen Ausschnitt, das „Fenster“, eines gegebenen [[Signal#Messtechnik|Signals]] berechnet. Das verwendete Fenster wird überlappend verschoben, d.&nbsp;h., wiederholt wird der letzte Wert aus dem betrachteten Ausschnitt gestrichen, der erste Wert nach dem Ausschnitt hinzugenommen und ein neuer Mittelwert berechnet. Für die Berechnung des Mittelwerts können die im Fenster vorkommenden Werte anschließend beliebig gewichtet werden.

Die Ergebnismenge der Mittelwerte ist für sich genommen unabhängig. Häufig wird sie jedoch in den Zusammenhang mit einer Position der Eingabemenge gebracht, die „hot spot“ genannt wird. Der hot spot kann im Bereich des Fensters liegen, muss es aber nicht. Bei Zeitreihen wird häufig der letzte Zeitpunkt als hot spot verwendet; in anderen Anwendungsfällen sind zentrierte Abbildungen üblich.

== Einfacher gleitender Durchschnitt ==
=== Berechnung und Zentrierung der Berechnung ===
[[Datei:MittelwertSignalvergleich3.svg|mini|400px|Zentrierte gleitende Durchschnitte der Breite 3 verglichen mit einem binomial gefalteten Signal gleicher Breite. Im Gegensatz zur Dämpfung bzw. Auslöschung der hohen Frequenzen zwischen <math>t=7</math> und <math>t=12</math> (schnelle Wechsel zwischen hohen und niedrigen Werten) wird die Signalphase vom gleitenden Durchschnitt invertiert, d.&nbsp;h., wo zuvor ein hoher Wert war, ist nun ein niedriger, und umgekehrt. Das Binomialfilter verursacht dagegen keine Phaseninversion.]]
Der einfache gleitende Durchschnitt ({{enS|''simple moving average'' (SMA)}}) <math>n</math>-ter Ordnung einer diskreten Zeitreihe <math>x(t)</math> ist die Folge der [[Arithmetischer Mittelwert|arithmetischen Mittelwerte]] von <math>n</math> aufeinanderfolgenden Datenpunkten. Da es sich um eine Zeitreihe handelt, liegt der ''hot spot'' auf dem letzten Zeitpunkt. Die nachfolgenden Ausführungen beziehen sich auf diesen Sonderfall.

:<math>
m_\text{MA}^{(n)}(t) = \frac{1}{n} \sum_{i=0}^{n-1} x(t-i)
</math>

In Analogie zu Filtern mit endlicher Impulsantwort wird <math>n+1</math> auch ''Ordnung'' genannt.

Solch ein gleitender Durchschnitt hat eine Verzögerung ([[Gruppenlaufzeit]]) von <math>\tau = \frac{n-1}{2}</math>, d.&nbsp;h. die gemittelten Werte „hinken“ um <math>\tau</math> Zeiteinheiten hinterher.

Diese Verzögerung kann korrigiert werden, indem man den gleitenden Durchschnitt um <math>\tau</math> verschiebt. Dies ist der sogenannte ''zentrierte Durchschnitt''. Dann sind jedoch keine Werte mehr für die ersten und letzten <math>\tau</math> Zeiteinheiten vorhanden. Ohne den Einsatz eines dynamischen Fensters kann diese Lücke nur geschlossen oder zumindest verkleinert werden, indem man eine geringere Ordnung, andere Gewichtungen oder einen [[Schätzfunktion|Schätzer]] verwendet.

Der zentrierte einfache gleitende Durchschnitt der Ordnung 3 ist also durch
:<math>
m_\text{MA}^{(3)}(t) = \frac{1}{3} \Big( x(t-2) + x(t-1) + x(t)  \Big)
</math>
gegeben.

Ein Beispiel für die Verwendung von derartigen gleitenden Durchschnitten sind die 38- bzw. 200-Tage-Durchschnittswerte von [[Börsenkurs]]en, die den gleitenden Durchschnitt der ''vergangenen'' <math>n</math> [[Handelstag|Börsentage]] eines [[Wertpapier]]kurses beschreiben.

Ein Gleitender-Mittelwert-Filter ist ein Tiefpass-Filter, jedoch werden einzelne Frequenzbereiche stärker oder schwächer gefiltert, und es kommt zu Signalverschiebungen ({{enS|[[Lag]]s}}).
[[Datei:Frequenzsweep-Grafisch.png|mini|Die Wirkung eines Glei&shy;tender-Mittel&shy;wert-Filters auf ein [[Chirp]]-Signal (oberes Drittel). Bis zum „Artefakt“ arbeitet der gleitende Durchschnitt (mittleres Drittel) gut als Tiefpass&shy;filter. Rechts davon jedoch wird das Signal wieder stärker durch&shy;gelassen, mal mit invertierter Phase, mal lediglich gedämpft. Zum Vergleich das Ergebnis nach Anwendung eines [[Binomialfilter|binomialen Filters]] (unteres Drittel).|398x398px]]

=== Gleitender Durchschnitt mit dynamischem Fenster ===

Der gleitende Durchschnitt <math>n</math>-ter Ordnung einer diskreten Menge <math>x(t)</math> mit <math>|x(t)|</math> Elementen hat <math>n-1</math> Wertepaare weniger als die Gesamtmenge. Da sehr stark oszillierende Graphen eine entsprechend hohe Ordnung zur ausreichenden Glättung der Ergebnismenge benötigen, kann dieser Effekt einen spürbaren Einfluss auf die Darstellung haben. Die resultierende Trendlinie liegt dann bei nicht zentriertem einfachen gleitenden Durchschnitt verstärkt versetzt und gibt somit eine zunehmend schlechtere Abstraktion der Originaldaten wieder.

Bei einer diskreten Menge kann der Versetzungseffekt durch die Verwendung eines dynamischen Mittelwertfensters vermieden werden. Bei diesem Verfahren wird jeder Mittelwert <math>m_\text{MA}^{(n)}(t)</math> aus den <math>n</math> Parametern <math>x(t-i)</math> bis <math>x(t+j)</math> berechnet. Die Fensterränder <math>i</math> und <math>j</math> leiten sich aus der jeweiligen relativen Position <math>p_\text{rel}(t)</math> ab, mit:

:<math>
p_\text{rel}(t) = \Big\lceil \frac{p_\text{abs}(t)}{|x(t)|} \cdot n \Big\rceil
</math>

<math>p_\text{abs}</math> ist dabei die (absolute) Position von <math>t</math> in <math>x(t)</math>. Das mit Hilfe von <math>i = p_\text{rel} - 1</math> und <math>j = n - p_\text{rel}</math> definierte dynamische Fenster zur Bildung der Mittelwerte deckt die volle Breite des ursprünglichen Graphen ab und abstrahiert die Originaldaten von Rundungsfehlern abgesehen versetzungsfrei. Geometrisch wird die [[Teilverhältnis|Position]] des Elements in der Gesamtmenge auf die Position des Elements in der Auswahlmenge abgebildet. Randelemente landen im entsprechenden Rand, das zentrale Element der Gesamtmenge liegt auch im Auswahlintervall mittig. Je mittiger das Element umso geringer der Versatz.

[[Datei:Gleitender Mittelwert auf diskreten Daten.png|alt=Gleitender Mittelwert auf diskreten Daten|mini|395x395px|Die Abbildung zeigt den klassischen gleitenden Mittelwertfilter (ohne Zentrierung) mit Versetzung nach rechts und dem Verlust von 25 Wertepaaren (rote Linie) im direkten grafischen Vergleich mit dem gleitenden Mittelwert unter Verwendung eines dynamischen Fensters (blaue Linie).]]

'''Beispiel'''

Das Wachstum der volkswirtschaftlichen [[Produktivität]] zwischen 1891 und 2012 schwankt sehr stark. Der Graph gebildet aus dem gleitenden Mittelwert der Ordnung <math>n = 25</math> glättet diese Kurve gut (rote Linie).

'''Berechnungsbeispiel für den ersten möglichen Wert von 1915'''

:<math>
m_\text{MA}^{(25)}(1915) = \frac{1}{n} \cdot \Big( x(t-24) + \ldots + x(t) \Big) = \frac{1}{25} \cdot \Big( x(1891) + \ldots + x(1915) \Big)
</math>

Die geglättete Kurve ist eine zwar aussagekräftige, aber deutlich versetzte Trendlinie. Zudem gehen dabei gut 20 % der Wertepaare verloren. Mit Hilfe eines dynamischen Fensters lassen sich die Werte über die gesamte Breite von <math>x(t)</math> berechnen. Das Ergebnis ist eine besonders im Inneren des Untersuchungszeitraums kaum versetzte Trendlinie (blaue Linie).

'''Berechnungsbeispiel für 1891, 1915 und 2012'''

Die Menge <math>x(t)</math> hat 122 Elemente: <math>|x(t)| = 122</math>, 1891 ist der Zeitindex des ersten Elements <math>p_\text{abs}(1891) = 1</math> und <math>n = 25</math>.

Damit folgen <math>p_\text{rel}(1891) = \Big\lceil \frac{1}{122} \cdot 25 \Big\rceil = 1</math> und <math>i = 0, j = 24</math> und es wird das Intervall <math>x(1891-0)</math> bis <math>x(1891+24)</math> betrachtet (die anderen Werte analog).
:<math>
\begin{align}
m_\text{MA}^{(25)}(1891)
&= \frac{1}{n} \cdot \Big( x(t-p_\text{rel}+1) + \ldots + x(t) + \ldots + x(t+n-p_\text{rel}) \Big) \\
&= \frac{1}{25} \cdot \Big( x(1891) + \ldots + x(1915) \Big)
\end{align}
</math>

:<math>
m_\text{MA}^{(25)}(1915) = \frac{1}{25} \cdot \Big( x(1911) + \ldots + x(1935) \Big)
</math>

:<math>
m_\text{MA}^{(25)}(2012) = \frac{1}{25} \cdot \Big( x(1988) + \ldots + x(2012) \Big)
</math>

=== Online-Berechnung ===

Durch die Überlappung bei der Berechnung des gleitenden Mittels für zwei aufeinanderfolgende Punkte wird die Summe für <math>n-1</math> Punkte doppelt berechnet. Um diesen [[Redundanz (Informationstheorie)|redundanten]] Aufwand zu reduzieren, existiert ein Online-[[Algorithmus]], der mit nur zwei Additionen und Multiplikationen pro Mittelwert auskommt:

:<math>
 \begin{align}
   m^{(n)}(t) &= \frac{1}{n} \sum_{i=1}^{n} x(t-n+i)  \\
              &= \frac{1}{n} \Big( \underbrace{ x(t-n) - x(t-n) }_{= 0} + \underbrace{ x(t-n+1) + x(t-n+2) + \dots + x(t) }_{ \sum_{i=1}^{n} x(t-n+i) } \Big)  \\
              &= \underbrace{ \frac{1}{n} \Big( x(t-n) + x(t-n+1) + x(t-n+2) + \dots + x(t-1) \Big) }_{= m^{(n)}(t-1)} - \frac{x(t-n)}{n} + \frac{x(t)}{n}  \\
              &= m^{(n)}(t-1) + \frac{x(t)}{n} - \frac{x(t-n)}{n}
 \end{align}
</math>

== Gewichteter gleitender Durchschnitt ==

Der gewichtete gleitende Durchschnitt <math>m(t)</math> der Ordnung <math>n</math> einer Zeitreihe <math>x(t)</math> ist – analog zum [[Filter mit endlicher Impulsantwort]] – definiert als:

:<math>
m_x^{(n)}(t-\tau) = \sum_{i=0}^{n-1} w_i \cdot x(t-i)
</math>

Dabei stellt <math>w_i</math> die Gewichtung der jeweiligen Datenpunkte dar (äquivalent zur Impulsantwort des Filters). Wenn <math>\tau > 0</math>, ist das Filter nicht kausal, sondern berücksichtigt zukünftige Werte bei der Durchschnittsbildung.
Die Summe über alle Gewichte muss 1 ergeben, da sonst noch eine Verstärkung (<math>\textstyle \sum_{i=0}^{n-1} w_i > 1</math>) oder Dämpfung (<math>\textstyle \sum_{i=0}^{n-1} w_i < 1</math>) hinzukommt.

Ein Beispiel ist das zentrierte Binomialfilter dritter Ordnung mit <math>w_{0}=\tfrac14</math>, <math>w_{1}=\tfrac12</math>, <math>w_{2}=\tfrac14</math> und <math>\tau = 1</math>:
:<math>
m_\text{BIN}(t-1) = \frac14 x(t) + \frac12 x(t-1) + \frac14 x(t-2)
</math>

== Zentrierter gleitender Durchschnitt mit spektralen Eigenschaften ==

[[Datei:Bin3ma3transf.svg|mini|Über&shy;tragungs&shy;funktion im Frequenz&shy;bereich des glei&shy;tenden Mittel&shy;werts der Breite 3 (<math>\varphi_{\text{MA}_3}</math>) im Vergleich mit der eines Binomial&shy;filters gleicher Breite (<math>\varphi_{\text{BIN}_3}</math>)]]

Bildet man den zentrierten gleitenden Durchschnitt <math>n</math>-ter Ordnung <math>\bar X^{(n)}</math> einer [[Schwach stationärer Prozess|schwach stationären]] Zeitreihe <math>X_t</math> mit [[Spektraldarstellung stationärer stochastischer Prozesse#Spektraldichte|Spektraldichte]] <math>f_X</math>, dann hat <math>\bar X^{(n)}</math> die gefilterte Spektraldichte
:<math>|\varphi(\omega)|^2\cdot f_X(\omega) </math>
mit der [[Übertragungsfunktion]]
:<math>|\varphi(\omega)|^2 := \frac{2\pi}n k_n(\omega)</math>,
wobei <math>k_n</math> den [[Fejér-Kern]] bezeichnet. An der grafischen Darstellung für <math>n=3</math> mit der Übertragungsfunktion <math>\varphi(\omega)_{\text{MA}_3} = (1 + 2\cos(\omega))/3</math> kann man die Tiefpass-Eigenschaft erkennen: Frequenzen nahe 0 werden nicht gedämpft. Andererseits zeigt dieses einfache Filter das übliche [[Sinc-Funktion|Antwortverhalten bei der Faltung mit einem Rechtecksignal]]. Bei einer Filterbreite von 3 werden die Frequenzen bis zum Punkt <math>\omega = \tfrac23\pi</math> zunehmend bis zur vollständigen Unterdrückung gedämpft. Frequenzen, die über diesen Punkt hinaus vorhanden sind, werden nicht etwa auch unterdrückt, sondern treten mit invertierter Phase auf.

Das kleinste [[Gauß-Filter|Binomialfilter]] mit ungerader Breite und mit den Gewichten <math>\tfrac14</math>, <math>\tfrac12</math>, <math>\tfrac14</math> ist ein Tiefpass-Filter mit <math>\varphi(\omega)_{\text{BIN}_3} = (1+ \cos(\omega))/2 \ge 0</math> für alle Frequenzen. Es dämpft die Frequenzen bis <math>\omega = \pi</math> zunehmend und mit konstanter [[Phasenverschiebung]].<ref>Jens-Peter Kreiß, Georg Neuhaus: ''Einführung in die Zeitreihenanalyse.'' Springer, 2006, ISBN 3-540-25628-8.</ref>

[[Datei:FrequenzSweepGlMWBinomial7.png|mini|Chirp-Signal bis zur [[Nyquist-Frequenz]] (oben) und mit gleitendem Durchschnitt (Breite: 7) gefaltete Variante (Mitte). Unten zum Vergleich ein mit einem Binomialfilter gefaltetes Signal (Breite: 7). Die Glättung mit dem gleitenden Durchschnitt sollte zu einem Signal führen, das die tiefen Frequenzen (links) unverändert enthält, die hohen Frequenzen (rechts) jedoch herausfiltert. Zwischen diesen Extremen wird zunehmend gedämpft. Der ungewichtete gleitende Durchschnitt (Mitte) erfüllt diese Aufgabe nur sehr unzureichend, das Binomialfilter dagegen erheblich besser.]]

== Linear gewichteter gleitender Durchschnitt ==

Ein linear gewichteter gleitender Durchschnitt (engl.: ''linear weighted moving average'' (LWMA, meist: WMA)) ordnet den Datenpunkten linear aufsteigende [[Gewichtung|Gewichte]] zu, d.&nbsp;h. je weiter die Werte in der Vergangenheit liegen, desto geringer ist ihr Einfluss:

:<math>
w_i = \frac{2i}{n (n+1)} \quad \forall i=1, \dots, n
</math>

:<math>
m_\text{LWMA}^{(n)}(t) = \frac{2}{n (n+1)} \sum_{i=1}^n i \cdot x(t-n+i)
</math>

== Exponentiell geglätteter Durchschnitt ==

{{Hauptartikel|Exponentielle Glättung}}

Der exponentiell geglättete Durchschnitt ordnet den Datenpunkten einer Zeitreihe exponentiell abnehmende Gewichte zu. Somit werden auch hier jüngere Datenpunkte stärker gewichtet als weiter zurückliegende, jedoch noch stärker als beim gewichteten gleitenden Durchschnitt.

:<math>
m_\text{EMA}^{(n)}(t) = \alpha \cdot x(t) + (1 - \alpha) \cdot m_\text{EMA}^{(n)}(t-1)
</math>

Da der exponentielle Durchschnitt nicht nur Werte aus der Zeitreihe, sondern auch vorangegangene Mittelwerte miteinbezieht, stellt er ein [[Filter mit unendlicher Impulsantwort]] dar. Ein entscheidender Vorteil ist seine wesentlich kürzere Verzögerung bei gleicher Glättung.

== Exponentiell gewichteter geglätteter Durchschnitt ==
Der exponentiell gewichtete geglättete Durchschnitt (EMWA) ist eine Verallgemeinerung des exponentiell geglätteten Durchschnitts, wobei Gewichte eingeführt werden.

== Siehe auch ==

* [[ARMA-Modell]]

== Literatur ==
* John G. Proakis, Dimitris G. Manolakis: ''Digital Signal Processing.'' 4. Auflage. Prentice Hall, 2007, ISBN 978-0-13-187374-2.

== Weblinks ==
{{Wikibooks|Statistik: Glättungsverfahren|Gleitende Durchschnitte}}

== Einzelnachweise ==
<references />

[[Kategorie:Mittelwert]]
[[Kategorie:Zeitreihenanalyse]]
[[Kategorie:Materialbedarfsermittlung]]