Generally when we implement Dupire's local volatility model, we follow the steps below:
- Calculate implied volatility from given historical data
- Fit the implied volatility skew. So we also know the corresponding call option prices.
- Use the Dupire's Formula to calculate the local volatility $$\sigma_{loc}^2(K,T) = \frac{\partial_TC(K,T) + rK\partial_KC(K,T)}{\frac{1}{2}K^2\partial_{KK}C(K,T)}$$.
- Use the SDE $$dS_t = r_tS_tdt + \sigma_{loc}(S_t, t) S_t dW_t$$ to do Monte Carlo simulation to get any call option price we want.
I have a few questions:
- What's the differences between the call option prices from step 2 and step 4?
- Using historical data, we can generate discrete points in the implied volatility surface. If we fix those points and regardless of the shape of the volatility skew, could we perfectly recover the historical call option prices using Dupire's Formula regardless of simulation path?
Is there a formal proof or a detailed explanation of the above two?