imported>Phzh: Form, typo

2025-09-18T22:32:43Z

Form, typo

Neue Seite

{{Belege fehlen}}
[[Datei:SQRTDiffusion.png|mini|hochkant=1.7|Drei Pfade von Wurzel-Diffusionsprozessen]]

Der '''Wurzel-Diffusionsprozess''' (daneben ist auch im deutschsprachigen Raum die englische Bezeichnung ''square-root diffusion'' gebräuchlich) ist ein [[stochastischer Prozess]], der über eine [[stochastische Differentialgleichung]] definiert ist. Seit diese Prozessklasse erstmals [[1951]] durch den [[Kroatien|kroatisch]]-[[Vereinigte Staaten|amerikanischen]] Mathematiker [[William Feller]] studiert wurde, hat sie zahlreiche Anwendungsfelder in verschiedenen Bereichen der [[Finanzmathematik]] gefunden.

== Definition und Parameter ==
[[Datei:CIR-Prozess1.png|mini|250px|Drei unabhängige Prozesse mit v0=1, theta=0.6, sigma=0.2, kappa=20]]

Seien reelle Parameter <math> \kappa \ge 0, \; \theta > 0, \sigma </math> und eine [[Wiener-Prozess|brownsche Bewegung]] <math> (W_t),\ t\ge 0</math> gegeben. Ein stochastischer Prozess <math> (X_t), \; t \ge 0 </math> heißt Wurzel-Diffusionsprozess mit jenen Parametern, wenn er die [[stochastische Differentialgleichung]]

:<math> \mathrm dX_t = \kappa (\theta-X_t)\mathrm dt + \sigma \sqrt{X_t}\mathrm dW_t </math>

löst. Zumeist wird mit <math> X_0 = x</math> eine [[Anfangsbedingung]] festgelegt, sodass aus obiger Differentialgleichung ein stochastisches [[Anfangswertproblem]] wird.

Ein Grund für die Beliebtheit dieses Prozesses liegt sicherlich darin, dass die Wirkungsweisen der einzelnen Parameter leicht interpretierbar sind (und somit gut als „Stellschrauben“ bei der Modellierung genutzt werden können):

* <math> \theta </math> ist das Gleichgewichtsniveau des Prozesses (engl. ''mean reversion level''). Ist <math> X_0 <\theta </math>, so ist der Driftterm <math> \kappa (\theta-X_t) \mathrm dt </math> positiv und „zieht“ den Prozess nach oben, andernfalls ist er negativ und <math>X</math> tendiert nach unten. <math> X_t </math> driftet also in Erwartung stets <math> \theta </math> entgegen.
* <math> \kappa </math> gibt die Stärke dieser Regulierungsfunktion an (engl. ''mean reversion speed''). Für <math> \kappa=0 </math> ist diese außer Kraft gesetzt und der Verlauf von <math>X</math> nur vom Diffusionsterm <math> \sqrt{X_t}\mathrm dW_t </math> abhängig. Je größer <math> \kappa </math> ist, umso „steifer“ ist <math>X</math> an <math> \theta </math> gebunden. Theoretisch wären auch negative Werte für <math> \kappa </math> denkbar, doch dann hätte der Prozess völlig andere Eigenschaften: oberhalb von <math> \theta </math> gestartet, würde er gegen <math> +\infty </math> divergieren, ansonsten würde er bald negativ – und damit aus dem eigenen Definitionsbereich herausfallen.
* Der Parameter <math> \sigma </math> ist die [[Volatilität]] des Prozesses. Er steuert, wie stark <math>X</math> den Schwankungen der Brownschen Bewegung ausgesetzt ist. Für große Werte von <math> \sigma </math> wird <math>X</math> über die Zeit sehr stark schwanken, bei <math> \sigma=0 </math> konvergiert <math>X</math> exponentiell gegen <math> \theta </math>.

== Eigenschaften ==
[[Datei:CIR-Prozess2.png|mini|250px|theta=1.4, sonst dieselben Werte wie oben]]
Obwohl obige Differentialgleichung keine geschlossene Lösung besitzt (d. h. es ist nicht möglich, den resultierenden Prozess <math>X</math> als Funktion von <math>t</math> und <math>W</math> auszudrücken), kann trotzdem viel über die Eigenschaften des Prozesses ausgesagt werden:

* Für <math> \kappa \ge 0 </math> ist <math>X</math> nie negativ, d. h. es gilt <math> \forall t \ge 0:\; P(X_t \ge 0)=1 </math>. Gilt zusätzlich die Stabilitätsbedingung <math> 2 \kappa \theta \ge \sigma^2 </math>, so ist <math>X</math> sogar fast sicher strikt positiv. Das liegt daran, dass für <math> X_t \to 0 </math> auch der Diffusionsterm gegen Null strebt und somit die positive Drift den Prozess wieder nach oben „zieht“.
* Gegeben sei ein vorausgegangener Wert <math>X_u</math> mit <math>u<t</math>, dann hat <math>X_t</math> eine [[Chi-Quadrat-Verteilung#Nichtzentrale Chi-Quadrat-Verteilung|Nichtzentrale Chi-Quadrat-Verteilung]], bei der der [[Nichtzentralitätsparameter]] von <math>X_u</math> abhängt.
* Langfristig konvergiert die Verteilung von <math>X_t</math> unabhängig vom Startpunkt gegen eine [[Gammaverteilung]] mit [[Erwartungswert]] <math> \theta </math> und [[Varianz (Stochastik)|Varianz]] <math> \frac{\sigma^2 \theta}{2 \kappa} </math>. Wird diese Verteilung als Startverteilung gewählt, hat <math>X</math> zu jedem Zeitpunkt dieselbe Verteilung.
* Der Wurzel-Diffusionsprozess ist ein [[Affiner Prozess]]

== Anwendung ==
=== Cox-Ingersoll-Ross-Modell ===
[[Datei:CIR-Prozess3.png|mini|250px|Hier ist sigma=0.4]]
In ihrer bahnbrechenden Arbeit ''A theory of the term structure of interest rates'' schlugen die amerikanischen Mathematiker Cox, Ingersoll und Ross 1985 den Wurzel-Diffusionsprozess erstmals als Modell für kurzfristige [[Zins]]sätze vor. Gleichzeitig präsentierten sie neue Ansätze zur Berechnung von Übergangswahrscheinlichkeiten des Prozesses. Ihr Zinsmodell wurde bald zu einem der Standardmodelle am Markt, sodass der Wurzel-Diffusionsprozess bald den Beinamen '''CIR-Prozess''' erhielt.

=== Heston-Modell für stochastische Volatilität ===
Im Jahr 1993 verlieh Steven L. Heston der »square-root diffusion« neue Beliebtheit, als er ein komplexes Kapitalmarktmodell entwickelte, indem er das Modell von [[Black-Scholes-Modell|Black und Scholes]] um eine stochastische Volatilität erweiterte, die bei Letzteren noch als konstant angenommen worden war. Diese wurde durch eine Wurzeldiffusion simuliert und konnte so erstmals über die brownsche Bewegung mit dem [[Börsenkurs]] [[Korrelation|korreliert]] werden: Das machte es möglich, das natürliche Phänomen (''leverage''), wonach fallende Kurse die Aufregung (die Volatilität) an den Märkten steigern, mathematisch zu erfassen. Das [[Heston-Modell]] wird heute als wichtigste Erweiterung zu Black-Scholes angesehen und verschaffte der Wurzel-Diffusion endgültig einen Stammplatz in den Lehrbüchern über stochastische Prozesse.

=== Stochastische Zeitverschiebung ===
[[Datei:CIR-Prozess4.png|mini|250px|Dieselben Werte wie ganz oben, aber theta=1 und kappa=50 (man beachte die Skala)]]
Eine weitere Möglichkeit, mit Hilfe des CIR-Prozesses allgemeine [[Lévy-Prozess]]e (eine weitaus breitere Klasse als die von Heston untersuchte Prozessklasse) mit stochastischer Volatilität zu versehen, basiert auf einer Idee von [[Wolfgang Döblin]] (dem Sohn [[Alfred Döblin]]s) aus den 30er-Jahren: Döblin hatte damals vorgeschlagen, die Volatilität eines Prozesses zu steuern, indem man die Geschwindigkeit der Zeit, in der der Prozess abläuft, kontrolliert. Viele Jahrzehnte später stellte sich heraus, dass der Wurzel-Diffusionsprozess bestens dazu geeignet ist, als vorgeschaltete „Zeitmaschine“ diese Geschwindigkeit zu beeinflussen. Ist zum Beispiel <math> (L_t),\; t \ge 0 </math> ein solcher Lévy-Prozess und <math>X</math> der oben definierte Prozess, so erhält man durch <math> Y_t := L_{X_t} , t \ge 0 </math> einen Prozess <math>Y</math> mit stochastischer Zeit.

== Simulation des Prozesses ==
In der Praxis hat man zwei Möglichkeiten, den CIR-Prozess zu simulieren, die beide ihre Vor- und Nachteile besitzen:
* Eine Möglichkeit besteht darin, die bekannte Übergangswahrscheinlichkeit (nichtzentrale Chi-Quadrat-Verteilung) auszunutzen und zu jedem gegebenen Punkt den jeweils nächsten (in beliebigem Abstand) direkt zu ziehen. Diese Methode bietet eine exakte (im Sinne von: erwartungstreue) Simulation des Prozesses, jedoch muss für jeden Schritt zur Ermittlung des Chi-Quadrat-Verteilten Updates eine [[Normalverteilung|normal-]] und eine [[Gammaverteilung|gammaverteilte]] [[Zufallsvariable]] gezogen werden. Letztere ist allerdings numerisch aufwändig zu ziehen, was diese Technik beispielsweise für [[Monte-Carlo-Simulation]]en ungeeignet erscheinen lässt.
* Die zweite Möglichkeit – nicht exakt, aber dafür schneller und für feinere Gitter geeignet – ist die Simulation durch ein einfaches [[Euler-Maruyama-Verfahren]] (auf diese Art sind die Beispielbilder zu diesem Artikel entstanden). Hierbei ist jedoch zu beachten, dass das Euler-Update in jedem Schritt mit positiver Wahrscheinlichkeit negativ werden kann (was beim echten Prozess nicht der Fall ist), was nicht nur numerische Probleme bereitet.

== Literatur ==
* {{Literatur
|Autor=Hull, John C.
|Titel=Options, Futures and Other Derivatives
|Verlag=Upper Saddle River, NJ: [[Prentice Hall]]
|Datum=2003
|ISBN=0-13-009056-5
|Sprache=en
|Online=https://archive.org/details/optionsfuturesot00hull_1}}
* {{Literatur
|Autor=Cox, J.C., J.E. Ingersoll and S.A. Ross
|Titel=A Theory of the Term Structure of Interest Rates
|Sammelwerk=[[Econometrica]]
|Band=53
|Nummer=2
|Datum=1985
|Seiten=385–407
|Sprache=en
|DOI=10.2307/1911242
|JSTOR=1911242}}
* {{Literatur
|Autor=Maghsoodi, Y.
|Titel=Solution of the extended CIR Term Structure and Bond Option Valuation
|Sammelwerk=Mathematical Finance
|Band=6
|Nummer=6
|Datum=1996
|Seiten=89–109
|Sprache=en
|DOI=10.1111/j.1467-9965.1996.tb00113.x}}
* {{Literatur
|Autor=Damiano Brigo, Fabio Mercurio
|Titel=Interest Rate Models — Theory and Practice with Smile, Inflation and Credit
|Auflage=2nd ed. 2006
|Verlag=Springer Verlag
|Datum=2001
|ISBN=3-540-22149-2
|Sprache=en}}
* {{Literatur
|Autor=Damiano Brigo, Fabio Mercurio
|Titel=A deterministic-shift extension of analytically tractable and time-homogeneous short rate models
|Sammelwerk=Finance & Stochastics
|Band=5
|Nummer=3
|Datum=2001
|Seiten=369–388
|Sprache=en
|Online=https://ideas.repec.org/a/spr/finsto/v5y2001i3p369-387.html
|DOI=10.1007/PL00013541}}
* [https://github.com/AlexandreMoulti/bachelier Open Source library implementing the CIR process in python]
* {{Literatur
|Autor=Giuseppe Orlando, Rosa Maria Mininni, Michele Bufalo
|Titel=Forecasting interest rates through Vasicek and CIR models: A partitioning approach
|Sammelwerk=Journal of Forecasting
|Band=39
|Nummer=4
|Datum=2020
|ISSN=1099-131X
|Seiten=569–579
|Sprache=en
|arXiv=1901.02246
|DOI=10.1002/for.2642}}

[[Kategorie:Stochastischer Prozess]]
[[Kategorie:Stochastische Differentialgleichung]]

Wurzel-Diffusionsprozess - Versionsgeschichte

imported>Phzh: Form, typo