Hedging and Pricing European-type, Early-Exercise and Discrete Barrier Options using Algorithm for the Convolution of Legendre Series

Resumo Técnico: Hedge e Precificação de Opções Europeias, de Exercício Antecipado e de Barreira Discreta utilizando um Algoritmo para a Convolução de Séries de Legendre

Definição do Problema
O artigo aborda o desafio de precificar e realizar o hedge de uma ampla variedade de derivativos financeiros — incluindo opções europeias, de exercício antecipado (americanas/bermudas) e de barreira com monitoramento discreto — sob processos de Lévy gerais. As técnicas numéricas existentes enfrentam limitações específicas:

• Métodos baseados em FFT (ex: CONV): Embora populares, muitas vezes exigem um grande número de termos de Fourier para alcançar precisão suficiente, tornando-os computacionalmente caros, particularmente para aproximar funções de densidade de probabilidade (PDFs) ou precificar opções com características não suaves.
• Métodos baseados em quadratura (ex: QUAD): Frequentemente dependem de PDFs de forma fechada ou requerem a inversão de funções características via FFT. Além disso, tipicamente produzem valores de pontos individuais em vez de curvas de precificação contínuas, e o cálculo das Gregas da opção (sensibilidades) é frequentemente limitado a diferenças finitas de primeira ordem, o que pode carecer de precisão.
• Métodos de expansão de série (ex: COS, SWIFT): Embora o método COS ofereça convergência exponencial, ele requer a seleção a priori de um domínio computacional; se o domínio for inadequado, os erros propagam-se recursivamente no tempo, levando a preços incorretos para opções de exercício antecipado. O método SWIFT aborda características de exercício antecipado, mas carece de uma prova teórica de convergência exponencial, uma vez que seu parâmetro de escala é escolhido heuristicamente. Adicionalmente, esses métodos frequentemente não abordam o cálculo das Gregas.

Metodologia
Os autores propõem o método CONLeg, que utiliza a convolução de séries de Legendre de suporte compacto para aproximar fórmulas de precificação e hedge de opções. A metodologia integra os seguintes componentes:

1. Convolução de Séries de Legendre: O algoritmo central, adaptado de Hale e Townsend (2014a), aproxima a convolução de duas funções (a função de payoff e a PDF refletida) usando séries de Legendre finitas. Esta abordagem evita o fenômeno de Gibbs associado às aproximações de FFT para funções não periódicas, dividindo o domínio de convolução em partes esquerda e direita ($h_L$ e $h_R$) e computando coeficientes via uma relação de recorrência de três termos.
2. Framework Chebfun: O método emprega o Chebfun, um sistema MATLAB de computação numérica com funções, para representar funções como séries de Chebyshev. Isso permite aproximações adaptativas de alta precisão (aprox. 15 dígitos).
3. Tratamento de PDFs Não Suaves: Para processos onde a PDF não é suave ou carece de uma forma fechada (ex: Variance Gamma, CGMY), o artigo utiliza:
   • Série de Fourier Complexa (CFS): Para representar a PDF quando uma forma fechada não está disponível.
   • Método Fourier–Padé: Para localizar singularidades na PDF e construir aproximações por partes precisas, resolvendo assim o fenômeno de Gibbs.
4. Formulação de Precificação e Gregas:
   • Precificação: O preço da opção é formulado como uma integral de convolução $V(x,K,t) \approx e^{-r(T-t)}K (f_M * g_N)(\tilde{x})$. Diferente dos métodos que computam valores pontuais, o CONLeg gera uma curva de precificação contínua (um polinômio por partes) sobre um intervalo de strikes ou preços à vista.
   • Hedge (Gregas): Como o preço é representado por uma série de Chebyshev/Legendre, as derivadas (Delta, Gamma, Vega) são computadas analiticamente através da diferenciação da série, evitando as imprecisões dos métodos de diferença finita.
5. Aplicação a Tipos de Opções:
   • Europeia/Bermuda/Americana: O método formula equações de precificação recursivas. Para opções americanas, utiliza extrapolação de Richardson em preços de opções bermudas.
   • Opções de Barreira: O método lida com o monitoramento discreto tratando os níveis de barreira como pontos de ruptura na aproximação polinomial por partes.
6. Seleção de Domínio: O intervalo computacional $[c, d]$ é determinado usando os cumulantes do processo de Lévy subjacente para garantir que o intervalo capture a maior parte da massa da PDF.

Principais Contribuições

• Geração de Curvas: A principal contribuição é a capacidade de formular curvas de precificação e de Gregas (funções de $S$ ou $K$) em vez de valores pontuais individuais. Isso é alcançado aproximando diretamente a integral de convolução.
• Abordagem Livre de Quadratura: O método oferece uma abordagem livre de quadratura ao aproveitar as propriedades de polinômios ortogonais (Legendre e Chebyshev) e o teorema da convolução, eliminando a necessidade de pesos de integração numérica.
• Robustez à Não-Suavidade: O método demonstra alta precisão mesmo quando a PDF neutra ao risco é não suave ou contém singularidades, desde que a técnica Fourier–Padé seja utilizada para localizar e tratar tais singularidades.
• Flexibilidade: A abordagem é aplicável a uma ampla gama de processos estocásticos (GBM, VG, NIG, CGMY), independentemente de existir uma PDF de forma fechada, pois pode trabalhar com funções características via representação CFS.
• Implementação: O artigo fornece uma implementação de prova de conceito usando Chebfun, destacando a conveniência de utilizar a computação de alto nível baseada em funções para modelagem financeira.

Resultados Numéricos
O artigo apresenta experimentos numéricos comparando o CONLeg contra métodos estabelecidos (COS, filter-COS, CONV, Lewis-FRFT, QUAD-CONV, SWIFT) sob vários modelos e tipos de opções:

• Opções Europeias: Sob o modelo Black-Scholes (GBM), o CONLeg atinge precisão de precisão de máquina ($R_\infty \approx 10^{-14}$) para preços, Delta e Gamma em menos de 0,1 segundo para 1000 pontos. Para opções de longo vencimento (T=50, 100), mantém alta precisão.
• PDFs de Curto Vencimento/Não Suaves: No modelo Variance Gamma (VG), com uma PDF de pico agudo e singularidade, o CONLeg supera o CONV, Lewis-FRFT e QUAD-CONV, alcançando precisão comparável aos métodos COS e filter-COS.
• Opções de Exercício Antecipado: Para opções bermudas e americanas sob os modelos NIG e CGMY, o CONLeg iguala ou supera ligeiramente a precisão do método COS dentro de tempos computacionais similares.
• Opções de Barreira: Para opções de barreira com monitoramento discreto (diário e mensal), o CONLeg apresenta desempenho comparável ao método SWIFT em termos de precisão e tempo de computação.

Significância e Alegações
Os autores alegam que o método CONLeg oferece um desempenho favorável ou comparável às técnicas existentes nos experimentos numéricos. Sua significância reside em:

1. Precisão: Proporciona alta precisão para opções dentro/fora do dinheiro e com vencimentos variados, mesmo sob PDFs não suaves.
2. Eficiência na Geração de Curvas: Gera eficientemente curvas inteiras de precificação e hedge, o que é vantajoso para calibração de modelos e gestão de risco, em vez de calcular pontos discretos.
3. Robustez Teórica: Evita as armadilhas de seleção de domínio do método COS para precificação recursiva e fornece uma maneira sistemática de lidar com singularidades via Fourier–Padé.
4. Praticidade: Os autores reconhecem que a implementação atual baseada em Chefun incurre em algum overhead devido aos procedimentos adaptativos, mas o método é uma "prova de conceito" que demonstra a viabilidade da convolução baseada em funções para modelagem financeira. Eles observam que uma versão não adaptativa seria necessária para uma comparação de eficiência justa quanto às taxas de convergência, mas os resultados atuais são promissores.

O artigo conclui afirmando que trabalhos futuros estenderão este método para opções dependentes do caminho sob processos de Lévy com mudança de tempo e modelos de volatilidade estocástica.