Valoración de opciones por diferencias finitas: signos, estabilidad, obstáculos y convergencia

Construya un valorador auditable por diferencias finitas con signos temporales correctos, operadores discretos, condiciones terminales, de frontera y de evento, complementariedad americana, amortiguación y convergencia del precio y las griegas.

Solo con fines educativos; no constituye asesoramiento de inversión. Invertir puede ocasionar pérdidas.

Respuesta directa

La valoración de opciones por diferencias finitas sustituye una ecuación continua de valoración por nodos de precio y tiempo. Una implementación válida exige más que una malla fina: deben concordar los signos del tiempo calendario y del tiempo hasta vencimiento, el operador espacial, las condiciones terminales y de frontera, los saltos por eventos, el esquema temporal, el solucionador del obstáculo americano, los residuos y los controles de convergencia.

Refinar la malla solo reduce el error numérico. No corrige un payoff, modelo estocástico, tratamiento de dividendos, derecho de liquidación, dato de mercado o supuesto de ejecución equivocado. Un solucionador puede converger con precisión a una respuesta errónea.

Derivar la malla antes de programarla

Fije el derecho, el payoff Phi(S), el estilo de ejercicio, la liquidación, los tipos, el coste de préstamo, los dividendos, el modelo de volatilidad y las coordenadas. Con tiempo calendario t, la PDE de Black-Scholes es V_t + L V = 0; con tiempo hasta vencimiento tau = T - t, es V_tau = L V, partiendo de V(S, 0) = Phi(S).
Elija una malla de precio o log-precio, el dominio truncado, la malla temporal y la alineación de spot, strike, barreras, dividendos y observaciones. Indique si N cuenta intervalos o nodos y pruebe S_max por separado de Delta S.
En la malla uniforme S_i = i Delta S, derive a_i = 0.5[sigma^2 i^2 - (r - q)i], b_i = -(sigma^2 i^2 + r) y c_i = 0.5[sigma^2 i^2 + (r - q)i], de modo que (L_h V)_i = a_i V_(i-1) + b_i V_i + c_i V_(i+1).
Defina exactamente el paso temporal: explícito V^(n+1) = (I + Delta tau L_h)V^n; implícito (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. Incorpore los términos de frontera en el lado derecho y nivel temporal correctos.
Compruebe globalmente signos de coeficientes, restricciones tipo CFL, positividad, monotonicidad, propiedades matriciales y oscilaciones. Los Euler explícito y hacia atrás suelen ser de primer orden temporal; el implícito es más disipativo, no universalmente menos exacto. La estabilidad lineal de Crank-Nicolson no garantiza monotonicidad ni resultados sin oscilaciones en los quiebres. Documente la amortiguación de Rannacher alrededor de cada discontinuidad donde se aplique.
Para ejercicio americano, resuelva un método discreto de complementariedad lineal identificado, como PSOR, penalización u operator splitting documentado. Exija V - Phi >= 0, A V - b >= 0 y (V - Phi) elementwise_times (A V - b) = 0; informe relajación, límite de iteraciones y residuo de complementariedad.
Interpole precio y griegas en spot; después refine dominio, espacio, tiempo, alineación de eventos y tolerancia del solucionador. Compare casos europeos con una fórmula y casos americanos con un árbol profundo o solucionador independiente; registre precio, Delta, Gamma, frontera de ejercicio, controles de no arbitraje y residuos.

Las fronteras pertenecen al contrato. Bajo supuestos estándar, una put europea puede usar V(0, tau) = K exp(-r tau) y tiende a cero cuando crece S, mientras que una put americana suele usar V(0, tau) = K. Una call europea con dividendos tiene otra asíntota de precio alto. Los dividendos discretos en efectivo requieren un salto e interpolación, como V(S, t-) = V(max(S - D, 0), t+), seguido en el orden correcto por cualquier decisión de ejercicio.

Las diferencias centradas pueden generar un coeficiente de convección negativo aunque el paso temporal sea pequeño; puede requerirse una discretización upwind, ajustada o en log-precio. Una PDE de difusión tampoco incorpora automáticamente saltos de precio. Los errores de modelo, parámetros, calibración, numéricos y de precio ejecutable deben permanecer separados.

Cuatro ejemplos resueltos

El signo temporal cambia con el reloj. De V_t + L V = 0 y tau = T - t, la regla de la cadena da V_tau = L V. Avanzar desde el payoff en tau = 0 hacia un tau mayor usa por tanto el operador positivo L; copiar el signo de retroceso del tiempo calendario resolvería otro problema.
Un nodo puede parecer seguro localmente. Sean sigma = 25%, r = 4%, q = 0, i = 95 y Delta tau = 0.0005. Entonces a_i = 280.13125, b_i = -564.1025 y c_i = 283.93125; los pesos explícitos son 0.140065625, 0.71794875 y 0.141965625, cuya suma es 0.99998. Los valores vecinos 6, 5, 4 se actualizan a 4.998. Para una put americana con obstáculo 5, una proyección explícita conserva 5; un sistema implícito o Crank-Nicolson todavía requiere su residuo LCP acoplado.
La restricción es global. Con los mismos datos en i = 200, los pesos explícitos son 0.623, -0.25002 y 0.627; el centro negativo destruye la monotonicidad aunque el nodo cercano a spot pareciera aceptable. Reducir el paso a Delta tau = 0.0002 produce 0.2492, 0.499992 y 0.2508, no negativos y con suma 0.999992.
Precio y griegas necesitan controles independientes. Para una put europea con S = $95, K = $100, T = 0.5, sigma = 25%, r = 4% y q = 0, Black-Scholes da d_1 = -0.088633257689, d_2 = -0.265409952985, precio $8.413258404601, Delta -0.535313311861 y Gamma 0.023662224614. Una malla debe converger por separado hacia los tres resultados antes de confiar en su prima americana o frontera libre.

Controles en siete pasos y modos de fallo

El tiempo calendario t y el tiempo hasta vencimiento tau pueden avanzar con signos incompatibles.
Pueden mezclarse unidades de días, años, tipos, volatilidad y dividendos.
El payoff, estilo de ejercicio, derecho de liquidación o nivel terminal pueden codificarse mal.
El rendimiento continuo y los dividendos discretos en efectivo pueden duplicarse u ordenarse mal.
S_max puede quedar demasiado cerca aunque Delta S parezca pequeño.
Las fórmulas de frontera inferior y superior pueden corresponder a otro derecho o supuesto de carry.
Dividendos, barreras, fixings u observaciones pueden caer entre nodos temporales e interpolarse mal.
El recuento de nodos e intervalos o las fórmulas derivadas de una malla no uniforme pueden tener un desfase de uno.
La deriva puede volver negativo un coeficiente espacial centrado y romper la monotonicidad.
Los pesos explícitos pueden incumplir estabilidad global o positividad lejos de spot.
Crank-Nicolson puede oscilar alrededor de quiebres del payoff, barreras o saltos por eventos.
Rannacher u otra amortiguación puede omitirse, aplicarse a eventos equivocados o quedar sin documentar.
Los residuos del sistema lineal pueden ser laxos pese a una superficie visualmente suave.
La relajación PSOR, tolerancia de parada, estimación inicial o límite de iteraciones pueden ser inestables.
Un máximo nodo a nodo aplicado una sola vez puede confundirse con una solución LCP implícita o Crank-Nicolson exacta.
La interpolación en spot puede sesgar valor, Delta, Gamma o frontera libre.
Gamma y las fronteras de ejercicio pueden seguir siendo ruidosas cuando el precio ya parece convergido.
El orden de convergencia observado o la extrapolación de Richardson pueden usarse fuera del régimen asintótico suave.
La difusión, interpolación de volatilidad local, saltos, parámetros o calibración pueden ser estructuralmente erróneos.
Un valor teórico puede presentarse como bid, ask, cobertura o arbitraje ejecutable.

Errores comunes

“Una malla más fina demuestra que el precio es correcto.” No demuestra ni el contrato ni el modelo.
“El implícito es intrínsecamente menos exacto que el explícito.” Ambos esquemas Euler suelen ser de primer orden temporal, con distinta estabilidad y comportamiento del error.
“Crank-Nicolson es incondicionalmente estable, por lo que es monótono y no oscila.” Son propiedades distintas.
“El valor americano es una solución sin restricción seguida de un máximo puntual.” Los sistemas implícitos acoplados requieren un método LCP o una aproximación cuantificada.
“Ajustar una cotización valida el solucionador y las griegas.” La calibración y el ruido de mercado pueden ocultar errores numéricos o estructurales.