A Sequential Quadratic Programming Method for Volatility Estimation in Option Pricing
27 Pages Posted: 14 May 2007
Date Written: March 29, 2006
Our goal is to identify the volatility function in Dupire's equation from given option prices. Following an optimal control approach in a Lagrangian framework, we propose a globalized sequential quadratic programming (SQP) algorithm with a modified Hessian - to ensure that every SQP step is a descent direction - and implement a line search strategy. In each level of the SQP method a linear-quadratic optimal control problem with box constraints is solved by a primal-dual active set strategy. This guarantees L^1 constraints for the volatility, in particular assuring its positivity. The proposed algorithm is founded on a thorough first- and second-order optimality analysis. We prove the existence of local optimal solutions and of a Lagrange multiplier associated with the inequality constraints. Furthermore, we prove a sufficient second-order optimality condition and present some numerical results underlining the good properties of the numerical scheme.
Keywords: Dupire equation, parameter identification, optimal control, optimality
JEL Classification: G13, C61
Suggested Citation: Suggested Citation