Brownian motion


Brownian motion
Zwei unabhängige Standard-Wiener-Prozesse

Ein Wiener-Prozess ist ein zeitstetiger stochastischer Prozess, der normalverteilte, unabhängige Zuwächse hat. Benannt wurde der Prozess, der auch als Brownsche Bewegung bekannt ist, nach dem amerikanischen Mathematiker Norbert Wiener. Seit der Einführung der Stochastischen Analysis (stochastische Integration, stochastische Differentialgleichungen) durch Itō Kiyoshi in den 1940er Jahren spielt der Wiener-Prozess die zentrale Rolle im Kalkül der zeitstetigen Stochastischen Prozesse und wird in zahllosen Gebieten der Natur- und Wirtschaftswissenschaften als Grundlage zur Simulation zufälliger Entwicklungen herangezogen.

Inhaltsverzeichnis

Definition

Ein stochastischer Prozess (W_t)_{t \in \mathbb{R}_{+}} auf dem Wahrscheinlichkeitsraum (Ω,F,P), der adaptiert an die Filtrierung (\mathcal{F}_t)_{t\in\mathbb{R}_{+}} ist, heißt (Standard-)Wiener-Prozess, wenn die vier folgenden Bedingungen gelten:

  1. W0 = 0 (P-fast sicher).
  2. Für gegebenes s mit t>s\geq0 sind alle Inkremente WtWs von \mathcal{F}_s stochastisch unabhängig. Der Wiener-Prozess hat also unabhängige Zuwächse.
  3. W_t-W_s \;\sim\; \mathcal{N} \left(0,t-s \right) \;\; \forall \; 0\leq s < t. Die Zuwächse sind also stationär und normalverteilt mit Mittelwert Null und Varianz t-s.
  4. Die einzelnen Pfade sind (P-)fast sicher stetig.

Der vierte Punkt kann auch aus der Definition insofern gestrichen werden, als dass sich mit dem Stetigkeitssatz von Kolmogorow−Centsow zeigen lässt, dass es unter den obigen Voraussetzungen immer eine fast sicher stetige Version des Prozesses gibt.

Umgekehrt lässt sich ein Wiener-Prozess (W_t)_{t \in \mathbb{R}_{+}} nach Lévy durch folgende zwei Eigenschaften charakterisieren:

  1. Wt ist ein stetiges lokales Martingal mit W0 = 0
  2. {W_t}^2-t ist ein Martingal

Eigenschaften

Einordnung

Eigenschaften der Pfade

Selbstähnlichkeiten

  • Auch das Negative eines Standard-Wiener-Prozesses, also (-W_t), \;t\ge0, ist ein Standard-Wiener-Prozess. Allgemeiner gilt auch das Reflexionsprinzip: Ein ab einer beliebigen Stoppzeit τ gespiegelter Wiener-Prozess ist wieder ein Wiener-Prozess. Der gespiegelte Prozess \hat{W}_t ist dabei wie folgt definiert: \hat{W}_t(\omega)=W_t(\omega) falls t\le \tau(\omega) und \hat{W}_t(\omega)=2W_{\tau(\omega)}(\omega)-W_t(\omega) falls t >\, \tau(\omega).
  • Der Wiener-Prozess ist selbstähnlich unter Streckung der Zeitachse, d.h. X_t=\tfrac{1}{\sqrt{\alpha}} \cdot W_{\alpha \cdot t} ist für jedes α > 0 ein Standard-Wiener-Prozess.
  • Inversion der Zeitachse: auch X_t=t W_{\frac{1}{t}} ist ein Standard-Wiener-Prozess
  • Verschiebung der Zeitachse: für jedes deterministische  s\ge 0 ist der stochastische Prozess Xt = Wt + sWs ebenfalls ein Wiener-Prozess; hier werden die Zuwächse vom Zeitpunkt s an betrachtet, d.h. W erfüllt die schwache Markow-Eigenschaft.

Generator

Eine Brownsche Bewegung auf der Kugel. Der Generator dieses Prozesses ist ½ mal der Laplace-Beltrami Operator auf einer Mannigfaltigkeit, hier einer Kugeloberfläche.

Der Generator eines eindimensionalen Wiener-Prozesses ist ½ mal der Laplace-Beltrami Operator. Damit lässt sich ein Wiener-Prozess auch auf anderen Mannigfaltigkeiten wie z.B. auf einer Kugel (siehe Bild) definieren.

Verallgemeinerter Wiener-Prozess

Ist (Wt) ein Standard-Wiener-Prozess, so nennt man den stochastischen Prozess

Xt = μt + σWt

Brownsche Bewegung mit Drift μ und Volatilität σ. Damit lassen sich auch stochastische Prozesse darstellen, die tendenziell eher fallen (μ < 0) oder tendenziell eher steigen (μ > 0). Dabei gilt

X_t -X_s \sim \mathcal{N}\left(\mu (t-s), \sigma^2 \left(t-s\right) \right).

Auch allgemeine Wiener-Prozesse sind Markow- und Lévy-Prozesse, aber die Martingaleigenschaft gilt nur noch in abgeschwächter Form:

Ist μ < 0, so ist Xt ein Supermartingal, ist μ > 0, so ist Xt ein Submartingal. Für μ = 0 ist Xt ein Martingal.

Der mehrdimensionale Fall

75 Schritte einer diskreten Annäherung an eine zweidimensionale Brownsche Bewegung. So ähnlich könnte sich auch das Partikel unter Browns Mikroskop bewegt haben.

Ein mehrdimensionaler stochastischer Prozess  W_t= (W_{1,t}, W_{2,t},\ldots W_{n,t})^T, \; t\ge 0 heißt n-dimensionaler (standard-)Wiener-Prozess oder n-dimensionale Brownsche Bewegung, falls die Koordinaten (Wi,t) unabhängige (Standard-)Wiener-Prozesse sind. Die Zuwächse Wi,tWi,s sind dann ebenfalls unabhängig und  \mathcal{N}(0,(t-s)I_n) -verteilt (n-dimensionale Normalverteilung), wobei In die Einheitsmatrix der Dimension n ist.

Der n-dimensionale Wiener-Prozess hat eine besonders schöne Eigenschaft, die ihn von den meisten anderen mehrdimensionalen Prozessen abhebt und die ihn für die Modellierung des Brownschen Partikels prädestiniert: Er ist invariant unter Drehungen der Koordinatenachsen. Das bedeutet, dass für jede orthogonale Matrix  Q \in \mathbb {R}^{n \times n} der gedrehte (oder gespiegelte) Prozess  X_t:=QW_t, \; t\ge 0 genau dieselbe Verteilung wie Wt besitzt.

Genau wie die eindimensionale Brownsche Bewegung kann man nun auch die n-dimensionale verallgemeinern: für jeden Vektor  \mu \in \mathbb{R}^n und jede Matrix  A \in \mathbb{R}^{n \times n} wird durch

 X_t := \mu t + A W_t, \; t \ge 0

eine Brownsche Bewegung mit Drift μ und Varianz AAT definiert. Dementsprechend gilt  X_t \sim \mathcal{N}(\mu t, t AA^T) . Hierbei können die einzelnen Koordinaten also auch miteinender korreliert sein.

Zusammenhang zu anderen stochastischen Prozessen

  • Ist  (X_t),\; t\ge 0 eine geometrische Brownsche Bewegung, so ist Wt: = ln(Xt) eine Brownsche Bewegung (mit Drift). Andererseits kann man aus jedem Wiener-Prozess  (X_t),\; t\ge 0 mit Drift μ und Volatilität σ durch  Y_t:=e^{X_t-\frac{\sigma^2 t}{2}} eine geometrische Brownsche Bewegung gewinnen.
  • Mit Hilfe des stochastischen Integralbegriffes von Itô lässt sich der Wiener-Prozess zum Itō-Prozess verallgemeinern.
  • Der symmetrische Random Walk kann als zeitdiskretes Pendant zum Wiener-Prozess angesehen werden, denn es gilt der folgende Konvergenzsatz: ist für  n \in \mathbb{N} der Random Walk  (R_t), \; t\in T auf dem diskreten Zeitgitter  T= \left\{ 0, \tfrac{1}{n}, \tfrac{2}{n}, \ldots \right\} so definiert, dass R0 = 0 gilt und R sich in jedem Zeitschritt mit Wahrscheinlichkeit ½ um  \sqrt{\tfrac{1}{n}} nach oben und mit Wahrscheinlichkeit ½ um \sqrt{\tfrac{1}{n}} nach unten bewegt, so konvergiert R für n \to \infty gegen einen standard-Wiener-Prozess (Invarianzprinzip von Donsker).
  • Ist  (W_t),\; t\ge 0 ein standard-Wiener-Prozess und T > 0, so ist  B_t:=W_t-\tfrac{t}{T}W_T ,\ 0\leq t\leq T eine Brownsche Brücke.

Geschichte des Wiener-Prozesses

Thorvald N. Thiele

1827 beobachtete der schottische Botaniker Robert Brown unter dem Mikroskop, wie Pflanzenpollen sich in einem Wassertropfen unregelmäßig hin- und herbewegten (daher der Name Brownsche Bewegung). Doch es waren nicht Naturwissenschaftler, die die Entwicklung des mathematischen Modells vorantrieben: 1880 beschrieb der Statistiker und Astronom Thorvald Nicolai Thiele (1838-1910) in Kopenhagen erstmals einen solchen „Prozess“ (die Theorie der stochastischen Prozesse war damals allerdings noch nicht entwickelt), als er wirtschaftliche Zeitreihen und die Verteilung von Residuen bei der Methode der kleinsten Quadrate studierte.

1900 griff der französische Mathematiker Louis Bachelier (1870-1946), ein Schüler Henri Poincarés, Thieles Idee auf, als er versuchte, die Kursbewegungen an der Pariser Börse zu analysieren. Beide Ansätze hatten letztendlich nur geringen Einfluss auf die zukünftige Entwicklung des Prozesses, zum Teil wohl aus dem Grunde, dass Finanzmathematik zu diesem Zeitpunkt eine untergeordnete Rolle in der Mathematik jener Zeit spielte. Heute jedoch gilt gerade die Finanzmathematik als Hauptanwendungsgebiet von Wiener-Prozessen. Dennoch bevorzugte zum Beispiel der Stochastiker William Feller die Bezeichnung Bachelier-Wiener-Prozess.

Louis Bachelier

Der Durchbruch kam jedoch, als Albert Einstein 1905 in seinem annus mirabilis [1], offenbar ohne Kenntnis von Bacheliers Arbeiten und unabhängig von ihm Marian Smoluchowski (1906), [2] den Wiener-Prozess in seiner heutigen Gestalt definierten. Einsteins Motivation war es, die Bewegung der brownschen Partikel durch die molekulare Struktur des Wassers zu erklären – ein Ansatz, der damals äußerst kontrovers war, heute aber unbestritten ist – und diese Erklärung mathematisch zu untermauern. Interessanterweise forderte er dabei eine weitere, physikalisch sinnvolle Eigenschaft, die Rektifizierbarkeit der Zufallspfade, für sein Modell nicht. Obwohl dies bedeutet, dass die Partikel in jeder Sekunde eine unendlich lange Strecke zurücklegen (was das gesamte Modell theoretisch disqualifiziert), bedeutete der einsteinsche Ansatz den Durchbruch sowohl für die molekulare Theorie als auch für den stochastischen Prozess.

Einen Beweis für die wahrscheinlichkeitstheoretische Existenz des Prozesses blieb Einstein allerdings schuldig. Dieser gelang erst 1923 dem US-Amerikanischen Mathematiker Norbert Wiener, der dabei neue Hilfsmittel von Lebesgue und Borel auf dem Gebiet der Maßtheorie ausnutzen konnte. Dennoch war sein Beweis so lang und kompliziert, dass ihn wohl nur eine handvoll Zeitgenossen verstehen konnten. Von Itō Kiyoshi ist überliefert, dass er einige seiner größten Fortschritte bei der Entwicklung des stochastischen Integrals bei dem Versuch erreichte, Wieners Arbeit nachzuvollziehen.

Letztendlich war es auch Itō, der dem Wiener-Prozess den Weg von der Physik in andere Wissenschaften ebnete: durch die von ihm aufgestellten stochastischen Differentialgleichungen konnte man die Brownsche Bewegung an mehr statistische Probleme anpassen. Bacheliers Ansatz scheiterte letztendlich daran, dass der Wiener-Prozess, unabhängig von seinem Startwert, im Laufe der Zeit fast sicher einmal negative Werte erreicht, was für Aktien unmöglich ist. Doch die durch eine stochastische Differentialgleichung abgeleitete geometrische Brownsche Bewegung löst dieses Problem und gilt seit der Entwicklung des berühmten Black-Scholes-Modells als Standard. Heute werden in praktisch allen Natur- und vielen Sozialwissenschaften brownsche Bewegungen und verwandte Prozesse als Hilfsmittel verwendet.

Simulation von Brownschen Pfaden

Um mit Hilfe von Zufallszahlen Pfade eines Wiener-Prozesses zu simulieren, stehen verschiedene Methoden zur Verfügung, die allesamt auf verschiedenen Eigenschaften des Prozesses aufbauen:

Einfacher Random Walk

Die einfachste Möglichkeit besteht darin, die oben erwähnte Konvergenz des einfachen Random Walk gegen einen Wiener-Prozess auszunutzen. Dazu muss man lediglich Bernoulli-verteilte Zufallsvariablen  B_1,B_2,B_3 \ldots simulieren, die untereinander unabhängig sind und jeweils mit 50-prozentiger Wahrscheinlichkeit die Werte 1 und −1 annehmen. Dann kann man zu einer vorgegebenen Schrittweite  \Delta t \ge 0 einen Wiener-Prozess an den Stellen 0, \;\Delta t,\; 2\,\Delta t,\; 3\,\Delta t \ldots durch

 W_{n\Delta t} \approx \sqrt{\Delta t}\sum_{i=1}^n B_i

approximieren. Der Vorteil dieser Methode liegt darin, dass nur sehr einfach herzustellende bernoulli-verteilte Zufallsvariablen benötigt werden. Allerdings handelt es sich nur um eine Approximation: Das Resultat ist kein Gauß-Prozess, sondern hat quasi binomialverteilte Zustände (genauer gesagt ist  \sum_{i=1}^n B_i + n binomial(n,0.5)-verteilt). Um die Normalverteilung hinreichend gut anzunähern, muss Δt deshalb sehr klein gewählt werden. Diese Methode ist deshalb nur zu empfehlen, wenn man den Prozess ohnehin auf einem sehr feinen Zeitgitter simulieren möchte.

Gaußscher Random Walk

Die folgende Methode ist dem einfachen Random Walk überlegen (sofern kein besonders feines Zeitgitter benötigt wird), da sie den Prozess exakt simuliert (d.h. die resultierenden Zustände stimmen in Verteilung mit denen eines Wiener-Prozesses überein):

 W_{n\Delta t} \approx \sqrt{\Delta t}\sum_{i=1}^n Z_i ,

wobei  Z_1, \; Z_2, \; Z_3\ldots unabhängige, standardnormalverteilte Zufallszahlen sind (beispielsweise erzeugt durch die Polar-Methode von Marsaglia). Diese als Gaußscher Random Walk bezeichnete Diskretisierung ist nur dann von Nachteil, wenn die vorhandenen normalverteilten Zufallsvariablen nicht von gleichmäßiger „Qualität“ sind. Wenn zum Beispiel Quasi-Zufallszahlen verwendet werden, weisen spät auftretende Zahlen bisweilen Abhängigkeitsstrukturen auf, die das Ergebnis verzerren können. In einem solchen Fall ist eine der folgenden Methoden vorzuziehen:

Brownsche Brücke

Die ersten fünf Halbierungsschritte der Brownschen Brücke. Die jeweils neu simulierte Iteration ist rot eingezeichnet

Diese Methode (die nur am Rande etwas mit dem gleichnamigen stochastischen Prozess zu tun hat) nutzt die Kovarianzstruktur des Wiener-Prozesses aus und legt ein höheres Gewicht auf frühe standardnormalverteilte Zufallsvariablen  Z_1, Z_2, \ldots .

Hier wird zuerst W1, welches normalverteilt mit Varianz 1 ist, durch W1 = Z1 simuliert. Nun wird das Intervall [0,1] schrittweise halbiert und folgender Schritt wiederholt:

 W_{\frac{1}{2}} ergibt sich als arithmetisches Mittel  \frac{1}{2}(W_0+W_1) plus eine weitere Normalverteilte Zufallsvariable, um die Varianz zu korrigieren. Also:

 W_{\frac{1}{2}}=\frac{1}{2}(W_0+W_1)+\frac{1}{2}Z_2

Analog:

\begin{align}
W_{\frac{1}{4}}&amp;amp;amp;=\frac{1}{2}(W_0+W_{\frac{1}{2}})+\frac{1}{\sqrt{8}}Z_3 \\
W_{\frac{3}{4}}&amp;amp;amp;=\frac{1}{2}(W_{\frac{1}{2}}+W_1)+\frac{1}{4}Z_4
\end{align}

und so weiter. Die Faktoren \tfrac{1}{2},\; \tfrac{1}{\sqrt{8}},\;\tfrac{1}{4}\ldots verringern sich dabei in jedem Halbierungsschritt um den Faktor \sqrt{2} und sorgen dafür, dass die Zustände die richtige Varianz erhalten.

Um einen Wiener-Prozess statt auf [0,1] auf ein beliebiges Intervall [0,a] auszuweiten, kann man nun die oben beschriebene Transformation  X_t=\sqrt{a} W_{\frac{t}{a}} anwenden; X ist dann ein Wiener-Prozess auf [0,a].

Spektralzerlegung

Bei der Spektralzerlegung wird der Wiener-Prozess in einer Art stochastischer Fourieranalyse als trigonometrische Polynome mit zufälligen Koeffizienten approximiert. Sind  Z_0, \;Z_1,\; Z_2\ldots unabhängig und standardnormalverteilt, so konvergiert die Reihe

S(t) = Z_0 \cdot t + \sum_{k=1}^\infty Z_k ~ \frac{\sqrt{2} \cdot \sin(k \pi t)}{k \pi}

gegen einen Wiener-Prozess. Diese Methode konvergiert bezüglich der  \mathcal{L}^2-Norm zwar mit maximaler Geschwindigkeit, beinhaltet aber im Gegensatz zur brownschen Brücke viele aufwändige trigonometrische Funktionsauswertungen. Daher findet sie, vor allem in der Monte-Carlo-Simulation, weniger oft Anwendung.

Geometrie

Die ein- und zwei-dimensionale Brownsche Bewegung ist rekurrent, in allen höheren Dimensionen ist sie transient. (Pólya: „Ein betrunkener Mann findet immer heim, ein betrunkener Vogel nicht.“) Siehe auch Markow-Kette.

Quellen

  1. Einstein, Albert: Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen.. In: Annalen der Physik. 17, 1905, S. 549-560. 
  2. Smoluchowski, M.: Zur kinetischen Theorie der Brownschen Molekularbewegung und der Suspensionen. In: Annalen der Physik. 21, 1906, S. 756-780. 

Literatur

  • Andrei N. Borodin, Paavo Salminen: Handbook of Brownian Motion - Facts and Formulae, Birkhäuser, Basel 2002, ISBN 3-7643-6705-9
  • Ioannis Karatzas, Steven E. Shreve: Brownian Motion and Stochastic Calculus (Graduate Texts in Mathematics), Springer, New York 1997, ISBN 0-3879-7655-8
  • Bernt Øksendal: Stochastic Differential Equations. An Introduction with Applications, Springer Berlin 2003, ISBN 3-540-04758-1
  • Philip E. Protter: Stochastic Integration and Differential Equations, Springer, Berlin 2003, ISBN 3-540-00313-4
  • John Michael Steele: Stochastic Calculus and Financial Applications, Springer, New York 2000, ISBN 0-387-95016-8

Siehe auch

Ornstein-Uhlenbeck-Prozess, Brownsche Brücke, Stochastische Differentialgleichung, Wurzel-Diffusionsprozess, geometrische Brownsche Bewegung


Wikimedia Foundation.

Schlagen Sie auch in anderen Wörterbüchern nach:

  • Brownian motion — Brownian Brown i*an, a. Pertaining to Dr. Robert Brown, who first demonstrated (about 1827) the commonness of the motion described below. [1913 Webster] {Brownian motion}, {Brownian movement}, the peculiar, rapid, vibratory movement exhibited by… …   The Collaborative International Dictionary of English

  • Brownian Motion —  Brownian Motion  Броуновское движение   Беспорядочное движение микроскопических, взвешенных в жидкости (или газе) частиц (броуновские частицы) твёрдого вещества (пылинки, крупинки взвеси, частички пыльцы растения и так далее), происходящее под… …   Толковый англо-русский словарь по нанотехнологии. - М.

  • Brownian motion — ► NOUN Physics ▪ the erratic movement of microscopic particles in a fluid, as a result of collisions with the surrounding molecules. ORIGIN named after the Scottish botanist Robert Brown (1773 1858) …   English terms dictionary

  • Brownian motion — This article is about the physical phenomenon; for the stochastic process, see Wiener process. For the sports team, see Brownian Motion (Ultimate). For the mobility model, see Random walk. Brownian motion (named after the botanist Robert Brown)… …   Wikipedia

  • Brownian motion — Any of various physical phenomena in which some quantity is constantly undergoing small, random fluctuations. It was named for Robert Brown, who was investigating the fertilization process of flowers in 1827 when he noticed a rapid oscillatory… …   Universalium

  • Brownian motion — Brown·ian motion .brau̇ nē ən n a random movement of microscopic particles suspended in liquids or gases resulting from the impact of molecules of the fluid surrounding the particles called also Brownian movement Brown brau̇n Robert (1773 1858)… …   Medical dictionary

  • Brownian motion — Brauno judesys statusas T sritis Standartizacija ir metrologija apibrėžtis Netvarkingas skystyje arba dujose pakibusios mikroskopinės dalelės judesys. atitikmenys: angl. Brownian motion; Brownian movement vok. Brownsche Bewegung, f rus.… …   Penkiakalbis aiškinamasis metrologijos terminų žodynas

  • Brownian motion — Brauno judesiai statusas T sritis chemija apibrėžtis Skystyje arba dujose pakibusių mikroskopinių dalelių judesiai. atitikmenys: angl. Brownian motion; Brownian movement rus. броуновские движения …   Chemijos terminų aiškinamasis žodynas

  • Brownian motion — Brauno judėjimas statusas T sritis fizika atitikmenys: angl. Brownian agitation; Brownian motion vok. Brownsche Bewegung, f rus. броуновское движение, n pranc. agitation brownienne, f; mouvement brownien, m; mouvement de Brown, m …   Fizikos terminų žodynas

  • Brownian motion — noun the random motion of small particles suspended in a gas or liquid • Syn: ↑Brownian movement, ↑pedesis • Hypernyms: ↑movement, ↑motion …   Useful english dictionary