Stochastische Integration

Die Theorie der stochastischen Integration befasst sich mit Integralen und Differentialgleichungen in der Stochastik. Sie verallgemeinert die Integralbegriffe von Henri Léon Lebesgue und Thomas Jean Stieltjes auf eine breitere Menge von Integratoren. Es sind stochastische Prozesse mit unendlicher Variation, insbesondere der Wiener-Prozess, als Integratoren zugelassen. Die Theorie der stochastischen Integration stellt dabei die Grundlage der stochastischen Analysis dar, deren Anwendungen sich zumeist mit der Untersuchung stochastischer Differentialgleichungen beschäftigen.

Integralbegriffe nach Itō und Stratonowitsch

Seien <math> (X_t),(Y_t),t\in [a,b] </math> zwei (nicht notwendigerweise unabhängige) reellwertige stochastische Prozesse auf einem gemeinsamen Wahrscheinlichkeitsraum <math> (\Omega,\mathcal{F},P) </math>. Als Itō-Integral (nach Itō Kiyoshi) von X nach Y über dem Intervall [a,b] bezeichnet man die Zufallsvariable

<math> I:=\int_a^b X_{t-}\,\mathrm dY_t:=\lim_{n \to \infty} \sum_{i=1}^n X_{a+(i-1)h} (Y_{a+ih}-Y_{a+(i-1)h}), \quad h=\frac{b-a}n.</math>

Das zugehörige Stratonowitsch-Integral (nach Ruslan Leontjewitsch Stratonowitsch) berechnet sich für dieselbe Wahl von <math>h</math> als

<math> S:=\int_a^b X_t\circ\mathrm dY_t:=\lim_{n \to \infty} \sum_{i=1}^n \frac{1}{2}(X_{a+(i-1)h}+X_{a+ih}) (Y_{a+ih}-Y_{a+(i-1)h}).</math>

Beim Itō-Integral wird der Integrand X also stets am Anfang des h-Intervalls ausgewertet, bei Stratonowitsch werden der Anfangs- und Endwert gemittelt. Bei gewöhnlichen (Riemann- oder Lebesgue-) Integralen von deterministischen (nicht zufälligen) und hinreichend glatten (beispielsweise stetigen) Funktionen hat dies keinen Einfluss auf das Ergebnis, doch im stochastischen Fall gilt: Sind X und Y nicht unabhängig, so kann das tatsächlich zu verschiedenen Werten führen (siehe Beispiel unten).

Als Klasse der möglichen Integratoren Y werden in der allgemeinsten Formulierung Semimartingale zugelassen, die Integranden X sind vorhersagbare Prozesse.

Beispiel

Sei <math> (W_t), t>0</math> ein (Standard-)Wiener-Prozess. Zu berechnen ist das Itō-Integral <math>\int_0^T W_t\,\mathrm dW_t</math>. Schreibt man der Kürze halber <math>B_i := W_{iT/n} , \Delta B_i :=B_{i+1}-B_i </math> und benutzt man die Identität

<math>B_{i+1}^2 -B_i^2=(B_{i+1}-B_i)^2 + 2 B_i (B_{i+1}-B_i),</math>

so erhält man aus obiger Integrationsvorschrift

<math>\begin{align}

 I &= \lim_{n \to \infty} \sum_{i=0}^{n-1} B_i (B_{i+1}-B_i)\\
   &= \lim_{n \to \infty} \left( \frac 12 \sum_{i=0}^{n-1}(B_{i+1}^2-B_i^2) -\frac 12 \sum_{i=0}^{n-1} (B_{i+1}-B_i)^2 \right)\\
   &= \frac 12 \lim_{n \to \infty} \sum_{i=0}^{n-1}(B_{i+1}^2-B_i^2) -\frac 12 \lim_{n \to \infty} \sum_{i=0}^{n-1} (\Delta B_i)^2\\
   &= \frac 12 \lim_{n \to \infty} \left(B_n^2-B_0^2\right) -\frac T2 \lim_{n \to \infty} \frac 1n \sum_{i=0}^{n-1} \left(\sqrt{\frac nT} \Delta B_i\right)^2.

\end{align}</math>

Benutzt man nun einerseits, dass <math>B_0=W_0=0, B_n=W_T</math> gilt, sowie andererseits die Eigenschaft, dass <math>\left(\sqrt{\frac nT} \Delta B_i\right)^2</math> i.i.d. <math>\chi^2</math>-verteilt ist (wegen der unabhängigen, normalverteilten Zuwächse der Brownschen Bewegung), so folgt mit dem Gesetz der großen Zahlen für den hinteren Grenzwert

<math>I=\frac{1}{2}W_T^2-\frac{T}{2}.</math>

Um das entsprechende Stratonowitsch-Integral zu berechnen, nutzt man die Stetigkeit der Brownschen Bewegung aus:

<math>\begin{align}

 S &= \lim_{n \to \infty} \sum_{i=0}^{n-1} \frac{1}{2}(B_{i+1}+B_i)(B_{i+1}-B_i)\\
   &= \lim_{n \to \infty} \sum_{i=0}^{n-1} \frac{1}{2}(B_{i+1}^2-B_i^2)\\
   &= \lim_{n \to \infty} \frac{1}{2}(B_{n}^2-B_0^2)\\
   &= \frac 12 W_T^2 

\end{align}</math> Itō- und Stratonowitsch-Integral über demselben Prozess führen also zu verschiedenen Ergebnissen, wobei das Stratonowitsch-Integral eher der intuitiven Ahnung aus der gewöhnlichen (deterministischen) Integralrechnung entspricht.

Martingaleigenschaft

Der bei weitem am häufigsten verwendete Integrator <math>Y </math> ist eine Brownsche Bewegung. Der entscheidende Vorteil, den das Stratonowitsch-Integral nicht hat und der letztendlich dazu führte, dass sich das Itō-Integral weitgehend als Standard durchgesetzt hat, ist die folgende Eigenschaft:

Sei <math>Y </math> ein Lévy-Prozess mit konstantem Erwartungswert, <math>X </math> eine nicht vorgreifende beschränkte Funktion von <math>Y </math> und <math>t </math> (d. h., für jedes <math>t>0</math> ist <math>X_t</math> messbar bezüglich der σ-Algebra <math>\sigma (Y_s; s<t)</math>, die von den Zufallsvariablen <math>Y_s,\,s<t</math> erzeugt wird), so ist der Prozess

<math>t \mapsto \int_0^t X_s\,\mathrm dY_s </math>

ein lokales Martingal bezüglich der natürlichen Filtrierung von <math>Y </math>. Unter zusätzlichen Beschränktheitsbedingungen ist der Integralprozess sogar ein Martingal.

Anwendung: Itō-Prozess

Ausgehend vom Itōschen Integralbegriff ist es nun möglich, eine breite Klasse von stochastischen Prozessen zu definieren: Demnach wird ein stochastischer Prozess <math>(X_t)</math> mit <math> t \geq 0 </math> Itō-Prozess genannt, wenn es eine Brownsche Bewegung <math>(W_t)</math> mit <math> t\geq 0</math> und stochastische Prozesse <math>(a_t(X_t, t))</math>, <math>(b_t(X_t, t))</math> gibt mit

<math>X_t = X_0 + \int_0^t a_s(X_s, s)\,\mathrm ds + \int_0^t b_s(X_s, s) \,\mathrm dW_s\,,</math>

wobei angenommen wird, dass die beiden Integrale existieren.<ref>Hui-Hsiung Kuo: Introduction to Stochastic Integration. Springer, 2006, ISBN 978-0387-28720-1, S. 102 (eingeschränkte Vorschau in der Google-Buchsuche).</ref> In Differentialschreibweise wird diese Gleichung als

<math>\mathrm{d}X_t = a_t(X_t, t)\, \mathrm{d}t + b_t(X_t, t)\, \mathrm{d}W_t</math>

notiert. Ein Itō-Prozess kann also als verallgemeinerter Wiener-Prozess mit zufälliger Drift und Volatilität angesehen werden.

Das Prädikat „<math>X</math> ist ein Itō-Prozess“ wird somit zu einem stochastischen Pendant zum Begriff der Differenzierbarkeit. Ausgehend hiervon wurden dann von Itō selbst die ersten stochastischen Differentialgleichungen definiert.

Hängen der Driftkoeffizient <math>a_t</math> und der Diffusionskoeffizient <math>b_t</math> nicht von der Zeit ab, so spricht man von Itō-Diffusion - hängen sie zusätzlich von der Zeit ab, so liegt dagegen ein allgemeinerer Itō-Prozess vor.

Durch zahlreiche Anwendungen in der mathematischen Modellierung, insbesondere in der Quantenphysik und der Finanzmathematik, hat sich der Itō-Kalkül inzwischen zu einem unverzichtbaren mathematischen Werkzeug entwickelt.

Siehe auch

• Diskretes stochastisches Integral

Literatur

• J. Jacod, A. Shiryaev: Limit theorems for stochastic processes. Springer, Berlin.
• P. Protter: Stochastic integrals and differential equations. Springer, Berlin.

Einzelnachweise

<references />