In der Mathematik, besonders der Wahrscheinlichkeitstheorie und Statistik, dient eine Übergangsmatrix (auch Prozessmatrix oder stochastische Matrix) dazu, die Übergangswahrscheinlichkeiten von (diskreten und kontinuierlichen) Markow-Ketten auszudrücken. Dadurch lassen sich künftige Entwicklungen vorausberechnen. In der Theorie der Markow-Ketten werden auch unendlichdimensionale Übergangsmatrizen definiert. In diesem Artikel werden jedoch nur Matrizen im Sinne der Linearen Algebra behandelt.
Eine Übergangsmatrix ist eine quadratische Matrix, deren Zeilen- oder Spaltensummen Eins betragen und deren Elemente zwischen Null und Eins liegen.
Prozessmatrizen dienen ebenfalls zur künftigen Berechnung dynamischer Entwicklungen. Im Gegensatz zu stochastischen Matrizen müssen sie jedoch keine Zeilen- bzw. Spaltensummen von 1 haben. Sie sind jedoch wie die stochastische Matrix quadratisch.
Weitere Unterscheidung
Eine Übergangsmatrix heißt zeilenstochastisch, wenn alle Einträge der Matrix zwischen 0 und 1 liegen und die Zeilensummen 1 ergeben. Eine Übergangsmatrix heißt spaltenstochastisch, wenn alle Einträge der Matrix zwischen 0 und 1 liegen und die Spaltensummen 1 ergeben. Eine Übergangsmatrix heißt doppelt-stochastisch, wenn sie sowohl zeilen- als auch spaltenstochastisch ist.
Äquivalent ist die folgende Definition: Eine Matrix heißt zeilen-(spalten-)stochastisch, wenn sie zeilen-(spalten-)weise aus Wahrscheinlichkeitsvektoren besteht.
Teilweise werden Matrizen mit Einträgen zwischen 0 und 1, deren Zeilensummen (bzw. Spaltensummen) kleiner als 1 sind, auch als substochastisch bezeichnet. In der Stochastik sind fast ausschließlich zeilenstochastische Matrizen gebräuchlich. Die Unterscheidung ist aber i. A. wenig gebräuchlich, da die Matrizen durch Transponierung ineinander übergehen.
Eigenschaften Eigenwerte und Eigenvektoren Den Eigenwerten und Eigenvektoren einer stochastischen Matrix kommt in der Stochastik eine besondere Rolle zu. Ist Eigenvektor zum Eigenwert , entspricht er einer stationären Verteilung der Markow-Kette (vgl. unten). Generell besitzt jede stochastische Matrix den Eigenwert 1. Ist z. B. zeilenstochastisch, so folgt mit der Zeilensummennorm, dass . Da der Spektralradius einer Matrix immer höchstens so groß wie ihre Norm ist, müssen alle Eigenwerte betragsmäßig kleiner oder gleich 1 sein. Ist nun ein Einsvektor (d. h. ein Vektor mit nur 1 als Einträgen), so gilt und 1 ist Eigenwert von . Der Beweis für spaltenstochastische Matrizen läuft analog, aber mit der Spaltensummennorm anstelle der Zeilensummennorm. Daraus folgt direkt, dass 1 auch immer betragsgrößter Eigenwert ist. Des Weiteren ist 1 auch immer ein halbeinfacher Eigenwert. Die Dimension des Eigenraumes lässt sich etwas schwerer berechnen. Mit dem Satz von Perron-Frobenius folgt: Ist die stochastische Matrix irreduzibel, so ist die Dimension des zum Eigenwert 1 gehörenden Eigenraumes gleich 1. Dies ist insbesondere der Fall, wenn die Einträge einer stochastischen Matrix echt größer als 0 sind.
Konvexität, Normen und Abgeschlossenheit
Die Menge der Übergangsmatrizen ist konvex. Sind also und zeilen- bzw. spaltenstochastische Matrizen, so ist wieder eine zeilen- bzw. spaltenstochastische Matrix für alle .
Direkt aus der Definition folgt, dass die Zeilensummennorm einer zeilenstochastischen Matrix 1 ist, genauso wie die Spaltensummennorm einer spaltenstochastischen Matrix.
Außerdem sind Übergangsmatrizen abgeschlossen bezüglich der Matrixmultiplikation: Sind spalten- bzw. zeilenstochastische Matrizen, so ist ebenfalls eine spalten- bzw. zeilenstochastische Matrix.
Beispiel
Das charakteristische Polynom einer -Übergangsmatrix lässt sich sehr leicht berechnen.
Mit der Spur und der Determinante gilt:
Aus der letzten Zeile ergibt sich, dass stets Eigenwert der Matrix P ist, unabhängig von der Wahl der Koeffizienten von P. Die anderen beiden Eigenwerte lassen sich dann über die p-q-Formel errechnen.
Anwendung zur Charakterisierung diskreter Markow-Ketten
Ist eine zeilenstochastische Matrix, so lässt sich damit auf folgende Weise eine zeitinvariante Markow-Kette mit endlichem Zustandsraum charakterisieren:
Die Einträge der Matrix sind genau die Übergangswahrscheinlichkeiten vom Zustand in den Zustand :
Der Wahrscheinlichkeitsvektor (ein Zeilenvektor) beschreibt den Zustand des Systems zum Zeitpunkt 0 und wird manchmal Startvektor genannt. Dabei beschreibt die Wahrscheinlichkeit zum Zeitpunkt 0 im Zustand zu sein. Die Wahrscheinlichkeiten zum Zeitpunkt 1 ergeben sich durch Multiplikation der Matrix von links mit dem Zeilenvektor :
Allgemein ergeben sich die Wahrscheinlichkeiten zum Zeitpunkt durch Multiplikation der Matrix von links mit dem Zeilenvektor :
Die Wahrscheinlichkeiten zu einem beliebigen Zeitpunkt in Abhängigkeit vom Startzustand sind daher
Eine besondere Rolle kommt den Linkseigenvektoren der Matrix zum Eigenwert zu, denn diese stellen die stationären Verteilungen der Markow-Kette dar.
Für spaltenstochastische Matrizen kann man analog vorgehen. Diese werden von rechts mit Spaltenvektoren multipliziert und man findet den gewöhnlichen Eigenvektor. Alternativ kann man auch die Matrix transponieren und das oben skizzierte Vorgehen nutzen.
Ein anwendungsorientiertes Beispiel für diese Verwendung von Übergangsmatrizen ist die Berechnung des PageRank mittels der Google-Matrix. Jeder Zustand entspricht dort einer Webseite im World Wide Web, die Übergangswahrscheinlichkeiten geben an, mit welcher Wahrscheinlichkeit ein Nutzer auf einen Link klickt. Die Grenzverteilung ist dann die relative Häufigkeit, mit welcher der Nutzer auf eine Webseite stößt, und damit ein Maß für die Wichtigkeit dieser Seite.
Auch die Rechtseigenvektoren einer Übergangsmatrix zum Eigenwert 1 spielen eine Rolle bei der Untersuchung von Markow-Ketten. Bei passender Normierung sind diese genau die Absorptionswahrscheinlichkeiten in einem absorbierenden Zustand.
Des Weiteren finden sich auch viele Eigenschaften einer Markow-Kette in der Übergangsmatrix wieder: Existiert eine Zahl , so dass ist, dann ist der Zustand j vom Zustand i aus erreichbar. Gilt außerdem auch für ein , so kommunizieren die Zustände i und j. Im Fall einer homogenen Markow-Kette mit endlichem Zustandsraum ist der Begriff der irreduziblen Markow-Kette äquivalent zur Irreduzibilität der Übergangsmatrix. Ist für ein , so ist die Markow-Kette irreduzibel und aperiodisch und konvergiert demnach gegen eine Grenzverteilung. Dieses Kriterium ist oft leichter zu überprüfen als die Irreduzibilität und Aperiodizität separat zu zeigen. Die Chapman-Kolmogorow-Gleichung lässt sich im Falle einer Übergangsmatrix als komponentenweise ausgeschriebene Matrixmultiplikation interpretieren.
Beispiele In den folgenden drei Beispiele werden spaltenstochastische Übergangsmatrizen verwendet. Dabei bezeichnet ein Eintrag in der -ten Zeile und -ten Spalte einer Übergangsmatrix Matrix die Übergangswahrscheinlichkeit vom Zustand in den Zustand . (Im Unterschied dazu bezeichnet bei einer Modellierung mit zeilenstochastischen Übergangsmatrizen ein Eintrag die Übergangswahrscheinlichkeit vom Zustand in den Zustand .) Würfelspiel Kniffel (Yahtzee) Beim Würfelspiel Kniffel (Yahtzee) werden Punkte für verschiedene Kategorien erspielt. Im Folgenden wird nur eine Spielrunde und die wertvollste Kategorie betrachtet, die ebenfalls Kniffel oder Yahtzee genannt wird.
Bei der Kategorie Kniffel geht es darum, mit 5 Würfeln innerhalb von höchstens 3 Würfen 5 gleiche Augenzahlen zu erzielen. Beim ersten Wurf werden alle 5 Würfel geworfen. Beim zweiten und beim dritten Wurf können Würfel beliebig ausgewählt werden. Diese werden erneut geworfen.
Um mit möglichst großer Wahrscheinlichkeit einen Kniffel zu erzielen, ist folgende naheliegende Strategie optimal:
Vor dem zweiten und vor dem dritten Wurf wird jeweils die oder eine der am häufigsten vorkommenden Augenzahlen bestimmt. Die Würfel, die diese Augenzahl zeigen, werden behalten. Alle anderen Würfel werden erneut geworfen.
Beispiele: Nach dem Wurf (1, 1, 1, 2, 2) werden also die drei Würfel mit der Augenzahl 1 behalten und die zwei Würfel mit der Augenzahl 2 erneut geworfen. Nach dem Wurf (4, 4, 5, 5, 6) werden die zwei Würfel mit der Augenzahl 4 behalten und die drei anderen Würfel erneut geworfen. Alternativ kann auch die Augenzahl 5 ausgewählt und die drei anderen Würfel erneut geworfen werden.
Die möglichen Variationen aus den Augenzahlen der 5 Würfel lassen sich in folgende 5 Klassen (Zustände) einteilen, wobei A, B, C, D, E jeweils verschiedene Augenzahlen darstellen: fünf verschiedene Augenzahlen (ABCDE) einmal oder zweimal 2 gleiche Augenzahlen (entweder AABCD oder AABBC) 3 gleiche Augenzahlen (entweder AAABC oder AAABB) 4 gleiche Augenzahlen (AAAAB) 5 gleiche Augenzahlen (AAAAA) Die Wahrscheinlichkeiten, mit dem ersten Wurf eine Variation aus einer dieser 5 Klassen zu erzielen, können in dieser Reihenfolge als Startvektor dargestellt werden:
Sie werden jeweils berechnet als Quotient aus der Anzahl der günstigen Variationen und der Anzahl 65 = 7776 aller Variationen. Weil jede Augenzahl mehrfach vorkommen kann, sind es Variationen mit Wiederholung. Die Wahrscheinlichkeiten, nach dem Erzielen einer bestimmten gegebenen Variationsklasse eine bestimmte Variationsklasse mit dem zweiten Wurf zu erzielen, sind hier die Übergangswahrscheinlichkeiten . Diese können als Übergangsmatrix
dargestellt werden. Zum Beispiel ist die Wahrscheinlichkeit, nach dem gegebenen ersten Wurf mit fünf verschiedenen Augenzahlen (ABCDE) im zweiten Wurf einmal oder zweimal 2 gleiche Augenzahlen (entweder AABCD oder AABBC) zu erzielen. Zur Berechnung dieser Wahrscheinlichkeiten sind manchmal Fallunterscheidungen notwendig, denn die häufigste Augenzahl nach dem zweiten Wurf kann eine andere sein als die häufigste Augenzahl nach dem ersten Wurf. Dafür wird die oben beschriebene optimale Strategie für das Erzielen eines Kniffels zugrunde gelegt.
Daraus ergibt sich der Vektor für die Wahrscheinlichkeiten nach zwei Würfen, indem die Übergangsmatrix mit dem Startvektor multipliziert wird:
Die gleichen Betrachtungen können für die Berechnung der Wahrscheinlichkeiten nach drei Würfen verwendet werden und man erhält den Vektor
Als Verallgemeinerung erhält man den Vektor für die Wahrscheinlichkeiten nach Würfen:
Das letzte Element von ist die Wahrscheinlichkeit für einen Kniffel nach Würfen. Die ausgerechneten Werte für die ersten 10 Würfe sind: Anzahl der WürfeWahrscheinlichkeit für einen Kniffel (Yahtzee)10,0007720,0126330,0460340,1005850,1705160,2490870,3305080,4104490,48601100,55553 Nach dem zehnten Wurf ist die Wahrscheinlichkeit, einen Kniffel zu erzielen, also zum ersten Mal größer als .DataGenetics: Yahtzee Probability
Die Ratte im Zimmer mini|Der Adjazenzgraph der stochastischen Matrix aus dem Beispiel und damit auch ein Übergangsgraph Peter besitzt eine Ratte. Ist die Ratte nicht im Käfig eingesperrt, so befindet sie sich entweder unter dem Schreibtisch (Zustand 3), hinter dem Schrank (Zustand 2) oder ist im Käfig, um zu fressen (Zustand 1). Die Ratte wechselt alle 5 Minuten ihren Ort. Ist sie gerade im Käfig, so bleibt sie mit Wahrscheinlichkeit 0,05 dort, mit Wahrscheinlichkeit 0,4 geht sie hinter den Schrank und mit Wahrscheinlichkeit 0,55 unter den Schreibtisch. Ist sie hinter dem Schrank, so bleibt sie mit Wahrscheinlichkeit 0,7 dort, mit Wahrscheinlichkeit 0,2 geht sie unter den Schreibtisch und mit Wahrscheinlichkeit 0,1 geht sie in den Käfig. Ist sie unter dem Schreibtisch, so bleibt sie mit Wahrscheinlichkeit 0,1 dort, mit Wahrscheinlichkeit 0,1 geht sie in den Käfig und mit Wahrscheinlichkeit 0,8 flüchtet sie hinter den Schrank. Das Verhalten der Ratte wird durch die spaltenstochastische Matrix beschrieben:
Peter lässt nun seine Ratte frei und will wissen, mit welcher Wahrscheinlichkeit sich die Ratte nach 20 Minuten im Käfig befindet. Der Startzustand des Systems ist