Precificação de opções por diferenças finitas: sinais, estabilidade, obstáculos e convergência

Construa um precificador auditável por diferenças finitas com sinais temporais corretos, operadores discretos, condições terminais, de contorno e de eventos, complementaridade americana, amortecimento e convergência de preço e gregas.

Somente para fins educacionais; não constitui recomendação nem aconselhamento de investimento. Investir pode causar perdas.

Resposta direta

A precificação de opções por diferenças finitas substitui uma equação contínua de avaliação por nós de preço e tempo. Uma implementação válida exige mais do que uma grade fina: sinais do tempo de calendário e do prazo até o vencimento, operador espacial, condições terminais e de contorno, saltos de eventos, esquema temporal, solucionador do obstáculo americano, resíduos e controles de convergência precisam ser coerentes.

Refinar a grade trata apenas o erro numérico. Isso não corrige payoff, modelo estocástico, tratamento de dividendos, direito de liquidação, dado de mercado ou premissa de execução errados. Um solucionador pode convergir precisamente para a resposta errada.

Derive a grade antes de programá-la

Fixe o direito, o payoff Phi(S), o estilo de exercício, a liquidação, as taxas, o custo de empréstimo, os dividendos, o modelo de volatilidade e as coordenadas. Com tempo de calendário t, a EDP de Black-Scholes é V_t + L V = 0; com prazo até o vencimento tau = T - t, é V_tau = L V, partindo de V(S, 0) = Phi(S).
Escolha uma grade de preço ou log-preço, domínio truncado, grade temporal e alinhamento de spot, strike, barreiras, dividendos e observações. Declare se N conta intervalos ou nós e teste S_max separadamente de Delta S.
Na grade 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) e 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 exatamente o passo 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. Insira termos de contorno no lado direito e na camada temporal corretos.
Verifique globalmente sinais dos coeficientes, restrições do tipo CFL, positividade, monotonicidade, propriedades matriciais e oscilações. Euler explícito e Euler para trás costumam ser de primeira ordem no tempo; o implícito é mais dissipativo, não universalmente menos preciso. A estabilidade linear de Crank-Nicolson não garante monotonicidade nem ausência de oscilações em quinas. Documente o amortecimento de Rannacher em torno de toda descontinuidade em que for usado.
Para exercício americano, resolva um método discreto de complementaridade linear identificado, como PSOR, penalização ou operator splitting documentado. Exija V - Phi >= 0, A V - b >= 0 e (V - Phi) elementwise_times (A V - b) = 0; informe relaxação, limite de iterações e resíduo de complementaridade.
Interpole preço e gregas no spot; depois refine domínio, espaço, tempo, alinhamento de eventos e tolerância do solucionador. Compare casos europeus com uma fórmula e casos americanos com uma árvore profunda ou solução independente; registre preço, Delta, Gamma, fronteira de exercício, testes de não arbitragem e resíduos.

As condições de contorno pertencem ao contrato. Sob premissas usuais, uma put europeia pode usar V(0, tau) = K exp(-r tau) e tende a zero conforme S cresce, enquanto uma put americana em geral usa V(0, tau) = K. Uma call europeia com dividendos tem outra assíntota de preço alto. Dividendos discretos em dinheiro exigem salto e interpolação, como V(S, t-) = V(max(S - D, 0), t+), seguidos na ordem correta por qualquer decisão de exercício.

Diferenças centradas podem produzir coeficiente de convecção negativo mesmo com passo temporal pequeno; pode ser necessária discretização upwind, ajustada ou em log-preço. Uma EDP de difusão também não incorpora saltos de preço automaticamente. Erros de modelo, parâmetro, calibração, numéricos e de preço executável devem permanecer separados.

Quatro exemplos resolvidos

O sinal temporal muda com o relógio. De V_t + L V = 0 e tau = T - t, a regra da cadeia fornece V_tau = L V. Avançar do payoff em tau = 0 para um tau maior usa, portanto, o operador positivo L; copiar o sinal regressivo do tempo de calendário resolveria outro problema.
Um nó pode parecer seguro localmente. Sejam sigma = 25%, r = 4%, q = 0, i = 95 e Delta tau = 0.0005. Então a_i = 280.13125, b_i = -564.1025 e c_i = 283.93125; os pesos explícitos são 0.140065625, 0.71794875 e 0.141965625, somando 0.99998. Valores vizinhos 6, 5, 4 são atualizados para 4.998. Para uma put americana com obstáculo 5, uma projeção explícita mantém 5; um sistema implícito ou Crank-Nicolson ainda precisa de seu resíduo LCP acoplado.
A restrição é global. Com os mesmos dados em i = 200, os pesos explícitos são 0.623, -0.25002 e 0.627; o centro negativo destrói a monotonicidade embora o nó perto do spot parecesse aceitável. Reduzir o passo para Delta tau = 0.0002 produz 0.2492, 0.499992 e 0.2508, todos não negativos e com soma 0.999992.
Preço e gregas exigem controles independentes. Para uma put europeia com S = $95, K = $100, T = 0.5, sigma = 25%, r = 4% e q = 0, Black-Scholes fornece d_1 = -0.088633257689, d_2 = -0.265409952985, preço $8.413258404601, Delta -0.535313311861 e Gamma 0.023662224614. Uma grade precisa convergir separadamente para os três resultados antes que seu prêmio americano ou sua fronteira livre mereçam confiança.

Controles em sete etapas e modos de falha

O tempo de calendário t e o prazo até o vencimento tau podem avançar com sinais incompatíveis.
Unidades de dia, ano, taxa, volatilidade e dividendo podem ser misturadas.
Payoff, estilo de exercício, direito de liquidação ou camada terminal podem ser codificados incorretamente.
Dividend yield contínuo e dividendos discretos em dinheiro podem ser duplicados ou ordenados incorretamente.
S_max pode ficar perto demais embora Delta S pareça pequeno.
Fórmulas de contorno inferior e superior podem pertencer ao direito ou às premissas de carry errados.
Dividendos, barreiras, fixings ou observações podem cair entre nós temporais e ser mal interpolados.
Contagens de nós e intervalos ou fórmulas derivativas em grade não uniforme podem apresentar erro de uma unidade.
O drift pode tornar negativo um coeficiente espacial centrado e eliminar a monotonicidade.
Pesos explícitos podem violar estabilidade global ou positividade longe do spot.
Crank-Nicolson pode oscilar perto de quinas do payoff, barreiras ou saltos de eventos.
Rannacher ou outro amortecimento pode ser omitido, aplicado nos eventos errados ou não documentado.
Resíduos do sistema linear podem estar frouxos apesar de uma superfície visualmente suave.
Relaxação PSOR, tolerância de parada, estimativa inicial ou limite de iterações podem ser instáveis.
Um máximo nó a nó aplicado uma só vez pode ser confundido com uma solução LCP implícita ou Crank-Nicolson exata.
A interpolação no spot pode enviesar valor, Delta, Gamma ou fronteira livre.
Gamma e fronteiras de exercício podem continuar ruidosos mesmo depois de o preço parecer convergido.
A ordem de convergência observada ou a extrapolação de Richardson pode ser usada fora do regime assintótico suave.
Difusão, interpolação de volatilidade local, saltos, parâmetros ou calibração podem estar estruturalmente errados.
Um valor teórico pode ser apresentado como bid, ask, hedge ou arbitragem executável.

Equívocos comuns

“Uma grade mais fina prova que o preço está correto.” Ela não prova nem o contrato nem o modelo.
“O implícito é intrinsecamente menos preciso que o explícito.” Ambos os esquemas Euler costumam ser de primeira ordem temporal, com estabilidade e erro diferentes.
“Crank-Nicolson é incondicionalmente estável, portanto monótono e livre de oscilações.” São propriedades distintas.
“O valor americano é uma solução sem restrições seguida de um máximo pontual.” Sistemas implícitos acoplados exigem método LCP ou aproximação quantificada.
“Ajustar uma cotação valida o solucionador e as gregas.” Calibração e ruído de mercado podem ocultar erros numéricos ou estruturais.