Stationäre Verteilung

Invariante Verteilung oder stationäre Verteilung ist ein Begriff aus der Theorie der Markow-Ketten. Diese sind spezielle stochastische Prozesse und damit Untersuchungsobjekte der Stochastik.

Eine Markow-Kette besitzt eine stationäre Verteilung genau dann, wenn es eine Startverteilung gibt, die sich im Zeitverlauf nicht ändert.

Stationäre Verteilungen sind nicht zu verwechseln mit der Stationarität eines stochastischen Prozesses oder mit stationären Übergangswahrscheinlichkeiten.

Gegeben sei eine homogene Markow-Kette ${\displaystyle (X_{n})_{n\in \mathbb {N} }}$ in diskreter Zeit mit höchstens abzählbarem Zustandsraum ${\displaystyle Z}$. Dann heißt eine Verteilung ${\displaystyle P_{\pi }}$ auf ${\displaystyle Z}$ stationäre Verteilung, wenn

${\displaystyle P_{\pi }(\{j\})=\sum _{i\in Z}P_{\pi }(\{i\})p(i,j)}$

für alle ${\displaystyle j\in Z}$ gilt, wobei ${\displaystyle p(i,j)=P(X_{n}=j|X_{n-1}=i)}$ die Übergangswahrscheinlichkeit von Zustand ${\displaystyle i}$ in den Zustand ${\displaystyle j}$ ist (die unabhängig vom Zeitpunkt ${\displaystyle n}$ ist). Im Falle eines endlichen Zustandsraumes entspricht dies

${\displaystyle \pi M=\pi }$,

wobei ${\displaystyle M}$ die Übergangsmatrix ist und ${\displaystyle \pi }$ ein Wahrscheinlichkeitsvektor als Zeilenvektor geschrieben. Damit sind die stationären Verteilungen in diesem Fall genau die Linkseigenvektoren der Übergangsmatrix zum Eigenwert 1, die bezüglich der Summennorm normiert wurden.

Im Allgemeinen müssen keine stationären Verteilungen existieren. Beispiel hierfür sind transiente Markow-Ketten. Diese besitzen nie stationäre Verteilungen. Umgekehrt lässt sich auch zeigen, dass irreduzible Markow-Ketten höchstens eine stationäre Verteilung besitzen. Für die Eindeutigkeit der stationären Verteilungen gelten folgende Aussagen:

• Eine irreduzible Markow-Kette besitzt genau dann eine stationäre Verteilung, wenn sie positiv rekurrent ist. Die Verteilung ist dann gegeben durch

${\displaystyle P_{\pi }(\{i\})={\frac {1}{\operatorname {E} (\tau _{ii})}}}$.

Hierbei ist ${\displaystyle \tau _{ii}}$ die Wiederkehrzeit in den Zustand ${\displaystyle i}$, wenn in diesem gestartet wurde.

• Im Falle eines endlichen Zustandsraumes sind Irreduzibilität der Markow-Kette und Irreduzibilität der Übergangsmatrix äquivalent. Daraus folgt aber sofort mit dem Satz von Perron-Frobenius, dass eine eindeutige invariante Verteilung (bzw. Linkseigenvektor) existiert und damit, dass die Markow-Kette positiv-rekurrent ist. Somit folgt hier aus Irreduzibilität positive Rekurrenz.
• Demnach hat eine irreduzible Markow-Kette mit endlichem Zustandsraum immer eine stationäre Verteilung. Diese entspricht genau dem normierten Linkseigenvektor zum Eigenwert 1 der Übergangsmatrix ${\displaystyle M}$ bzw. dem normierten Eigenvektor der transponierten Übergangsmatrix ${\displaystyle M^{T}}$ zum Eigenwert 1.
• Erfüllt eine Verteilung die Detailed-Balance-Gleichung, so ist diese Verteilung eine stationäre Verteilung.

Eine irreduzible, positiv rekurrente Markow-Kette konvergiert genau dann gegen eine stationäre Verteilung, wenn sie aperiodisch ist. Konvergenz bedeutet hier, dass

${\displaystyle \lim _{n\to \infty }P(X_{n}=i)=P_{\pi }(\{i\})}$

für jede Startverteilung von ${\displaystyle X_{0}}$ und jeden Zustand ${\displaystyle i\in Z}$ gilt.

Ist der Zustandsraum endlich, so konvergieren dann die Zeilen von ${\displaystyle M^{n}}$ gegen die stationäre Verteilung.

Bei endlichen Zustandsräumen findet sich oft das Konvergenzkriterium, dass ein ${\displaystyle k\in \mathbb {N} }$ existieren muss, sodass für die Übergangsmatrix ${\displaystyle M}$ gilt, dass ${\displaystyle M^{k}>0}$ ist. Dies entspricht der Überprüfung der Matrix auf Aperiodizität und Irreduzibilität.

Verzichtet man auf die Aperiodizität, so lässt sich folgende Aussage zeigen: Ist eine Markow-Kette irreduzibel und rekurrent, so ist

${\displaystyle \lim _{N\to \infty }{\frac {1}{N}}\sum _{n=0}^{N-1}P(X_{n}=i)={\frac {1}{\operatorname {E} (\tau _{ii})}}=P_{\pi }(\{i\})}$.

Der Mittelwert der Eintrittswahrscheinlichkeiten konvergiert also komponentenweise gegen die stationäre Verteilung.

Allgemeiner gilt: Ist die Markow-Kette nicht irreduzibel, so zerfällt sie in mehrere Teilmengen von Zuständen, die miteinander kommunizieren und alle dieselbe Periode ${\displaystyle d}$ besitzen. Sei ${\displaystyle K(i)}$ so eine Menge mit einem beliebigen, aber fixierten Zustand ${\displaystyle i}$ aus dieser Menge. Dann lässt sich jeder Zustand ${\displaystyle j}$ aus dieser Menge in einer eindeutigen Zahl ${\displaystyle k}$ von Schritten von ${\displaystyle i}$ aus erreichen. Ist die Teilmenge ${\displaystyle T(i)}$ nun rekurrent, so gilt

${\displaystyle \lim _{n\to \infty }p_{ij}^{nd+k}={\frac {d}{\operatorname {E} (\tau _{jj})}}}$.

Für einige Anwendungen ist vor allem interessant, wie schnell die stationäre Verteilung erreicht wird. Es lässt sich zeigen, dass wenn ${\displaystyle 1=\lambda _{1}>|\lambda _{2}|\geq |\lambda _{3}|\dots }$ für ${\displaystyle M}$ gilt und alle Eigenwerte einfach sind, die folgende Abschätzung gilt:

${\displaystyle \Vert v_{0}M^{n}-P_{\pi }\Vert \leq C|\lambda _{2}|^{n}}$

für eine beliebige Startverteilung ${\displaystyle v_{0}}$. Wichtigster Einflussfaktor auf die Konvergenzgeschwindigkeit ist also der betragsmäßig zweitgrößte Eigenwert.

Es lassen sich noch vergleichbare Aussagen für schwächere Voraussetzungen an die Übergangsmatrix zeigen, dabei müssen aber Korrekturterme für die Jordanblöcke eingeführt werden.

Betrachten wir die Markow-Kette mit der folgenden Übergangsmatrix:

${\displaystyle M={\begin{bmatrix}1&0&0&0&0\\0{,}1&0&0{,}8&0{,}1&0\\0&0{,}5&0&0&0{,}5\\0&0&0&0&1\\0&0&0&1&0\end{bmatrix}}}$

Diese Markow-Kette hat zwei absorbierende Mengen: ${\displaystyle A=\{1\}}$ und ${\displaystyle B=\{4,5\}}$. Da diese Zustände nicht mehr verlassen werden können, haben sie einen Einfluss auf die Existenz der stationären Verteilungen. Dies zeigt sich auch in den normierten Eigenvektoren. Diese sind ${\displaystyle \pi _{1}=(1,0,0,0,0)}$ und ${\displaystyle \pi _{2}=(0,0,0,0{,}5,0{,}5)}$. Somit ist hier die stationäre Verteilung nicht eindeutig.

Betrachten wir die Markow-Kette mit dem rechts dargestellten Übergangsgraph. Der Einfachheit halber setzen wir alle Übergangswahrscheinlichkeiten auf 0,5. Die Markow-Kette ist irreduzibel, da man sich von jedem Zustand in maximal zwei Schritten in jeden anderen Zustand bewegen kann. Sie ist aber auch periodisch, da eine Rückkehr zum Startpunkt nur zu geraden Zeitpunkten möglich ist. Die ebenfalls irreduzible Übergangsmatrix ist dann

${\displaystyle M={\begin{bmatrix}0&0{,}5&0&0{,}5\\0{,}5&0&0{,}5&0\\0&0{,}5&0&0{,}5\\0{,}5&0&0{,}5&0\end{bmatrix}}}$

Der Satz von Perron-Frobenius garantiert die Eindeutigkeit des Eigenvektors. Da die Matrix zusätzlich doppelt-stochastisch ist, hat sie die stationäre Verteilung ${\displaystyle \pi ={\frac {1}{4}}(1,1,1,1)}$. Die Matrixpotenzen konvergieren aber nicht, insbesondere ist ${\displaystyle M^{2n+1}=M}$ und

${\displaystyle M^{2n}=M^{2}={\begin{bmatrix}0{,}5&0&0{,}5&0\\0&0{,}5&0&0{,}5\\0{,}5&0&0{,}5&0\\0&0{,}5&0&0{,}5\end{bmatrix}}}$

Betrachten wir aber nun die Mittelwerte, so konvergieren diese gegen die entsprechende Komponente der stationären Verteilung:

${\displaystyle \lim _{N\to \infty }{\frac {1}{N}}\sum _{n=0}^{N-1}P(X_{n}=1|X_{0}=1)=\lim _{N\to \infty }{\frac {1}{N}}\sum _{n=0}^{N-1}p_{11}^{n}\rightarrow {\frac {1}{4}}{\text{ für }}N\to \infty }$

Dies folgt hier mithilfe der entsprechenden Einträge (im Beispiel die erste Zeile und Spalte) der obigen Übergangsmatrix.

Betrachte die rechts dargestellte Markow-Kette mit den Übergangswahrscheinlichkeiten wie in der Übergangsmatrix angegeben:

${\displaystyle M={\begin{bmatrix}0{,}5&0&0{,}5\\0{,}1&0{,}5&0{,}4\\0&0{,}5&0{,}5\end{bmatrix}}}$

Es gilt dann

${\displaystyle M^{2}={\frac {1}{100}}{\begin{bmatrix}25&25&50\\10&45&45\\5&50&45\end{bmatrix}}>0}$

Somit ist die Markow-Kette irreduzibel (und damit auch positiv rekurrent), aperiodisch und konvergiert demnach gegen eine stationäre Verteilung. Ein Eigenvektor von ${\displaystyle M^{T}}$ zum Eigenwert 1 ist ${\displaystyle v=(1,5,5)}$, Normierung auf 1 bzgl. der Summennorm liefert als eindeutige invariante Verteilung

${\displaystyle \pi ={\frac {1}{11}}(1,5,5)=(0.0909,0.4545,0.4545)}$.

Berechnet man die Matrixpotenzen, so stimmen bei ${\displaystyle M^{6}}$ schon zwei Nachkommastellen mit der exakten Lösung überein, bei ${\displaystyle M^{12}}$ schon mehr als vier Nachkommastellen. Umgekehrt kann man aus der als Linkseigenvektor berechneten stationären Verteilung bei Konvergenz sofort den Grenzwert der Matrixpotenzen angeben, da diese zeilenweise gegen die stationäre Verteilung konvergieren:

${\displaystyle \lim _{n\to \infty }M^{n}={\frac {1}{11}}\cdot {\begin{bmatrix}1&5&5\\1&5&5\\1&5&5\end{bmatrix}}}$

Aus der stationären Verteilung kann man auch die erwartete Rückkehrzeit berechnen, diese ist genau der Kehrwert der entsprechenden Komponente der Verteilung. Somit ist hier die durchschnittliche Zeit beim Start in 1 bis zur Rückkehr

${\displaystyle \operatorname {E} (\tau _{11})={\frac {1}{\pi _{1}}}=11}$.

Die Spielfigur Pac-Man frisst in einem Labyrinth kleine Kugeln und Kirschen und wird dabei von Gespenstern verfolgt. Der Einfachheit halber ist die Spielwelt in diesem Beispiel ein kleines ${\displaystyle 3\times 3}$-Gitter und die Monster bewegen sich rein zufällig. Jedes horizontal und vertikal angrenzende Spielfeld ist mit gleicher Wahrscheinlichkeit der nächste Aufenthaltsort des Gespensts, mit Ausnahme eines Geheimgangs zwischen den Zuständen ${\displaystyle 2}$ und ${\displaystyle 8}$ (siehe nebenstehenden Übergangsgraphen).

Der Zustandsraum lautet ${\displaystyle Z=\{1,2,3,4,5,6,7,8,9\}}$. In der nun folgenden Übergangsmatrix ${\displaystyle P}$ wurden Einträge mit Wahrscheinlichkeit ${\displaystyle 0}$ entfernt, um eine bessere Übersichtlichkeit zu erhalten:

${\displaystyle P={\begin{pmatrix}&{\frac {1}{2}}&&{\frac {1}{2}}&&&&&\\{\frac {1}{4}}&&{\frac {1}{4}}&&{\frac {1}{4}}&&&{\frac {1}{4}}&\\&{\frac {1}{2}}&&&&{\frac {1}{2}}&&&\\{\frac {1}{3}}&&&&{\frac {1}{3}}&&{\frac {1}{3}}&&\\&{\frac {1}{4}}&&{\frac {1}{4}}&&{\frac {1}{4}}&&{\frac {1}{4}}&\\&&{\frac {1}{3}}&&{\frac {1}{3}}&&&&{\frac {1}{3}}\\&&&{\frac {1}{2}}&&&&{\frac {1}{2}}&\\&{\frac {1}{4}}&&&{\frac {1}{4}}&&{\frac {1}{4}}&&{\frac {1}{4}}\\&&&&&{\frac {1}{2}}&&{\frac {1}{2}}&\end{pmatrix}}}$

Diese Markow-Kette ist irreduzibel, da sich die Gespenster in endlicher Zeit von jedem beliebigen Zustand in jeden beliebigen Zustand begeben können. Dank des Geheimgangs sind hierfür nur maximal drei Zustandswechsel nötig. Durch den Geheimgang ist die Markow-Kette auch aperiodisch, weil die Monster sowohl in einer geraden als auch in einer ungeraden Anzahl an Zustandswechseln von jedem beliebigen Zustand in jeden beliebigen Zustand wechseln können. Ohne den Geheimgang wäre die Markow-Kette periodisch, weil dann ein Übergang von einem geraden in einen geraden Zustand bzw. von einem ungeraden in einen ungeraden Zustand nur in einer geraden Anzahl von Zustandswechseln möglich wäre.

Wegen der Irreduzibilität und Aperiodizität gibt es genau eine stabile Gleichgewichtsverteilung, welche die Markow-Kette nach einer unendlich langen Zeit annimmt. Die Aufenthaltswahrscheinlichkeiten für die einzelnen Zustände ändern sich nach langer Zeit (fast) nicht mehr. Die stationäre Verteilung lässt sich naiv bestimmen, indem in die Gleichung

${\displaystyle \pi _{n}=\pi _{0}\cdot P^{n}}$

für eine beliebige Startverteilung ${\displaystyle \pi _{0}}$ ein großes ${\displaystyle n}$ eingesetzt wird, weil die Matrixpotenzen wie im obigen Beispiel konvergieren. Um eine analytische Lösung zu berechnen, ist das lineare Gleichungssystem

${\displaystyle \pi =\pi \cdot P}$

nach ${\displaystyle \pi }$ aufzulösen, unter der Nebenbedingung einer Zeilensumme von ${\displaystyle 1}$. Das Einsetzen der naiven Lösung in diese Gleichung dient als Kontrolle. Die obige Gleichung ist äquivalent zu

${\displaystyle (P^{T}-I)\cdot \pi ^{T}=0}$.

Die Übergangsmatrix wird demnach transponiert und die Einheitsmatrix subtrahiert. Der gesuchte Vektor der Zustandswahrscheinlichkeiten ist nun ein Spaltenvektor. Die Lösung des linearen Gleichungssystems ergibt:

${\displaystyle \pi =(7.7,15.4,7.7,11.5,15.4,11.5,7.7,15.4,7.7)\ \%}$

Die Gespenster halten sich demnach am häufigsten in der Mitte auf, weniger oft am Rand und am seltensten in der Ecke. Eine Ausnahme bilden die Randzustände ${\displaystyle 2}$ und ${\displaystyle 8}$, die aufgrund des Geheimwegs durchschnittlich genauso oft besucht werden wie das zentrale Spielfeld.

Betrachten wir eine Markow-Kette auf dem Zustandsraum ${\displaystyle Z=\mathbb {N} _{0}}$ mit Übergangswahrscheinlichkeiten

${\displaystyle p_{ij}={\begin{cases}{\frac {1}{i+1}},&{\text{ falls }}j=i+1,\\1-{\frac {1}{i+1}},&{\text{ falls }}j=0,\\0&{\text{ sonst.}}\end{cases}}}$

Mit einer gewissen Wahrscheinlichkeit steigt man also zu einer Zahl eins höher auf. Falls dies nicht geschieht, startet man wieder am Nullpunkt. Es gilt:

1. Alle Zustände kommunizieren miteinander, da jeder Zustand die 0 erreichen kann und die 0 jeden Zustand erreicht. Demnach ist die Markow-Kette irreduzibel.
2. Die Rückkehrzeiten in die 0 sind ${\displaystyle \{2,3,4,\dotsc \}}$ und demnach hat die 0 die Periode 1, ist also aperiodisch. Aufgrund der Irreduzibilität ist dann auch die ganze Markow-kette aperiodisch.
3. Die erwartete Rückkehrzeit zu 0 ist gegeben durch:

${\displaystyle \operatorname {E} (\tau _{00})=\sum _{k=2}^{\infty }{\frac {k}{(k-1)!}}\cdot (1-{\frac {1}{k}})=e}$

Somit ist die 0 positiv rekurrent und aufgrund der Irreduzibilität der Markow-Kette auch die ganze Kette positiv rekurrent.

Die Kette konvergiert also gegen eine von der Startverteilung unabhängige invariante Verteilung. Da wir bereits wissen, dass ${\displaystyle P_{\pi }(\{0\})={\frac {1}{e}}}$, können wir die Definition als Rekursionsgleichung nutzen und folgern, dass ${\displaystyle P_{\pi }(\{i\})={\frac {1}{i!e}}}$ gilt. Dass die Berechnung der stationären Verteilung direkt möglich ist, ist aber die Ausnahme.

Stationäre Verteilungen haben zahlreiche Anwendungen in der Praxis. Stellt man sich die Markow-Kette als zufällig durch einen Graphen wandernden Punkt vor, so ist die i-te Komponente der stationäre Verteilung genau die relative Wahrscheinlichkeit, ihn im Zustand i anzutreffen. Ein Beispiel hierfür ist das Ehrenfest-Modell. Es modelliert die Diffusion von Molekülen durch eine Membran. Der stationäre Zustand ist dann genau die Konzentration, wenn sich ein Diffusionsgleichgewicht eingestellt hat. Ein anderes Beispiel ist die Berechnung des PageRanks mittels des Zufalls-Surfer-Modells. Hier modelliert die Markow-Kette das Verhalten eines Internetnutzers: Mit einer bestimmten Wahrscheinlichkeit wählt er einen der Links auf der Homepage, auf der er sich befindet, aus, andernfalls wählt er eine andere Homepage per Browser aus. Die stationäre Verteilung ist dann die relative Wahrscheinlichkeit, dass der Zufallssurfer auf eine bestimmte Seite trifft. Dies ist dann ein Maß für die Wichtigkeit dieser Seite und auch der normierte PageRank dieser Seite.

Des Weiteren spielen stationäre Verteilungen eine wichtige Rolle bei Markov-Chain-Monte-Carlo-Verfahren. Sie sind genau die Verteilungen, für die eine Stichprobe erstellt werden soll.

• Ulrich Krengel: Einführung in die Wahrscheinlichkeitstheorie und Statistik. 8. Auflage, Vieweg, 2005. ISBN 978-3-8348-0063-3.
• Hans-Otto Georgii: Stochastik: Einführung in die Wahrscheinlichkeitstheorie und Statisti. 4. Auflage, de Gruyter, 2009. ISBN 978-3-11-021526-7.
• Christian Hesse: Angewandte Wahrscheinlichkeitstheorie. Eine fundierte Einführung mit über 500 realitätsnahen Beispielen und Aufgaben, Vieweg, Braunschweig/Wiesbaden 2003, ISBN 978-3-528-03183-1.
• Achim Klenke: Wahrscheinlichkeitstheorie. 2. Auflage. Springer-Verlag, Berlin/Heidelberg 2008, ISBN 978-3-540-76317-8.