# A Shannon Wavelet Method for Pricing Forward Starting Options under the Double Exponential Jump Framework with Two-Factor Stochastic Volatilities

Authors: Peiyi Li
Journal: Mathematical Problems in Engineering
 (2022)

---

## Abstract

This paper provides an efficient pricing method for forward starting options based on Shannon wavelet expansions. Specifically, we derive the pricing results under a more realistic stock model that incorporates the double exponential jump, stochastic jump intensity, and two-factor stochastic volatilities to capture various features observed in financial markets. We obtain the characteristic function related to the payoff function; then, the options can be well evaluated by the Shannon wavelet method. Numerical experiments show that this method is fast and accurate compared to the Monte Carlo simulation. Finally, we study the influences of changing some important parameters to further illustrate the robustness and stability of the model.


---

## Body

## 1. Introduction

A forward starting option belongs to the family of path-dependent options, it is purchased and paid for now but becomes active later at a starting datet0 with an expiration date T further in the future. Forward starting options are often used by companies as employee stock options to motivate staff. In addition, insurance companies often choose forward starting options as hedging tools to reduce the risk contained in life insurance products.A closed-form solution for the price of forward starting options was first solved by Rubinstein [1] based on the Black–Scholes (BS) framework [2]. However, many studies have shown that the BS model fails to capture some key features in financial markets, such as nonconstant volatility of stock price. Therefore, many researchers have introduced stochastic volatility in price dynamics to explain volatility smile. So far, the most common models are proposed by Heston [3], Stein and Stein [4], and Hull and White [5]. Based on this, Kruse and Nögel [6], Amerio [7], and Haastrecht and Pelsser [8] considered pricing forward starting options by adding stochastic volatility to financial markets. However, few researchers such as Da Fonseca et al. [9] have shown that single-factor stochastic volatility models do not perform well when trying to capture the long-term structure of implied volatility. Then, Christoffersen et al. [10] proposed the double Heston model, allowing the volatility dynamics to be driven by two uncorrelated stochastic processes to provide a better fit to the empirical implied volatility.In addition to stochastic volatility, extreme price movements exist in underlying price dynamics, especially when financial crises happen. Merton [11] first incorporated log-normal jumps in asset price. Kou and Wang [12] has proposed another model assuming that the jump size follows an asymmetric double exponential distribution, it can explain extreme movements in asset price and capture the leptokurtic feature as well. Moreover, prior empirical research [13] revealed a low level of correlation between stock volatility and jump intensity; then, Huang et al. [14–16] assumed in their work that the stochastic volatility and jump intensity are governed by separate processes. Therefore, in this paper, we consider combining double stochastic volatilities, asymmetric double exponential jump with stochastic intensity in the modeling framework to valuate forward starting options.While the above stochastic structures work well in effectively capturing various market features and obtaining more accurate option prices as a result, how to find a closed-form solution under these complicated models becomes a demanding task. The numerical pricing methods for models with stochastic volatility can be classified into three main groups [14, 17]: the Monte Carlo simulation, numerical integration methods, and partial (integro-) differential equation methods. In general, the price of European-style options is determined by the discounted expectation of payoff function. Calculation of the expectation requires knowledge about the probability density function (PDF) of stock price. Although the PDF may not be available in most sophisticated price processes, the Feynman–Kac theorem connects the numerical solution of a partial (integro-) differential equation to the expectation of payoff function. Then, we can calculate the characteristic function for relevant price processes, which is the same as the Fourier transform of the PDF.The integration methods are used when the characteristic function of the stock price is known, and a well-known example is the fast Fourier transform (FFT) method [18]. However, the convergence of FFT depends on a damping factor, so it is often unstable. Based on Fourier-cosine series expansion, Fang and Oosterlee [17] proposed a more efficient alternative called the COS method. It can achieve an exponential convergence and has been widely used in the valuation of options [19–22]. One of the drawbacks of this method is that cosine expansions are formed on a global basis and exhibit periodicity behavior, and errors may accumulate near the domain boundaries.Ortiz-Gracia and Oosterlee [23] proposed another valuation method for plain vanilla options called the Shannon Wavelet Inverse Fourier Technique (SWIFT). It performs well with more flexibility and accuracy as wavelet functions can capture local features better than cosine series. Another main improvement of the SWIFT method is that it does not require integral truncation ranges because the local wavelet basis can adaptively indicate if the error is under a defined tolerance.The main goal of this paper is to provide an efficient method to price forward starting options. The contributions of this paper are mainly two-fold. Firstly, we utilize the SWIFT method and derive a closed-form solution. Some numerical integration methods we mentioned above have been used to price forward starting options. Kruse and Nögel [6] did the pioneering research, and they applied the FFT method under the Heston model. Zhang and Geng [20] applied the COS method to price forward starting options under the double Heston model. In this paper, we choose the SWIFT method to valuate forward starting options, due to the established accuracy and robustness of Shannon wavelets in option pricing, as demonstrated in a wealth of existing literature [14, 16, 24, 25]. Secondly, this paper proposes a more realistic model that incorporates two-factor stochastic volatilities, asymmetric double exponential jump with stochastic jump intensity to capture various features observed in financial markets. It is an extension of the work by Christoffersen et al. [10], Kou and Wang [12] and Huang et al. [16].The rest of the paper is arranged as follows: Section2 sets up the model of stock price dynamics. In Section 3, we derive the characteristic function of the log-asset price. Section 4 first briefly introduces the Shannon wavelet method; based on this, we derive a formula for pricing European-style forward starting options. Numerical results and sensitivity analysis are given in Section 5, and we conclude the paper in Section 6.


## 2. The Model

We assume thatΩ,ℱt,Q is a complete probability space where four Brownian motions B1st, B2st, B1vt, B2vt, and Bλt0≤t≤T are defined, ℱt is a filtration, and Q is a risk-neutral measure. The stock price St, two volatilities v1t and v2t, and jump intensity λt are specified as follows:(1)dSt=r−λtmStdt+∑j=12vjtStdBjst+yt−1StdNt,dv1t=κ1θ1−v1tdt+ε1v1tdB1vt,dv2t=κ2θ2−v2tdt+ε2v2tdB2vt,dλt=κλθλ−λtdt+ελλtdBλt,where r is a constant interest rate and κj, θjj=1,2, and θλ represent the mean-reverting speeds, long-term volatilities, and long-term jump intensity, respectively. The stochastic processes vjtj=1,2 and λt satisfy the Feller conditions [26], i.e., 2κjθj≥εj2j=1,2,λ to make the relevant process remain positive. B1st and B1vt, and B2st and B2vt are two pairs of Brownian motions with correlation coefficients ρ1 and ρ2, respectively.Nt is a Poisson process with stochastic jump intensity λt and jump size yt=eY, where Y is a random variable and follows an asymmetric double exponential distribution with the density function fYy=p1/η1e−1/η1y1y≥0+q1/η2e1/η2y1y〈0,0〈η1〈1,η2〉0, where p,q≥0,p+q=1 are up-move and down-move probabilities and m represents average jump amplitude, m:=EeY−1=p/1−η1+q/1+η2. We further suppose that the processes B1st, B2st, B1vt, B2vt, and Bλt are all independent of Nt and Y.Remark 1.
Model equation (1) contains some special models as follows:(1)
The geometric Brownian motion by settingY=κj=θj=εj=0j=1,2,λ;(2)
The Heston model withY=κj=θj=εj=0j=2,λ;(3)
The double Heston model whenY=κλ=θλ=ελ=0;(4)
The double exponential jump model withκj=θj=εj=0j=1,2.


## 3. Derivation of Characteristic Function

In this paper, we consider the forward starting options with starting datet0 and expiration date T0≤t0≤T. Let K denote the strike price. In this section, we first obtain the characteristic function at time T for log-price Xt:=nSt, then derive the characteristic function of lnST/St0−lnK.Lemma 1.
Assume that the underlying asset follows equation (1), then under the risk-neutral measure Q, the conditional characteristic function φ· is given by(2)φu,t=EQeiuXT|ℱt=expiuXt+Au,t+∑j=12Bju,tvjt+Cu,tλt,where(3)Au,t=−riuT−t+κ1θ1ε12α+T−t−2lnα+eζ1T−t−α−2ζ1+κ2θ2ε22β+T−t−2lnβ+eζ2T−t−β−2ζ2+κλθλελ2γ+T−t−2lnγ+eζ3T−t−γ−2ζ3,B1u,t=−ui+ueζ1T−t−1α+eζ1T−t−α−,B2u,t=−ui+ueζ2T−t−1β+eζ2T−t−β−,Cu,t=2Λu−miueζ3T−t−1γ+eζ3T−t−γ−,with(4)α±=κ1−ρ1ε1iu±ζ1,β±=κ2−ρ2ε2iu±ζ2,γ±=κλ±ζ3,ζ1=κ1−ρ1ε1iu2+ε12ui+u,ζ2=κ2−ρ2ε2iu2+ε22ui+u,ζ3=κλ2+2ελ2miu−Λu,Λu=p1−iuη1+q1+iuη2−1.Proof.
See Appendix A.Lemma 2.
Given that the processVt follows Cox–Ingersoll–Ross model:(5)dVt=κθ−vtdt+εvtdBt,the expectation of expiuvt is given by(6)EQexpiuvt|ℱ0=1−iuε21−e−κt2κ−2κθε2expiue−κtv01−iuε21−e−κt/2κ.Proof.
See [27].Theorem 1.
Suppose that the stock price follows (1), the characteristic function of lnST/St00≤t0≤T is given by(7)Φu=expAu,t0+∑j=12Bj′+C′,where(8)Bj′=−2κjθjεj2lnhj+Bju,t0e−κjt0vj0hj,j=1,2,hj=1−Bju,t0εj21−e−κjt02κj,j=1,2,C′=−2κλθλελ2lnhλ+Cu,t0e−κλt0λ0hλ.Proof.
See Appendix B.


## 4. Pricing Method

### 4.1. Shannon Wavelet Series Expansion

Letℒ2R denote the space of square integrable functions. We start with a family of closed subspaces Vjj∈Z of ℒ2R with the following properties:(9)i⋯⊂V−2⊂V−1⊂V0⊂V1⊂V2⋯,ii∪j∈ZVj¯=ℒ2R,∩j∈ZVj=0,iiifx∈Vj⇔f2x∈Vj+1.Definition 1.
A functionϕ∈ℒ2R is said to generate a multi-resolution analysis (MRA) if it generates a sequence of closed subspaces Vjj∈Z that satisfy (i), (ii), and (iii), and the set ϕj,kx:=2j/2ϕ2jx−kk∈Z forms an orthonormal basis of Vj.Definition 2.
A functionϕ∈ℒ2R is called a basic scaling function or father wavelet if it generates a MRA.
We define a projectionPm:ℒ2R⟶Vm.
Observe that at a chosen level of resolutionm, every function f∈ℒ2R can be approximated:(10)fx≈Pmfx=∑k∈ZCm,kϕm,kx,where Pmf⋅ converges to f· as m⟶∞, cm,k=∫−∞+∞fxϕm,kx¯dx are called scaling coefficients.
The scaling functions of Shannon wavelets inVm is computed as(11)ϕm,kx=2m/2ϕ2mx−k=2m/2sinc2mx−k=2m/2sin2mx−kπ2mx−kπ.


### 4.2. Pricing Method for Forward Starting Options

Here, we consider the final payoff function of forward starting options asαST/St0−K+, with α=1 for call and α=−1 for put, where t0 is the starting time, T is the maturity time, and K denotes the strike price. Then, the price of forward starting options can be expressed as(12)vt0,T=e−rTEQαSTSt0−K+.Letδ=lnST/St0−lnK; then, the option price equation (12) becomes(13)vt0,T=e−rTK∫−∞+∞αeδ−1+fδdδ,where fδ represents the probability density function (PDF) of δ.Based on the wavelets theory discussed above,fx can be expressed through Shannon wavelet expansions:(14)fx≈Pmfx=∑k∈Zcm,kϕm,kx≈∑k=k1k2cm,kϕm,kx,the second approximately equal sign holds because if fx tends to zero at infinite, cm,k satisfy limn⟶∞cm,k=0 (see [23] for details).Substituting equations (13) into (14) and exchanging the summation and integration,(15)vt0,T=e−rTK∫−∞+∞αeδ−1+∑k=k1k2cm,kϕm,kδdδ=e−rTK∑k=k1k2cm,kVm,k,where Vm,k=∫−∞+∞αeδ−1+ϕm,kδdδ.


#### 4.2.1. Coefficients Computation

Proposition 1.
(16)sinct=∏j=1+∞cosπt2j.Proposition 2.
(17)∏j=1Jcosπt/2j=12J−1∑j=12J−1cos2j−12Jπt.Proof.
The results of Propositions1 and 2 can be easily derived by two-fold duplication formula sin2a=2cosasina and product-to-sum formula 2cosacosb=cosa+b+cosa−b, respectively.
According to Propositions1 and 2, sinc(t) can be approximated as(18)sinct≈sinc∗t≔12J−1∑j=12J−1cos2j−12Jπt.
Then, we replaceϕ in (11) with (18),(19)cm,k=∫−∞+∞fsϕm,ks¯ds≈cm,k∗=2m/22J−1∑j=12J−1∫−∞+∞fscos2j−12Jπ2ms−kds.
Consider the Fourier transform offx:(20)f^u=∫−∞+∞eiuxfxdx,and Ref^x=∫−∞+∞fxcosuxdx, where Re· denotes the real part of argument, the approximate computation of the scaling coefficients cm,k can be obtained:(21)cm,k≈cm,k∗=2m/22J−1∑j=12J−1Ref^2j−1π2m2Jeikπ2j−1/2J.


#### 4.2.2. Pricing Formula

For formula (15), we truncated the integration range to a finite domain Dm=k1/2m,k2/2m without significant loss (see [23]), then(22)Vm,k≈Vm,k∗=∫Dmαeδ−1+ϕm,kδdδ=K2m/2∫Dm∩0,+∞eδ−1sinc2mδ−kdδ,for callK2m/2∫Dm∩−∞,01−eδsinc2mδ−kdδ,for put..If we define(23)∏1,kp,q=∫pqexcosCj2mx−kdx,∏2,kp,q=∫pqcosCj2mx−kdx,with Cj=2j−1/2Jπ; through a straight forward calculation, the above integrals can be expressed as(24)∏1,kp,q=Cj2m1+Cj2m2eqsinCj2mq−k−eqsinCj2mp−k+1Cj2meqcosCj2mq−k−epcosCj2mp−k,∏2,kp,q=1Cj2msinCj2mq−k−sinCj2mp−k.Thereafter,(25)Vm,k∗=∑j=12J−1∏1,kmaxk1,02m,k22m−∏2,kmaxk1,02m,k22m,ifk20,0,ifk2≤0,for a call option, and(26)Vm,k∗=∑j=12J−1∏1,kk12m,mink2,02m−∏2,kk12m,mink2,02m,ifk10,0,ifk1≥0,for a put.Based on what we have discussed above, we can obtain the formula on the price on forward starting options:(27)vt0,T=e−rTK∑k=k1k2cm,kVm,k,with cm,k approximated by (21) and Vm,k approximated by equations (25) and (26).


## 4.1. Shannon Wavelet Series Expansion

Letℒ2R denote the space of square integrable functions. We start with a family of closed subspaces Vjj∈Z of ℒ2R with the following properties:(9)i⋯⊂V−2⊂V−1⊂V0⊂V1⊂V2⋯,ii∪j∈ZVj¯=ℒ2R,∩j∈ZVj=0,iiifx∈Vj⇔f2x∈Vj+1.Definition 1.
A functionϕ∈ℒ2R is said to generate a multi-resolution analysis (MRA) if it generates a sequence of closed subspaces Vjj∈Z that satisfy (i), (ii), and (iii), and the set ϕj,kx:=2j/2ϕ2jx−kk∈Z forms an orthonormal basis of Vj.Definition 2.
A functionϕ∈ℒ2R is called a basic scaling function or father wavelet if it generates a MRA.
We define a projectionPm:ℒ2R⟶Vm.
Observe that at a chosen level of resolutionm, every function f∈ℒ2R can be approximated:(10)fx≈Pmfx=∑k∈ZCm,kϕm,kx,where Pmf⋅ converges to f· as m⟶∞, cm,k=∫−∞+∞fxϕm,kx¯dx are called scaling coefficients.
The scaling functions of Shannon wavelets inVm is computed as(11)ϕm,kx=2m/2ϕ2mx−k=2m/2sinc2mx−k=2m/2sin2mx−kπ2mx−kπ.


## 4.2. Pricing Method for Forward Starting Options

Here, we consider the final payoff function of forward starting options asαST/St0−K+, with α=1 for call and α=−1 for put, where t0 is the starting time, T is the maturity time, and K denotes the strike price. Then, the price of forward starting options can be expressed as(12)vt0,T=e−rTEQαSTSt0−K+.Letδ=lnST/St0−lnK; then, the option price equation (12) becomes(13)vt0,T=e−rTK∫−∞+∞αeδ−1+fδdδ,where fδ represents the probability density function (PDF) of δ.Based on the wavelets theory discussed above,fx can be expressed through Shannon wavelet expansions:(14)fx≈Pmfx=∑k∈Zcm,kϕm,kx≈∑k=k1k2cm,kϕm,kx,the second approximately equal sign holds because if fx tends to zero at infinite, cm,k satisfy limn⟶∞cm,k=0 (see [23] for details).Substituting equations (13) into (14) and exchanging the summation and integration,(15)vt0,T=e−rTK∫−∞+∞αeδ−1+∑k=k1k2cm,kϕm,kδdδ=e−rTK∑k=k1k2cm,kVm,k,where Vm,k=∫−∞+∞αeδ−1+ϕm,kδdδ.


### 4.2.1. Coefficients Computation

Proposition 1.
(16)sinct=∏j=1+∞cosπt2j.Proposition 2.
(17)∏j=1Jcosπt/2j=12J−1∑j=12J−1cos2j−12Jπt.Proof.
The results of Propositions1 and 2 can be easily derived by two-fold duplication formula sin2a=2cosasina and product-to-sum formula 2cosacosb=cosa+b+cosa−b, respectively.
According to Propositions1 and 2, sinc(t) can be approximated as(18)sinct≈sinc∗t≔12J−1∑j=12J−1cos2j−12Jπt.
Then, we replaceϕ in (11) with (18),(19)cm,k=∫−∞+∞fsϕm,ks¯ds≈cm,k∗=2m/22J−1∑j=12J−1∫−∞+∞fscos2j−12Jπ2ms−kds.
Consider the Fourier transform offx:(20)f^u=∫−∞+∞eiuxfxdx,and Ref^x=∫−∞+∞fxcosuxdx, where Re· denotes the real part of argument, the approximate computation of the scaling coefficients cm,k can be obtained:(21)cm,k≈cm,k∗=2m/22J−1∑j=12J−1Ref^2j−1π2m2Jeikπ2j−1/2J.


### 4.2.2. Pricing Formula

For formula (15), we truncated the integration range to a finite domain Dm=k1/2m,k2/2m without significant loss (see [23]), then(22)Vm,k≈Vm,k∗=∫Dmαeδ−1+ϕm,kδdδ=K2m/2∫Dm∩0,+∞eδ−1sinc2mδ−kdδ,for callK2m/2∫Dm∩−∞,01−eδsinc2mδ−kdδ,for put..If we define(23)∏1,kp,q=∫pqexcosCj2mx−kdx,∏2,kp,q=∫pqcosCj2mx−kdx,with Cj=2j−1/2Jπ; through a straight forward calculation, the above integrals can be expressed as(24)∏1,kp,q=Cj2m1+Cj2m2eqsinCj2mq−k−eqsinCj2mp−k+1Cj2meqcosCj2mq−k−epcosCj2mp−k,∏2,kp,q=1Cj2msinCj2mq−k−sinCj2mp−k.Thereafter,(25)Vm,k∗=∑j=12J−1∏1,kmaxk1,02m,k22m−∏2,kmaxk1,02m,k22m,ifk20,0,ifk2≤0,for a call option, and(26)Vm,k∗=∑j=12J−1∏1,kk12m,mink2,02m−∏2,kk12m,mink2,02m,ifk10,0,ifk1≥0,for a put.Based on what we have discussed above, we can obtain the formula on the price on forward starting options:(27)vt0,T=e−rTK∑k=k1k2cm,kVm,k,with cm,k approximated by (21) and Vm,k approximated by equations (25) and (26).


## 4.2.1. Coefficients Computation

Proposition 1.
(16)sinct=∏j=1+∞cosπt2j.Proposition 2.
(17)∏j=1Jcosπt/2j=12J−1∑j=12J−1cos2j−12Jπt.Proof.
The results of Propositions1 and 2 can be easily derived by two-fold duplication formula sin2a=2cosasina and product-to-sum formula 2cosacosb=cosa+b+cosa−b, respectively.
According to Propositions1 and 2, sinc(t) can be approximated as(18)sinct≈sinc∗t≔12J−1∑j=12J−1cos2j−12Jπt.
Then, we replaceϕ in (11) with (18),(19)cm,k=∫−∞+∞fsϕm,ks¯ds≈cm,k∗=2m/22J−1∑j=12J−1∫−∞+∞fscos2j−12Jπ2ms−kds.
Consider the Fourier transform offx:(20)f^u=∫−∞+∞eiuxfxdx,and Ref^x=∫−∞+∞fxcosuxdx, where Re· denotes the real part of argument, the approximate computation of the scaling coefficients cm,k can be obtained:(21)cm,k≈cm,k∗=2m/22J−1∑j=12J−1Ref^2j−1π2m2Jeikπ2j−1/2J.


## 4.2.2. Pricing Formula

For formula (15), we truncated the integration range to a finite domain Dm=k1/2m,k2/2m without significant loss (see [23]), then(22)Vm,k≈Vm,k∗=∫Dmαeδ−1+ϕm,kδdδ=K2m/2∫Dm∩0,+∞eδ−1sinc2mδ−kdδ,for callK2m/2∫Dm∩−∞,01−eδsinc2mδ−kdδ,for put..If we define(23)∏1,kp,q=∫pqexcosCj2mx−kdx,∏2,kp,q=∫pqcosCj2mx−kdx,with Cj=2j−1/2Jπ; through a straight forward calculation, the above integrals can be expressed as(24)∏1,kp,q=Cj2m1+Cj2m2eqsinCj2mq−k−eqsinCj2mp−k+1Cj2meqcosCj2mq−k−epcosCj2mp−k,∏2,kp,q=1Cj2msinCj2mq−k−sinCj2mp−k.Thereafter,(25)Vm,k∗=∑j=12J−1∏1,kmaxk1,02m,k22m−∏2,kmaxk1,02m,k22m,ifk20,0,ifk2≤0,for a call option, and(26)Vm,k∗=∑j=12J−1∏1,kk12m,mink2,02m−∏2,kk12m,mink2,02m,ifk10,0,ifk1≥0,for a put.Based on what we have discussed above, we can obtain the formula on the price on forward starting options:(27)vt0,T=e−rTK∑k=k1k2cm,kVm,k,with cm,k approximated by (21) and Vm,k approximated by equations (25) and (26).


## 5. Numerical Results

In this section, some numerical experiments are performed by utilizing the SWIFT method to price forward starting options. We first discuss the convergence of errors and find an appropriate level of resolutionm. Then, we display some numerical results to show the performance of the SWIFT method. Finally, we explore the price sensitivity on changing model parameters.


### 5.1. Error Convergence

We choose the appropriate level of resolutionm by two error convergence analyses. Firstly, we follow the line of [16] to compute the tail mass that the characteristic function Φ not recovered, which can be approximated by(28)14πT2mπΦ−2mπ+Φ2mπ.Figure1 shows the loss of the tail mass for characteristic function across a range of m.Figure 1
Error of approximation for characteristic function.In addition, Huang et al. [16] and Ortiz-Gracia and Oosterlee [23] have discussed that one advantage of utilizing Shannon wavelet is that once the scaling coefficients cm,k are derived, and the area under the approximated density function can be easily calculated as(29)A=12mcm,k12+∑k1kk2cm,k+cm,k22.Hence, we further investigate the lost area under probability density function (PDF) through equation (29). Figure 2 presents the error when approximating the density function with Shannon wavelet expansions.Figure 2
Error of approximation for density function.We can observe from Figure2 that the differences is negligible for m≥5. Meanwhile, as can be seen in Figure 1, m=5 is an appropriate level of approximation within a loss tolerance level of TOL=10−90. Therefore, we choose m=5 for numerical experiments.


### 5.2. Numerical Analysis

Here, following Fang [17], the approximate truncated range can be calculated asμ1−Lμ2+μ4,μ1+Lμ2+μ4, where L=10 and μn is defined as the nth derivative of the underlying characteristic function, i.e., μn=1/in∂nlnΦu/∂un∣u=0. Besides, J≔log2πN with N≔maxk1,k2.We use the pricing formula derived above to do numerical experiments and analyse the speed of computation and the accuracy and the stability of the SWIFT method.The values of parameters are listed in Table1 for all numerical examples. All of these numerical examples were performed in Python 3.7. Also, the computer we used equips an Intel Core i7 CPU with a 2.2 GHz processor.Table 1
Values of parameters.
Parametersrκ1θ1ε1ρ1κ2θ2ρ2Value0.1210.60.2−0.560.4−0.5Parametersε2κλθλελpqt0TValue0.950.060.50.40.60.52Here, we use three different approaches (SWIFT, COS, and Monte Carlo simulation) to calculate the price of forward starting call options. For Monte Carlo method, we use 100,000 numbers of simulations with 200 numbers of time steps. The COS method has been proved to be highly efficient and accurate in a wealth of literature [17, 19, 22, 28]. Monte Carlo simulation is a typical numerical method in the domain of option pricing, and it can be flexibly used for various exotic options and thus becomes one of the most common approaches in practice. Therefore in this paper, we choose the results of the COS method as the benchmark and compare the performance between SWIFT and Monte Carlo methods.Table2 displays the resulting prices and average CPU time of the three methods. The pricing results show that the price differences between the SWIFT method and the COS method are negligible compared to the price differences between Monte Carlo simulation and COS, which means that the SWIFT method is more accurate than Monte Carlo simulation. For Monte Carlo method, we use 100,000 numbers of simulations with 200 numbers of time steps, and it takes more than 100 seconds to calculate. The SWIFT and COS make big improvement in computation speed, and the average CPU time they need is significantly less than Monte Carlo.Table 2
Prices of forward starting options.
KSwiftCosMonte Carlo0.80.4811220.4811340.4753530.850.4628780.4628980.4556580.90.445620.4456490.4507000.950.4292770.4293160.42448710.4137840.4138330.4037351.050.3990820.3991430.3866371.10.3851170.3851910.3968951.150.3718380.3719270.3656901.20.3592010.3593060.356996Average CPU time (s)0.1390410.060573105.449163Then, we examine the pricing error of SWIFT and Monte Carlo with different expiration timeT=1,1.5,2 and different strike prices K=0.8,0.9,1,1.1,1.2. Table 3 shows the result. The error of using the SWIFT method are much lower than those computed using Monte Carlo simulation for all expiration time and strike prices. It demonstrates that the SWIFT method is more stable than Monte Carlo simulation.Table 3
Differences of the option prices.
K0.80.911.11.2T = 1Swift3.94E − 054.44E − 059.68E − 061.38E − 053.99E − 05Monte Carlo9.21E − 031.60E − 027.73E − 036.34E − 034.99E − 03T = 1.5Swift6.81E − 053.22E − 051.07E − 056.15E − 062.29E − 05Monte Carlo4.21E − 031.53E − 027.70E − 039.21E − 036.58E − 03T = 2Swift1.20E − 052.90E − 054.90E − 057.40E − 051.05E − 04Monte Carlo5.78E − 035.05E − 031.01E − 021.17E − 022.31E − 03


### 5.3. Sensitivity Analysis

We will conduct sensitivity analyses to study the impact of some parameters on pricing forward starting options. These analyses can provide evidence for robustness and stability of the SWIFT method under different market conditions. We mainly study the sensitivity of option prices derived by the SWIFT method to changes in the following parameters: (1) (Figure3) instantaneous volatility of v1t: ε1, (2) (Figure 4) the speeds of mean-reverting: κ2, (3) (Figure 5) long-term jump intensity: θλ, and (4) (Figure 6) reciprocal of the mean value for the negative jumps: η2.Figure 3
Prices with respect toε1.Figure 4
Prices with respect toκ2.Figure 5
Prices with respect toθλ.Figure 6
Prices with respect toη2.


## 5.1. Error Convergence

We choose the appropriate level of resolutionm by two error convergence analyses. Firstly, we follow the line of [16] to compute the tail mass that the characteristic function Φ not recovered, which can be approximated by(28)14πT2mπΦ−2mπ+Φ2mπ.Figure1 shows the loss of the tail mass for characteristic function across a range of m.Figure 1
Error of approximation for characteristic function.In addition, Huang et al. [16] and Ortiz-Gracia and Oosterlee [23] have discussed that one advantage of utilizing Shannon wavelet is that once the scaling coefficients cm,k are derived, and the area under the approximated density function can be easily calculated as(29)A=12mcm,k12+∑k1kk2cm,k+cm,k22.Hence, we further investigate the lost area under probability density function (PDF) through equation (29). Figure 2 presents the error when approximating the density function with Shannon wavelet expansions.Figure 2
Error of approximation for density function.We can observe from Figure2 that the differences is negligible for m≥5. Meanwhile, as can be seen in Figure 1, m=5 is an appropriate level of approximation within a loss tolerance level of TOL=10−90. Therefore, we choose m=5 for numerical experiments.


## 5.2. Numerical Analysis

Here, following Fang [17], the approximate truncated range can be calculated asμ1−Lμ2+μ4,μ1+Lμ2+μ4, where L=10 and μn is defined as the nth derivative of the underlying characteristic function, i.e., μn=1/in∂nlnΦu/∂un∣u=0. Besides, J≔log2πN with N≔maxk1,k2.We use the pricing formula derived above to do numerical experiments and analyse the speed of computation and the accuracy and the stability of the SWIFT method.The values of parameters are listed in Table1 for all numerical examples. All of these numerical examples were performed in Python 3.7. Also, the computer we used equips an Intel Core i7 CPU with a 2.2 GHz processor.Table 1
Values of parameters.
Parametersrκ1θ1ε1ρ1κ2θ2ρ2Value0.1210.60.2−0.560.4−0.5Parametersε2κλθλελpqt0TValue0.950.060.50.40.60.52Here, we use three different approaches (SWIFT, COS, and Monte Carlo simulation) to calculate the price of forward starting call options. For Monte Carlo method, we use 100,000 numbers of simulations with 200 numbers of time steps. The COS method has been proved to be highly efficient and accurate in a wealth of literature [17, 19, 22, 28]. Monte Carlo simulation is a typical numerical method in the domain of option pricing, and it can be flexibly used for various exotic options and thus becomes one of the most common approaches in practice. Therefore in this paper, we choose the results of the COS method as the benchmark and compare the performance between SWIFT and Monte Carlo methods.Table2 displays the resulting prices and average CPU time of the three methods. The pricing results show that the price differences between the SWIFT method and the COS method are negligible compared to the price differences between Monte Carlo simulation and COS, which means that the SWIFT method is more accurate than Monte Carlo simulation. For Monte Carlo method, we use 100,000 numbers of simulations with 200 numbers of time steps, and it takes more than 100 seconds to calculate. The SWIFT and COS make big improvement in computation speed, and the average CPU time they need is significantly less than Monte Carlo.Table 2
Prices of forward starting options.
KSwiftCosMonte Carlo0.80.4811220.4811340.4753530.850.4628780.4628980.4556580.90.445620.4456490.4507000.950.4292770.4293160.42448710.4137840.4138330.4037351.050.3990820.3991430.3866371.10.3851170.3851910.3968951.150.3718380.3719270.3656901.20.3592010.3593060.356996Average CPU time (s)0.1390410.060573105.449163Then, we examine the pricing error of SWIFT and Monte Carlo with different expiration timeT=1,1.5,2 and different strike prices K=0.8,0.9,1,1.1,1.2. Table 3 shows the result. The error of using the SWIFT method are much lower than those computed using Monte Carlo simulation for all expiration time and strike prices. It demonstrates that the SWIFT method is more stable than Monte Carlo simulation.Table 3
Differences of the option prices.
K0.80.911.11.2T = 1Swift3.94E − 054.44E − 059.68E − 061.38E − 053.99E − 05Monte Carlo9.21E − 031.60E − 027.73E − 036.34E − 034.99E − 03T = 1.5Swift6.81E − 053.22E − 051.07E − 056.15E − 062.29E − 05Monte Carlo4.21E − 031.53E − 027.70E − 039.21E − 036.58E − 03T = 2Swift1.20E − 052.90E − 054.90E − 057.40E − 051.05E − 04Monte Carlo5.78E − 035.05E − 031.01E − 021.17E − 022.31E − 03


## 5.3. Sensitivity Analysis

We will conduct sensitivity analyses to study the impact of some parameters on pricing forward starting options. These analyses can provide evidence for robustness and stability of the SWIFT method under different market conditions. We mainly study the sensitivity of option prices derived by the SWIFT method to changes in the following parameters: (1) (Figure3) instantaneous volatility of v1t: ε1, (2) (Figure 4) the speeds of mean-reverting: κ2, (3) (Figure 5) long-term jump intensity: θλ, and (4) (Figure 6) reciprocal of the mean value for the negative jumps: η2.Figure 3
Prices with respect toε1.Figure 4
Prices with respect toκ2.Figure 5
Prices with respect toθλ.Figure 6
Prices with respect toη2.


## 6. Conclusion

In this paper, we use a novel valuation method to price forward starting options. We consider a model combining the asymmetric exponential jump, stochastic jump intensity, and two-factor stochastic volatilities to capture various features observed in financial markets and derive a more realistic pricing framework. We calculate the characteristic function related to the final payoff by applying Feynman–Kac formula, and then solve the differential equations. After deriving the characteristic function, the SWIFT method is applied to compute the pricing results of forward starting options. Numerical experiments are performed to show the efficiency compared to the Monte Carlo method. Finally, we investigate the impact of changing model parameters to the resulting option values, thus proving the robustness and stability of our model.



---

*Source: 4234627-2022-06-06.xml*