Zufallsfeld

Ein Zufallsfeld, auch zufälliges Feld[1], engl. random field, wird benötigt, wenn man zufallsbeeinflusste Phänomene im Raum modellieren will, z. B. den Kohlendioxidgehalt in der Atmosphäre in Ballungsräumen, oder die Niederschlagsmenge in verschiedenen Regionen Deutschlands.

Mathematische Definition

Ein Zufallsfeld ist eine Familie ( Y ( x , ) ) x X {\displaystyle (Y(x,\cdot ))_{x\in {\mathcal {X}}}} mit X R d {\displaystyle {\mathcal {X}}\subseteq \mathbb {R} ^{d}} von reellwertigen Zufallsvariablen (Zufallsgrößen) auf einem Wahrscheinlichkeitsraum ( Ω , Σ , P ) {\displaystyle (\Omega ,\Sigma ,P)} . Dabei ist d {\displaystyle d} eine natürliche Zahl.

Im Fall d = 1 {\displaystyle d=1} spricht man von einem stochastischen Prozess, dann spielt x {\displaystyle x} häufig die Rolle des Zeitparameters und wird dann in der Regel mit t {\displaystyle t} bezeichnet. In einem solchen Fall bezeichnet die Indexmenge X {\displaystyle {\mathcal {X}}} eine Menge von Zeitpunkten. In den Fällen d = 2 {\displaystyle d=2} und d = 3 {\displaystyle d=3} bezeichnet x {\displaystyle x} häufig die Koordinaten eines Ortes und die Indexmenge X {\displaystyle {\mathcal {X}}} ist eine Menge von Orten. Für ein realisiertes ω Ω {\displaystyle \omega \in \Omega } ist dann x Y ( x , ω ) = y ( x ) {\displaystyle x\mapsto Y(x,\omega )=y(x)} für x X {\displaystyle x\in {\mathcal {X}}} eine Realisierung des Zufallsfeldes, die auch Trajektorie oder Pfade heißt. Häufig interessiert man sich z. B. für die Wahrscheinlichkeit, dass eine Trajektorie einen bestimmten Wert übersteigt, siehe z. B.[2] (Anwendungsbeispiel: Hochwasserschutz)

Trend, Kovarianz, Stationarität, Isotropie

Die Erwartungswertfunktion

E Y ( x , ) = m ( x ) , x X {\displaystyle EY(x,\cdot )=m(x),\quad x\in {\mathcal {X}}}

wird als Trend und die Zweite-Moment-Funktion

E [ Y ( x , ) m ( x ) ] [ Y ( z , ) m ( z ) ] = r ( x , z ) , x , z X {\displaystyle E[Y(x,\cdot )-m(x)][Y(z,\cdot )-m(z)]=r(x,z),\quad x,z\in {\mathcal {X}}}

als Kovarianzfunktion des Zufallsfeldes bezeichnet. Zufallsfelder, für die Trend und Kovarianzfunktion existieren und endlich sind, heißen Zufallsfelder 2. Ordnung. Zufallsfelder mit konstantem Trend und verschiebungsinvarianter Kovarianzfunktion, d. h.

m ( x ) = m , r ( x , z ) = r ( x z ) {\displaystyle m(x)=m,\quad r(x,z)=r(x-z)}

nennt man stationär im weiteren Sinn. Stationarität im engeren Sinn erfordert Verschiebungsinvarianz nicht nur der ersten beiden Momente, sondern aller endlichdimensionalen Verteilungen des zufälligen Feldes. Ist die Kovarianzfunktion rotationsinvariant, d. h.

r ( x , z ) = r ( x z ) , {\displaystyle r(x,z)=r(\|x-z\|),\quad \|\cdot \|} Euklidischer Abstand,

dann nennt man das Zufallsfeld isotrop.

Vorhersage von Werten des Zufallsfeldes

Hat man das Zufallsfeld an den Orten x 1 , . . . , x n {\displaystyle x_{1},...,x_{n}} beobachtet mit den Resultaten y 1 , . . . , y n {\displaystyle y_{1},...,y_{n}} , so kann man daraus eine Vorhersage Y ^ ( x ) {\displaystyle {\hat {Y}}(x)} des Zufallsfeldes an einer nichtbeobachteten Stelle x {\displaystyle x} konstruieren. Die beste Vorhersage, die den mittleren quadratischen Fehler minimiert, ist die bedingte Erwartung von Y ( x , ) {\displaystyle Y(x,\cdot )} , gegeben die Beobachtungen y 1 , . . . , y n {\displaystyle y_{1},...,y_{n}} , d. h.

Y ^ ( x ) = E [ Y ( x , ) | Y 1 = y 1 , . . . , Y n = y n ] {\displaystyle {\hat {Y}}(x)=E[Y(x,\cdot )|Y_{1}=y_{1},...,Y_{n}=y_{n}]} .

Diese Vorhersage lässt sich ohne weitere Verteilungsannahmen nicht berechnen. Bei bekannter Kovarianzfunktion des Zufallsfeldes lässt sich jedoch mit wenig Aufwand die beste lineare Vorhersage berechnen.

Anwendung in der Geostatistik

In der Geostatistik wird in der Regel anstatt der Kovarianzfunktion in äquivalenter Weise das Semivariogramm

V ( x , z ) = 1 2 E [ Y ( x , ) Y ( z , ) ] 2 {\displaystyle V(x,z)={\frac {1}{2}}E[Y(x,\cdot )-Y(z,\cdot )]^{2}} ,

d. h. der halbe (semi) Erwartungswert der quadratischen Differenz Y ( x ) Y ( z ) {\displaystyle Y(x)-Y(z)} , benutzt. Die beste lineare Vorhersage heißt in geostatistischer Terminologie Kriging, siehe z. B.[3] Die Dimension d {\displaystyle d} des Zufallsfeldes ist hier in der Regel auf natürliche Weise gegeben, z. B. d = 2 {\displaystyle d=2} für die Oberflächentemperatur eines Sees, d = 3 {\displaystyle d=3} für Probleme der Lagerstättenerkundung im Bergbau, d = 4 {\displaystyle d=4} für raum-zeitliche Phänomene in der Meteorologie.

Literatur

  • R. Adler, J. Taylor: Random Fields and Geometry. Springer, New York 2007.
  • N. Cressie: Statistics for Spatial Data. World Scientific, Singapore 2007.
  • P. H. Müller (Hrsg.): Lexikon der Stochastik – Wahrscheinlichkeitsrechnung und mathematische Statistik. 5. Auflage. Akademie-Verlag, Berlin 1991, ISBN 978-3-05-500608-1, Zufälliges Feld, S. 513–516. 
  • E. Vanmarcke: Random Fields: Analysis and Synthesis. World Scientific, Singapore 2010.

Einzelnachweise

  1. P. H. Müller (Hrsg.): Lexikon der Stochastik – Wahrscheinlichkeitsrechnung und mathematische Statistik. 5. Auflage. Akademie-Verlag, Berlin 1991, ISBN 978-3-05-500608-1, Zufälliges Feld, S. 513–516. 
  2. R. J. Adler: The Geometry of Random Fields. Wiley, Chichester/ New York/ Toronto 1981.
  3. J. P. Chiles, P. Delfiner: Geostatistics: Modelling Spatial Uncertainty. Wiley, New York 1999.