Numerical and Analytical Methods for Variational Inequalities and Related Problems 2013
View this Special IssueResearch Article  Open Access
Fast Fourier Transform Based Power Option Pricing with Stochastic Interest Rate, Volatility, and Jump Intensity
Abstract
Firstly, we present a more general and realistic doubleexponential jump model with stochastic volatility, interest rate, and jump intensity. Using FeynmanKac formula, we obtain a partial integrodifferential equation (PIDE), with respect to the moment generating function of log underlying asset price, which exists an affine solution. Then, we employ the fast Fourier Transform (FFT) method to obtain the approximate numerical solution of a power option which is conveniently designed with different risks or prices. Finally, we find the FFT method to compute that our option price has better stability, higher accuracy, and faster speed, compared to Monte Carlo approach.
1. Introduction
In financial markets, option pricing problem becomes important for market practitioners ranging from individuals to financial institutes. The classic methodology for option pricing is developed by Black and Scholes [1]. In their model, the stock price follows a geometric Brownian motion and the option pricing formula in no arbitrage. However, a large number of empirical examples show that the initial assumption of their model generates high pricing deviations between theory and practice, for example the volatility clustering in [2] and the incompleteness of markets in [3]. Due to those drawbacks, some market practitioners begin to propose new underlying asset processes to fit empirical facts and develop new option types to control and reduce risks.
Doubleexponential jump model with stochastic volatility, interest rate, and jump intensity could overcome these drawbacks of the BlackScholes model in [2–8]. For example, SantaClara and Yan [4, 5] find that innovations to two risks, respectively, denoted by diffusion volatility process and jump intensity process in their model, have affected the expected return in the stock market. The empirical results show that those two processes are largely uncorrelated and do not support models that make jump intensity vary with the level of diffusive volatility. Based on the results of their theory, we believe that two components of risk are more conducive to the optimization of our model. An excellent contribution of our model is developing the model of SantaClara and Yan in [4, 5], combing stochastic volatility, interest rate, and meanreverting jump intensity processes. Thus, our model better corresponds with the real market than other constant jump intensity models.
As for new option types, compared with other exotic options, power option with considerable operability is simpler and more convenient in practical financial market [9]. Unfortunately, research on power option pricing is literally nil according to our best knowledge. Therefore, in this paper, we define a power option and creatively study and explain how the power coefficient affects the valuation of option.
However, it is difficult to get an exact solution expression of option price. Several numerical methods have been proposed to solve it, including finite difference methods by Zhang and Wang in [10, 11], Monte Carlo simulation by Green in [12] and Glasserman in [13] and the FFT method by Carr and Madan in [14]. However, the finite difference method and Monte Carlo simulation are difficult to be applied in option pricing because they require substantially more computing time [14–17]. Compared to those approaches, the FFT is more effective and simpler to compute in [14–17]. Therefore, in this paper, we employ the FFT to compute the valuation of power potion.
In conclusion, the prominent contribution in our paper is modeling a more general and realistic model with stochastic interest rate, volatility, and jump intensity. Defining the power option, discussing the effects of power coefficient, and using the FFT to obtain its numerical results are another three contributions.
This paper proceeds as follows. In Section 2, we build our model which links doubleexponential jump, stochastic interest rate, volatility, and jump intensity. The characteristic function is derived in Section 3. In Section 4, we discuss the approximate solutions for power option pricing by the FFT. Numerical examples and discussions of the power coefficient are shown in Section 5. Conclusions are in Section 6.
2. Our Model
Let be a probability space where is the filtration generated by the Brownian motion and jump process at time , and is a riskneutral probability. The underlying asset price at time is given by where , , , and are Brownian motions; the correlation coefficient of and is , and other motions are independent on each other; is a Poisson process with stochastic intensity , and jump size is a random variable; is average jump amplitude; is the stochastic interest rate and is stochastic volatility of the underlying asset return; the parameters , , and and the meanreverting rates , , and are positive constants; the constants , , and are the longterm interest rate, volatility, and jump intensity, respectively.
We suppose that jump size has an asymmetric double exponential distribution , where , are probability of upmove jump and downmove jump, respectively. and are mean of positive jumps and negative jumps, respectively. Therefore, .
3. The Characteristic Function
Given the log asset price process in (1), it is possible to obtain the characteristic function of at time . Under the riskneutral measure , firstly, we define an explicit expression of the moment generating function (MGF) of at time , and the complexvalued characteristic function is given by .
From (3), we find that can be interpreted as a contingent claim that pays off at time . Therefore, it is possible to get the MGF and then obtain the characteristic function in following theorem.
Theorem 1. The characteristic function of under the measure is given by where
Proof. From (1), FeynmanKac formula gives the following PIDE for the characteristic function,
We conjecture that the solution of PDE (6) is
with initial conditions , , , , , and .
We consider the integral term in (6):
where .
Substituting (7) and (8) into (6) yields
If we denote , then (9) leads to a system of seven ordinary differential equations:
We consider the second ODE (11), in which the general solution can be derived for the Riccati equation. Making the substitution,
Hence,
The ODE (11) has a general solution of the form
where
According to (15) and initial condition , and are two constants to be determined from boundary conditions,
We obtain that the solutions of (21) is
Therefore, the exact solution of (11) is
Substituting (23) into (14) and using initial conditions , we have
The OEDs for , , , and are solved by analogy. Hence
Using , it is easy to get Theorem 1.
4. Power Option Pricing Using the Fast Fourier Transform
4.1. Power Option
Under the riskneutral measure , the price of the power call option at time with strike price and maturity date is given by
As a European option, the valuation of power at initial time is
Remark 2. In particularly, when , the power option degenerates into the vanilla option.
4.2. Fast Fourier Transform
Let denote the log of the strike price , and let be the desired value of an initial power call option with strike price . The riskneutral density of the price at time is . From (27), the initial value of power call option is related to the riskneutral density by
Here we employ the price of a zero coupon bond at time to express the stochastic discounted process where Then (28) can be written as
Note that tends to as tends to , and hence the call pricing function (31) is not square integrable. According to Wong and Lo [15], we now consider the modified call price function defined by The Fourier transform of is defined by Then the inverse transform of is given by Then Using the Trapezoid rule for the integral in (35), the value of is approximated as where
The FFT returns values of and for a regular spacing size of where is a power of 2. The value for is which corresponds to log strike prices ranging from to , where .
Substituting (37) into (36) yields here we denote .
With Simpson's rule weightings, the value of European call option is where
5. Numerical Examples and Discussion
In this section, we present several numerical experiments to compare numerical differences between the FFT and Monte Carlo simulation and discuss the effects of power coefficient. Now, we perform some Monte Carlo simulations for the endowment dynamics by using the discrete scheme of system (1) as follows: where , , , and ; is the time interval; , , and are samples from the standard normal distribution with correlation coefficient 0; , if , , is the th jump interval which obeys an exponential distribution with parameter .
The option prices are calculated by Monte Carlo simulation using
In our numerical examples by the FFT, we take and which lead to the log strike spacing of and the damping coefficient . Other parameters are denoted as , , , , , , , , , , , , , , , , , , , , and . The results are presented in Tables 1, 2, and 3.



Numerical examples presented that the FFT is faster than Monte Carlo method with different parameters. To different strike prices, the FFT takes around approximate 340 seconds to produce 4098 option prices in three numerical experiments. However, Monte Carlo technique takes around 1260 seconds. On the level of the accuracy, comparing with Monte Carlo, we show that the percentage differences of the FFT are less than 0.1% in Tables 1, 2, and 3. Therefore, the FFT is more accurate and faster than Monte Carlo.
To examine effects of the power coefficient , we calculate valuations with , , and . Firstly, we can see that with the rising, the gaps between and become greater in Figures 1, 2, and 3. Secondly, accordoing to different , we find that valuations of option are increasing with the growth of . It is because the fluctuations and risks of the underlying asset price are enhancing with the growth of , and in order to compensate risk, the option seller need gets higher price of option. It means that the seller can design different power options with different risks or prices by defining different the power coefficients. Therefore, power option is one of convenient options which are designed easily.
6. Conclusion
In this paper, we study a more general and realistic double exponentialjump model which links the stochastic interest rate, volatility, and jump intensity. If we decide the constant interest rate, volatility, or jump intensity, our model can be degenerated into many previous models, for example, Black and Scholes’ model in [1], SantaClara and Yan’ model in [4, 5], Carr and Madan’s model in [14], Zhang and Wang’s model in [10, 11], and Zhang and Wang’s model in [17].
Moreover, we obtain the characteristic function of the log asset price, which is important to not only take the FFT in [14–17] but also derive the option pricing formula, for example, in [18]. Then, in our paper, a power option is defined and its valuations are computed with the FFT by characteristic function. Through numerical examples, comparing the results produced from the FFT and the Monte Carlo simulations of power option valuation, we find that the FFT is approximate with Monte Carlo approach but has better stability, higher accuracy, and faster speed in pricing power option. In addition, we show that the seller can design different power options with different risks or prices by defining different the power coefficients.
Acknowledgment
This work is supported by the Fundamental Research Funds for the Central Universities (JBK130401).
References
 F. Black and M. Scholes, “The pricing of options and corporate liabilities,” The Journal of Political Economy, vol. 81, no. 3, pp. 637–654, 1973. View at: Google Scholar
 K. Chen and S. H. Poon, “Variance swap premium under stochastic volatility and selfexciting jumps,” SSRN, 2200172, 2013. View at: Google Scholar
 R. Cont and P. Tankov, Financial Modelling with Jump Processes, Financial Mathematics, Chapman & Hall, Boca Raton, Fla, USA, 2004. View at: MathSciNet
 P. SantaClara and S. Yan, “Crashes, volatility, and the equity premium: Lessons from S&P 500 options,” Review of Economics and Statistics, vol. 92, no. 2, pp. 435–451, 2010. View at: Publisher Site  Google Scholar
 P. SantaClara and S. Yan, “Jump and volatility risk and risk premia: a new model and lessons from S&P 500 options,” National Bureau of Economic Research, no. w10912, 2004. View at: Google Scholar
 A. Amin, “Hybrid equity, FX, and interest rate models with stochastic volatility and jump diffusion,” FX, and Interest Rate Models with Stochastic Volatility and Jump Diffusion, 2010. View at: Google Scholar
 C. Chang, C. D. Fuh, and S. K. Lin, “A tale of two regimes: theory and empirical evidence for a Markovmodulated jump diffusion model of equity returns and derivative pricing implications,” Journal of Banking & Finance, vol. 37, no. 8, pp. 3204–3217, 2013. View at: Google Scholar
 J. P. F. Charpin and M. Cummins, “Fast Fourier transform option pricing: efficient approximation methods under multifactor stochastic volatility and jumps,” Topics in Numerical Methods for Finance, vol. 19, pp. 115–137, 2012. View at: Google Scholar
 P. G. Zhang, “An introduction to exotic options,” European Financial Management, vol. 1, no. 1, pp. 87–95, 1995. View at: Google Scholar
 K. Zhang and S. Wang, “Pricing options under jump diffusion processes with fitted finite volume method,” Applied Mathematics and Computation, vol. 201, no. 12, pp. 398–413, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. Zhang and S. Wang, “A computational scheme for options under jump diffusion processes,” International Journal of Numerical Analysis and Modeling, vol. 6, no. 1, pp. 110–123, 2009. View at: Google Scholar  MathSciNet
 P. J. Green, “Reversible jump Markov chain Monte Carlo computation and Bayesian model determination,” Biometrika, vol. 82, no. 4, pp. 711–732, 1995. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Glasserman, Monte Carlo Methods in Financial Engineering, vol. 53, Springer, New York, NY, USA, 2004. View at: MathSciNet
 P. Carr and D. Madan, “Option valuation using the fast Fourier transform,” Journal of Computational Finance, vol. 2, no. 4, pp. 61–73, 1999. View at: Google Scholar
 H. Y. Wong and Y. W. Lo, “Option pricing with mean reversion and stochastic volatility,” European Journal of Operational Research, vol. 197, no. 1, pp. 179–187, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Pillay and J. G. O'Hara, “FFT based option pricing under a mean reverting process with stochastic volatility and jumps,” Journal of Computational and Applied Mathematics, vol. 235, no. 12, pp. 3378–3384, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Zhang and L. Wang, “Fast Fourier transform option pricing with stochastic interest rate, stochastic volatility and double jumps,” Applied Mathematics and Computation, vol. 219, no. 23, pp. 10928–10933, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 A. L. Lewis, Option Valuation under Stochastic Volatility, Finance Press, Newport Beach, Calif, USA, 2000. View at: MathSciNet
Copyright
Copyright © 2013 Jiexiang Huang et al. 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.