Stochastic Systems: Modeling, Analysis, Synthesis, Control, and their Applications to EngineeringView this Special Issue
Research Article | Open Access
Su-mei Zhang, Li-he Wang, "A Fast Fourier Transform Technique for Pricing European Options with Stochastic Volatility and Jump Risk", Mathematical Problems in Engineering, vol. 2012, Article ID 761637, 17 pages, 2012. https://doi.org/10.1155/2012/761637
A Fast Fourier Transform Technique for Pricing European Options with Stochastic Volatility and Jump Risk
We consider European options pricing with double jumps and stochastic volatility. We derived closed-form solutions for European call options in a double exponential jump-diffusion model with stochastic volatility (SVDEJD). We developed fast and accurate numerical solutions by using fast Fourier transform (FFT) technique. We compared the density of our model with those of other models, including the Black-Scholes model and the double exponential jump-diffusion model. At last, we analyzed several effects on option prices under the proposed model. Simulations show that the SVDEJD model is suitable for modelling the long-time real-market changes and stock returns are negatively correlated with volatility. The model and the proposed option pricing method are useful for empirical analysis of asset returns and managing the corporate credit risks.
The classical Black-Scholes (BS) model  has long been known to result in systematically biased option valuation. By adding jumps to the archetypal price process with Gaussian innovations Merton  is able to partly explain the observed deviations from the benchmark model which are characterized by fat tail and excess kurtosis in the returns distribution. For an overview of “stylized facts” on asset returns see Cont . Statistical properties of implied volatilities are summarized in Cont et al. . In the sequel also other authors develop more realistic models, for example, the pure jump models of Eberlein and Keller , Madan et al. , and Duffie et al. , stochastic volatility models of Steven , and stochastic volatility model with normal jumps of Bates  and Keppo et al. . The double exponential jump-diffusion (DEJD) model, recently proposed by Kou , generates a highly skewed and leptokurtic distribution and is capable of matching key features of stock and index returns. Moreover, the DEJD model leads to tractable pricing formulas for exotic and path dependent options . Accordingly, the DEJD model has gained wide acceptance. However, the DEJD model cannot capture the volatility clustering effects, which can be captured by stochastic volatility models . Jump-diffusion models and the stochastic volatility model complement each other: the stochastic volatility model can incorporate dependent structures better, while the DEJD model has better analytical tractability, especially for path-dependent options. Since allowing interest rates to be stochastic does not improve pricing performance any further , the model that combines stochastic volatility and double exponential jump-diffusion (SVDEJD) may be more reasonable.
In the BS setting, the probability measure has a well-known analytic form , but, under stochastic volatility, it can only be obtained numerically. Monte Carlo simulation and the finite difference method are usually used to value the options. But, the two techniques require substantially more computing time and thus are difficult to be applied in real option pricing. Recently, being fast, accurate, and easy to implement, Fourier transforms have been widely used in valuing financial derivatives, for example, Carr and Madan  propose Fourier transforms with respect to log-strike price; Duffie et al.  offer a comprehensive survey that the Fourier methods are applicable to a wide range of stochastic processes; Carr and Wu  apply the transforms to time-changed Lévy processes and the class of generalized affine models. Hurd and Zhou  express the spread option payoff in terms of the gamma function and FFT technique. For an overview of option pricing using Fourier transforms, see Schmelzle .
The current paper extends the study of option pricing under the DEJD model in three ways. First, we propose a model which combines the double jumps and stochastic volatility. Second, using the martingale method, Fourier transform formula, and Feynman-Kac theorem, we obtain a closed-form solution for European call options pricing under the proposed model. Third, we obtain fast and accurate numerical solutions for European call options pricing by FFT technique.
The rest of the paper is organized as follows. Section 2 develops the underlying pricing model. Section 3 derives a closed-form solution for European call options pricing under the proposed model. Section 4 provides approximation solutions for European call options pricing by FFT technique. Section 5 numerically compares the density of the solutions to the alternative models and analyzes several effects on potion prices. Section 6 concludes. Applied program codes in Matlab package are presented in the appendix.
2. The Model
We consider an arbitrage-free, frictionless financial market where only riskless asset and risky asset are traded continuously up to a fixed horizon date . Let be a complete probability space with a filtration satisfying the usual conditions, that is, the filtration is continuous on the right and contains all -null sets. Suppose , are both standard Brownian motion, which is adapted, and has correlation with .
Let represent the price for a stock or a stock portfolio. Generally, instantaneous variance of asset returns in financial markets shows randomness; thus, the continuous part of the price process, defined as , is where is risk-free rate and is nonnegative constant, and suppose , which can be set equal to 1 without any loss of generality. The size of the diffusion component is determined by , which represents, absent of any jump occurring, the level of (stochastic) return variance attributable to diffusion variations. For tractability, let follow a square-root process: where nonnegative constants , and , respectively, reflect the speed of adjustment, the long-run mean, and the variation coefficient of , and suppose .
It has been suggested from extensive empirical studies that markets tend to have both overreaction and underreaction to various good or bad news. One may interpret the jump part of the model as the market response to outside news. Good or bad news arrives according to a Poisson process, and the asset price changes in response according to the jump size distribution. According to Kou , the jumps in the log-price are modeled as a sequence of i.i.d. nonnegative random variables that occur at times determined by an independent Poisson process with constant intensity such that has an asymmetric double exponential distribution with the density where denotes the indicator function, so equals 1 if , but 0 otherwise. are up-move jump and down-move jump, respectively. Except for which has correlation with , all sources of randomness, , and , are assumed to be independent.
Because of jumps and stochastic volatility, the risk-neutral probability measure is not unique. Following Naik and Lee  and Kou , by using the rational expectations argument with a HARA-type utility function for the representative agent, one can choose a particular risk-neutral measure so that the equilibrium price of an option is given by the expectation under this risk-neutral measure of the discounted option payoff. Throughout this paper, we assume that there exists a martingale probability measure being equivalent to . Let be the sum of all the jumps which occur up to and including time , , we have Obviously, is a -martingale. Finally, the price process is defined as
Remark 2.1. The model contains most existing models as special cases. For example, we obtain (1) the BS model by setting and ; (2) the SV model by setting ; (3) the DEJD model by setting .
Let , where is standard Brownian motion that is adapted, and independent of , , and random variables . From It’s formula, we have
3. A Closed-Form Solution of European Option Pricing
In this section, we derive closed-form solution of a European call option pricing under the SVDEJD model. For a European put option, we can obtain easily corresponding result by the put-call parity . For this purpose, we need the following results.
Lemma 3.1. Supposing the variance process follows (2.2) and are any complex, one has where
Proof. Let . Because of the affine structure of the variance process (2.2), we obtain that has a solution of the following form: From the Feynman-Kac formula, is the solution of the following backward Parabolic partial differential equation with the Cauchy problem: Putting (3.3) in (3.4), we have Solving (3.5), we can obtain the result of Lemma 3.1.
Lemma 3.2. Supposing the asset price follows (2.6) and is any complex, one has where
Proof. Let . Because is independent of , and, we have From (2.2) and (2.3), we have Because is a standard Brownian motion, we have Then, Let , and . From Lemma 3.1, we have From (3.8), (3.9), (3.10) and (3.12), we can obtain the required Lemma 3.2.
Lemma 3.3. Suppose is the characteristic function of ; then where
Theorem 3.4. Let denote the log of the strike price , , and the desired value of a -maturity call option with strike . Assume that, under , the underlying nondividend-paying stock price and its components are given by (2.1)–(2.5), is the characteristic function of , is the density of ; then the initial call value is written as where represents real part.
Proof. From the risk-neutral theory, we have Introducing a change of measure from to by a Radon-Nikodym derivative, we get With this new measure , the Fourier transform of is defined as Because of the no-arbitrage condition, we can obtain From the Fourier transform formula, the probability density for our model is given by Hence, Changing the order of integration, we have From (3.17), (3.20), and (3.23), we can obtain the required Theorem 3.4.
Remark 3.5. In (3.16), tends to not zero as goes to . Hence, is not (absolutely integrable) and a Fourier transform does not exist.
4. Fast Fourier Transform for European Option Pricing
Since the integrand in (3.16) is singular at the required evaluation point , the FFT cannot be applied directly to evaluate the integrals we mentioned above. Therefore, instead of solving for the risk-neutral exercise probabilities of finishing in-the-Money (ITM), Carr and Madan  introduce a new technique with the key idea to calculate the Fourier transform of a modified call option price with respect to the logarithmic strike price. With this specification and a FFT routine, a whole range of option prices can be obtained within a single Fourier inversion. In this section, we develop the numerical solutions of the prices by using the idea of Carr and Madan .
4.1. Fourier Transform of ITM and at-the-Money (ATM) Option Prices
This method is viable when is chosen in a way that the damped option price is well behaved. Damping the option price with makes it integrable for the negative axis . On the other hand, for the option prices increase by the exponential , which influences the integrability for the positive axis. A sufficient condition of to be integrable for both sides (square integrability) is given by being finite, that is, Thus we need , which is equivalent to Therefore, is well behaved when the moments of order of the underlying asset exist and are finite. If not all moments of exist, this will impose an upper bound on . We find that one quarter of this upper bound serves as a good choice for .
Using the trapezoid rule for the integral on the right-hand side of (4.1) and setting , an approximation for is The FFT returns values of , and we employ a regular spacing of size so that our values for are This gives us log-strike levels ranging from to , where
In order to apply FFT we define
To obtain an accurate integration with larger values of , we incorporate Simpson’s rule weightings into our summation. From (4.1)–(4.7)and Simpson’s rule weightings, we obtain ATM and ITM call value as where is the Kronecker delta function that is unity for and zero otherwise. The summation in (4.8) is an exact application of the FFT.
4.2. Fourier Transform of out-of-the-Money (OTM) Option Prices
In the previous section call values are calculated by an exponential function to obtain square integrable function whose Fourier transform is an analytic function of the characteristic function of the log-price. But, for very short maturities, the call value approaches its non analytic intrinsic value causing the integrand in the Fourier inversion to be high oscillate, and therefore difficult to integrate numerically. We introduce an alternative approach that works with time values only, which is quite similar to the previous approach. But in this case the call price is obtained via the Fourier transform of a modified time value, where the modification involves a hyperbolic sine function instead of an exponential function.
Let denote the time value of an OTM option, that is, for we have the put price for and for we have the call price. Scaling for simplicity, is defined by where is the risk-neutral density of the log-price . Let be the Fourier transform of : By considering a damping function , the time value of an option follows a Fourier inversion: where .
The use of the FFT for calculating OTM option prices is similar to (4.8). The only differences are that they replace the multiplication by with a division by and the function call to is replaced by a function call to .
5. Simulation Studies
In this section, to compare across the BS, DEJD, and SVDEJD models, we analyze the probability densities of these models. Then, we analyze mainly the impact of and volatility of volatility on option pricing under the SVDEJD model. For our FFT methods, we used points in our quadrature, implying a log-strike spacing of , which is adequate for practice. For the choice of the dampening coefficient in the transform of the modified call price, we used a value of . For the modified time value, we used . Other parameter values used in the computation are listed in Table 1. (We have used analytic moments to set plausible parameter values for the model. For a formal econometric estimator, one could use these moments to develop a generalized method of moments estimator within the framework of Hall and Inoue .)
5.1. Probability Densities under Alternative Models
We compare the probability densities of the SVDEJD model, the BS model, and the DEJD model to verify the rationality of our model. Suppose is the characteristic function of and the probability density of our model. From FFT algorithm, can be approximated by The density has the mean and variance given by
Figure 1 shows the figures of the probability density , compared with the normal density with the same mean and variance given by (5.2). The first figure compares the overall shapes of the densities of the SVDEJD model and the BS model, the second one details the shapes around the peak areas, and the last one shows the right tail. From Figure 1, we can see that the leptokurtic and skewness feature of the density of our model is quite evident. Moreover, additional numerical plots suggest that the feature of skewness becomes more significant if increases, which is impossible for the DEJD model.
We also compare the short-term and long-term densities of the SVDEJD model, the BS model, and the DEJD one. Figure 2 shows their densities under months and years. From Figure 2, we can see that the SVDEJD model and the DEJD model generate virtually identical densities for short-term options, with a slight departure occurring between the two densities in the upper tail. This means that differential pricing performance between the SVDEJD model and the DEJD model is unlikely to occur when they are applied to price short-term OTM puts and that only when they are applied to deep ITM puts (and deep OTM calls) can differences be observed between these models. Yet, compared to the BS model density, the densities of the two models are distinctly different: they all have leptokurtic and skewness feature. Therefore, the two models can potentially correct the BS model’s tendency to underprice deep OTM puts and overprice deep OTM calls. The long-term density curves in Figure 2 still show significantly different pricing structures between the BS and its two alternatives. But, more importantly, the densities of the SVDEJD model and the DEJD model also exhibit different shapes now. The SVDEJD density has higher peak and assigns more weight to both the entire lower tail and the far upper tail, but less weight to those payoffs than the DEJD.
Our simulation studies have demonstrated that the SVDEJD model has better performance than the DEJD one on pricing long-term options, while both the DEJD model and the SVDEJD model have better performance than the BS model.
5.2. Effects of the Main Parameter on Option Values
In Table 2, we use the SVDEJD model to examine the effects of volatility of volatility , the correlation coefficient , exercise price , and maturity time on option values. We analyze the prices of three-month call options and two-year call options. To examine the effect of the negative correlation coefficient, we have calculated the model with , , and . The prices for three-month call options associated with volatility of volatility and are relatively close. With , the largest price difference is 0.0149; with , the largest price difference is only 0.0008; with , the largest price difference is 0.0213. However, the difference is significantly larger when longer-time horizons such as two-year call options are valued. With , the largest price difference is an increase on 0.0149 to 0.0545 and the effect of volatility of volatility is an increase for ITM calls and a decrease for OTM calls. With , the largest price difference is increase of 0.0008 to 0.01 and the effect of volatility of volatility is a small decreas that is negligible for option values. With , the largest price difference is an increase of 0.0213 to 0.0633 and the effect of volatility of volatility is a decrease for ITM and ATM calls and an increase for OTM calls. The correlation parameter has several effects depending on the relation between the strike price and the current stock price. A negative tends to produce higher values for ITM calls and lower values for OTM money calls.
We have also compared the model with the BS model, which can be interpreted as a first-order approximation with no jumps, and . A common practice is to set the implied volatility in the BS model so that the model matches the price for the option with a strike price closest to the current stock price. For some comparisons not reported here, we have set the implied volatility in the BS model so that it matches the price generated by the stochastic volatility model for an ATM option. The BS implied volatility is very close to the expected volatility under the risk-neutral distribution when short-term options are valued. When longer-term options are used, there is a significant difference between the BS implied volatility and the expected volatility. As an approximation, the BS model tends to undervalue ITM calls and overvalue OTM calls.
The SVDEJD model incorporates several important features of stock returns. We derive a closed-form solution for European call options in the model by using the martingale method, Fourier inversion transform formula, and Feynman-Kac theorem. Using FFT, we obtain fast and accurate numerical solution to European option under the model. The comparison of densities of the alternative models shows that the SVDEJD model has better pricing performance on long-time options. An analysis of the model reveals that volatility of volatility and the correlation coefficient have significant impact on option values, especially long-time option, stock returns are negatively correlated with volatility, and these negative correlations are important for option valuation.
A.1. Matlab Codes for ITM and ATM Options Pricing by FFT
function CV =inSVDexpJ(ata1, ata2, lamta, sigma, thetav, alphav,rho, sigmav, r, p, s0, v0, strike, T)x0 = log(s0)alpha = 2.55N = 4096c = 600eta = c/Nb =pi/etau = [0:N-1]*etalamda = 2*b/Nposition = (log(strike) + b)/lamda + 1v = u - (alpha+1)*ik=p*ata1/(ata1-1)+(1-p)*ata2/(ata2+1)-1l=p*ata1./(ata1-i*v)+(1-p)*ata2./(ata2+i*v)m=sqrt((alphav-i*v*rho*sigma*sigmav).∧2+i*v.*(1-i*v)*(sigma * sigmav)∧2)n=2*m+(alphav-m-i*v*rho*sigma*sigmav).*(1-exp(-m*T))A=(2*m./n).∧(2*thetav/sigmav∧2)B=i*v*x0+(thetav*(alphav-m)*T)/sigmav∧2-(i*v*rho*sigma*thetav*T)/sigmav…+lamta*T*(l-i*v*k-1)+i*v*r*TC=(i*v.*(i*v-1)*sigma∧2.*(1-exp(-m*T)))./ncharFunc=A.*exp(C*v0+B)ModifiedCharFunc = charFunc*exp(-r*T)./(alpha∧2+alpha - u.∧2 + i*(2*alpha +1)*u)SimpsonW = 1/3*(3 + (−1).∧[1:N] - [1, zeros(1,N-1)])FftFunc = exp(i*b*u).*ModifiedCharFunc*eta.*SimpsonWpayoff = real(fft(FftFunc))CallValueM = (exp(-log(strike)*alpha))’*payoff/piformat shortCV= CallValueM(round(position)).
A.2. Matlab Codes for OTM Options Pricing by FFT
function CV=outSVDexpJ(ata1,ata2,lamta,sigma,thetav,alphav, rho,sigmav,r,p,s0,v0,strike,T)x0 = log(s0)alpha = 1.55N=4096c = 600eta = c/Nb =pi/etau = [0:N-1]*etalamda = 2*b/Nposition = (log(strike) + b)/lamda + 1w1 = u-i*alphaw2 = u+i*alphav1 = u-i*alpha -iv2 = u+i*alpha -ik=p*ata1/(ata1-1)+(1-p)*ata2/(ata2+1)-1l1=p*ata1./(ata1-i*v1)+(1-p)*ata2./(ata2+i*v1)m1=sqrt((alphav-i*v1*rho*sigma*sigmav).∧2+i*v1.*(1-i*v1)*(sigma*sigmav)∧2)n1=2*m1+(alphav-m1-i*v1*rho*sigma*sigmav).*(1-exp(-m1*T))A1=(2*m1./n1).∧(2*thetav/sigmav∧2)B1=i*v1*x0+(thetav*(alphav-m1)*T)/sigmav∧2-(i*v1*rho*sigma*thetav*T)/sigmav…+lamta*T*(l1-i*v1*k-1)+i*v1*r*TC1=(i*v1.*(i*v1-1)*sigma∧2.*(1-exp(-m1*T)))./n1charFunc1=A1.*exp(C1*v0+B1)ModifiedCharFunc1 = exp(-r*T)*(1./(1+i*w1)…-exp(r*T)./(i*w1) - charFunc1./(w1.∧2 - i*w1))l2=p*ata1./(ata1-i*v2)+(1-p)*ata2./(ata2+i*v2)m2=sqrt((alphav-i*v2*rho*sigma*sigmav).∧2+i*v2.*(1-i*v2)*(sigma*sigmav)∧2)n2=2*m2+(alphav-m2-i*v2*rho*sigma*sigmav).*(1-exp(-m2*T))A2=(2*m2./n2).∧(2*thetav/sigmav∧2)B2=i*v2*x0+(thetav*(alphav-m2)*T)/sigmav∧2-(i*v2*rho*sigma*thetav*T)/sigmav…+lamta*T*(l2-i*v2*k-1)+i*v2*r*TC2=(i*v2.*(i*v2-1)*sigma∧2.*(1-exp(-m2*T)))./n2charFunc2=A2.*exp(C2*v0+B2)ModifiedCharFunc2 = exp(-r*T)*(1./(1+i*w2)- exp(r*T)./(i*w2)…- charFunc2./(w2.∧2 - i*w2))ModifiedCharFuncCombo = (ModifiedCharFunc1 - ModifiedCharFunc2)/2SimpsonW = 1/3*(3 + (−1).∧[1:N] - [1, zeros(1,N-1)])FftFunc = exp(i*b*u).*ModifiedCharFuncCombo*eta.*SimpsonWpayoff = real(fft(FftFunc))CallValueM = payoff/pi/sinh(alpha*log(strike))format shortCV= CallValueM(round(position)).
This work is supported by the National Natural Science Foundation of China (Grant no. 11171266) and the Science Plan Foundation of the Education Bureau of Shanxi Province (nos. 11JK0491 and 11JK0493). The authors are grateful to the reviewers as well as the editor for their valuable comments and suggestions, which led to a substantial improvement of the paper.
- F. Black and M. S. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, no. 3, pp. 637–654, 1973.
- R. C. Merton, “Option pricing when underlying stock returns are discontinuous,” Journal of Financial Economics, vol. 3, no. 1-2, pp. 125–144, 1976.
- R. Cont, “Empirical properties of asset returns: stylized facts and statistical issues,” Quantitative Finance, vol. 1, no. 2, pp. 223–236, 2001.
- R. Cont, J. da Fonseca, and V. Durrleman, “Stochastic models of implied volatility surfaces,” Economic Notes, vol. 31, no. 2, pp. 361–377, 2002.
- E. Eberlein and U. Keller, “Hyperbolic distributions in finance,” Bernoulli, vol. 1, no. 3, pp. 281–299, 1995.
- D. B. Madan, P. P. Carr, and E. C. Chang, “The variance gamma process and option pricing,” European Finance Review, vol. 2, no. 1, pp. 79–105, 1998.
- D. Duffie, J. Pan, and K. Singleton, “Transform analysis and asset pricing for affine jump-diffusions,” Econometrica, vol. 68, no. 6, pp. 1343–1376, 2000.
- L. H. Steven, “A closed-form solution for options with stochastic volatility with applications to bond and currency options,” Review of Financial Studies, vol. 6, no. 2, pp. 327–343, 1993.
- D. S. Bates, “Jumps and stochastic volatility: exchange rate processes implicit in deutsche mark options,” Review of Financial Studies, vol. 9, no. 1, pp. 69–107, 1996.
- J. Keppo, X. Meng, S. Shive, and M. Sullivan, “Modelling and hedging options under stochastic pricing parameters,” in Proceedings of the Industrial and Operations Engineering, University of Michigan at Ann Arbor, 2003.
- S. G. Kou, “A jump-diffusion model for option pricing,” Management Science, vol. 48, no. 8, pp. 1086–1101, 2002.
- S. G. Kou and H. Wang, “Option pricing under a double exponential jump diffusion model,” Management Science, vol. 50, no. 9, pp. 1178–1192, 2004.
- J. R. Birge and V. Linetsky, Handbooks in Operations Research and Management Science, vol. 15, North-Holland, 2008.
- G. Bakshi, C. Cao, and Z. Chen, “Pricing and hedging long-term options,” Journal of Econometrics, vol. 94, no. 1-2, pp. 277–318, 2000.
- T. W. Epps, Pricing Derivative Securities, World Scientific, Singapore, 2nd edition, 2007.
- P. P. Carr and D. B. Madan, “Option valuation using the fast Fourier transform,” Journal of Computational Finance, vol. 2, no. 4, pp. 61–73, 1999.
- P. P. Carr and L. Wu, “Time-changed Lévy processes and option pricing,” Journal of Financial Economics, vol. 71, no. 1, pp. 113–141, 2004.
- T. R. Hurd and Z. Zhou, “A Fourier transform method for spread option pricing,” SIAM Journal on Financial Mathematics, vol. 1, pp. 142–157, 2010.
- M. Schmelzle, “Option pricing formulae using Fourier transform: theory and application,” 2010, http://pfadintegral.com.
- V. Naik and M. Lee, “General equilibrium pricing of options on the market portfolio with discontinuous returns,” Review of Financial Studies, vol. 3, no. 4, pp. 493–521, 1990.
- A. R. Hall and A. Inoue, “The large sample behaviour of the generalized method of moments estimator in misspecified models,” Journal of Econometrics, vol. 114, no. 2, pp. 361–394, 2003.
Copyright © 2012 Su-mei Zhang and Li-he Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.