Stochastische Differentialgleichung

Der Begriff der stochastischen Differentialgleichung (Abkürzung SDGL oder englisch SDE für stochastic differential equation) ist in der Mathematik eine Verallgemeinerung des Begriffs der gewöhnlichen Differentialgleichung auf stochastische Prozesse. Stochastische Differentialgleichungen werden in zahlreichen Anwendungen eingesetzt, um zeitabhängige Vorgänge zu modellieren, die neben deterministischen Einflüssen zusätzlich stochastischen Störfaktoren (Rauschen) ausgesetzt sind.

Die mathematische Formulierung des Problems stellte die Mathematiker vor große Probleme, und so wurde die formale Theorie der stochastischen Differentialgleichungen erst in den 1940er Jahren durch den japanischen Mathematiker Itō Kiyoshi formuliert. Gemeinsam mit der stochastischen Integration begründet die Theorie der stochastischen Differentialgleichungen die stochastische Analysis.

Stochastische Differentialgleichungen sind im Allgemeinen keine Differentialgleichungen und zu unterscheiden von den zufälligen Differentialgleichungen (RDE für {{Modul:Vorlage:lang}} Modul:Vorlage:lang:103: attempt to index field 'wikibase' (a nil value)). Letztere besitzen differenzierbare Pfade.

Von der Differential- zur Integralgleichung

Genau wie bei deterministischen Funktionen möchte man auch bei stochastischen Prozessen den Zusammenhang zwischen dem Wert der Funktion und ihrer momentanen Änderung (ihrer Ableitung) in einer Gleichung formulieren. Was im einen Fall zu einer gewöhnlichen Differentialgleichung führt, ist im anderen Fall problematisch, da viele stochastische Prozesse, wie beispielsweise der Wiener-Prozess, nirgends differenzierbar sind.

Jedoch lässt sich eine gewöhnliche Differentialgleichung

immer auch äquivalent als Integralgleichung

schreiben, die ohne explizite Erwähnung der Ableitung auskommt. Bei stochastischen Differentialgleichungen geht man nun den umgekehrten Weg, d. h., man definiert den Begriff mit Hilfe der zugehörigen Integralgleichung.

Die Formulierung

Gegeben sei ein filtrierter Wahrscheinlichkeitsraum <math>(\Omega,\mathcal{F},P,\mathbb{F}=(\mathcal{F}_t)_{t\geq 0})</math>, der die üblichen Bedingungen erfüllt (d. h. <math>\mathbb{F}</math> ist rechtsstetig und augmentiert), und mit <math>\operatorname{Lin}(\mathbb{R}^n,\R^d)</math> sei die Gesamtheit aller linearen Abbildungen von <math>\R^n</math> nach <math>\R^d</math> bezeichnet.

Itōsche Differentialgleichung

Seien zwei stetige Funktionen

<math>a\colon\mathbb{R}^d \times \mathbb{R}_{+} \to \mathbb{R}^d</math>

<math>\sigma\colon\R^d \times \R_{+} \to \operatorname{Lin}(\R^n,\R^d)</math>

sowie eine <math>\mathbb{F}</math>-adaptierte <math>n</math>-dimensionale brownsche Bewegung <math> (W_t)_{t \geq 0}</math> gegeben. Die dazugehörige stochastische Integralgleichung

<math> X_t=X_0 + \int_0^t a(X_s,s) \,{\rm d}s + \int_0^t \sigma(X_s,s)\,{\rm d}W_s </math>

wird durch Einführung der Differentialschreibweise

<math> {\rm d} X_t=a(X_t,t)\,{\rm d}t+\sigma(X_t,t)\,{\rm d}W_t</math>

zur stochastischen Differentialgleichung. Durch die Substitution <math>Y_t=(t,W_t)^T</math> lässt sich die Gleichung auch verkürzt aufschreiben. Das erste Integral ist als Lebesgue-Integral und das zweite als Itō-Integral zu lesen. Zu gegebenen Funktionen <math>a</math> und <math>\sigma</math> (auch als Drift und Diffusionskoeffizient bezeichnet) und einer brownschen Bewegung <math>(W_t)</math> wird hier also ein <math>d</math>-dimensionaler Prozess <math>(X_t)</math> gesucht, der die obige Integralgleichung erfüllt. Dieser Prozess ist dann eine Lösung der obigen SDGL. Um den infinitesimalen Generator zu berechnen, wendet man die Itō-Formel an und integriert dann.

Auf Mannigfaltigkeiten

Stochastische Differentialgleichungen können auch allgemeiner auf Mannigfaltigkeiten betrachtet werden, diese sind Untersuchungsgegenstand der stochastischen Differentialgeometrie.

Eine SDGL auf einer Mannigfaltigkeit <math>M</math> ist ein Paar <math>(A,Z)</math>, wobei

<math>Z=(Z_t)_{t\in\mathbb{R}_+}</math> ein stetiges Semimartingal auf einem endlichdimensionalen <math>\R</math>-Vektorraum <math>E</math> ist.
<math>A:M\times E\to TM</math> ein Homomorphismus von Vektorbündeln über <math>M</math>

<math>A:(x,e)\mapsto A(x)e</math>

ist, wobei <math>A(x):E\to TM</math> eine lineare Abbildung bezeichnet.

Die stochastische Differentialgleichung <math>(A,Z)</math> notieren wir als Fisk-Stratonowitsch-Integral

Existenz und Eindeutigkeit

Ist <math>A</math> eine beliebige, auf demselben Wahrscheinlichkeitsraum wie <math>W</math> definierte Zufallsvariable, so wird aus der obigen SDGL durch Hinzufügen der Bedingung <math> X_0=A </math> fast sicher ein stochastisches Anfangswertproblem als Pendant zum Anfangswertproblem für gewöhnliche Differentialgleichungen.

Auch zum Existenz- und Eindeutigkeitssatz von Picard und Lindelöf findet sich hier eine Entsprechung: wenn die folgenden drei Eigenschaften erfüllt sind:

<math> A \in L^2(P) </math>, d. h., <math>A</math> hat endliche Varianz.
Lipschitz-Bedingung: Es gibt eine Konstante <math>K \geq 0</math>, sodass für alle <math>x,y \in \R</math> und alle <math>t \geq 0</math> gilt

<math>|a(x,t)-a(y,t)|+|b(x,t)-b(y,t)| \le K|x-y|</math>.

Lineare Beschränktheit: Es gibt eine Konstante <math>C \geq 0</math>, sodass für alle <math>x \in \R</math> und alle <math>t \geq 0</math> gilt

<math>|a(x,t)| + |b(x,t)| \leq C(1+|x|)</math>.

Dann besitzt das Anfangswertproblem eine (bis auf fast sichere Gleichheit) eindeutige Lösung <math>X</math>, die zudem zu jedem Zeitpunkt <math>t</math> endliche Varianz besitzt.

Allgemeine Situation: lokale Lipschitz-Bedingung und Maximallösungen

Wir betrachten die allgemeine Form einer stochastischen Differentialgleichung

<math>\mathrm{d}Y_t=\alpha(t,Y_t)\mathrm{d}X_t</math>

wobei

<math>X</math> ein stetiges Semimartingal in <math>\R^n</math> und <math>Y</math> ein stetiges Semimartingal in <math>\R^d</math> ist,
<math>\alpha:\R_{+}\times U \to \operatorname{Lin}(\R^{n};\R^{d})</math> ist eine Abbildung von einer nichtleeren offenen Menge <math>U\subset \R^d</math>, wobei <math>\operatorname{Lin}(\R^{n};\R^{d})</math> der Raum aller linearen Abbildungen von <math>\R^{n}</math> nach <math>\R^{d})</math> ist.

Ob die Lösung der Gleichung explodiert oder nicht, hängt von der Wahl der Funktion <math>\alpha</math> ab. Deshalb führen wir nun eine lokale Lipschitz-Bedingung für <math>\alpha</math> ein. Für <math>t\geq 0</math>, eine kompakte Menge <math>K\subset U</math> und eine Konstante <math>L(t,K)</math> sei

<math>|\alpha(s,y)-\alpha(s,x)|\leq L(t,K)|y-x|,\quad x,y\in K,\;0\leq s\leq t,</math>

wobei <math>|\cdot|</math> die euklidische Norm bezeichnet. Diese Bedingung garantiert die Existenz und Eindeutigkeit einer Maximallösung.

Sei nun <math>\alpha</math> stetig und erfülle die oben genannte lokale Lipschitz-Bedingung, weiter sei <math>F:\Omega\to U</math> eine Initialbedingung, das heißt eine messbare Funktion bezüglich der Initial-σ-Algebra <math>\mathcal{F}_0</math>. Sei <math>\zeta:\Omega\to \overline{\R}_{+}</math> ein vorhersehbare Stoppzeit mit <math>\zeta>0</math> fast sicher. Ein <math>U</math>-wertiges Semimartingal <math>(Y_t)_{t<\zeta}</math> heißt Maximallösung von

<math>dY_t=\alpha(t,Y)dX_t,\quad Y_0=F</math>

mit Lebenszeit <math>\zeta</math>, falls

für eine (und somit für alle) ankündigende Stoppzeiten <math>\zeta_n\nearrow\zeta</math> der gestoppte Prozess <math>Y^{\zeta_n}</math> eine Lösung der gestoppten stochastischen Differentialgleichung

<math>dY=\alpha(t,Y)dX^{\zeta_n}</math>

ist,

auf <math>\{\zeta <\infty\}</math> fast sicher <math>Y_{t}\to\partial U</math> mit <math>t\to \zeta</math> gilt.<ref>Wolfgang Hackenbroch und Anton Thalmaier: Stochastische Analysis: Eine Einführung in die Theorie der stetigen Semimartingale. Hrsg.: Vieweg+Teubner Verlag Wiesbaden. ISBN 978-3-519-02229-9, S. 297–299. </ref>

<math>\zeta</math> nennt man auch Explosionszeit.

Beispiele

Die SDGL für die geometrische brownsche Bewegung lautet <math> {\rm d}S_t= r S_t {\rm d}t + \sigma S_t {\rm d}W_t </math>. Sie wird beispielsweise im Black-Scholes-Modell zur Beschreibung von Aktienkursen verwendet.
Die SDGL für einen Ornstein-Uhlenbeck-Prozess ist <math>\mathrm dX_t=\theta(\mu-X_t) \mathrm dt + \sigma \mathrm dW_t</math>. Sie wird unter anderem im Vasicek-Modell zur finanzmathematischen Modellierung von Zinssätzen über den Momentanzins verwendet.
Die SDGL für den Wurzel-Diffusionsprozess nach William Feller lautet <math> {\rm d}X_t = \kappa (\theta - X_t){\rm d}t+ \sigma \sqrt{X_t}{\rm d}W_t </math>
Kardar-Parisi-Zhang-Gleichung
Dysons brownsche Bewegung
Tsirelsons stochastische Differentialgleichung

Lösen von stochastischen Differentialgleichungen und Simulation der Lösungen

Genau wie bei deterministischen gibt es auch bei stochastischen Differentialgleichungen keinen allgemeinen Ansatz zur Ermittlung der Lösung. In manchen Fällen (wie bei der oben erwähnten Black-Scholes-SDGL, deren Lösung eine geometrische brownsche Bewegung ist) ist es auch hier möglich, die Lösung zu „erraten“ und durch Ableiten zu verifizieren (wobei das Differenzieren hier mit Hilfe der Itō-Formel erfolgt).

Die Lösungen einer stochastischen Differentialgleichung sind wiederum in starke und schwache Lösungen unterteilt. Der Unterschied liegt darin, dass bei einer starken Lösung mit Initialwert <math>x_0</math> der Wahrscheinlichkeitsraum und die brownsche Bewegung schon a priori gegeben sind, bei der schwachen Lösung können diese selber gewählt werden und es muss nur <math>P\circ X^{-1}_0=\mu</math> für ein gegebenes <math>\mu</math> gelten. Der Satz von Yamada-Watanabe liefert einen Zusammenhang zwischen beiden Begriffen.

In den meisten Fällen, die in der Praxis auftauchen, wie zum Beispiel auch im Fall des Wurzel-Diffusionsprozesses, ist jedoch keine geschlossene Form der Lösung zu erreichen. Doch ist man zumeist auch nur daran interessiert, Zufallspfade der entsprechenden Lösung zu simulieren. Dies kann approximativ durch numerische Diskretisierungsverfahren erreicht werden, etwa durch das Euler-Maruyama-Schema (das dem expliziten Euler-Verfahren für gewöhnliche Differentialgleichungen nachempfunden ist) oder das Milstein-Verfahren.

Stochastische Delay-Differentialgleichungen

Bei einer stochastischen Delay-Differentialgleichung (SDDE, stochastic delay differential equation) hängt der zukünftige Zuwachs nicht nur von dem derzeitigen Zustand, sondern auch von den Zuständen in einem davorliegenden beschränkten Zeitintervall ab. Existenz und Eindeutigkeit sind unter ähnlichen Bedingungen wie in „normalen“ SDGLs gegeben. Seien

<math>f\colon[0,\infty)\times C([-r,0],\mathbb{R}^d)\to\mathbb{R}^d</math>,

<math>g\colon[0,\infty)\times C([-r,0], \mathbb{R}^d)\to\mathbb{R}^{d\times m}</math>

stetig, <math>r>0</math> und <math>W</math> sei eine m-dimensionale Brownsche Bewegung. Dann ist eine stochastische Delay-Differentialgleichung eine Gleichung der Form

wobei <math>X_t(s):=X(t+s) \ \forall s\in[-r,0]</math>

Die dazugehörige Differentialschreibweise lautet dann

<math>{\rm d}X(t)=f_t(X_t)\,{\rm d}t+g_t(X_t)\,{\rm d}W(t)</math>.

Siehe auch

Langevin-Gleichung

Literatur

Bernt Øksendal: Stochastic Differential Equations. An Introduction with Applications. 6. Auflage. Springer, Berlin 2003, ISBN 3-540-04758-1.
Philip E. Protter: Stochastic Integration and Differential Equations. Springer, Berlin 2003, ISBN 3-540-00313-4.

Einzelnachweise