Optionsbewertung mit finiten Differenzen: Vorzeichen, Stabilität, Hindernisse und Konvergenz

Entwickeln Sie einen prüfbaren Finite-Differenzen-Löser mit korrekten Zeitvorzeichen, diskreten Operatoren, End-, Rand- und Ereignisbedingungen, amerikanischer Komplementarität, Dämpfung sowie Preis- und Greeks-Konvergenz.

Nur zu Bildungszwecken; stellt keine Anlageberatung oder Anlageempfehlung dar. Anlagen können zu Verlusten führen.

Direkte Antwort

Die Optionsbewertung mit finiten Differenzen ersetzt eine kontinuierliche Bewertungsgleichung durch Preis- und Zeitknoten. Eine belastbare Implementierung benötigt mehr als ein feines Gitter: Vorzeichen von Kalenderzeit und Restlaufzeit, räumlicher Operator, End- und Randbedingungen, Ereignissprünge, Zeitschema, Löser für das amerikanische Hindernis, Residuen und Konvergenzkontrollen müssen zusammenpassen.

Gitterverfeinerung reduziert nur numerische Fehler. Sie korrigiert weder falsche Auszahlungen noch ein falsches stochastisches Modell, eine falsche Dividendenbehandlung, einen falschen Abwicklungsanspruch, Marktdaten oder Ausführungsannahmen. Ein Löser kann präzise zur falschen Antwort konvergieren.

Das Gitter vor der Programmierung herleiten

Legen Sie Anspruch, Auszahlung Phi(S), Ausübungsart, Abwicklung, Zinssätze, Leihe, Dividenden, Volatilitätsmodell und Koordinaten fest. Mit Kalenderzeit t lautet die Black-Scholes-PDE V_t + L V = 0; mit Restlaufzeit tau = T - t lautet sie V_tau = L V, ausgehend von V(S, 0) = Phi(S).
Wählen Sie ein Preis- oder Logpreisgitter, den abgeschnittenen Bereich, das Zeitgitter und die Ausrichtung von Spot, Basispreis, Barrieren, Dividenden und Beobachtungen. Geben Sie an, ob N Intervalle oder Knoten zählt, und testen Sie S_max getrennt von Delta S.
Leiten Sie auf dem gleichmäßigen Gitter S_i = i Delta S die Koeffizienten a_i = 0.5[sigma^2 i^2 - (r - q)i], b_i = -(sigma^2 i^2 + r) und c_i = 0.5[sigma^2 i^2 + (r - q)i] her, sodass (L_h V)_i = a_i V_(i-1) + b_i V_i + c_i V_(i+1).
Definieren Sie den Zeitschritt exakt: explizit V^(n+1) = (I + Delta tau L_h)V^n; implizit (I - Delta tau L_h)V^(n+1) = V^n; Crank-Nicolson (I - 0.5 Delta tau L_h)V^(n+1) = (I + 0.5 Delta tau L_h)V^n. Fügen Sie Randterme auf der richtigen rechten Seite und Zeitebene ein.
Prüfen Sie global Koeffizientenvorzeichen, CFL-artige Beschränkungen, Positivität, Monotonie, Matrixeigenschaften und Schwingungen. Expliziter und rückwärts gerichteter Euler sind normalerweise erster Ordnung in der Zeit; implizit ist stärker dissipativ, nicht allgemein ungenauer. Lineare Stabilität von Crank-Nicolson garantiert weder Monotonie noch schwingungsfreie Ergebnisse an Knicken. Dokumentieren Sie jede verwendete Rannacher-Dämpfung um Diskontinuitäten.
Lösen Sie bei amerikanischer Ausübung ein benanntes diskretes lineares Komplementaritätsverfahren, etwa PSOR, Penalty oder dokumentiertes Operator Splitting. Fordern Sie V - Phi >= 0, A V - b >= 0 und (V - Phi) elementwise_times (A V - b) = 0; berichten Sie Relaxation, Iterationsgrenze und Komplementaritätsresiduum.
Interpolieren Sie Preis und Greeks am Spot und verfeinern Sie danach Bereich, Raum, Zeit, Ereignisausrichtung und Lösertoleranz. Vergleichen Sie europäische Fälle mit einer Formel und amerikanische mit einem tiefen Baum oder unabhängigen Löser; protokollieren Sie Preis, Delta, Gamma, Ausübungsgrenze, Arbitrageprüfungen und Residuen.

Randbedingungen gehören zum Vertrag. Unter Standardannahmen kann ein europäischer Put V(0, tau) = K exp(-r tau) verwenden und für wachsendes S gegen null gehen, während ein amerikanischer Put häufig V(0, tau) = K verwendet. Ein europäischer Call mit Dividenden hat eine andere Hochpreis-Asymptote. Diskrete Bardividenden erfordern einen Sprung samt Interpolation, etwa V(S, t-) = V(max(S - D, 0), t+), gefolgt von einer etwaigen Ausübungsentscheidung in der richtigen Reihenfolge.

Zentrale Differenzen können auch bei kleinem Zeitschritt einen negativen Konvektionskoeffizienten erzeugen; ein Upwind-, angepasstes oder Logpreisverfahren kann nötig sein. Eine Diffusions-PDE erhält auch nicht automatisch Kurssprünge. Modell-, Parameter-, Kalibrierungs-, numerische und ausführbare Preisfehler sind getrennt zu halten.

Vier durchgerechnete Beispiele

Das Zeitvorzeichen ändert sich mit der Uhr. Aus V_t + L V = 0 und tau = T - t folgt nach der Kettenregel V_tau = L V. Vom Payoff bei tau = 0 zu größerem tau wird daher mit dem positiven Operator L fortgeschritten; das Rückwärtsvorzeichen der Kalenderzeit würde ein anderes Problem lösen.
Ein Knoten kann lokal unauffällig wirken. Für sigma = 25%, r = 4%, q = 0, i = 95 und Delta tau = 0.0005 ergeben sich a_i = 280.13125, b_i = -564.1025 und c_i = 283.93125; die expliziten Gewichte sind 0.140065625, 0.71794875 und 0.141965625 mit Summe 0.99998. Nachbarwerte 6, 5, 4 werden zu 4.998. Bei einem amerikanischen Put mit Hindernis 5 behält eine explizite Projektion 5; ein implizites oder Crank-Nicolson-System benötigt weiterhin sein gekoppeltes LCP-Residuum.
Die Beschränkung gilt global. Mit denselben Daten bei i = 200 lauten die expliziten Gewichte 0.623, -0.25002 und 0.627; die negative Mitte zerstört die Monotonie, obwohl der spotnahe Knoten akzeptabel erschien. Eine Verringerung auf Delta tau = 0.0002 ergibt 0.2492, 0.499992 und 0.2508; alle sind nicht negativ und summieren sich zu 0.999992.
Preis und Greeks brauchen unabhängige Kontrollen. Für einen europäischen Put mit S = $95, K = $100, T = 0.5, sigma = 25%, r = 4% und q = 0 liefert Black-Scholes d_1 = -0.088633257689, d_2 = -0.265409952985, Preis $8.413258404601, Delta -0.535313311861 und Gamma 0.023662224614. Ein Gitter muss separat zu allen drei Ergebnissen konvergieren, bevor man seinem amerikanischen Aufschlag oder seiner freien Grenze vertraut.

Siebenstufige Kontrollen und Fehlerquellen

Kalenderzeit t und Restlaufzeit tau können mit widersprüchlichen Vorzeichen fortgeschrieben werden.
Einheiten für Tage, Jahre, Zins, Volatilität und Dividende können vermischt werden.
Auszahlung, Ausübungsart, Abwicklungsanspruch oder Endschicht können falsch codiert sein.
Kontinuierliche Rendite und diskrete Bardividenden können doppelt gezählt oder falsch angeordnet werden.
S_max kann zu nah liegen, obwohl Delta S klein erscheint.
Untere und obere Randformeln können zum falschen Anspruch oder Carry gehören.
Dividenden, Barrieren, Fixings oder Beobachtungen können zwischen Zeitknoten liegen und schlecht interpoliert werden.
Knoten- und Intervallzählung oder Ableitungsformeln nicht gleichmäßiger Gitter können um eins verschoben sein.
Drift kann einen zentralen Raumkoeffizienten negativ machen und die Monotonie brechen.
Explizite Gewichte können fern vom Spot globale Stabilität oder Positivität verletzen.
Crank-Nicolson kann an Auszahlungsknicken, Barrieren oder Ereignissprüngen schwingen.
Rannacher- oder andere Dämpfung kann fehlen, an falschen Ereignissen erfolgen oder undokumentiert bleiben.
Residuen des linearen Systems können trotz optisch glatter Fläche zu groß sein.
PSOR-Relaxation, Abbruchtoleranz, Startwert oder Iterationsgrenze können instabil sein.
Ein einmaliges knotenweises Maximum kann mit einer exakten impliziten oder Crank-Nicolson-LCP-Lösung verwechselt werden.
Interpolation am Spot kann Wert, Delta, Gamma oder freie Grenze verzerren.
Gamma und Ausübungsgrenzen können noch rauschen, obwohl der Preis konvergiert erscheint.
Beobachtete Konvergenzordnung oder Richardson-Extrapolation kann außerhalb des glatten asymptotischen Bereichs verwendet werden.
Diffusion, lokale Volatilitätsinterpolation, Sprünge, Parameter oder Kalibrierung können strukturell falsch sein.
Ein theoretischer Wert kann als ausführbarer Bid, Ask, Hedge oder Arbitrage dargestellt werden.

Häufige Fehlannahmen

„Ein feineres Gitter beweist, dass der Preis stimmt.“ Es beweist weder Vertrag noch Modell.
„Implizit ist grundsätzlich ungenauer als explizit.“ Beide Euler-Verfahren sind meist erster Ordnung in der Zeit, mit unterschiedlicher Stabilität und Fehlercharakteristik.
„Crank-Nicolson ist unbedingt stabil, also monoton und schwingungsfrei.“ Das sind getrennte Eigenschaften.
„Amerikanischer Wert ist eine unbeschränkte Lösung mit anschließendem punktweisen Maximum.“ Gekoppelte implizite Systeme verlangen ein LCP-Verfahren oder eine quantifizierte Näherung.
„Die Übereinstimmung mit einem Kurs validiert Löser und Greeks.“ Kalibrierung und Marktrauschen können numerische oder strukturelle Fehler verdecken.