Cuando calibramos nuestra volatilidad local, necesitamos crear una superficie de volatilidad implícita suavizada y sin arbitraje (por ejemplo SSVI o Sabr).
Imaginemos la superficie de volatilidad implícita suavizada como una caja negra donde proporcionamos una huelga y un vencimiento y obtenemos nuestra volatilidad implícita de vuelta.
Ahora hay una fórmula bien conocida que puede dar las vols locales como función de las vols implícitas si asumimos que tenemos disponible una grilla continua de vols implícitas:
$\sigma_{\mathrm{Dup}}(T,K)^2 = \frac{ \frac{\partial w}{\partial T} }{1 - \frac{y}{w} \frac{\partial w}{\partial y}+ \frac{1}{4}\left( - \frac{1}{4} - \frac{1}{w} + \frac{y^2}{w^2} \right) \left(\frac{\partial w}{\partial y}\right)^2 + \frac{1}{2}\frac{\partial^2 w}{\partial y^2} }$
Dada la superficie suavizada y libre de arbitraje de volatilidad implícita, podemos calcular las derivadas anteriores numéricamente y podemos extraer cualquier vol local que deseemos.
Pero aún veo que uno de los enfoques muy populares (algo similar presentado en el documento "Interpolación de Volatilidad" de Jesper Andreasen y Brian Huge) donde intentan resolver un problema de minimización mal planteado con alguna restricción de suavidad en una grilla de EDP combinada con la EDP de Dupire hacia adelante, por ejemplo. También asumen alguna estructura de la volatilidad local (constante por partes en el tiempo, lineal en la dimensión de la huelga) para poder aplicar este enfoque de FDM.
Esto parece ser más complicado y también, en términos de modelado, la constante por partes en las dimensiones de tiempo se siente muy restrictiva y mucho menos precisa (por ejemplo, la interpolación lineal en la dimensión de la huelga).
Siento que me estoy perdiendo algo. ¿Quizás los resultados del primer método no son lo suficientemente suaves? ¿Tal vez hay otras desventajas que me estoy perdiendo?