Time Discretized Variational Iteration Method for the Stochastic Volatility Process with Jumps ()
1. Introduction
The stochastic volatility process with jumps (SVPJ) is governed by the stochastic differential equation (SDE) [1] [2] [3]
(1.1)
where parameters are detailed in Table 1, and
being the Poision process.
The SVPJ was introduced by Andrew Matytsin in 1999 [4]. Special cases of the SVPJ are the Heston’s stochastic volatility and Merton’s jump diffusion models. The classical Black-Scholes model is the foundation of the Merton’s model which has the advantage of replicating asset returns that are discontinuous. The jumps, however, do not depend on the diffusion. Bate [5] was first to merge the above two powerful models in order to ascertain the risk neutral mean for the variances and asset value. Thus, Bate’s stochastic volatility with jump is obtained by setting
, with Brownian motion being correlated as
(1.2)
Eraker et al., [6] following the Bate’s procedure developed a model introducing jumps in volatility as given below:
(1.3)
The authors considered the strategy called the “likelihood—base estimation” to provide estimates of jump times, parameters and spot volatility [7] [8].
Over the years, many works by mathematicians on the numerical evaluation of these models have been on the increase. However, the variational iteration method in a discretized setting has not been implemented in literature in the evaluation of these models.
Thus, the objective of this paper is to seek the numerical solution of the SVPJ through a discretized variational iteration method (DVIM) to obtain sample paths for the state variable and variance process at different timesteps and replications
in order to estimate the expected jump times. In this approach, the standard variational iteration method as proposed by He [9] [10] [11] [12] is being discretized in time within the time interval
with equidistant step size
in order to estimate the expected jump times.
The method is quite robust in its analysis as it discretizes the state and variance parameters in a stochastic jump process. The method is easy to implement without making any weak assumptions or perturbation to produce convergent solutions that depict the degree of randomness in the area of option price in real time.
The parameters
and
are related by [13] [14]
(1.4)
2. Materials and Methods
Discretized Variational Iteration Method
Early numerical studies coupled with theoretical investigations by Wright [15], Clements and Anderson [16] and others [17] [18] [19] [20] [21] revealed that not every heuristic discrete time approximation of a stochastic differential equation converges absolutely to the solution process as the maximum step size tends to zero. As revealed, one cannot simply employ deterministic solvers of ordinary differential equations for stochastic differential equations. Thus, the need for a careful and systematic investigation through the time discretized variational iteration method for the solution process of (1.1) cannot be over-emphasized.
Now, we consider a discretization
with
(2.1)
over the time interval
with
(2.2)
being the equidistant step size.
An Ito process
with initial condition
, has the differential form
, (2.3)
is called the Ito stochastic differential equation [22], where
and
are the drift and diffusion coefficients respectively.
Now, the mathematical formulation of the proposed method for the SVPJ process is as follows:
We construct a correction functional for the SVPJ process as
(2.4)
subject to the initial conditions
(2.5)
with
, (2.6)
(2.7)
and
, (2.8)
for
.
The
and
are the Lagrange multipliers which can be obtained through the variational theory. Here,
and
are restricted variables. We take the variation
on both sides of (2.4) as follows:
(2.9)
Solving (2.9), we have,
and
.
. (2.10)
The scheme (2.4) only produces approximations at the times of discretization. For intermediate iterations, values from previous discretization point could be used. It should be noted that
in (2.7) are normally distributed random variables
with means and variances given as [23] [24] [25]
(2.11)
3. Sample Paths Approximation
Usually, solutions of stochastic differential equations are difficult to achieve, so often times simulations are used to illustrate the trajectories of the solution process at the different time discretization for
. If the solution is known explicitly, then, mathematical software can be used to compute the absolute error using the absolute error criterion, which is given as [26]
and
(3.1)
which quantify the sample path closeness at the extreme of the time interval
. To this end, the discretized variational iteration method for the SVPJ model corresponding to N different simulations of the Wiener and Ito processes are implemented.
4. Numerical Simulations
Here, we use the scheme (2.4) for
to construct an SVPJ process and estimate the expected jump times with the given variable parameters in line with [7] [8] as follows:
The sample paths as shown in Figure 1(a) and Figure 1(b) respectively, generates 10 replications with 30 timesteps for the state variable and variance process. The average computed expected jump times for the state variable and variance process is 20.40344029.
Similarly, using different parameters
and
, with every other parameters remains the same as in the previous case. The sample paths for the state variable and variance process are shown in Figure 2(a) and Figure 2(b), respectively. It generates 10 replications with 30 timesteps. The average computed expected jump times is 52.40607869.
5. Discussion of Results
The DVIM has been applied successfully for the SVPJ process. Results were generated and presented graphically via MAPLE 18 software. The sample paths were generated with 10 replications and 30 timesteps to quantify the closeness at the extreme of the time interval
. It was observed that the average computed expected jump times, which is 20.40344029, for the state variable and variance process is moderated by the parameters
and
. For instance, at each
Figure 1. (a) Sample path for the state variable with 10 replications and 30 timesteps at different iterates; (b) Sample path for the variance process with 10 replications and 30 timestep at different iterates.
Figure 2. (a) Sample path for the state variable with 10 replications and 30 timesteps at different iterates. (b) Sample path for the variance process with 10 replications and 30 timestep at different iterates.
iterate for
, with parameters
and
, the stochastic jumps for the state variable are less cluttered compared to the variance process with maximum volatility cluttering as shown in Figure 1(a) and Figure 1(b). The value 20.40344029, in terms of option pricing implies that volatility cluttering and mean reversion between equity returns experiences less uphill movement in price assert. In contrast, at each iterate for
, with parameters when
and
, the state variable cluttered maximally compared to the variance process with less volatility cluttering as shown in Figure 2(a) and Figure 2(b), respectively. The average computed expected jump time is 52.40607869. In terms of option pricing, the value 52.40607869 suggests a better bargain compared to the value 20.40344029 due to the fact that it yields less volatility rate. We tend to point out that the negativity of the correlation parameter does not imply the curves in the figures must appear negative (see [1] [3] for details).
6. Conclusion
This paper focused on the numerical application of the SVPJ model. We have successfully discretized the variational iteration method for the purpose of generating sample paths with different replications and timesteps for computing the average expected jumps times experience in the volatility cluttering. Results obtained are satisfactory and are in line with those available in literature [13] [14] and very relevant in the area of option pricing and Forex trading.