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.

Geschichte

Schon Norbert Wiener untersuchte Integrale von deterministischen Integranden $f(t)$ bezüglich der brownschen Bewegung der Form

$\int f(t)\mathrm {d} _{t}X(t,\omega )$

und mehrdimensionale stochastische Integrale dieser Form. Itō Kiyoshi verallgemeinerte diese Resultate und die moderne Theorie der stochastischen Integration baut im Wesentlichen auf seiner Arbeit auf. 2000 wurde ein versiegelter Umschlag von Wolfgang Döblin aus dem Jahre 1940 geöffnet. Darin befanden sich Resultate über die stochastische Integration, die er Itō Kiyoshi vorwegnahm. Döblin verstarb allerdings im selben Jahr, weshalb die Arbeit unentdeckt blieb.

Stochastische Integration

Es existieren verschiedene stochastische Integralbegriffe.

Generell muss, um klassische stochastische Integrale zu konstruieren (Wiener, Itō, Stratonowitsch), der Integrand $X$ gewisse Kriterien der Messbarkeit und Integrierbarkeit erfüllen. Sei hier ${\mathcal {M}}_{0,\operatorname {loc} }^{c}$ der Raum der ${\mathcal {F}}_{t}$-adaptierten, stetigen lokalen Martingale $M=(M_{t})_{t\geq 0}$ mit $M_{0}=0$. Für $M\in {\mathcal {M}}_{0,\operatorname {loc} }^{c}$ mit $\mathbb {E} [\langle M\rangle _{\infty }]<\infty$ definiert man den L²-Hilbert-Raum der Äquivalenzklassen von ${\mathcal {L}}([0,\infty )\times \Omega ,{\mathcal {B}}([0,\infty ))\otimes {\mathcal {F}}_{\infty },\mu _{M})$, wobei $\mu _{M}$ für $\Xi \in {\mathcal {B}}([0,\infty ))\otimes {\mathcal {F}}_{\infty }$ durch

$\mu _{M}(\Xi )=\mathbb {E} \left[\int _{0}^{\infty }1_{\Xi }(s,\omega )\mathrm {d} \langle M\rangle _{s}(\omega )\right]$

definiert ist (die Norm wird durch $\mu _{M}$ induziert). Die richtige Wahl der Integranden sind die $L^{2}(\mu _{M})$-integrierbaren progressiv-messbaren $X$. Möchte man allgemeiner gegen nicht-stetige Semimartingale integrieren, dann muss man die Klasse der Integranden auf vorhersagbare Prozesse beschränken.

Integralbegriff nach Wiener

Sei $C$ der klassische Wiener-Raum ausgestattet mit dem Wiener-Maß $\mu _{W}$. Sei $F$ ein Funktional auf $C$, dann nennt man

$\int _{C}F(x)\mathrm {d} \mu _{W}$

Wiener-Integral. Allgemein werden Integrale einer deterministischen Funktion bezüglich eines Wiener-Prozesses so bezeichnet.

Integralbegriffe nach Itō und Stratonowitsch

Itō-Integral

Das Itō-Integral ist zunächst für Semimartingale $Y$ und für elementare vorhersagbare Prozesse definiert, d. h. für (an eine Filtration $({\mathcal {F}}_{t})_{t}$ adaptierte) stochastische Prozesse $H$ der Form

$H_{t}=h_{0}1_{\{0\}}(t)+\sum \limits _{i=0}^{n-1}h_{i}1_{(t_{i},t_{i+1}]}(t),\quad 0=t_{0}<...<t_{n},\quad n\in \mathbb {N} ,\quad h_{i}{\text{ }}{\mathcal {F}}_{t_{i}}{\text{-messbar}},$

durch

$I_{Y}(H):=\sum \limits _{i=0}^{n-1}h_{i}\left(Y_{t_{i}}-Y_{t_{i+1}}\right)$

Die elementaren Prozesse können alternativ auch allgemeiner mit Stoppzeiten anstelle von deterministischen Zeitpunkten $t_{i}$ definiert werden. Sei $\mathbf {L}$ der Raum der adaptierten Càglàd-Prozesse und $\mathbf {S}$ der Raum der elementaren vorhersagbaren Prozesse. Wir nennen die Topologie, welche durch die gleichmäßige Konvergenz auf kompakten Mengen in Wahrscheinlichkeit erzeugt wird, die UCP-Topologie (UCP für englisch uniformly on compact in probability). Man kann nun zeigen, dass $\mathbf {S}$ in der UCP-Topologie dicht in $\mathbf {L}$ liegt. Damit lässt sich das stochastische Integral (als lineare Abbildung $H\longmapsto I_{Y}(H)$) auf $\mathbf {L}$ fortsetzen. Konkret: Das Ito-Integral eines Prozesses $X\in \mathbf {L}$ ist also definiert als der Grenzwert

$\int _{0}^{t}X_{s}\mathrm {d} Y_{s}:=\int _{0}^{t}X_{s}^{(n)}\mathrm {d} Y_{s}\quad {\text{in UCP}}$

für jede Folge von Prozessen $(X^{(n)})_{n}\subseteq \mathbf {S}$, die gegen $X$ konvergieren (bzgl. der UCP-Topologie). Die Definition ist in der Tat unabhängig von der gewählten Folge.

In der allgemeinsten Formulierung werden als Integratoren Semimartingale $Y$ und als Integranden vorhersagbare Prozesse $X$ zugelassen (die zusätzlich gewisse Integrierbarkeitsbedingungen erfüllen). Sind die Integratoren $Y$ zusätzlich stetig, genügt es für die Integranden $X$ progressiv-messbar (und in $L(Y)$) zu sein.

Als Folge der (abstrakten) Konstruktion des Integrals erhält man folgenden anschaulicheren Zusammenhang:

Seien $Y$ ein Semimartingal und $X$ ein adaptierter Càdlàg- oder Càglàd-Prozess. Dann gilt für jede Folge reeller Zahlen $\left(T_{n}\right)_{n}$ mit $\lim _{n\to \infty }T_{n}=\infty$ und für jede Folge $\left(\pi _{n}\right)_{n}$ von Partitionen des Intervalls $[0,T_{n}]$ mit $\lim _{n\to \infty }\sup _{k}|t_{k+1}^{n}-t_{k}^{n}|=0$ die Konvergenz

$\sup _{t\in [0,T]}{\Big |}\int _{0}^{t}X_{s-}\,\mathrm {d} Y_{s}-\sum _{i}X_{t_{i}^{n}}(Y_{t_{i+1}^{n}\land t}-Y_{t_{i}^{n}\land t}){\Big |}\xrightarrow {n\to \infty } 0$

in Wahrscheinlichkeit $\forall \,T>0$, wobei $X_{s-}=\lim _{u\to s,u<s}X_{u}$ (die linkstetige Version von $X$) und $t\land s=\min(t,s)$. Dies lässt sich auch kompakter schreiben als

$\lim _{n\to \infty }\sum _{i}X_{t_{i}^{n}}(Y_{t_{i+1}^{n}\land t}-Y_{t_{i}^{n}\land t})=\int _{0}^{t}X_{s-}\,\mathrm {d} Y_{s}\quad {\text{in UCP}}$

Die Aussage gilt sogar allgemeiner für Folgen von random partitions tending to the identity, was aber mehr Notation für die Definition des Begriffs erfordert.

Stratonowitsch-Integral

→ Hauptartikel: Stratonowitsch-Integral

Für ein Semimartingal $Y$ und ein adaptierter Càdlàg-Prozess $X$, sodass die quadratische Kovariation $\left[X,Y\right]$ existiert, kann man das Stratonowitsch-Integral oder Fisk-Stratonowitsch-Integral (nach Ruslan Leontjewitsch Stratonowitsch und Donald Fisk) definiert durch

${\begin{aligned}\int _{0}^{t}X_{s-}\circ \mathrm {d} Y_{s}&=\int _{0}^{t}X_{s-}\mathrm {d} Y_{s}+{\frac {1}{2}}\left[X,Y\right]_{t}-{\frac {1}{2}}\sum \limits _{0\leq s\leq t}\Delta Y_{s}\Delta X_{s}\\&=\int _{0}^{t}X_{s-}\mathrm {d} Y_{s}+{\frac {1}{2}}\left[X,Y\right]_{t}^{c}\end{aligned}}$

wobei $\int _{0}^{t}X_{s-}\mathrm {d} Y_{s}$ das Itō-Integral und $\Delta Y_{s}=Y_{s}-Y_{s-}$ die Sprungstelle von $Y$ an der Stelle $s$ sind.

Daraus folgt ähnlich wie beim Ito-Integral eine anschaulichere Darstellung des Integrals:

Seien $Y$ und $X$ wie in der obigen Definition und gelte zusätzlich, dass $Y$ und $X$ keine Sprünge zum gleichen Zeitpunkt haben, d. h. $\sum \limits _{0\leq t}\Delta Y_{t}\Delta X_{t}=0$. Dann gilt für jede Folge reeller Zahlen $\left(T_{n}\right)_{n}$ mit $\lim _{n\to \infty }T_{n}=\infty$ und für jede Folge $\left(\pi _{n}\right)_{n}$ von Partitionen des Intervalls $[0,T_{n}]$ mit $\lim _{n\to \infty }\sup _{k}|t_{k+1}^{n}-t_{k}^{n}|=0$ die Konvergenz

$\lim _{n\to \infty }\sum _{i}{\frac {1}{2}}(X_{t_{i}^{n}}+X_{t_{i+1}^{n}})(Y_{t_{i+1}^{n}\land t}-Y_{t_{i}^{n}\land t})=\int _{0}^{t}X_{s-}\,\mathrm {d} Y_{s}\quad {\text{in UCP}}$.

Auch hier gilt die Aussage noch allgemeiner für sequences of random partitions tending to the identity.

Vergleich der Integrale

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).

Eine Brownsche Bewegung $B_{s}$ und das Integral von $B_{s}\,\mathrm {d} B_{s}$

Verallgemeinerungen

Integralbegriff nach Ogawa

→ Hauptartikel: Ogawa-Integral

Der Integralbegriff ist für nicht-adaptierte Integranden. Man bildet eine Zufallsreihe mit Hilfe eines orthonormalen Systems im $L^{2}$-Hilbertraum und lässt diese dann gegen das Ogawa-Integral konvergieren. Der entsprechende Kalkül wird nicht-kausales Kalkül genannt.

Integralbegriff nach Marcus

Eine Verallgemeinerung des Fisk-Stratonowitsch-Integrals auf allgemeine Semimartingale mit Jumps ist das Marcus-Integral. Stochastische Differentialgleichungen mit diesem Integralbegriff nennt man vom Marcus-Typ. Marcus entwickelte ein Kalkül, welches auf dem Kalkül von McShane basiert.

Integralbegriff nach Hitsuda-Skorochod

→ Hauptartikel: Skorochod-Integral

Eine Erweiterung des Itō-Integrals auf nicht-adaptierte Prozesse ist das Hitsuda-Skorochod-Integral. Das Integral ist ein Spezialfall des adjungierten Operators des Ableitungsoperator der Malliavin-Ableitung. Im Falle der Integrierbarkeit bezüglich der brownschen Bewegung und der Adaptierbarkeit des Integranden erhält man gerade das Itō-Integral. Alternativ lässt sich das Integral auch über die Wiener-Chaos-Zerlegung definieren.

Integralbegriff nach Walsh

Das Walsh-Integral ist ein Integral bezüglich eines Martingal-Maßes um stochastische partielle Differentialgleichungen zu studieren. Das Integral wurde von John B. Walsh eingeführt. Von Robert C. Dalang existiert eine Erweiterung für distributionelle Integranden.

Beispiele

• Sei $(W_{t}),t>0$ ein (Standard-)Wiener-Prozess. Trivialerweise $\int _{a}^{b}\mathrm {d} W_{t}=W_{b}-W_{a}$ für $a,b\geq 0$.
• Sei $(W_{t}),t>0$ ein (Standard-)Wiener-Prozess. Zu berechnen ist das Itō-Integral $\int _{0}^{T}W_{t}\,\mathrm {d} W_{t}$. Schreibt man der Kürze halber $B_{i}:=W_{iT/n},\Delta B_{i}:=B_{i+1}-B_{i}$ und benutzt man die Identität

$B_{i+1}^{2}-B_{i}^{2}=(B_{i+1}-B_{i})^{2}+2B_{i}(B_{i+1}-B_{i}),$

so erhält man aus obiger Integrationsvorschrift

${\begin{aligned}I&=\lim _{n\to \infty }\sum _{i=0}^{n-1}B_{i}(B_{i+1}-B_{i})\\&=\lim _{n\to \infty }\left({\frac {1}{2}}\sum _{i=0}^{n-1}(B_{i+1}^{2}-B_{i}^{2})-{\frac {1}{2}}\sum _{i=0}^{n-1}(B_{i+1}-B_{i})^{2}\right)\\&={\frac {1}{2}}\lim _{n\to \infty }\sum _{i=0}^{n-1}(B_{i+1}^{2}-B_{i}^{2})-{\frac {1}{2}}\lim _{n\to \infty }\sum _{i=0}^{n-1}(\Delta B_{i})^{2}\\&={\frac {1}{2}}\lim _{n\to \infty }\left(B_{n}^{2}-B_{0}^{2}\right)-{\frac {T}{2}}\lim _{n\to \infty }{\frac {1}{n}}\sum _{i=0}^{n-1}\left({\sqrt {\frac {n}{T}}}\Delta B_{i}\right)^{2}.\end{aligned}}$

Benutzt man nun einerseits, dass $B_{0}=W_{0}=0,B_{n}=W_{T}$ gilt, sowie andererseits die Eigenschaft, dass $\left({\sqrt {\frac {n}{T}}}\Delta B_{i}\right)^{2}$ i.i.d. $\chi ^{2}$-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

$I={\frac {1}{2}}W_{T}^{2}-{\frac {T}{2}}.$

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

${\begin{aligned}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 {1}{2}}W_{T}^{2}\end{aligned}}$

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 $Y$ 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 $Y$ ein Lévy-Prozess mit konstantem Erwartungswert, $X$ eine nicht vorgreifende beschränkte Funktion von $Y$ und $t$ (d. h., für jedes $t>0$ ist $X_{t}$ messbar bezüglich der σ-Algebra $\sigma (Y_{s};s<t)$, die von den Zufallsvariablen $Y_{s},\,s<t$ erzeugt wird), so ist der Prozess

$t\mapsto \int _{0}^{t}X_{s}\,\mathrm {d} Y_{s}$

ein lokales Martingal bezüglich der natürlichen Filtrierung von $Y$. 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 $(X_{t})$ mit $t\geq 0$ Itō-Prozess genannt, wenn es eine Brownsche Bewegung $(W_{t})$ mit $t\geq 0$ und stochastische Prozesse $(a_{t}(X_{t},t))$, $(b_{t}(X_{t},t))$ gibt mit

$X_{t}=X_{0}+\int _{0}^{t}a_{s}(X_{s},s)\,\mathrm {d} s+\int _{0}^{t}b_{s}(X_{s},s)\,\mathrm {d} W_{s}\,,$

wobei angenommen wird, dass die beiden Integrale existieren. In Differentialschreibweise wird diese Gleichung als

$\mathrm {d} X_{t}=a_{t}(X_{t},t)\,\mathrm {d} t+b_{t}(X_{t},t)\,\mathrm {d} W_{t}$

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

Das Prädikat „$X$ 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 $a_{t}$ und der Diffusionskoeffizient $b_{t}$ 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 statistischen Physik und der Finanzmathematik, hat sich der Itō-Kalkül inzwischen zu einem unverzichtbaren mathematischen Werkzeug entwickelt.

Siehe auch

• Diskretes stochastisches Integral
• Euler-Maruyama-Verfahren

Literatur

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