Valutazione delle opzioni alle differenze finite: segni, stabilità, ostacoli e convergenza

Costruire un valutatore auditabile alle differenze finite con segni temporali corretti, operatori discreti, condizioni terminali, al contorno e di evento, complementarità americana, smorzamento e convergenza di prezzo e greche.

Solo a scopo educativo; non costituisce consulenza né raccomandazione d’investimento. Gli investimenti possono comportare perdite.

Risposta diretta

La valutazione delle opzioni alle differenze finite sostituisce un’equazione continua di valutazione con nodi di prezzo e tempo. Un’implementazione valida richiede più di una griglia fitta: segni del tempo di calendario e del tempo a scadenza, operatore spaziale, condizioni terminali e al contorno, salti di evento, schema temporale, solutore dell’ostacolo americano, residui e controlli di convergenza devono essere coerenti.

Il raffinamento della griglia riduce soltanto l’errore numerico. Non corregge payoff, modello stocastico, trattamento dei dividendi, diritto di regolamento, dato di mercato o ipotesi di esecuzione errati. Un solutore può convergere con precisione verso la risposta sbagliata.

Derivare la griglia prima di programmarla

Definire diritto, payoff Phi(S), stile di esercizio, regolamento, tassi, costo del prestito, dividendi, modello di volatilità e coordinate. Con il tempo di calendario t, la PDE di Black-Scholes è V_t + L V = 0; con il tempo a scadenza tau = T - t, è V_tau = L V, partendo da V(S, 0) = Phi(S).
Scegliere griglia di prezzo o log-prezzo, dominio troncato, griglia temporale e allineamento di spot, strike, barriere, dividendi e osservazioni. Dichiarare se N conta intervalli o nodi e testare S_max separatamente da Delta S.
Sulla griglia uniforme S_i = i Delta S, derivare a_i = 0.5[sigma^2 i^2 - (r - q)i], b_i = -(sigma^2 i^2 + r) e c_i = 0.5[sigma^2 i^2 + (r - q)i], così che (L_h V)_i = a_i V_(i-1) + b_i V_i + c_i V_(i+1).
Definire esattamente il passo temporale: esplicito V^(n+1) = (I + Delta tau L_h)V^n; implicito (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. Inserire i termini al contorno nel lato destro e nel livello temporale corretti.
Verificare globalmente segni dei coefficienti, restrizioni tipo CFL, positività, monotonicità, proprietà della matrice e oscillazioni. Euler esplicito ed Euler all’indietro sono normalmente del primo ordine nel tempo; l’implicito è più dissipativo, non universalmente meno accurato. La stabilità lineare di Crank-Nicolson non garantisce monotonicità né risultati privi di oscillazioni presso le discontinuità. Documentare lo smorzamento di Rannacher intorno a ogni discontinuità in cui viene usato.
Per l’esercizio americano, risolvere un metodo discreto di complementarità lineare identificato, come PSOR, penalizzazione oppure operator splitting documentato. Imporre V - Phi >= 0, A V - b >= 0 e (V - Phi) elementwise_times (A V - b) = 0; riportare rilassamento, limite di iterazioni e residuo di complementarità.
Interpolare prezzo e greche allo spot, quindi raffinare dominio, spazio, tempo, allineamento degli eventi e tolleranza del solutore. Confrontare casi europei con una formula e casi americani con un albero profondo o un solutore indipendente; registrare prezzo, Delta, Gamma, frontiera di esercizio, controlli di non arbitraggio e residui.

Le condizioni al contorno appartengono al contratto. In ipotesi standard, una put europea può usare V(0, tau) = K exp(-r tau) e tende a zero al crescere di S, mentre una put americana usa comunemente V(0, tau) = K. Una call europea con dividendi ha un’altra asintotica ad alto prezzo. I dividendi discreti in contanti richiedono salto e interpolazione, come V(S, t-) = V(max(S - D, 0), t+), seguiti nel corretto ordine dall’eventuale decisione di esercizio.

Le differenze centrate possono produrre un coefficiente convettivo negativo anche con un passo temporale piccolo; può servire una discretizzazione upwind, fitted o in log-prezzo. Una PDE di diffusione non acquisisce automaticamente salti di prezzo. Errori di modello, parametro, calibrazione, numerici e di prezzo eseguibile devono restare separati.

Quattro esempi svolti

Il segno temporale cambia con l’orologio. Da V_t + L V = 0 e tau = T - t, la regola della catena fornisce V_tau = L V. Procedere dal payoff a tau = 0 verso un tau maggiore usa quindi l’operatore positivo L; copiare il segno all’indietro del tempo di calendario risolverebbe un problema diverso.
Un nodo può sembrare localmente sicuro. Siano sigma = 25%, r = 4%, q = 0, i = 95 e Delta tau = 0.0005. Si ottengono a_i = 280.13125, b_i = -564.1025 e c_i = 283.93125; i pesi espliciti sono 0.140065625, 0.71794875 e 0.141965625, con somma 0.99998. I valori adiacenti 6, 5, 4 diventano 4.998. Per una put americana con ostacolo 5, una proiezione esplicita mantiene 5; un sistema implicito o Crank-Nicolson richiede ancora il proprio residuo LCP accoppiato.
La restrizione è globale. Con gli stessi input a i = 200, i pesi espliciti sono 0.623, -0.25002 e 0.627; il centro negativo distrugge la monotonicità anche se il nodo vicino allo spot sembrava accettabile. Ridurre il passo a Delta tau = 0.0002 produce 0.2492, 0.499992 e 0.2508, non negativi e con somma 0.999992.
Prezzo e greche richiedono controlli indipendenti. Per una put europea con S = $95, K = $100, T = 0.5, sigma = 25%, r = 4% e q = 0, Black-Scholes fornisce d_1 = -0.088633257689, d_2 = -0.265409952985, prezzo $8.413258404601, Delta -0.535313311861 e Gamma 0.023662224614. Una griglia deve convergere separatamente verso tutti e tre i risultati prima che il premio americano o la frontiera libera siano affidabili.

Controlli in sette fasi e modalità di errore

Il tempo di calendario t e il tempo a scadenza tau possono avanzare con segni incoerenti.
Unità di giorno, anno, tasso, volatilità e dividendo possono essere mescolate.
Payoff, stile di esercizio, diritto di regolamento o livello terminale possono essere codificati male.
Dividend yield continuo e dividendi discreti in contanti possono essere duplicati o ordinati male.
S_max può essere troppo vicino anche se Delta S sembra piccolo.
Le formule al contorno inferiore e superiore possono appartenere al diritto o alle ipotesi di carry sbagliati.
Dividendi, barriere, fixing o osservazioni possono cadere tra nodi temporali ed essere interpolati male.
Il conteggio di nodi e intervalli o le formule delle derivate su griglia non uniforme possono essere sfasati di uno.
Il drift può rendere negativo un coefficiente spaziale centrato e compromettere la monotonicità.
I pesi espliciti possono violare stabilità globale o positività lontano dallo spot.
Crank-Nicolson può oscillare intorno a discontinuità del payoff, barriere o salti di evento.
Lo smorzamento Rannacher o altro può essere omesso, applicato agli eventi errati o non documentato.
I residui del sistema lineare possono essere ampi nonostante una superficie visivamente regolare.
Rilassamento PSOR, tolleranza di arresto, stima iniziale o limite di iterazioni possono essere instabili.
Un massimo nodo per nodo applicato una sola volta può essere scambiato per una soluzione LCP implicita o Crank-Nicolson esatta.
L’interpolazione allo spot può distorcere valore, Delta, Gamma o frontiera libera.
Gamma e frontiere di esercizio possono restare rumorosi anche quando il prezzo sembra convergente.
L’ordine di convergenza osservato o l’estrapolazione di Richardson può essere usato fuori dal regime asintotico regolare.
Diffusione, interpolazione della volatilità locale, salti, parametri o calibrazione possono essere strutturalmente errati.
Un valore teorico può essere presentato come bid, ask, copertura o arbitraggio eseguibile.

Idee errate comuni

“Una griglia più fitta dimostra che il prezzo è corretto.” Non dimostra né il contratto né il modello.
“L’implicito è intrinsecamente meno accurato dell’esplicito.” Entrambi gli schemi Euler sono normalmente del primo ordine temporale, con stabilità ed errore diversi.
“Crank-Nicolson è incondizionatamente stabile, quindi monotono e privo di oscillazioni.” Sono proprietà distinte.
“Il valore americano è una soluzione senza vincoli seguita da un massimo puntuale.” I sistemi impliciti accoppiati richiedono un metodo LCP o un’approssimazione quantificata.
“La corrispondenza con una quotazione valida il solutore e le greche.” Calibrazione e rumore di mercato possono nascondere errori numerici o strutturali.