Next Article in Journal
A Look at Generalized Degenerate Bernoulli and Euler Matrices
Previous Article in Journal
Enhancing Precision in Large-Scale Data Analysis: An Innovative Robust Imputation Algorithm for Managing Outliers and Missing Values
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pricing Variance Swaps under MRG Model with Regime-Switching: Discrete Observations Case

1
School of Information Science and Technology, Jinan University, Guangzhou 510632, China
2
Shenzhen Campus, Jinan University, Shenzhen 518053, China
*
Authors to whom correspondence should be addressed.
Mathematics 2023, 11(12), 2730; https://doi.org/10.3390/math11122730
Submission received: 25 May 2023 / Revised: 12 June 2023 / Accepted: 13 June 2023 / Published: 16 June 2023
(This article belongs to the Section Financial Mathematics)

Abstract

:
In this paper, we creatively price the discretely sampled variance swaps under the mean-reverting Gaussian model (MRG model in short) with regime-switching asymmetric double exponential jump diffusion. We extend the traditional MRG model by further considering the trend of the financial market as well as a sudden and unexpected event of the market. This new model is meaningful because it uses observable Markov chains that represent market states to adjust its parameters, which helps capture the movement of the market and fluctuations in asset prices. By utilizing the characteristic function and the conditional transition characteristic function, we obtain analytical solutions for pricing formulae. Note that this is our first effort to provide the analytical solution for the ordinary differential equations satisfied by the Feynman–Kac theorem. To achieve this, we have developed a new methodology in Proposition 2 that involves dividing the sampling interval into more detailed switching and non-switching intervals. One significant advantage of our closed-form solution is its high computational accuracy and efficiency. Subsequent semi-Monte Carlo simulations will provide specific validation results.

1. Introduction

With the globalization of economies and intensifying competition in financial markets, more and more investors and scholars have realized the importance of hedging against volatility risk. Against this background, volatility derivatives have emerged as instruments that hedge against this volatility risk. Their most prominent feature is that they provide direct exposure to asset volatility without continuous delta-hedging [1]. Among the various volatility derivatives, variance or volatility swaps are considered the main and essential ones. Essentially, they are forward contracts for future realized volatility or variance. The long position of a volatility or variance swap pays a fixed delivery price at expiry and receives the floating amounts of annualized realized volatility or variance, whereas the short position is just the opposite [2]. In other words, the payoff of a variance swap at maturity T is v 2 ( T ) K × I , where v 2 ( T ) is the realized annual variance of the return, K is the strike price such that the initial value of the variance swap equal to zero for both long and short positions, and I is the notional amount of the swap in dollars per annualized volatility point squared.
In the literature, there are roughly two types of approaches for pricing volatility swaps or variance swaps: static replication methods and model-based stochastic methods. Static replication means that a variance swap could be replicated theoretically by using a portfolio of standard options and without specifying the volatility function. For example, Demeterfi [1] produced some of the most influential pioneering work that demonstrated how to use a standard option portfolio to replicate variance swaps. Obviously, the outstanding advantage of static replication is that it does not require a specific model to solve pricing problems [3], and this method has been widely used for the VIX index since 2003 [4]. However, according to Sahalia et al. [5], this method has drawbacks that make it unsuitable for the real market. These include the inability to promptly reflect jumping variance swap rates, as well as the need for constant adjustment in options. Therefore, we consider model-based stochastic methods in this paper.
The study of variance swaps pricing in stochastic models can be broadly classified into two approaches. The first method entails pricing under stochastic volatility with continuously sampled realized variance. It provides an effective way to reduce the computational difficulty, and recent work can be found in [6,7,8,9,10]. However, this approach is an idealization that assumes continuous sampling time. According to the research conducted by Zhu and Lian [11], using this method to estimate realized variance may lead to significant inaccuracies in variance swaps with longer tenors or less frequent sampling rates when compared to discretely sampled variance swaps traded in the OTC market. Fortunately, the second method was developed to overcome this limitation, which entails pricing the realized variance obtained directly from discrete sampling. Broadie and Jain [12] used discretely sampled data and various stochastic models to examine the fair strike prices of variance swaps and volatility swaps. They also compared their pricing formulae with a continuously sampled situation and found that as observation frequency increases towards infinity, the discrete price converges to the continuous price. Rujivan and Zhu [13] proposed a simplified analytical approach for pricing discretely sampled variance swaps when the realized variance is defined as the squared log return of the underlying price. Zhang [14] obtained an analytical solution for variance swaps under the MRG stochastic volatility model. Lian [15] first explained the conditions for discrete volatility swaps to exhibit convergence behavior. Recent work in this field can be seen in [16,17,18].
However, the regime-switching model is rarely considered in most of the literature cited above. It is worth noting that the long-term trends in the financial market, which are influenced by structural changes in macroeconomic conditions, can have a significant impact on prices. That is, incorporating market trends into regime considerations can help convert market volatility uncertainty into regime switching to some extent, and it can also better assist investors in capturing market risks and opportunities [19]. As Elliott and Lian [20] pointed out, different market states (good or bad economy) can be reflected by different regimes, and the parameters of the regime-switching model are modulated by an observable Markov chain. In other words, the regime-switching model can be described as a Markov-modulated model. Because the market is generally incomplete, the Esscher transform has been extensively studied as an effective way to determine the unique equivalent martingale measure (see [21,22,23]). To the best of our knowledge, although Markov-modulated Heston models have been widely applied and studied for pricing options [4,24,25,26,27], the pricing of discrete-sampling variance swaps under an Markov-modulated MRG model has not yet been investigated. Compared to the Heston model and its reformulated version, the MRG volatility model has the following highlighted strengths. First, it has the mean-reverting property that ensures the autocorrelation of volatility [28] when volatility clustering occurs. Second, it has a relatively explicit solution for pricing. Therefore, in this paper, we attempt to study the pricing of discrete-sampling variance swaps under the Markov-modulated MRG model for the first time. This research has significant meaning and value. It can enrich the existing theory of regime-switching models and help market participants better understand market trends.
At the same time, there is a part of the regime-switching model that has not been well explored, which is the underlying asset process with jump diffusion. Sudden and unexpected events in real-world financial markets can result in significant fluctuations in asset prices, which cannot be disregarded, especially when pricing volatility derivatives. Many researchers have examined the jump effect with numerical tests and have proven that it can increase the degree of matching to the reality of financial markets [29]. In particular, when the jump diffusion follows the asymmetric double exponential distribution proposed by Kou [30], the model can help capture market-implied volatility smiles and market jump behavior [31]. However, little work has been done on pricing variance swaps under a regime-switching model with asymmetric double exponential jump diffusion, especially the Markov-modulated MRG model.
Motivated by the reasons stated above, we introduce a novel model in this paper—a Markov-modulated MRG model with asymmetric double exponential jump diffusion—to price discretely sampled variance swaps. To our knowledge, this is the first attempt to apply this new framework to pricing discretely sampled variance swaps. Hopefully, this attempt to price variance swaps under this new model could provide some thinking for other researchers and fill the research gap. By utilizing the Feynman–Kac theorem and transition probability function, we can greatly simplify the pricing process and theoretically derive the general expression of the characteristic function. This in turn allows us to solve the pricing problem more effectively. The main contributions of this paper can be summarized as follows: First, one potential issue for the regime-switching model that incorporates various economic factors is that the state space of the Markov chain can become quite large. This can lead to significant computational complexity in pricing and may pose challenges in practice. For practical purposes, it is significant to reduce the additional burden caused by the Markov chain. In this paper, we creatively obtain analytic pricing formulae that can reduce the computational burden and give the analytical solution of the ordinary differential equations satisfied by the Feynman–Kac theorem. To achieve this, we propose a new method of dividing the sampling time into more detailed switching intervals and non-switching intervals. Our Proposition 2 presents a specific expression for this solution. Second, we have made a creative improvement to the semi-Monte Carlo simulation method that verifies the accuracy and efficiency of our pricing formulae. Specifically, we use independent and identically exponentially distributed random numbers to generate repeated simulation paths randomly, as opposed to the traditional method of generating repeated simulation paths. Hopefully, it can be viewed as a step forward for the semi-Monte Carlo simulation in the regime-switching field.

2. Model Setup

In this section, we will formally introduce the MRG model with regime-switching jump diffusion and the computation methods for calculating realized variance. Additionally, we will discuss the specific details of the regime-switching framework.

2.1. The Regime-Switching Frame and Markov Chain

Regime-switching models have generated considerable interest for their applications to a variety of financial problems. As the trend of the financial market significantly influences the movement of individual stocks, it is essential to enable the key parameters of the stock to react to general market fluctuations. The market regime is one possible way to reflect the state of the economy, the mood of investors, and other economically realistic factors. If the Markov chain is observable, then the state can be interpreted as the market regime [20,25,32]. In other words, we can let parameters of the stochastic modes be modulated by the Markov chain, then the market factors can be considered in the model. It can help researchers capture the market movement by representing different market regimes.
Here, we follow the notation introduced in [4]: X = { X ( t ) } is a continuous-time finite-state X : = { e 1 , e 2 , , e N } observable Markov chain on a complete probability ( Ω , F , P ) , where P presents the real-world probability measure, X represents N different regimes of the financial market and the ith state e i = { 0 , , 1 i , , 0 } . Moreover, the dynamic of X is given by
d X ( t ) = Q ( t ) X ( t ) d t + d M ( t ) ,
where { M ( t ) } is an R N -valued martingale increment process with respect to the natural filtration generated by X, and Q ( t ) = ( q i j ) i , j = 1 , 2 , N is a transition matrix. The element q i j presents the transition rate from state e i to e j , satisfying q i j > 0 , i j and q i i = j i q i j (it shows that the row sum of matrix Q ( t ) is always zero).

2.2. The MRG Model with Regime-Switching Jump Diffusion

Under the Markov-modulated MRG model with double exponential jump diffusion, the dynamic of the stock price process S ( t ) with a short-term volatility process v ( t ) is given by
d S ( t ) S ( t ) = μ X ( t ) d t + v ( t ) d W 1 ( t ) + d i = 1 N t ( C i 1 ) , d v ( t ) = k θ X ( t ) v ( t ) d t + σ d W 2 ( t ) .
All variables here are set on ( Ω , F , P ) unless otherwise specified, where P presents the real-world probability measure. Assume · , · is the inner product and for every i = μ , θ , λ , r , the elements in i : = ( i 1 , i 2 , , i N ) are non-negative.
In the first equation of (2), μ X ( t ) : = μ , X ( t ) is the appreciation rate; { W 1 ( t ) } is a standard Brownian motion; { N t } is a Poisson process with rate λ X ( t ) : = λ , X ( t ) ; and { C i } is a sequence of independent identically distributed (i.i.d) nonnegative random variables such that Y i = ln ( C i ) has an asymmetric double exponential distribution with the density
f ( y ) = p η 1 e η 1 y I y 0 + q η 2 e η 2 y I y < 0 , η 1 > 1 , η 2 > 0 ,
where p , q > 0 , p + q = 1 represent the probabilities of upward and downward jumps.
In the second equation of (2), θ X ( t ) : = θ , X ( t ) is the long-term average value of volatility. To simplify our calculation and deduce an analytic solution for pricing variance swaps, we regard the volatility σ and the mean recovery speed parameter k as constants; { W 2 ( t ) } is another standard Brownian such that d W 1 ( t ) d W 2 ( t ) = ρ d t , ρ ( 1 , 1 ) .
Let the dynamic of the bank account process B ( t ) follow
d B ( t ) = r X ( t ) B ( t ) d t , B ( 0 ) = 1
where r X ( t ) : = r , X ( t ) is the instantaneous market interest rate.
For simplicity, we also let F t S , F t V , and F t X be the natural filtrations generated by W 1 ( t ) , W 2 ( t ) , and Markov chain X ( t ) , respectively.
Remark 1. 
This is the new model that we propose. We take the appreciation rate, long-term average value of volatility, and instantaneous market interest rate to be Markov-modulated in a creative way. In other words, the values of these parameters are not fixed compared to the traditional MRG model. Additionally, we add an asymmetric double exponential jump-diffusion to the asset process and also let the jump intensity be Markov-modulated parameters.

2.3. Variance Swaps

The realized variance in variance swaps is typically calculated through discrete sampling over the period [ 0 , T ] in financial markets, where T represents maturity time. Two common methods for calculating realized variance are actual-return realized variance
v d 1 2 ( T ) = A F N i = 1 N S ( t i ) S ( t i 1 ) S ( t i 1 ) 2 × 100 2 ,
and log-return realized variance
v d 2 2 ( T ) = A F N i = 1 N ln 2 S ( t i ) S ( t i 1 ) × 100 2 ,
where S t k , k = 0 , , N denotes the kth observation of the stock price at time t k , which implies that the total number of observations is N. In addition, A F is an annualized factor with a range from 4 to 252, which can be formed as A F = 1 t k t k 1 = N T .

3. Change of Numeraire

3.1. Change of Risk-Neutral Measure Q ^

Because the market modeling by the observable Markov chain is incomplete in general [20], the need to seek for “proper” measure is badly needed. After plenty of studies, researchers have found that the Esscher transform is the tool that satisfies the requirement. Until now, there has been mature technology to calculate the equivalent martingale measure via the Esscher transform (see [4,22,33] for more details). Here, we adopt the assumption that the continuous and jump components of the underlying asset can be respectively interpreted as systematic and idiosyncratic risk [34]. In other words, the jump component cannot affect the pricing. Following the setting in [4] and the method used in [23], we can quickly derive that W ˜ i ( t ) , i = 1 , 2 is a standard Brownian with respect to risk-neutral measure Q ^ , where
d W ˜ 1 ( t ) = d W 1 ( t ) + μ X ( t ) r X ( t ) v ( t ) d t , d W ˜ 2 ( t ) = ρ d W ˜ 1 ( t ) + 1 ρ d W 1 ( t ) .
We then substitute (7) into (2), and the dynamic process of S ( t ) and v ( t ) under the equivalent martingale measure Q ^ can be rewritten as
d S ( t ) S ( t ) = r X ( t ) d t + v ( t ) d W ˜ 1 ( t ) + d i = 1 N ˜ t ( C i 1 ) , d v ( t ) = k θ ˜ X ( t ) v ( t ) d t + σ d W ˜ 2 ( t ) ,
where
θ ˜ X ( t ) = θ X ( t ) + μ X ( t ) r X ( t ) v ( t ) · ρ σ k , λ ˜ X ( t ) : = λ X ( t ) p + q .

3.2. Change of T-Forward Measure Q ^ T

The value of a variance swap at time t under the risk-neutral measure Q ^ is
V i ( X , t ) = E Q ^ e t T r X ( t ) d t v i 2 ( T ) K i × I | F t S F t V F t X , i = d 1 , d 2 .
It is worth noting that the interest rate r X ( t ) under a regime-switching frame is not a fixed constant in the usual sense, and it will be influenced by Markov chain X ( t ) . Therefore, we need to consider the dependence between discount factor e t T r X ( t ) d t and variance swap payment, which will greatly increase the difficulty of solving pricing problems; for example, the calculation of the fair strike prices will no longer be a simple calculation of
K i = E Q ^ v i 2 ( T ) F 0 X F 0 S F 0 V , i = d 1 , d 2 ,
but
K i = E Q ^ e t T r X ( t ) d t v i 2 ( T ) F 0 S F 0 V F 0 X E Q ^ e t T r X ( t ) d t F 0 S F 0 V F 0 X , i = d 1 , d 2 .
Fortunately, we find that this problem will be solved by changing the numeraire to a zero-coupon bond and derive the T-forward measure Q ^ T . According to Theorem 9.2.2 and Definition 9.4.1 in Shreve [35], we can change the risk-neutral measure Q ^ into Q ^ T , that is, the numeraire will be changed from the original money market account to a zero-coupon bond, and the T-forward measure Q ^ T can be defined as
Q ^ T ( A ) = 1 B ( 0 , T , X ) A e 0 T r X ( t ) d t d Q ^ , A F .
Let B ( t , T , X ) be a zero-coupon bond that pays one unit at maturity T, then it is given by [25,36,37]
B ( t , T , X ) = E Q ^ e t T r X ( t ) d t | F t S F t V F t X = exp { B ( T t ) } , X I e ,
where I e denotes N-dimensional unit vector, operator B : = d i a g r Q , d i a g r denotes the matrix with vector r on its diagonal, and Q denotes the transposition of Q. Therefore, in this paper, the volatility vector of zero-coupon bonds is zero. In other words, discounting the underlying assets with it will not affect its volatility vector. Now, we can present (9) as
V i ( X , t ) = E Q ^ e t T r X ( t ) d t v i 2 ( T ) K i × I F t S F t V F t X = B ( t , T , X ) E Q ^ T v i 2 ( T ) K i × I F t S F t V F t X ,
and obtain the fair strike price
K i = E Q ^ T v i 2 ( T ) F 0 X F 0 S F 0 V , i = d 1 , d 2 .

4. Pricing Variance Swaps

In this section, we will take the distribution characteristics of the new model as a breakthrough, and derive the pricing formulae of variance swaps under discrete sampling. Here, we employ the characteristic function method to accomplish this goal.
Now we are in a position to formally introduce the characteristic function method. Let y ( T ) = ln S ( T + Δ ) ln S ( T ) , and assume its characteristic function is given by
f ϕ ; t , T , Δ , v ( t ) = E Q ^ T [ e ϕ y ( T ) F t X F t S F t V ] .
It should be noted that the imaginary number i = 1 is absorbed in the parameter ϕ in (14), and the result will not be changed without the explicit expression of the imaginary number [20]. Similar to the method used in [4], we first consider the conditional characteristic function under the given Markov chain information F T + Δ X . When the information F T + Δ X is given, the Markov-modulated parameters θ ˜ X ( t ) , λ ˜ X ( t ) , and r X ( t ) are deterministic functions, which means that the information available at time T + Δ is sufficient to evaluate these parameters. Namely, we can simplify (14) as
f ϕ ; t , T , Δ , v ( t ) = E Q ^ T f ϕ ; t , T , Δ , v ( t ) F T + Δ X F t X .

4.1. Conditional Characteristic Function

In this section, we will focus on calculating the concrete expression of conditional characteristic function f ϕ ; t , T , Δ , v ( t ) F T + Δ X and give the following lemma and proposition.
Lemma 1. 
If the volatility process v ( t ) satisfies the dynamics (8), then it has an explicit expression as follows:
v ( t ) = e k t v ( 0 ) + k 0 t θ ˜ X ( u ) e k ( t u ) d u + σ 0 t e k ( t u ) d W ˜ 2 ( u ) .
Proof. 
Let
f ( t , b ) : = e k t v ( 0 ) + σ e k t b + k 0 t θ ˜ X ( u ) e k ( t u ) d u ,
and
B ( t ) : = 0 t e k u d W ˜ 2 ( u ) .
For the Ito–Doeblin formula, we shall need the following the partial derivatives of f ( t , b ) :
f t ( t , b ) = k θ ˜ X ( t ) d t k f ( t , b ) , f b ( t , b ) = σ e k t , f b b ( t , b ) = 0 .
We will also need the partial derivative of B ( t ) , which is d B ( t ) = e k t d W ˜ 2 ( u ) . Then, the Ito–Doeblin formula states that
d v ( t ) = d f ( t , B ( t ) ) = k θ ˜ X ( t ) k f ( t , b ) d t + σ d W ˜ 2 ( u ) ,
which shows that the dynamics (8) satisfies the expression (16).    □
Proposition 1. 
Given the information F T + Δ X , if the stochastic process follows (8) and the parameters satisfy 1 2 σ ^ τ 2 > C ( ϕ ; T ) , then the characteristic function is
f ( ϕ ; t , T , Δ , v ( t ) F T + Δ X ) = L 1 2 ( ϕ ; T , τ ) exp N ( ϕ ; T , τ ) + S ( ϕ ; T , τ ) e k τ v ( t ) + R ( ϕ ; T , τ ) e 2 k τ v 2 ( t ) .
with
L ( ϕ ; T , τ ) = 1 2 σ ^ τ 2 C ( ϕ ; T ) , N ( ϕ ; T , τ ) = F ( ϕ ; T ) + M 2 ( ϕ ; T , τ ) μ ^ 2 ( t , T ) L ( ϕ ; T , τ ) 2 σ ^ τ 2 L ( ϕ ; T , τ ) , S ( ϕ ; T , τ ) = M ( ϕ ; T ) μ ^ ( t , T ) L ( ϕ ; T , τ ) σ ^ τ 2 L ( ϕ ; T , τ ) , R ( ϕ ; T , τ ) = L ( ϕ ; T , τ ) + 1 2 σ ^ τ 2 L ( ϕ ; T , τ ) , M ( ϕ ; T , τ ) = μ ^ ( t , T ) + σ ^ τ 2 E ( ϕ ; T ) μ ^ τ = e k τ v ( t ) + μ ^ ( t , T ) , σ ^ τ 2 = σ 2 2 k ( 1 e 2 k τ ) , μ ^ ( t , T ) = k t T θ ˜ X ( u ) e k ( T u ) d u , τ = T t .
where i ( ϕ ; T ) , i = C , E , F satisfies the following ordinary differential equations:
d F d t + 1 2 σ 2 E 2 + k θ ˜ E + σ 2 C + r ˜ ϕ + λ ˜ E Q ^ T [ e ϕ Y 1 ] = 0 , F ( ϕ ; T + Δ ) = 0 , d E d t + 2 k θ ˜ C + ( 2 σ 2 C k + ρ σ ϕ ) E = 0 , E ( ϕ ; T + Δ ) = 0 , d C d t + 2 σ 2 C 2 + ( 2 ρ σ ϕ 2 k ) C + 1 2 ϕ ( ϕ 1 ) = 0 , C ( ϕ ; T + Δ ) = 0 .
Proof. 
The proof of this proposition is presented in Appendix A.    □

4.2. Characteristic Function

In this section, we will consider the effect of Markov chain X ( t ) on the characteristic function. Using Formula (17) in the previous section, we can derive that
f ϕ ; t , T , Δ t , v ( t ) = E Q ^ T f ϕ ; t , T , Δ , v ( t ) F T + Δ X F t X = E Q ^ T L 1 2 ( ϕ ; T , τ ) exp R ( ϕ ; T , τ ) e 2 k τ v 2 ( t ) + S ( ϕ ; T , τ ) e k τ v ( t ) + N ( ϕ ; T , τ ) | F t X = L 1 2 ( ϕ ; T , τ ) exp R ( ϕ ; T , τ ) e 2 k τ v 2 ( t ) + S ( ϕ ; T , τ ) e k τ v ( t ) · E Q ^ T exp N ( ϕ ; T , τ ) | F t X ,
which means the core calculation of characteristic function is E Q ^ T exp N ( ϕ ; T , τ ) | F t X .
Note that N ( ϕ ; T ) is defined by (18) and i ( ϕ ; T ) , i = C , E , F satisfies the ordinary differential equations in (19), then we can rewrite N ( ϕ ; T , τ ) as
N ( ϕ ; T , τ ) = F ( ϕ ; T ) + σ ^ τ 2 2 L ( ϕ ; T , τ ) E 2 ( ϕ ; T ) + μ ^ 2 ( t , T ) 2 σ ^ τ 2 L ( ϕ ; T , τ ) E ( ϕ ; T ) + μ ^ 2 ( t , T ) 2 σ ^ τ 2 μ ^ 2 ( t , T ) L ( ϕ ; T , τ ) 1 .
In addition, if we try to derive the partial derivatives of (21), then
N ( ϕ ; T , τ ) t = E 2 ( ϕ ; T ) · a t + E ( ϕ ; T ) · b t + c t · d + d t · c
where
a : = σ ^ τ 2 2 L ( ϕ ; T , τ ) , b : = μ ^ 2 ( t , T ) 2 σ ^ τ 2 L ( ϕ ; T , τ ) ,
c : = μ ^ 2 ( t , T ) 2 σ ^ τ 2 , d : = μ ^ 2 ( t , T ) L ( ϕ ; T , τ ) 1 .
Therefore, the partial derivatives of a , b , c , d satisfy
a t = σ ^ 2 e 2 k τ 1 2 σ ^ τ 2 C ( ϕ ; T ) 2 σ ^ 2 C ( ϕ ; T ) e 2 k τ σ ^ τ 2 2 1 2 σ ^ τ 2 C ( ϕ ; T ) 2 , b t = μ ^ ( t , T ) k θ ˜ X ( t ) e k ( T t ) σ ^ τ 2 1 2 σ ^ τ 2 C ( ϕ ; T ) σ ^ τ 4 1 2 σ ^ τ 2 C ( ϕ ; T ) 2 μ ^ 2 ( t , T ) 2 σ ^ 2 e 2 k τ 1 2 σ ^ τ 2 C ( ϕ ; T ) + 4 σ ^ τ 2 C ( ϕ ; T ) e 2 k τ σ ^ τ 2 4 σ ^ τ 4 1 2 σ ^ τ 2 C ( ϕ ; T ) 2 , c t = 4 μ ^ ( t , T ) σ ^ τ 2 k θ ˜ X ( t ) e k ( T t ) + 2 σ ^ 2 e 2 k τ μ ^ 2 ( t , T ) 4 σ ^ τ 4 , d t = 2 μ ^ ( t , T ) k θ ˜ X ( t ) e k ( T t ) 1 2 σ ^ τ 2 C ( ϕ ; T ) 2 σ ^ 2 e 2 k τ μ ^ 2 ( t , T ) C ( ϕ ; T ) 1 2 σ ^ τ 2 C ( ϕ ; T ) 2 .
In other words, E Q ^ T exp N ( ϕ ; T , τ ) F t X can be rewritten as
E Q ^ T exp N ( ϕ ; T , τ ) | F t X = E Q ^ T exp 0 t N ( ϕ ; T , τ ) t d t | F t X = E Q ^ T exp 0 t E 2 ( ϕ ; T ) · a t + E ( ϕ ; T ) · b t + c t · d + d t · c d t | F t X .
According to the Proposition 3.2 in [20], if Q denotes the transposition of Q, then
E Q ^ T exp ϕ 1 t T + Δ h , X ( s ) u ( s ) d s | F t X = Φ ( t , T ; h ) X ( t ) , I e ,
where
Φ ( t , T ; h ) = exp t T + Δ Q + ϕ 1 u ( s ) d i a g [ h ] d s ,
and
I e = ( 1 , 1 , 1 ) R N .
Therefore, we can calculate (24) in theory if ϕ 1 = 1 and h , X ( t ) = N ( ϕ ; T , τ t . So far, the theoretical derivation of the characteristic function has been completed in this paper. In the next section, we will first give the analytical solution of the ordinary differential equations in (19) that ensure high calculation precision.

4.3. Derivation of ODEs

In this section, we will focus on deriving the analytical solutions of (19) in Proposition 1. These solutions can greatly reduce errors caused by directly simulating numerical solutions of differential equations and improve computational efficiency.
Proposition 2. 
For any T t T + Δ , if function i ( ϕ ; t ) , i = C , E , F follows the ordinary differential equation
d F d t + 1 2 σ 2 E 2 + k θ ˜ E + σ 2 C + r ˜ ϕ + λ ˜ E Q ^ T [ e ϕ Y 1 ] = 0 , F ( ϕ ; T + Δ ) = 0 , d E d t + 2 k θ ˜ C + ( 2 σ 2 C k + ρ σ ϕ ) E = 0 , E ( ϕ ; T + Δ ) = 0 , d C d t + 2 σ 2 C 2 + ( 2 ρ σ ϕ 2 k ) C + 1 2 ϕ ( ϕ 1 ) = 0 , C ( ϕ ; T + Δ ) = 0 ,
then we can let w i denote the times of regime switching in [ T , T + Δ ] , and obtain the following analytical solutions of i ( ϕ ; t ) , i = C , E and F ( ϕ ; T ) :
C ( ϕ ; t ) = d b 4 σ 2 · u ( t ) 1 g · u ( t ) 1 , E ( ϕ ; t ) = D ( ϕ ; t ) · E ˜ ( ϕ ; t , θ ˜ ) , F ( ϕ ; T ) = F 1 ( ϕ ; T ) + F 2 ( ϕ ; T ) + F 3 ( ϕ ; T )
where
D ( ϕ ; t ) = a g · u ( t ) g · u ( t ) 1 · exp d 2 t , E ˜ ( ϕ ; t , θ ˜ ) = h i + E 0 ( ϕ ; w i + 1 , θ ˜ w i ) , h i = j = i + 1 z E 0 ( ϕ ; w j + 1 , θ ˜ w j ) , A ( w ) = u 1 2 ( w ) + u 1 2 ( w ) , E 0 ( ϕ ; w , θ ˜ ) = 2 θ ˜ d · exp d 2 ( T + Δ ) A ( w ) A ( t ) , b = 2 ρ σ ϕ 2 k , d = b 2 4 σ 2 ϕ ( ϕ 1 ) , a = k ( b + d ) 2 σ 2 , g = b d b + d , x = 1 2 d ( T + Δ t ) , α = arctanh g 1 g + 1 , y = α + x u ( t ) = exp d ( T + Δ t ) , F 1 ( ϕ ; T ) = i = 1 z σ 2 h i 2 a 2 g 4 d · exp d ( T + Δ ) · coth ( y w i + 1 ) coth ( y w i ) + σ 2 h i 2 · a 2 g θ ˜ X ( w i ) d · exp d ( T + Δ ) · A ( w i + 1 ) d · coth ( y w i + 1 ) coth ( y w i ) + g + 1 d g · 1 sinh ( y w i + 1 ) 1 sinh ( y w i ) g + 1 d g · ln tanh ( w i w i 2 2 ) tanh ( w i + 1 w i + 1 2 2 ) + 8 · a g σ θ ˜ X ( w i ) d 2 · A 2 ( w i + 1 ) 2 d ( g + 1 ) · 1 ( g 1 ) + ( g + 1 ) tanh ( x w i + 1 ) 1 ( g 1 ) + ( g + 1 ) tanh ( x w i ) 1 8 g 2 d · d g 2 + 1 · w i + 1 w i + ( g + 1 ) 2 · 1 tanh ( y w i + 1 ) 1 tanh ( y w i ) + 2 ( g 2 1 ) · ln sinh ( y w i + 1 ) sinh ( y w i ) + 2 A ( w i + 1 ) 4 g 3 / 2 d · ( g 1 ) · ln tanh ( w i w i 2 2 ) tanh ( w i + 1 w i + 1 2 2 ) + ( g + 1 ) · 1 sinh ( y w i + 1 ) 1 sinh ( y w i ) , F 2 ( ϕ ; T ) = i = 1 z k θ ˜ X ( w i ) h i · a g d · exp 1 2 d ( T + Δ ) · ln g · u ( w i + 1 ) + 1 g · u ( w i + 1 ) 1 ln g · u ( w i ) + 1 g · u ( w i ) 1 + k θ ˜ X ( w i ) · 2 a θ ˜ X ( w i ) d 2 · A ( w i + 1 ) · g · ln g · u ( w i + 1 ) + 1 g · u ( w i + 1 ) 1 + g + 1 · ln g · u ( w i + 1 ) 1 + g d · w i + 1 g + 1 · ln g · u ( w i ) 1 g d · w i A ( w i + 1 ) · g · ln g · u ( w i ) + 1 g · u ( w i ) 1 , F 3 ( ϕ ; T ) = i = 1 z ϕ r ˜ X ( w i ) + λ ˜ X ( w i ) ( p η 1 η 1 ϕ + q η 2 η 2 + ϕ 1 ) · w i + 1 w i + σ 2 · G ( ϕ ; w i + 1 ) G ( ϕ ; w i ) , G ( ϕ ; t ) = d + b 4 σ 2 · t + 1 2 σ 2 · ln 1 1 g · u ( t ) .
Proof. 
The proof of this proposition is presented in Appendix B.    □
Remark 2. 
The basic idea of this proposition is to simplify the numerical solution process of the differential equation. Considering that the Markov-modulated parameters remain constant before the next regime-switching moment, this proposition creatively divides the integration interval of the ordinary differential equation into the regime-switching interval and non-regime-switching one. Therefore, one can solve the simple function definite integral within these small intervals.

4.4. Pricing Formulae

In this section, we will use the relationship between the characteristic function and the expectation to derive the pricing formulae. Recall that the fair strike price K i , i = d 1 , d 2 of variance swaps should satisfy V ( X , 0 ) = 0 . In other words, we can obtain that
K d 1 = E Q ^ T A F N i = 1 N S ( t i ) S ( t i 1 ) S ( t i 1 ) 2 × 100 2 | F 0 X F 0 S F 0 V = A F N i = 1 N E Q ^ T e y ( t i ) 1 2 | F 0 X F 0 S F 0 V × 100 2 = A F N i = 1 N E Q ^ T e 2 y ( t i ) 2 e y ( t i ) + 1 | F 0 X F 0 S F 0 V × 100 2 = A F N i = 1 N f ( 2 ; 0 , t i , Δ t , v ( 0 ) ) 2 f ( 1 ; 0 , t i , Δ t , v ( 0 ) ) + 1 × 100 2 .
Similarly, we can compute
K d 2 = E Q ^ T A F N i = 1 N ln 2 S ( t i ) S ( t i 1 ) × 100 2 | F 0 X F 0 S F 0 V = A F N i = 1 N E Q ^ ln 2 S ( t i ) S ( t i 1 ) | F 0 X F 0 S F 0 V × 100 2 = A F N i = 1 N E Q ^ T y 2 ( t i ) F 0 X F 0 S F 0 V × 100 2 = A F N i = 1 N f ( 2 ) ( 0 ; 0 , t i , Δ t , v ( 0 ) ) × 100 2 ,
where f ( 2 ) ( 0 ; 0 , t i , Δ t , v ( 0 ) ) = f ( 2 ) ( ϕ ; 0 , t i , Δ t , v ( 0 ) ) ϕ 2 .

5. Numerical Tests and Discussion

In this section, we will use semi-Monte Carlo simulation to verify the accuracy and efficiency of our pricing formulae. It provides direct numerical evidence of the comparison between the Monte Carlo simulation of the model and our formulae. We will also examine the rationality of the promotion term in the new model, as it will help readers understand the effects of these realistic economic factors. In addition, we will analyze the sensitivity of (28) and (29) by studying a 1 % change of model parameters. This analysis not only helps to assess the stability and reliability of the model, but also allows us to identify the parameters that exert the greatest influence on its results.
We assume a typical market with only two regimes, namely, the market state switches between regime 1 and regime 2. These market states can be interpreted as bull markets (a good economy) and bear markets (a bad economy), respectively. Thus, the state space in X ( t ) is X : = { e 1 , e 2 } , and the transition rate matrix is given by
Q = 0.1 0.1 0.4 0.4 .
Moreover, we will adopt the values of parameters in Table 1 that are used in [4,14].

5.1. Semi-Monte Carlo Stimulation

In this section, we will test the accuracy and efficiency of our pricing Formulae (28) and (29) by comparing the results generated by semi-Monte Carlo simulation. This comparison more directly and effectively demonstrates the rationality of the new model and the advantages of the pricing formulae.
We utilize Euler–Maruyama discretization for our model (8) to Monte Carlo simulation:
M ( t i + 1 ) = M ( t i ) + r X ( t i ) 1 2 v 2 ( t i ) d t + v ( t i ) d t W ˜ 1 ( t i ) + Y i , S ( t i + 1 ) = exp M ( t i + 1 ) , v ( t i + 1 ) = k θ ˜ X ( t i ) v ( t i ) d t + σ d t ρ W ˜ 1 ( t i ) + 1 ρ 2 W ( t i )
where W ˜ 1 ( t i ) and W ( t i ) are two independent standard normal random variables. If we divide the sampling interval [ 0 , T ] into n sub-intervals, then the time length of each part satisfies d t = T / n .
Compared to traditional Monte Carlo simulation, semi-Monte Carlo simulation is an effective method for improving the simulation efficiency of regime-switching models [4,20,38]. In traditional Monte Carlo simulation, randomness exists not only in the simulation of regime-switching paths but also in the simulation of asset price and volatility processes. This means that two identical regime paths can have different calculation results. In contrast, the randomness of semi-Monte Carlo simulation only exists in the regime-switching path. Therefore, if the paths are the same, the fair price of the simulation must be the same. Given that semi-Monte Carlo simulations require fewer computation paths than Monte Carlo simulations, this paper leverages independent and identically exponentially distributed random variables to effectively generate simulation paths and capture the necessary randomness. Specifically, we will adopt the following four steps to implement the semi-Monte Carlo simulation process:
Step 1: Generate N = 20,000 independent identically exponentially distributed (i.i.d) random variables with parameter β e 1 . These variables represent the time until the first regime-switching event occurs for each of the N sample paths. Based on the initial market regime e 1 and the probability of remaining in this regime q 11 = 0.1 , we can estimate the number of sample paths that do not experience a regime switch before the maturity time T = 1 as N 0 e 0.1 · T · N 0.90 · N . Note that the actual number of such paths may differ due to randomness in the generation of the i.i.d random variables. Therefore, we only need to run the following computation for ( N N 0 ) + 1 sample paths, where N N 0 is the number of paths that do experience a regime switch before maturity, and + 1 represents the one path that does not experience any regime switches.
Step 2: For every sample path i = 1 , 2 , , ( N N 0 ) + 1 , calculate the conditional characteristic function (17) according to Propositions 1 and 2. When there is no switch in the regime switching path, it is worth noting that the difference between the analytical solutions of the differential equation system (Proposition 2) only depends on the interval length and the initial regime at the starting point (keeping in mind that there may have been a switch before the path interval, and the initial regime may be different). Additionally, the difference in the initial regime is proportional. Therefore, before starting step 2, we assume that there are no regime switches in the whole process. Then, we calculate and save the results under different sampling frequencies. After entering the loop simulation, for any path, we divide it into intervals with and without switching based on the sampling frequency. In the subsequent calculations, we call and adjust these results when there is no switching in the initial regime.
Step 3: Utilize the conditional characteristic function in step 2 to calculate the fair strike price K d j i , i = 1 , 2 , , ( N N 0 ) + 1 , j = 1 , 2 .
Step 4: Calculate the mean of the fair strike price K d j i as the ultimate fair strike price K d j , j = 1 , 2 , that is
K d j = N 0 · K d j ( n o n r e g i m e s w i t c h i n g ) + N N 0 · K d j ( r e g i m e s w i t c h i n g ) N .
Moreover, we provide an algorithm (Algorithm 1) that calculates the fair strike price for a variance swap. Here, we only take K 1 as an example.
Figure 1 shows images of the fair strike price obtained from semi-Monte Carlo simulation of Equation (28) and those calculated by Monte Carlo simulation under different observation frequencies ranging from 4 to 252. Actually, the observation frequencies represent the number of trading days in a year, with 4 and 252 representing quarterly and daily observations, respectively.
Algorithm 1: Calculation of the strike price
Mathematics 11 02730 i001
It is evident that our solutions match the simulations well, which effectively verifies the accuracy of our pricing formulae. Furthermore, our pricing formula has higher computational efficiency, with a simulation time of 11.30401 s. In contrast, the Monte Carlo simulation time for the model requires 3.962049 min (using the simulation running time of K d 1 as an example).

5.2. Regime-Switching Effect

In this section, we will compare the fair strike price K d i , i = 1 , 2 computed from the regime-switching hybrid model and the non-regime-switching hybrid model, as shown in Figure 2. Here, the K d i , i = 1 , 2 without regime switching mainly refer to the situation where the model parameters are kept constant at the initial regime value throughout the entire contract period T. For example, if the initial regime is bearish, then the model parameters during the maturity T are all the values under regime 2.
Figure 2 shows that the price K d i , i = 1 , 2 can be significantly impacted by regime switching. More specifically, the regime tends to lower K d i , i = 1 , 2 to some level when it switches from a bull to a bear regime. On the other hand, when the regime switches from a bear to a bull regime, the price tends to rise somewhat. At the same time, we also observe that K d i , i = 1 , 2 are relatively higher when the initial regime is in a bull state (regime 1) compared to when the market regime is in a bear state (regime 2), which is consistent with the findings in [4] and real financial markets. Particularly, when the observation frequencies rises from quarterly to daily, we keep track of the price under regime 1 (see Table 2).
Furthermore, we can see from Table 2 that K d 1 is considerably higher than K d 2 , and regime switching can result in large variations in K d 2 . This is mostly due to the fact that (6) prioritizes describing the relative volatility of asset prices, whereas (5) focuses more on the intensity of asset price swings.

5.3. Jump Diffusion Effect

In this section, we will investigate the impact of jump diffusion by comparing the fair strike prices K d i , i = 1 , 2 with different jump intensities λ ˜ X ( t ) . Unless otherwise specified, we assume that there is no jump effect in the stock price process when λ ˜ X ( t ) = 0 . Since the major concern in this section is not the regime-switching effect, we choose to explore the jump diffusion effect under regime 1. Figure 3 shows K d i , i = 1 , 2 at different observation frequencies as the jump intensity λ ˜ X ( t ) increases from 0 to 0.4.
From Figure 3, we can observe that at each observation frequency, the jump in S ( t ) leads to higher prices, and as the jump intensity λ ˜ X ( t ) increases, the prices further increase. One possible explanation for this is that jump diffusion may lead to changes in the underlying asset prices, resulting in higher realized variances and higher prices.

5.4. Sensitivity Analysis

In this section, we will examine the sensitivity of our model parameters by computing the fair strike prices K d i , i = 1 , 2 on the 1 % change in it. As regime switching and observation frequency are not the main concern in this section, we consider the sensitivity of parameters under daily sampling and the initial state of regime 1. Table 3 shows the percentage change in K d i , i = 1 , 2 caused by changes in model parameters.
From Table 3, we can see that the three parameters—long-term volatility mean θ ˜ X ( t ) , initial volatility v ( 0 ) , and mean reversion speed of volatility k—have a significant impact on the volatility of fair execution prices K d i , i = 1 , 2 . Furthermore, changes in k have the most significant impact on the volatility of prices. This is consistent with the real financial market, where asset price volatility is primarily related to the stochastic volatility process of the asset price. Additionally, we find that the initial stock price S ( 0 ) has no impact on the volatility of fair execution price K d i . i = 1 , 2 .

5.5. Comparison of the New Model with the Traditional MRG Model

In this section, we will consider the fair strike price under different models and examine how the inclusion of more realistic factors influences price. In other words, we will compare the fair strike price under the new model illustrated in (8) with prices derived from the traditional MRG model.
Both the Heston model and the MRG model are classic models in the stochastic field, and it would be comprehensive to compare our new model with both of them. However, as pointed out by Zhang [14], the connection between the MRG model and the Heston model can be made clear under a specific assumption—that the long-term average value of volatility θ ˜ is equal to zero. Therefore, in this section, we will focus on comparing our new model with the traditional MRG model only. Recall that the traditional MRG models assume a continuous price process with stochastic volatility, and all the parameters in the models are constant. To maintain this setting, we suppose that the parameters in the traditional models are all set in regime 1, and non-regime switching occurs in maturity.
It should be noted that MRG model may not adequately capture the long-term trends in the financial market that are influenced by structural changes in macroeconomic conditions. However, empirical finance literature extensively documents the correlation between the financial market and the macroeconomic environment. Moreover, sudden and unexpected events in real-world financial markets can result in significant fluctuations in asset prices, which cannot be disregarded, especially when pricing volatility derivatives. Therefore, adding the assumption in our new model can make it more comprehensive and meaningful. Figure 4 shows the strike prices obtained from Zhang’s [14] formula and ours. (Here, we only take K d 1 as an example to show the difference.)
From Figure 4, we can see that the price under our new model is cheaper than the one in the traditional MRG model. This can be explained as follows. First, both regime switching and jump diffusion contribute to the variation of the strike price, as discussed in the previous subsection. Second, this can be interpreted by the assumption that the parameters in the traditional MRG model are all set in regime 1 (good economy), which is larger than the ones in regime 2 (bad economy). In other words, a bearish market can drag down the price to some extent. Moreover, the small value of the long-term volatility mean may lead to a lower jump probability, meaning that the regime-switching effect has a greater effect on the price than the jump-diffusion effect. The sensitivity analysis in Table 3 can provide some evidence.

6. Conclusions

In this paper, we investigate the fair strike prices for discretely sampled variance swaps under the Markov-modulated MRG model with asymmetric double exponential jump diffusion. Compared to recent work in the variance swaps pricing field, the major contribution of this paper is to propose an effective approach to obtain a general pricing formula and improve the computational efficiency of semi-Monte Carlo simulation. The pricing problem can be reduced to an expectation calculation under the risk-neutral T-forward measure. We obtain analytical solutions for the pricing problem using the characteristic function and the conditional transition characteristic function. Additionally, we carry out a series of semi-Monte Carlo simulations to verify the accuracy and effectiveness of the pricing formulae. Our experiments aim to explore the influence of regime-switching and jump-diffusion terms, as well as analyze the sensitivity of pricing formulae parameters to the fair strike price. We find that when regime switching occurs, either the fair strike price for good economics will decrease or the fair strike price for bad economics will increase. When there is a jump in market asset prices, the diffusion of the jump will accelerate changes in stock prices, resulting in higher fair execution prices. In addition, we also find that the stochastic fluctuation parameters of the propagation model are more sensitive to price fluctuations.
Furthermore, we would like to emphasize the practical applications of our model. The pricing method proposed in our proposition is generalizable and provides a new way to accurately compute Markov-modulated variables. This approach can also be extended to other regime-switching stochastic models. For example, it could simplify the pricing formulae presented in Elliott [25] and Yang [4] without the need for overly complex expressions. Our attempt to price variance swaps under this new model can provide valuable insights and numerical evidence for financial institutions and investors to better manage and control market volatility risk, ultimately reducing trading losses.
However, it is a pity that we do not deeply consider the effect of a long contract life where the interest rate can change stochastically. In other words, if the interest rate is assumed to follow a specific stochastic process instead of a deterministic one, the ultimate risk-neutral measure and pricing formulae may differ. It would be worthwhile to explore how to adjust the pricing approach in such cases for future studies. Moreover, it would be worthwhile to explore reasonable methods for calibrating model parameter values, which we lack in this paper and hope that future research can address.

Author Contributions

Conceptualization, C.W.; Methodology, A.Z. and J.W.; Software, J.W.; Writing—original draft, A.Z.; Supervision, C.W. All authors have read and agreed to the published version of the manuscript.

Funding

The project is supported by the Fundamental Research Funds for the Central Universities (23IJKY05) and the National Nature Science Foundation of China (11601306).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Proposition 1

In this appendix, a brief proof of Proposition 1 will be given. Because the effect of Markov chain X ( t ) can be temporarily ignored for given information F T + Δ X , we can utilize the tower law of expectation to simplify the conditional characteristic function as
f ϕ ; t , T , Δ t , v ( t ) F T + Δ X = E Q ^ T e ϕ y ( T ) F t S F t V F T + Δ X = E Q ^ T E Q ^ T e ϕ y ( T ) F T S F T V F T + Δ X F t S F t V F T + Δ X = E Q ^ T E Q ^ T e ϕ y ( T ) F T S F T V F t S F t V .
Therefore, we can solve the outer expectation of the inner expectation E Q ^ T e ϕ y ( T ) F T S F T V to compute (A1).
We begin by solving the inner expectation E Q ^ T e ϕ y ( T ) F T S F T V . To keep the notation as simple as possible, the subscript of i X ( t ) , i = θ ˜ , r , λ ˜ , μ ˜ will be omitted here. Let
U ϕ ; t , X , v ( t ) = E Q ^ T e ϕ y ( T ) F t S F t V , T t T + Δ ,
where X = ln S ( t ) ln S ( T ) . According to the Feynman–Kac Theorem, function (A2) follows the partial differential equation
U t + ( r ˜ 1 2 v 2 ) U X + k ( θ ˜ v ) U v + 1 2 v 2 2 U X 2 + 1 2 σ 2 2 U v 2 + ρ σ v 2 U v X + λ ˜ E Q ^ T [ U ( ϕ ; t , X + Y , v ) U ( ϕ ; t , X , v ) ] = 0 ,
with
U ( ϕ ; t = T + Δ , X , v ) = e ϕ X .
As the statement in [14], we can assume that the solution of this partial Equation (A3) is
U ( ϕ ; t , X , v ) = e F ( ϕ ; t ) + E ( ϕ ; t ) v + C ( ϕ ; t ) v 2 + ϕ X .
Substituting (A5) into Equations (A3) and (A4), we can obtain the following ODEs:
d F d t + 1 2 σ 2 E 2 + k θ ˜ E + σ 2 C + r ˜ ϕ + λ ˜ E Q ^ T [ e ϕ Y 1 ] = 0 , F ( ϕ ; T + Δ ) = 0 , d E d t + 2 k θ ˜ C + ( 2 σ 2 C k + ρ σ ϕ ) E = 0 , E ( ϕ ; T + Δ ) = 0 , d C d t + 2 σ 2 C 2 + ( 2 ρ σ ϕ 2 k ) C + 1 2 ϕ ( ϕ 1 ) = 0 , C ( ϕ ; T + Δ ) = 0 .
Note that X = ln S ( T ) ln S ( T ) = 0 at time t = T , then the inner expectation follows:
E Q ^ T e ϕ y ( T ) F T S F T V = e F ( ϕ ; T ) + E ( ϕ ; T ) v ( T ) + C ( ϕ ; T ) v 2 ( T ) .
We next aim to compute the outer expectation. Substituting the expression of (A6) into (A1), the outer expectation follows:
f ϕ ; t , T , Δ t , v ( t ) F T + Δ X = E Q ^ T e F ( ϕ ; T ) + E ( ϕ ; T ) v ( T ) + C ( ϕ ; T ) v 2 ( T ) F t S F t V F T + Δ X .
Adopting the solution procedure used by [11,16], we can solve it by utilizing
f ϕ ; t , T , Δ , v ( t ) F T + Δ X = f v ( T ) p v ( T ) | v ( t ) d v ( T ) ,
where f ( v ) = e F ( ϕ ; T ) + E ( ϕ ; T ) v + C ( ϕ ; T ) v 2 and p v ( T ) | v ( t ) presents the transition density from state ( v ( t ) , t ) to state ( v ( T ) , T ) . As Lemma 1 holds, the volatility process v ( T ) at time T can be
v ( T ) = e k T v ( 0 ) + k 0 T θ ˜ X ( u ) e k ( T u ) d u + σ 0 T e k ( T u ) d W ˜ 2 ( u ) = e k ( T t ) v ( t ) + k t T θ ˜ X ( u ) e k ( T u ) d u + σ t T e k ( T u ) d W ˜ 2 ( u ) ,
which implies that conditional expectaion μ ^ τ is given by
μ ^ τ : = E Q ^ T [ v ( T ) | v ( t ) ] = e k τ v ( t ) + k t T θ ˜ X ( u ) e k ( T u ) d u : = e k τ v ( t ) + μ ^ ( t , T )
and the conditional variance σ ^ τ 2 follows
σ ^ τ 2 : = V a r [ v ( T ) | v ( t ) ] = σ 2 2 k ( 1 e 2 k τ ) .
Therefore, we can argue that the transition density p v ( T ) | v ( t ) must be
p v ( T ) | v ( t ) = 1 2 π σ ^ τ e 1 2 σ ^ τ 2 v ( T ) μ ^ τ 2 , τ = T t .
Substituting (A9) into (A8), we then obtain that
f ( ϕ ; t , T , Δ , v ( t ) F T + Δ X ) = 1 2 π σ ^ τ exp 1 2 σ ^ τ 2 v ( T ) μ ^ τ 2 + F ( ϕ ; T ) + E ( ϕ ; T ) v ( T ) + C ( ϕ ; T ) v 2 ( T ) d v ( T ) = 1 2 π σ ^ τ exp C ( ϕ ; T ) 1 2 σ ^ τ 2 v 2 ( T ) + E ( ϕ ; T ) + μ ^ τ σ ^ τ 2 v ( T ) + F ( ϕ ; T ) μ ^ τ 2 2 σ ^ τ 2 d v ( T ) = 1 2 π σ ^ τ exp 1 2 σ ^ τ 2 1 2 σ ^ τ 2 C ( ϕ ; T ) v 2 ( T ) 2 σ ^ τ 2 E ( ϕ ; T ) + μ ^ τ v ( T ) + μ ^ τ 2 2 σ ^ τ 2 F ( ϕ ; T ) d v ( T ) = 1 L ( ϕ ; T , τ ) L ( ϕ ; T , τ ) 2 π σ ^ τ exp L ( ϕ ; T , τ ) 2 σ ^ τ 2 v ( T ) σ ^ τ 2 E ( ϕ ; T ) + μ ^ τ L ( ϕ ; T , τ ) 2 + μ ^ τ 2 2 σ ^ τ 2 F ( ϕ ; T ) L ( ϕ ; T , τ ) σ ^ τ 2 E ( ϕ ; T ) + μ ^ τ 2 L 2 ( ϕ ; T , τ ) d v ( T ) .
To work out this integral (A10), we recall that the probability density function of the standard normal distribution has the characteristics
e a ( x b ) 2 d x = 2 π + 1 2 π e 1 2 2 a ( x b ) 2 d x = π a + 1 2 π e 1 2 y 2 d y = π a , a > 0 ,
and for any a 0 ,
e a ( x b ) 2 d x = e a y 2 d y = 2 0 e a y 2 d y 2 0 1 d y = .
Thus we use this fact to argue that, for 1 2 σ ^ τ 2 > C ( ϕ ; T ) , the conditional function (A10) can be derived below:
f ( ϕ ; t , T , Δ , v ( t ) F T + Δ X ) = L 1 2 ( ϕ ; T , τ ) exp L ( ϕ ; T , τ ) 2 σ ^ τ 2 μ ^ τ 2 L ( ϕ ; T , τ ) 2 σ ^ τ 2 F ( ϕ ; T ) L ( ϕ ; T , τ ) σ ^ τ 2 E ( ϕ ; T ) + μ ^ τ 2 L 2 ( ϕ ; T , τ ) = L 1 2 ( ϕ ; T , τ ) exp F ( ϕ ; T ) + σ ^ τ 2 E ( ϕ ; T ) + μ ^ τ 2 L ( ϕ ; T , τ ) μ ^ τ 2 2 σ ^ τ 2 L ( ϕ ; T , τ ) = L 1 2 ( ϕ ; T , τ ) exp F ( ϕ ; T ) + σ ^ τ 2 E ( ϕ ; T ) + e k τ v ( t ) + μ ^ ( t , T ) 2 2 σ ^ τ 2 L ( ϕ ; T , τ ) e k τ v ( t ) + μ ^ ( t , T ) 2 2 σ ^ τ 2 = L 1 2 ( ϕ ; T , τ ) exp N ( ϕ ; T , τ ) + S ( ϕ ; T , τ ) e k τ v ( t ) + R ( ϕ ; T , τ ) e 2 k τ v 2 ( t ) .
where the parameters are defined in (18), and for 1 2 σ ^ τ 2 C ( ϕ ; T ) , the conditional function (A10) tends to infinity.

Appendix B. Proof of Proposition 2

In this appendix, a brief proof of Proposition 2 will be given. As shown in Appendix A of [4], the differential equation
d C d t + 2 σ 2 C 2 + ( 2 ρ σ ϕ 2 k ) C + 1 2 ϕ ( ϕ 1 ) = 0 , C ( ϕ ; T + Δ ) = 0 .
There exists an analytical solution
C ( ϕ ; t ) = d b 4 σ 2 · u ( t ) 1 g · u ( t ) 1 ,
and all the parameters here are given by (27). Let
d G ( ϕ ; t ) d t = C ( ϕ ; t ) , a = k ( b + d ) 2 σ 2 .
We then rewrite d E d t + 2 k θ ˜ C + ( 2 σ 2 C k + ρ σ ϕ ) E = 0 as
d exp 2 σ 2 G ( ϕ ; t ) + b 2 t E ( ϕ ; t ) d t = 2 k θ ˜ C ( ϕ ; t ) exp 2 σ 2 G ( ϕ ; t ) + b 2 t .
To solve this, we first compute
G ( ϕ ; t ) = d b 4 σ 2 · 1 g 0 t 1 + 1 g g · u ( s ) 1 d t = d + b 4 σ 2 · t + ( 1 g ) · 0 t 1 g · u ( s ) 1 d t = d + b 4 σ 2 · t d + b 4 σ 2 · 1 g d · 0 t g d u ( s ) g · u ( s ) g · u ( s ) 1 = d + b 4 σ 2 · t + 1 2 σ 2 · 0 t 1 g · u ( s ) g · u ( s ) 1 d ( g · u ( s ) ) = d + b 4 σ 2 · t + 1 2 σ 2 · ln 1 1 g · u ( t ) .
Then substitute (A12) into (A11) to obtain
d exp 2 σ 2 G ( ϕ ; t ) + b 2 t E ( ϕ ; t ) d t = 2 k θ ˜ · d b 4 σ 2 · u ( t ) 1 g · u ( t ) 1 · exp d 2 t · 1 1 g · u ( t ) = 2 k θ ˜ · d b 4 σ 2 · exp d 2 t · u ( t ) 1 g · u ( t ) = a θ ˜ · exp d 2 t · 1 exp d ( T + Δ t ) .
Solving (A13), we have
E ( ϕ ; t ) = a g · u ( t ) g · u ( t ) 1 · t T + Δ θ ˜ exp d 2 ( s t ) · 1 u ( s ) 1 d t .
Note that there are some Markov-modulated parameters θ ˜ , r , λ ˜ , μ ˜ in functions E ( ϕ ; t ) and F ( ϕ ; t ) that can lead to “bad” expressions. We first consider dividing the sampling time interval [ T , T + Δ ] in more detail. If we denote the time when regime switching happens as w 1 , w 2 , , w z , then
[ T , T + Δ ] = i = 0 z [ w i , w i + 1 ) , i = T , 1 , 2 , , z , T + Δ ,
which yields that the Markov-modulated parameters are fixed constants in state e i , i = 1 , 2 , , N . Therefore, there exists i such that t [ w i , w i + 1 ) ,
E ˜ ( ϕ ; t , θ ˜ ) : = t T + Δ θ ˜ X ( s ) exp d 2 s · exp d ( T + Δ s ) 1 d s = E ˜ ( ϕ ; w i + 1 ) + θ ˜ X ( w i ) t w i + 1 exp d 2 s · exp d ( T + Δ s ) 1 d s = E ˜ ( ϕ ; w i + 1 ) + θ ˜ X ( w i ) t w i + 1 exp d ( T + Δ ) + d 2 s exp d 2 s d s = E ˜ ( ϕ ; w i + 1 , θ ˜ ) + 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) A ( w ) A ( t ) .
In particular, for i = z , E ˜ ( ϕ ; w z + 1 ) = E ( ϕ ; Δ ) = 0 . In other words, (A15) follows that
E ˜ ( ϕ ; t , θ ˜ ) = E ˜ ( ϕ ; w i , θ ˜ ) + E 0 ( ϕ ; w i + 1 , θ ˜ X ( w i ) ) = j = i + 1 z E 0 ( ϕ ; w j + 1 , θ ˜ w j ) + E 0 ( ϕ ; w i + 1 , θ ˜ X ( w i ) ) = h i + E 0 ( ϕ ; w i + 1 , θ ˜ w i ) .
Then the solution (A14) of function E ( ϕ ; t ) can ultimately simplify as
E ( ϕ ; t ) : = D ( ϕ ; t ) · E ˜ ( ϕ ; t , θ ˜ ) ,
where
h i = j = i + 1 z E 0 ( ϕ ; w j + 1 , θ ˜ w j ) , A ( w ) = u 1 2 ( w ) + u 1 2 ( w ) , D ( ϕ ; t ) = a g · u ( t ) g · u ( t ) 1 · exp d 2 t , E 0 ( ϕ ; w , θ ˜ ) = 2 θ ˜ d · exp d 2 ( T + Δ ) A ( w ) A ( t ) .
Now we are in a position to solve
d F d t + 1 2 σ 2 E 2 + k θ ˜ E + σ 2 C + r ˜ ϕ + λ ˜ E Q ^ T [ e ϕ Y 1 ] = 0 , F ( ϕ ; T + Δ ) = 0 .
According to (18), the computation of conditional characteristic function only requires F ( ϕ ; T ) . Therefore, we will focus on the calculation of F ( ϕ ; T ) by following three parts:
F ( ϕ ; T ) = 1 2 σ 2 T T + Δ E 2 ( ϕ ; t ) d t + k T T + Δ θ ˜ E ( ϕ ; t ) d t + σ 2 T T + Δ C ( ϕ ; t ) d t + ϕ T T + Δ r ˜ d t + E Q ^ T [ e ϕ Y 1 ] T T + Δ λ ˜ d t = i = 1 z 1 2 σ 2 w i w i + 1 E 2 ( ϕ ; t ) d t + i = 1 z k w i w i + 1 θ ˜ X ( t ) E ( ϕ ; t ) d t + i = 1 z σ 2 w i w i + 1 C ( ϕ ; t ) d t + ϕ w i w i + 1 r ˜ d t + E Q ^ T [ e ϕ Y 1 ] w i w i + 1 λ ˜ d t : = F 1 ( ϕ ; T ) + F 2 ( ϕ ; T ) + F 3 ( ϕ ; T ) ,
where
F 1 ( ϕ ; T ) = i = 1 z 1 2 σ 2 w i w i + 1 E 2 ( ϕ ; t ) d t , F 2 ( ϕ ; T ) = i = 1 z k w i w i + 1 θ ˜ X ( t ) E ( ϕ ; t ) d t , F 3 ( ϕ ; T ) = i = 1 z σ 2 w i w i + 1 C ( ϕ ; t ) d t + ϕ w i w i + 1 r ˜ d t + E Q ^ T [ e ϕ Y 1 ] w i w i + 1 λ ˜ d t .
We begin by deriving the formula for F 2 ( ϕ ; T ) . In fact,
k w i w i + 1 θ ˜ X ( t ) E ( ϕ ; t ) d t = k θ ˜ X ( w i ) w i w i + 1 D ( ϕ ; t ) · E ˜ ( ϕ ; t , θ ˜ ) d t = k θ ˜ X ( w i ) w i w i + 1 D ( ϕ ; t ) · h i + E 0 ( ϕ ; w i + 1 , θ ˜ ) d t = k θ ˜ X ( w i ) h i w i w i + 1 D ( ϕ ; t ) d t + k θ ˜ X ( w i ) w i w i + 1 D ( ϕ ; t ) · E 0 ( ϕ ; w i + 1 , θ ˜ ) d t ,
which shows that we need to derive the two integrals of (A19). Note that
w i w i + 1 D ( ϕ ; t ) · E 0 ( ϕ ; w i + 1 , θ ˜ ) d t = w i w i + 1 a g · u ( t ) g · u ( t ) 1 · exp d 2 t · 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) A ( w i + 1 ) A ( t ) d t = 2 a g θ ˜ X ( w i ) d · w i w i + 1 u ( t ) g · u ( t ) 1 · u 1 2 ( t ) A ( w i + 1 ) A ( t ) d t = 2 a g θ ˜ X ( w i ) d · A ( w i + 1 ) w i w i + 1 u 1 2 ( t ) g · u ( t ) 1 d t w i w i + 1 u 1 2 ( t ) g · u ( t ) 1 · A ( t ) d t = 2 a g θ ˜ X ( w i ) d · A ( w i + 1 ) · 2 d g · w i w i + 1 1 g · u ( t ) 2 1 · d g · u ( t ) g + 1 g · w i w i + 1 1 g · u ( t ) 1 d t 1 g · w i w i + 1 d t = 2 a θ ˜ X ( w i ) d 2 · A ( w i + 1 ) · g · ln g · u ( w i + 1 ) + 1 g · u ( w i + 1 ) 1 + g + 1 · ln g · u ( w i + 1 ) 1 + g d · w i + 1 A ( w i + 1 ) · g · ln g · u ( w i ) + 1 g · u ( w i ) 1 g + 1 · ln g · u ( w i ) 1 g d · w i .
and
w i w i + 1 D ( ϕ ; t ) d t = a · w i w i + 1 g · u ( t ) g · u ( t ) 1 · exp d 2 t d t = a g · exp 1 2 d ( T + Δ ) · w i w i + 1 u ( t ) g · u ( t ) 1 · u 1 2 ( t ) d t = a g · exp 1 2 d ( T + Δ ) · 2 d g · w i w i + 1 1 g · u ( t ) 2 1 · d g · u ( t ) = a g d · exp 1 2 d ( T + Δ ) · ln g · u ( w i + 1 ) + 1 g · u ( w i + 1 ) 1 ln g · u ( w i ) + 1 g · u ( w i ) 1 .
Then we can accomplish the computation of F 2 ( ϕ ; T ) by substituting (A21) and (A20) into (A19).
We next aim to derive the formula for F 3 ( ϕ ; T ) . It is easy to see that
σ 2 w i w i + 1 C ( ϕ ; t ) d t + ϕ w i w i + 1 r ˜ d t + E Q ^ T [ e ϕ Y 1 ] w i w i + 1 λ ˜ d t = σ 2 · G ( ϕ ; w i + 1 ) G ( ϕ ; w i ) + ϕ r ˜ X ( w i ) + λ ˜ X ( w i ) E Q ^ T [ e ϕ Y 1 ] · w i + 1 w i
where G ( ϕ ; t ) follows (A12). Since Y follows an asymmetric double exponential distribution (3), then the expectation satisfies
E Q ^ T e ϕ Y 1 F t X = + e ϕ y · p η 1 e η 1 y I y 0 + q η 2 e η 2 y I y < 0 d y 1 = p η 1 · 0 + e ( ϕ η 1 ) y d y + q η 2 · 0 e ( ϕ + η 2 ) y d y 1 = p η 1 η 1 ϕ + q η 2 η 2 + ϕ 1 .
Therefore, we can accomplish the computation of F 3 ( ϕ ; T ) by substituting (A23) into (A22).
We finally try to utilize hyperbolic function to derive the formula for F 1 ( ϕ ; T ) . There are three integrals in F 1 ( ϕ ; T ) :
1 2 σ 2 w i w i + 1 E 2 ( ϕ ; t ) d t = w i w i + 1 D 2 ( ϕ ; t ) · E 2 ( ϕ ; w i + 1 , θ ˜ ) d t = 1 2 σ 2 w i w i + 1 D 2 ( ϕ ; t ) · h i + E 0 ( ϕ ; w i + 1 , θ ˜ ) 2 d t = 1 2 σ 2 h i 2 · w i w i + 1 D 2 ( ϕ ; t ) d t + σ 2 h i · w i w i + 1 D 2 ( ϕ ; t ) · E 0 ( ϕ ; w i + 1 , θ ˜ ) d t + 1 2 σ 2 · w i w i + 1 D 2 ( ϕ ; t ) · E 0 2 ( ϕ ; w i + 1 , θ ˜ ) d t .
Let
x : = 1 2 d ( T + Δ t ) , α : = arctanh g 1 g + 1 , y : = α + x .
Then we can obtain that cosh ( α ) = g + 1 4 g and sinh ( α ) = g 1 4 g . In this sense, we can naturally derive
D ( ϕ ; t ) · E 0 ( ϕ ; w i + 1 , θ ˜ ) = a g · u ( t ) g · u ( t ) 1 · exp d 2 t · 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) · A ( w i + 1 ) A ( t ) = 2 a g θ ˜ X ( w i ) d · u ( t ) g · u ( t ) 1 · u 1 2 ( t ) · A ( w i + 1 ) 2 cosh ( x ) = 2 a g θ ˜ X ( w i ) d · A ( w i + 1 ) 2 cosh ( x ) g · u 1 2 ( t ) u 1 2 ( t ) = 2 a g θ ˜ X ( w i ) d · A ( w i + 1 ) 2 cosh ( x ) ( g 1 ) cosh ( x ) + ( g + 1 ) sinh ( x ) = 4 a g θ ˜ X ( w i ) d · A ( w i + 1 ) A ( w i + 1 ) 2 2 sech ( x ) 1 ( g 1 ) + ( g + 1 ) tanh ( x ) ,
and in the time interval [ w i , w i + 1 ) ,
w i + 1 w i D 2 ( ϕ ; t ) · E 0 2 ( ϕ ; w i + 1 , θ ˜ ) d t = 4 a g θ ˜ X ( w i ) d 2 · { w i + 1 w i A ( w i + 1 ) A ( w i + 1 ) 2 2 2 sech 2 ( x ) ( g 1 ) + ( g + 1 ) tanh ( x ) 2 d t w i + 1 w i 1 A ( w i + 1 ) sech ( x ) ( g 1 ) + ( g + 1 ) tanh ( x ) 2 d t } .
Further, if we denote x w j : = 1 2 d ( T + Δ w j ) , j = i , i + 1 , then
w i + 1 w i A ( w i + 1 ) A ( w i + 1 ) 2 2 2 sech 2 ( x ) ( g 1 ) + ( g + 1 ) tanh ( x ) 2 d t = 2 d · w i + 1 w i A ( w i + 1 ) A ( w i + 1 ) 2 2 2 ( g 1 ) + ( g + 1 ) tanh ( x ) 2 d tanh ( x ) = A 2 ( w i + 1 ) 2 d · 1 ( g + 1 ) ( g 1 ) + ( g + 1 ) tanh ( x w i + 1 ) 1 ( g + 1 ) ( g 1 ) + ( g + 1 ) tanh ( x w i ) ,
and
w i + 1 w i 1 A ( w i + 1 ) sech ( x ) ( g 1 ) + ( g + 1 ) tanh ( x ) 2 d t = w i + 1 w i cosh 2 ( x ) A ( w i + 1 ) cosh ( x ) ( g 1 ) cosh ( x ) + ( g + 1 ) sinh ( x ) 2 d t = 1 4 g · w i + 1 w i cosh 2 ( x ) A ( w i + 1 ) cosh ( x ) sinh ( α ) cosh ( x ) + cosh ( α ) sinh ( x ) d t = 1 4 g · w i + 1 w i cosh 2 ( y α ) A ( w i + 1 ) cosh ( y α ) sinh 2 ( y ) d t = u 1 + u 2 .
where
u 1 = 1 4 g 2 · w i + 1 w i ( g + 1 ) 2 cosh 2 ( y ) 2 ( g 2 1 ) cosh ( y ) sinh ( y ) + ( g 1 ) 2 sinh 2 ( y ) sinh 2 ( y ) d t = g + 1 4 g 2 · w i + 1 w i 1 + csch 2 ( y ) d t 2 ( g 2 1 ) 4 g 2 · w i + 1 w i cosh ( y ) sinh ( y ) d t + g 1 4 g 2 · w i + 1 w i d t = 1 8 g 2 d · d g 2 + 1 · w i + 1 w i + ( g + 1 ) 2 · 1 tanh ( y w i + 1 ) 1 tanh ( y w i ) + 2 ( g 2 1 ) · ln sinh ( y w i + 1 ) sinh ( y w i ) , u 2 = 1 4 g · A ( w i + 1 ) 4 g · w i + 1 w i csch ( y ) d t 1 4 g · A ( w i + 1 ) 4 g · w i + 1 w i cosh ( y ) sinh 2 ( y ) d t = 2 A ( w i + 1 ) d 4 g 3 / 2 · ( g 1 ) · ln tanh y w i + 1 y w i + 1 2 2 tanh y w i y w i 2 2 + ( g + 1 ) · 1 sinh ( y w i + 1 ) 1 sinh ( y w i ) .
Therefore, we can substitute (A26) and (A27) into (A25) and obtain that
w i + 1 w i D 2 ( ϕ ; t ) · E 0 2 ( ϕ ; w i + 1 , θ ˜ ) d t = 4 a g θ ˜ X ( w i ) d 2 · { A 2 ( w i + 1 ) 2 d ( g + 1 ) · 1 ( g 1 ) + ( g + 1 ) tanh ( x ) 1 8 g 2 d · d g 2 + 1 · t + ( g + 1 ) 2 tanh ( y ) + 2 ( g 2 1 ) · ln [ sinh ( y ) ] + 2 A ( w i + 1 ) 4 g 3 / 2 d · ( g 1 ) · ln | tanh y 2 | + ( g + 1 ) · 1 sinh ( y ) } | w i w i + 1 .
Similarly, we can simplify
D 2 ( ϕ ; t ) · E 0 ( ϕ ; w i + 1 , θ ˜ ) = a 2 g 2 · u 2 ( t ) g · u ( t ) 1 2 · exp d t · 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) · A ( w i + 1 ) 2 cosh ( x ) = 2 a 2 g 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) · u 2 ( t ) g · u ( t ) 1 2 · u 1 ( t ) · A ( w i + 1 ) 2 cosh ( x ) = 2 a 2 g 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) · A ( w i + 1 ) 2 cosh ( x ) ( g 1 ) cosh ( x ) + ( g + 1 ) sinh ( x ) 2 = 2 a 2 g 2 θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) · 1 4 g · A ( w i + 1 ) 2 cosh ( x ) sinh ( α ) cosh ( x ) + cosh ( α ) sinh ( x ) 2 = a 2 g θ ˜ X ( w i ) 2 d · exp d 2 ( T + Δ ) · A ( w i + 1 ) 2 cosh ( y α ) sinh 2 ( y ) .
Note that
w i + 1 w i A ( w i + 1 ) sinh 2 ( y ) d t = 2 A ( w i + 1 ) d · coth ( y w i + 1 ) coth ( y w i ) ,
and
w i + 1 w i cosh ( y α ) sinh 2 ( y ) d t = w i + 1 w i cosh ( y ) cosh ( α ) sinh ( α ) sinh ( y ) sinh 2 ( y ) d t = g + 1 d g · w i + 1 w i 1 sinh 2 ( y ) d sinh ( y ) + g 1 d g · w i + 1 w i 1 sinh ( y ) d y = g + 1 d g · 1 sinh ( y w i + 1 ) 1 sinh ( y w i ) + g 1 d g · ln tanh y w i + 1 y w i + 1 2 2 tanh y w i y w i 2 2 .
then the integral of (A29) in time interval [ w i , w i + 1 ) follows
w i + 1 w i D 2 ( ϕ ; t ) · E 0 ( ϕ ; w i + 1 , θ ˜ ) d t = a 2 g θ ˜ X ( w i ) d · exp d 2 ( T + Δ ) · A ( w i + 1 ) d · coth ( y ) + g + 1 d g · 1 sinh ( y ) g 1 d g · ln | tanh y 2 | | w i w i + 1 .
Up to now, the only integral in (A24) that we do not compute is w i + 1 w i D 2 ( ϕ ; t ) d t . Since
D 2 ( ϕ ; t ) = a 2 g 2 · u 2 ( t ) g · u ( t ) 1 2 · exp d t = a 2 g 2 · exp d ( T + Δ ) · u 2 ( t ) g · u ( t ) 1 2 · u 1 ( t ) = a 2 g 2 · exp d ( T + Δ ) · 1 ( g 1 ) cosh ( x ) + ( g + 1 ) sinh ( x ) 2 = a 2 g 4 · exp d ( T + Δ ) · 1 sinh ( α ) cosh ( x ) + cosh ( α ) sinh ( x ) 2 = a 2 g 4 · exp d ( T + Δ ) · 1 sinh 2 ( y ) ,
then the integral of (A31) in time interval [ w i , w i + 1 ) follows
w i + 1 w i D 2 ( ϕ ; t ) d t = a 2 g 4 · exp d ( T + Δ ) · w i + 1 w i 1 sinh ( y ) 2 d t = a 2 g 2 d · exp d ( T + Δ ) · coth ( y w i + 1 ) coth ( y w i ) .
Now we can ultimately derive F 1 ( ϕ ; T ) by substituting (A28) and (A30) into (A24).

References

  1. Demeterfi, K.; Derman, E.; Kamal, M.; Zou, J. More than you ever wanted to know about volatility swaps. Goldman Sachs Quant. Strateg. Res. Notes 1999, 41, 1–56. [Google Scholar]
  2. Zhu, S.-P.; Lian, G.-H. Analytically pricing volatility swaps under stochastic volatility. J. Comput. Appl. Math. 2015, 288, 332–340. [Google Scholar] [CrossRef]
  3. Carr, P.; Lee, R. Volatility derivatives. Annu. Rev. Financ. Econ. 2009, 1, 319–339. [Google Scholar] [CrossRef] [Green Version]
  4. Yang, Y.; Liu, S.; Yonghong Wu, Y.; Wiwatanapataphee, B. Pricing of volatility derivatives in a heston–cir model with markov-modulated jump diffusion. J. Comput. Appl. Math. 2021, 393, 113–277. [Google Scholar] [CrossRef]
  5. Ait-Sahalia, Y.; Karaman, M.; Mancini, L. The Term Structure of Variance Swaps and Risk Premia; Swiss Finance Institute Research Paper; Swiss Finance Institute: Zürich, Switzerland, 2018; pp. 18–37. [Google Scholar]
  6. Heston, S.L.; Nandi, S. Derivatives on Volatility: Some Simple Solutions Based on Observables; Federal Reserve Bank of Atlanta WP: Atlanta, GA, USA, 2000. [Google Scholar]
  7. Sepp, A. Pricing options on realized variance in the heston model with jumps in returns and volatility. J. Comput. Financ. 2008, 11, 33–70. [Google Scholar] [CrossRef] [Green Version]
  8. Carr, P.; Lee, R. Hedging variance options on continuous semimartingales. Financ. Stochastics 2010, 14, 179–207. [Google Scholar] [CrossRef]
  9. Swishchuk, A. Modeling of variance and volatility swaps for financial markets with stochastic volatilities. WILMOTT Mag. 2004, 2, 64–72. [Google Scholar]
  10. Broadie, M.; Jain, A. Pricing and hedging volatility derivatives. J. Deriv. 2008, 15, 7–24. [Google Scholar] [CrossRef] [Green Version]
  11. Zhu, S.-P.; Lian, G.-H. A closed-form exact solution for pricing variance swaps with stochastic volatility. Math. Financ. Int. J. Math. Stat. Financ. Econ. 2011, 21, 233–256. [Google Scholar]
  12. Broadie, M.; Jain, A. The effect of jumps and discrete sampling on volatility and variance swaps. Int. J. Theor. Appl. Financ. 2008, 11, 761–797. [Google Scholar] [CrossRef]
  13. Rujivan, S.; Zhu, S.-P. A simplified analytical approach for pricing discretely-sampled variance swaps with stochastic volatility. Appl. Math. Lett. 2012, 25, 1644–1650. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, L.-W. A closed-form pricing formula for variance swaps with mean-reverting gaussian volatility. ANZIAM J. 2014, 55, 362–382. [Google Scholar]
  15. Lian, G.; Chiarella, C.; Kalev, P.S. Volatility swaps and volatility options on discretely sampled realized variance. J. Econ. Dyn. Control 2014, 47, 239–262. [Google Scholar] [CrossRef]
  16. Cao, J.-P.; Fang, Y.-B. An analytical approach for variance swaps with an ornstein-uhlenbeck process. ANZIAM J. 2017, 59, 83–102. [Google Scholar]
  17. Mao, C.; Liu, G.; Wang, Y. A closed-form pricing formula for log-return variance swaps under stochastic volatility and stochastic interest rate. Mathematics 2021, 10, 5. [Google Scholar] [CrossRef]
  18. Zhu, S.P.; Lian, G.H. Pricing forward-start variance swaps with stochastic volatility. Appl. Math. Comput. 2015, 250, 920–933. [Google Scholar] [CrossRef]
  19. Maheu, J.M.; McCurdy, T.H. Identifying bull and bear markets in stock returns. J. Bus. Econ. Stat. 2000, 18, 100–112. [Google Scholar]
  20. Elliott, R.J.; Lian, G.-H. Pricing variance and volatility swaps in a stochastic volatility model with regime switching: Discrete observations case. Quant. Financ. 2013, 13, 687–698. [Google Scholar] [CrossRef]
  21. Bo, W.; Ruili, S. The memms for markov-modulated gbms. Chin. J. Appl. Probab. Stat. 2013, 29, 179–187. [Google Scholar]
  22. Elliott, R.J.; Chan, L.; Siu, T.K. Option pricing and esscher transform under regime switching. Ann. Financ. 2005, 1, 423. [Google Scholar] [CrossRef]
  23. Liao, J.; Shu, H.; Wei, C. Pricing power options with a generalized jump diffusion. Commun.-Stat.-Theory Methods 2017, 46, 11026–11046. [Google Scholar] [CrossRef]
  24. Bo, L.; Wang, Y.; Yang, X. Markov-modulated jump–diffusions for currency option pricing. Insur. Math. Econ. 2010, 46, 461–469. [Google Scholar] [CrossRef]
  25. Elliott, R.J.; Siu, T.K.; Chan, L. Pricing volatility swaps under heston’s stochastic volatility model with regime switching. Appl. Math. Financ. 2007, 14, 41–62. [Google Scholar] [CrossRef]
  26. Lin, S.; He, X.-J. Analytically pricing variance and volatility swaps with stochastic volatility, stochastic equilibrium level and regime switching. Expert Syst. Appl. 2023, 217, 119592. [Google Scholar] [CrossRef]
  27. Hamilton, J.D. Regime switching models. In Macroeconometrics and Time Series Analysis; Palgrave Macmillan: London, UK, 2010; pp. 202–209. [Google Scholar]
  28. Gatheral, J. The Volatility Surface: A Practitioner’s Guide; John Wiley & Sons: New York, NY, USA, 2011. [Google Scholar]
  29. Tankov, P. Financial Modelling with Jump Processes; Chapman and Hall/CRC: Boca Raton, FL, USA, 2003. [Google Scholar]
  30. Kou, S.G. A jump-diffusion model for option pricing. Manag. Sci. 2002, 48, 1086–1101. [Google Scholar] [CrossRef] [Green Version]
  31. Huang, J.; Zhu, W.; Ruan, X. Option pricing using the fast fourier transform under the double exponential jump model with stochastic volatility and stochastic intensity. J. Comput. Appl. Math. 2014, 263, 152–159. [Google Scholar] [CrossRef]
  32. Vo, M.T. Regime-switching stochastic volatility: Evidence from the crude oil market. Energy Econ. 2009, 31, 779–788. [Google Scholar] [CrossRef]
  33. Yang, H. Esscher transform. Wiley Statsref Stat. Ref. Online 2014. [Google Scholar] [CrossRef]
  34. Merton, R.C. Option pricing when underlying stock returns are discontinuous. J. Financ. Econ. 1976, 3, 125–144. [Google Scholar] [CrossRef] [Green Version]
  35. Shreve, S.E. Stochastic Calculus for Finance II: Continuous-Time Models; Springer: Berlin/Heidelberg, Germany, 2004; Volume 11. [Google Scholar]
  36. Buffington, J.; Elliott, R.J. American options with regime switching. Financ. Stochastics 2002, 5, 497–514. [Google Scholar] [CrossRef]
  37. Buffington, J.; Elliott, R.J. Regime switching and european options. In Stochastic Theory and Control: Proceedings of a Workshop Held in Lawrence, Kansas; Springer: Berlin/Heidelberg, Germany, 2002; pp. 73–82. [Google Scholar]
  38. Liu, S.; Yang, Y.; Zhang, H.; Wu, Y. Variance swap pricing under markov-modulated jump-diffusion model. Discret. Dyn. Nat. Soc. 2021, 2021, 9814605. [Google Scholar] [CrossRef]
Figure 1. Comparison of the fair strike price obtained from model Monte Carlo simulation and our semi-Monte Carlo simulation. (a) K d 1 with regime 1. (b) K d 2 with regime 1.
Figure 1. Comparison of the fair strike price obtained from model Monte Carlo simulation and our semi-Monte Carlo simulation. (a) K d 1 with regime 1. (b) K d 2 with regime 1.
Mathematics 11 02730 g001aMathematics 11 02730 g001b
Figure 2. Comparison of the fair strike price obtained from regime-switching or non-regime-switching models. (a,b) K d i , i = 1 , 2 with regime 1. (c,d) K d i , i = 1 , 2 with regime 2.
Figure 2. Comparison of the fair strike price obtained from regime-switching or non-regime-switching models. (a,b) K d i , i = 1 , 2 with regime 1. (c,d) K d i , i = 1 , 2 with regime 2.
Mathematics 11 02730 g002
Figure 3. Comparison of the fair strike price obtained from various jump intensities. (a) K d 1 with regime 1. (b) K d 2 with regime 1.
Figure 3. Comparison of the fair strike price obtained from various jump intensities. (a) K d 1 with regime 1. (b) K d 2 with regime 1.
Mathematics 11 02730 g003
Figure 4. Comparison of the fair strike price obtained under the traditional MRG model with our new model.
Figure 4. Comparison of the fair strike price obtained under the traditional MRG model with our new model.
Mathematics 11 02730 g004
Table 1. Model parameters.
Table 1. Model parameters.
NotationParametersRegime 1Regime 2
TMaturity1.001.00
v ( 0 ) Initial volatility0.220.22
S ( 0 ) Initial stock price1.001.00
θ ˜ X ( t ) Long-term volatility mean 0.075 0.040
kMean-reverting speed2.002.00
σ Volatility of volatility0.100.10
r X ( t ) Interest rate0.080.06
ρ Correlation coefficient−0.40−0.40
λ ˜ X ( t ) Jump intensity0.300.20
pProbabilities of upward jumps0.200.20
qProbabilities of downward jumps0.800.80
η 1 Inverse mean one25.0025.00
η 2 Inverse mean two50.0050.00
Table 2. Fair strike price under various observation frequencies (regime 1).
Table 2. Fair strike price under various observation frequencies (regime 1).
Observation
Frequency
K d 1
(Regime-
Switching)
K d 2
(Regime-
Switching)
K d 1
(Non-Regime-
Switching)
K d 2
(Non-Regime-
Switching)
Quarterly (N = 4)684.6565655.9717685.9228661.9446
Monthly (N = 12)670.1038650.0076663.4953657.1487
Fortnightly (N = 26)657.2872647.8999657.8372656.8177
Weekly (N = 52)651.1446646.8401655.0941654.8623
Daily (N = 252)645.4533645.5934652.4541654.7870
Table 3. Model parameters sensitivity analysis.
Table 3. Model parameters sensitivity analysis.
ParametersValuePercent Change in K d 1 Percent Change in K d 2
v ( 0 ) 0.220.712377800.71249860
S ( 0 ) 1.000.000000000.00000000
θ ˜ X ( t ) 0.075 1.2309270001.23111800
k2.000.104577900.10459840
σ 0.100.057274340.05774045
r X ( t ) 0.080.0013876120.00043574
ρ −0.40−0.000314550.00015732
λ ˜ X ( t ) 0.300.006068600.00587247
p0.200.002702670.00222682
η 1 25.00−0.00725323−0.00598776
η 2 50.00−0.00538579−0.00594781
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zou, A.; Wang, J.; Wu, C. Pricing Variance Swaps under MRG Model with Regime-Switching: Discrete Observations Case. Mathematics 2023, 11, 2730. https://doi.org/10.3390/math11122730

AMA Style

Zou A, Wang J, Wu C. Pricing Variance Swaps under MRG Model with Regime-Switching: Discrete Observations Case. Mathematics. 2023; 11(12):2730. https://doi.org/10.3390/math11122730

Chicago/Turabian Style

Zou, Anqi, Jiajie Wang, and Chiye Wu. 2023. "Pricing Variance Swaps under MRG Model with Regime-Switching: Discrete Observations Case" Mathematics 11, no. 12: 2730. https://doi.org/10.3390/math11122730

APA Style

Zou, A., Wang, J., & Wu, C. (2023). Pricing Variance Swaps under MRG Model with Regime-Switching: Discrete Observations Case. Mathematics, 11(12), 2730. https://doi.org/10.3390/math11122730

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop