Next Article in Journal
Conformable Fractional Martingales and Some Convergence Theorems
Next Article in Special Issue
An Analysis and Comparison of Multi-Factor Asset Pricing Model Performance during Pandemic Situations in Developed and Emerging Markets
Previous Article in Journal
Similitude for the Dynamic Characteristics of Dual-Rotor System with Bolted Joints
Previous Article in Special Issue
Deep Learning Models for Predicting Monthly TAIEX to Support Making Decisions in Index Futures Trading
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Closed-Form Pricing Formula for Log-Return Variance Swaps under Stochastic Volatility and Stochastic Interest Rate

College of Mathematics, Harbin Normal University, Harbin 150000, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(1), 5; https://doi.org/10.3390/math10010005
Submission received: 24 November 2021 / Revised: 14 December 2021 / Accepted: 15 December 2021 / Published: 21 December 2021
(This article belongs to the Special Issue Mathematics and Financial Economics)

Abstract

:
At present, the study concerning pricing variance swaps under CIR the (Cox–Ingersoll–Ross)–Heston hybrid model has achieved many results; however, due to the instantaneous interest rate and instantaneous volatility in the model following the Feller square root process, only a semi-closed solution can be obtained by solving PDEs. This paper presents a simplified approach to price log-return variance swaps under the CIR–Heston hybrid model. Compared with Cao’s work, an important feature of our approach is that there is no need to solve complex PDEs; a closed-form solution is obtained by applying the martingale theory and It o ^ ’s lemma. The closed-form solution is significant because it can achieve accurate pricing and no longer takes time to adjust parameters by numerical method. Another significant feature of this paper is that the impact of sampling frequency on pricing formula is analyzed; then the closed-form solution can be extended to an approximate formula. The price curves of the closed-form solution and the approximate solution are presented by numerical simulation. When the sampling frequency is large enough, the two curves almost coincide, which means that our approximate formula is simple and reliable.

1. Introduction

Since the break of the global financial crisis in 2008, with the sharp rise and fall of the stock market, financial markets have shown high volatility and risk. Almost all financial markets in the world have experienced volatility risk; some local financial problems are likely to spread and become a serious global problem with the globalization of economies. Under this background, trading or hedging this risk has become more and more important. The management of volatility risk has attracted extensive attention from financial practitioners. Both individuals and financial institutions are very interested in trading volatility/variance to obtain returns under high risk or effectively hedge volatility risk. Volatility derivatives have gradually become an important financial instrument. Carr and Lee [1] pointed out that the trading volume of volatility derivatives in the current financial market is increasing significantly. Among many volatility derivatives, volatility/variance swaps are the most popular. They are essentially a kind of forward contract. The swap parties agree to swap the volatility/variance (realized volatility/variance) generated by the price of the underlying asset during the contract period with the fixed value agreed in advance on the maturity date. In other words, the holder of the contract actually exchanges the future uncertain volatility/variance with the fixed volatility/variance value, so as to avoid the volatility risk.
In the financial market, traders prefer trading volatility swaps directly. However, due to the convexity of the square root process, the pricing problem of volatility swaps is very complex. Financial practitioners try to weaken the condition and apply variance to describe the volatility risk. The advantage of variance swaps lies in its good mathematical properties (such as additivity), which enables scholars to systematically and accurately study variance swaps and provide investors with accurate volatility risk exposure. Moreover, the pricing of volatility swaps is often approached by the results of variance swaps. Therefore, the pricing of variance swaps has attracted more attention from financial practitioners. Many scholars study variance swaps by establishing stochastic volatility models to describe the price process of underlying assets [2,3,4,5]. Swishchuk [6,7] defines and studies the single factor and multi factor stochastic volatility model with delay, they presented the pricing formula of variance swaps under the corresponding model. When the volatility process of the underlying asset satisfies the BN-S (Barndorff-Nielsen and Shephard) model, Habtemicael and Sengupta [8] studied the pricing problem of variance and volatility swaps under the assumption of no arbitrage. Issaka and Sengupta [9] obtained various results related to variance swaps under the BN-S model. They also introduced a price-weighted index modulated by market variance and verified its impact on pricing. In addition, Sengupta et al. [10] designed a combination of futures and variance swaps under the BN-S model, which can optimally hedge oil price risk.
The above literature is carried out under the assumption of continuous sampling. However, Jarrow et al. [11] pointed out that the results obtained under continuous sampling can only be seen as an approximation of the actual financial market. For more accurate pricing, scholars tried to solve the pricing problem of volatility and variance swaps under the assumption of discrete sampling. Refs. [12,13,14] show that the pricing of derivatives under discrete sampling assumption is more complex than continuous sampling assumption. The results of variance swaps with discrete sampling are usually obtained based on the stochastic volatility model. Little and Pant [15] proposed a finite difference method, by reducing dimension and introducing new variables, they solved the pricing problem of variance swaps under Black-Scholes framework. A path dependent option pricing simulation method under Heston [16] stochastic volatility model is proposed in [17]. Zhu and Lian [18,19] derived the analytical solution of variance swaps based on Heston [16] stochastic volatility model. Similarly, Zhang [20] obtained the analytical solution of variance swaps under the MRG (mean-reversing Gaussian) stochastic volatility model. It is worth noting that if only a single factor stochastic volatility model is considered, the results can not truly reflect the actual financial market. Therefore, scholars introduce other volatility factors to the model to better adapt to the data in the financial market. Jump process is often added into the underlying asset price process. Jain [21] proves that the introduction of jump diffusion process will affect the result of variance swaps. Liu [22] obtained a analytical solution of variance swaps based on the Hawkes jump diffusion process. In the mean time, the stochastic interest rate is also embedded into the underlying asset price process. Kim [23] proves that the introduction of stochastic interest rate can indeed get better results. Cao et al. [24] combined the Heston [16] stochastic volatility model with Cox-Ingersoll-Ross (CIR) [25] stochastic interest rate model, using the dimension reduction technology in Little and Pant [15] and characteristic function, they obtained a semi-closed solution of variance swaps under the CIR–Heston hybrid model. Cao et al. [26] further extended the partial coefficient correlation model to the full correlation structure model, and presented a semi-closed solution through the derivation of characteristic function. Based on the work of [20,24], Zhao [27] introduced the Vasicek interest rate process to the MRG(mean-reverting Gaussian) model to obtain a closed-form solution. Recently, Xu et al. [28] obtained the pricing formula of variance swaps with the liquidity risk of the underlying assets. As to the pricing of variance swaps under Markov modulation model, readers can refer to [29,30].
In recent years, the pricing problem of log-return variance swaps has also made great progress. Zhu and Lian [31] obtained an exact solution of the PDE system based on the Heston’s [16] two-factor model. Bernard and Cui [32] solved the pricing problem of log-return variance swaps with discrete sampling under three different stochastic volatility models. Based on Heston’s [16] two factor stochastic volatility model, Zhu and Lian [33] proposed a general method of forward start variance swaps with discrete sampling. By using the forward characteristic function, two analytical solutions of forward start variance swaps are derived. Furthermore, there are some studies on pricing log-return variance swaps under other models [34,35]. As to the national policies, in order to deal with some possible risks in the future, which can adjust the net reserve value and promote sustainable economic growth, readers can refer to [36].
To our best knowledge, the pricing problem of variance swaps based on the CIR–Heston hybrid model has been deeply studied in Cao [24,26] by PDEs. However, solving PDEs is a very cumbersome process and only a semi-closed solution can be obtained. In contrast, this paper study the problem from the perspective of stochastic analysis without solving PDEs. By applying the martingale theory and It o ^ ’s lemma, the solving process can be greatly simplified and a closed-form solution of log-return variance swaps is obtained. The significant advantage of our closed-form solution is there is no need to spend time adjusting parameters by numerical simulation. In addition, one limitation of solving the pricing problem by PDEs is not convenient to analyze the influence of sampling frequency on the pricing formula. Another important feature of this paper is that the closed-form solution can be further extended to a simple and reliable approximate formula. The numerical simulation results show that there is little difference between the exact solution and the approximate solution when the sampling frequency is large enough.
The rest of this paper is organized as follows: In Section 2, we give a description of CIR–Heston hybrid model and obtain a closed-form formula for discretely sampled variance swaps. In Section 3, using It o ^ ’s lemma and the martingale theory, we derive an approximate formula in the case of large sampling. The numerical analysis is shown to illustrate our main results in Section 4. In Section 5, we conclude the results.

2. Pricing Variance Swaps and Our Model

In this section, we introduce the related concepts of variance swaps and CIR–Heston bybrid model, and we also show the solution of variance swaps valuation in detail.

2.1. CIR–Heston Hybrid Model

Let ( Ω , F t , Q ) is a probability measurable space, the price of the underlying asset S t , the instantaneous interest rate r t and the instantaneous stochastic volatility v t , can be described by the following CIR–Heston hybrid model:
d S t = r t S t d t + v t S t d W t S d v t = k ( θ v t ) d t + σ v t d W t v d r t = h ( η r t ) d t + ξ r t d W t r
θ is the long-term mean of volatility, σ is the volatility of volatility, η is the long-term mean of the instantaneous interest rate, ξ is the volatility of the interest rate. k and h are the mean-reverting speed parameter of v t and r t , respectively. W t S , W t v , W t r are three one-dimensional Brownian motions. We make the following assumptions:
(1)
d W t S d W t v = ρ d t , d W t S d W t v = 0 , d W t v d W t r = 0 , and 1 ρ 1 .
(2)
To ensure the value of r t and v t are always positive, we set 2 k θ σ 2 , 2 ξ h η 2 .
(3)
All the parameters are denoted under the risk-neutral measure Q .

2.2. Variance Swaps

We firstly give a introduction for the measurement of the realized variance of the underlying asset, which is expressed as the sum of squares of log-returns during the contract period in [15]:
σ R 2 = A F N i = 1 N ln 2 ( S t i S t i 1 ) × 100 2
where S t i is the price of the underlying asset at time t i . Suppose S t has N observations in the period [0,T]: 0 t 1 < t 2 < t 3 < t N T . 100 2 is the parameter used to adjust the result. A F is the annualized factor converting the expression to an annualized variance, related to sampling frequency. If the sampling frequency is every trading day, then A F = 252 , which means that there are 252 trading days in one year; If the sampling frequency is every week, then A F = 52 . If every month, then A F = 12 , Typically, assuming equally spaced discrete observations and take A F = 1 Δ t = N T .
Variance swaps is a kind of forward contract about the variance of the future return of the underlying asset. The long position of the contract pay a fixed value on the maturity date and receive the amount of the realized variance of the underlying asset during the contract period, while the short position just the opposite. Suppose the price of an underlying asset is S t with the maturity date T, the return of the long position on the maturity date is:
( σ R 2 K ) × L
where K is the annualized delivery price of the variance swaps; L is the notional amount of the swap in dollars per annualized variance point, also known as nominal principal.

2.3. Measure Transformation

According to the risk-neutral pricing principle, the value of variance swaps contract at time t is the present value of future expected return, thus:
V t = E t Q [ e t T r s d s ( σ R 2 K ) × L ]
Due to the contract is fair to both parties, therefore, the price of log-return variance swaps contract is actually the value of K that ensures V t = 0 at the initial time, thus K = E 0 Q [ e 0 T r s d s σ R 2 ] . It is quite difficult to directly calculate the value of K because of the stochastic interest rate, we need to change the expression (4) from the risk-neutral measure Q to the T-forward measure Q T . According to Brigo and Mercurio [37], the value of zero coupon bonds at time t = 0 is given by E 0 Q [ e 0 T r s d s ] = P ( 0 , T ) . Then, Equation (4) can be written as
V t = E t Q [ e t T r s d s ( σ R 2 K ) × L ] = P ( 0 , T ) E 0 T [ ( σ R 2 K ) × L ]
where E 0 T [ · ] is the expectation under the T-forward measure Q T at time t = 0 . Hence, the price of log-return variance swaps contract K should satisfies K = E 0 T ( σ R 2 ) . Further, according to the mathematical expression of the realized variance in (2):
E 0 T [ σ R 2 ] = E 0 T [ 1 T i = 1 N ln 2 ( S t i S t i 1 ) ] × 100 2 = 1 T E 0 T [ i = 1 N ln 2 ( S t i S t i 1 ) ] × 100 2
the process of computing the value of K can be reduced to compute the following N conditional expectation:
E 0 T [ ln 2 ( S t i S t i 1 ) ] , i = 1 , , N
Obviously, for continue solving, we also need to change SDEs in (1) from the risk-neutral measure Q to the T-forward measure Q T . By using Cholesky decomposition, we can obtain:
d S t S t d v t d r t = r t k ( θ v t ) h ( η r t ) d t + × C × d W ^ t S d W ^ t v d W ^ t r
where
= v t 0 0 0 σ v t 0 0 0 ξ r t , C = 1 0 0 ρ 1 ρ 2 0 0 0 1
and
C C T = 1 ρ 0 ρ 1 0 0 0 1 , d W ˜ t S d W ˜ t v d W ˜ t r = C × d W ^ t S d W ^ t v d W ^ t r
According to Brigo and Mercuio [37] and Cao [24], the expression of SDEs in (1) under the T-forward measure Q T can be represented as:
d S t = r t S t d t + v t S t d W ^ t S d v t = k ( θ v t ) d t + σ v t ( ρ d W ^ t S + 1 ρ 2 d W ^ t v ) d r t = ( h η [ h + B ( t , T ) η 2 ] r t ) d t + ξ r t d W ^ t r
where
B ( t , T ) = 2 ( e ( T t ) h 2 + 2 η 2 1 ) 2 h 2 + 2 η 2 + ( h + h 2 + 2 η 2 ) ( e ( T t ) h 2 + 2 η 2 1 )
and W ^ t S , W ^ t v , W ^ t r are three independent Brownian motions.
Note that the form of d r t in (8) is complex, in order to make the solution process more clear, let
h * = h + B ( t , T ) η 2 , η * = h η h + B ( t , T ) η 2
then we can change (8) to the following system under Q T
d S t = r t S t d t + v t S t d W ^ t S d v t = k ( θ v t ) d t + σ v t ( ρ d W t S ^ + 1 ρ 2 d W ^ t v ) d r t = h * ( η * r t ) d t + ξ r t d W ^ t r

2.4. Pricing Formula for Variance Swaps

We present a two-step approach to compute the value of K. First of all, the filtration satisfies that F 0 F t i 1 F t i , in terms of the tower property of conditional expectation:
E 0 T [ ln 2 ( S t i S t i 1 ) ] = E 0 T [ E t i 1 T [ ln 2 ( S t i S t i 1 ) ] ]
where E t i 1 T [ · ] is defined as the conditional expectation with respect to F t i 1 under Q T .
Hence, the first step is to compute E t i 1 T [ · ] , S t i 1 , r t i 1 and v t i 1 are F t i 1 measurable.
Set X t = ln ( S t ) , based on (8), X t satisfies the following dynamics:
d X t = ( r t 1 2 v t ) d t + v t d W ^ t S .
To compute the conditional moment of X t , we need to calculate the conditional moments of r t and v t firstly, which are shown in the following theorem.
Theorem 1.
Suppose that S t follows the dynamics described in (12), for any 0 t i 1 t , we have:
E t i 1 T ( r t ) = η * + ( r t i 1 η * ) e h * ( t t i 1 )
E t i 1 T ( v t ) = θ + ( v t i 1 θ ) e k ( t t i 1 )
V a r t i 1 ( r t ) = ξ 2 [ η * 2 h * ( 1 e 2 h * ( t t i 1 ) ) + r t i 1 θ h * ( e h * ( t t i 1 ) e 2 h * ( t t i 1 ) ) ]
V a r t i 1 ( v t ) = σ 2 [ θ 2 k ( 1 e 2 k ( t t i 1 ) ) + v t i 1 θ k ( e k ( t t i 1 ) e 2 k ( t t i 1 ) ) ]
the conditional moments are F t i 1 measurable, thus: E t i 1 T = E T ( · | F t i 1 ) , V a r t i 1 = V a r ( · | F t i 1 ) .
Proof. 
The proof of this theorem is left in Appendix A. ☐
In order to obtain the expectation of ln 2 ( S t i S t i 1 ) , Cao [24,26] use the characteristic function and the generalized Fourier transform to solve PDEs, however, the computation process is cumbersome. We present a simpler approach. Considering the definition of variance
E t i 1 T [ ln 2 ( S t i S t i 1 ) ] = E t i 1 T [ ln ( S t i S t i 1 ) ] 2 + V a r t i 1 [ ln ( S t i S t i 1 ) ]
we transform the problem into calculating the expectation and variance of log-return with respect to F t i 1 . The result are presented by Theorem 2.
Theorem 2.
If S t follows the dynamics described in (10). Let Δ t = t i t i 1 , then we have
E t i 1 T [ ln ( S t i S t i 1 ) ] = P 1 ( Δ t ) + P 2 ( Δ ) ( r t i 1 η * ) + P 3 ( Δ ) ( v t i 1 θ )
V a r t i 1 [ ln ( S t i S t i 1 ) ] = Q 1 ( Δ t ) + Q 2 ( Δ t ) ( r t i 1 η * ) + Q 3 ( Δ t ) ( v t i 1 θ )
where S t i 1 , r t i 1 , v t i 1 are F t i 1 measurable, and
P 1 ( Δ t ) = ( η * 1 2 θ ) Δ t , P 2 ( Δ t ) = 1 h * ( 1 e ( h * ) 2 Δ t ) , P 3 ( Δ t ) = 1 2 k ( 1 e k Δ t )
Q 1 ( Δ t ) = ( ξ 2 ( h * ) 2 η * + θ σ 2 ( 1 ρ 2 ) 4 k 2 + θ θ ρ σ k + σ 2 ρ 2 θ 4 k 2 ) Δ t ( θ σ 2 ( 1 ρ 2 ) 2 k 3 + θ ρ σ k 2 + σ 2 ρ 2 θ 2 k 3 ) ( 1 e k Δ t ) + ( θ σ 2 ( 1 ρ 2 ) 8 k 3 + θ σ 2 ρ 2 8 k 3 ) ( 1 e 2 k Δ t ) 2 ξ 2 ( h * ) 3 η * ( 1 e h * Δ t ) + ξ 2 2 ( h * ) 3 η * ( 1 e 2 h * Δ t ) σ 2 ρ 2 2 k 2 e k Δ t Δ t Q 2 ( Δ t ) = 2 ξ 2 ( h * ) 2 ( Δ t e h * Δ t ) + ξ 2 ( h * ) 3 ( 1 e 2 h * Δ t ) Q 3 ( Δ t ) = ( σ 2 ( 1 ρ 2 ) 2 k 2 + ρ σ k ) ( Δ t e k Δ t ) + ( σ 2 ρ 2 4 k 3 + 1 k σ ρ k 2 ) ( 1 e k Δ t ) + σ 2 ( 1 ρ 2 ) 4 k 3 ( 1 e 2 k Δ t ) + σ 2 ρ 2 4 k 3 ( e k Δ t e 2 k Δ t )
Proof. 
We show the details of the derivation in Appendix B. ☐
Square (18) on both sides,
E t i 1 T [ ln ( S t i S t i 1 ) ] 2 = P 1 2 ( Δ t ) + P 2 2 ( r t i 1 η * ) + P 3 2 ( Δ t ) ( v t i 1 θ ) 2 + 2 P 1 ( Δ t ) P 2 ( Δ t ) ( r t i 1 η * ) + 2 P 1 ( Δ t ) P 3 ( Δ t ) ( v t i 1 θ ) + 2 P 2 ( Δ t ) P 3 ( Δ t ) ( r t i 1 η * ) ( v t i 1 θ )
then plug (19) and (22) into (17) and the expression of inner expectation under the information set F t i 1 is given by
E t i 1 T [ ln 2 ( S t i S t i 1 ) ] = D 1 ( Δ t ) + D 2 ( Δ t ) ( r t i 1 η * ) + D 3 ( Δ t ) ( v t i 1 θ ) + D 4 ( Δ t ) ( r t i 1 η * ) 2 + D 5 ( Δ t ) ( v t i 1 θ ) 2 + D 6 ( Δ t ) ( r t i 1 η * ) ( v t i 1 θ )
where:
D 1 ( Δ t ) = Q 1 ( Δ t ) + P 1 2 ( Δ t ) D 2 ( Δ t ) = 2 P 1 ( Δ t ) P 2 ( Δ t ) + Q 2 ( Δ t ) D 3 ( Δ t ) = 2 P 1 ( Δ t ) P 3 ( Δ t ) + Q 3 ( Δ t ) D 4 ( Δ t ) = P 2 2 ( Δ t ) D 5 ( Δ t ) = P 3 2 ( Δ t ) D 6 ( Δ t ) = 2 P 2 ( Δ t ) P 3 ( Δ t )
Next, the second step is to compute the outer expectation. Note that D k ( Δ t ) , k = 1 , , 6 are all constants, so only the expectations of r s and v s under the information set F 0 need to be computed. The result can be derived by the following theorem.
Theorem 3.
For any i = 1 , 2 , 3 , , N , v 0 , r 0 > 0 , the price of log-return variance swaps in discrete sampling at the initial moment can be written as
E 0 T [ ln 2 ( S t i S t i 1 ) ] = D ˜ 1 ( Δ t ) + D ˜ 2 ( Δ t ) ( r 0 η * ) + D ˜ 3 ( Δ t ) ( v 0 θ ) + D ˜ 4 ( Δ t ) ( r 0 η * ) 2 + D ˜ 5 ( Δ t ) ( v 0 θ ) 2 + D ˜ 6 ( Δ t ) ( r 0 η * ) ( v 0 θ )
where Δ t = t i t i 1 , and
D 1 ˜ ( Δ t ) = D 1 ( Δ t ) + D 4 ( Δ t ) η * ξ 2 2 h * ( 1 e 2 h * t i 1 ) + D 5 ( Δ t ) θ σ 2 2 k ( 1 e 2 k t i 1 ) D 2 ˜ ( Δ t ) = e h * t i 1 D 2 ( Δ t ) + D 4 ( Δ t ) σ 2 h * ( e h * t i 1 e 2 h * t i 1 ) D 3 ˜ ( Δ t ) = e k t i 1 D 3 ( Δ t ) + D 5 ( Δ t ) σ 2 k ( e k t i 1 e 2 k t i 1 ) ) D 4 ˜ ( Δ t ) = e 2 h * t i 1 D 4 ( Δ t ) D 5 ˜ ( Δ t ) = e 2 k t i 1 D 5 ( Δ t ) D 6 ˜ ( Δ t ) = e h * t i 1 e k t i 1 D 6 ( Δ t )
Proof. 
Refer to Appendix C. ☐
Now, we obtain the strike price K for the variance swaps contract at the initial time as
K = E 0 T [ σ R 2 ] = 100 2 T i = 1 N E 0 T [ ln 2 ( S t i S t i 1 ) ]
= 100 2 T 1 ( r 0 η * ) ( v 0 θ ) ( r 0 η * ) 2 ( v 0 θ ) 2 ( r 0 η * ) ( v 0 θ ) T i = 1 N D ˜ 1 ( Δ , t i 1 ) i = 1 N D ˜ 2 ( Δ , t i 1 ) i = 1 N D ˜ 3 ( Δ , t i 1 ) i = 1 N D ˜ 4 ( Δ , t i 1 ) i = 1 N D ˜ 5 ( Δ , t i 1 ) i = 1 N D ˜ 6 ( Δ , t i 1 )
Remark 1.
The closed-form solution of the log-return variance swaps consists of constant term, first-order term ( r 0 η * ) , ( v 0 θ ) and second-order term ( r 0 η * ) 2 , ( v 0 θ ) 2 , ( r 0 η * ) ( v 0 θ ) . The expression in this form is very advantageous. Compared with Cao [24,26], our formula can be intuitively observed that the parameter θ will have great impact on the results. Moreover, if r 0 , v 0 and η * , θ are relatively close respectively, ( r 0 η * ) 2 , ( v 0 θ ) 2 , ( r 0 η * ) ( v 0 θ ) will tend to be 0, the second terms will have little impact on the results.
Remark 2.
Another contribution of our method is to consider the influence of sampling frequency. Roughly speaking, when the sampling frequency is large enough, the higher-order term about Δ t will tend to 0, then the exact solution can be extended to a simple and reliable approximate formula. The detailed derivation process and results will be shown in Section 3.

3. The Approximate Formula

As mentioned above, when the sampling frequency is large enough, the first-order term of Δ t will play a major role in the exact formula. Then the pricing formula can be simplified in theory. However, if we take Taylor expansion to the higher-order directly, the calculation is very complex [22]. The feature of our approach is to note the fact that both v t and r t are bounded functions, then we deduce a simple and reliable approximate formula by applying the integral mean value theorem and It o ^ ’s lemma. We show the result in the form of Theorem 4.
Theorem 4.
For any i = 1 , 2 , 3 , , N , v 0 , r 0 > 0 , if the sampling frequency is large enough, N , Δ t 0 , the price of log-return variance swaps in discrete sampling at the initial moment can be written as
K = E 0 T [ σ R 2 ( 0 , N , T ) ] = 100 2 T i = 1 N E 0 T [ ln 2 ( S t i S t i 1 ) ] = 100 2 × [ θ + 1 T i = 1 N ( v 0 θ ) e k t i 1 Δ t ) ] + O ( 1 N )
Proof. 
Integrating both sides to (12) from t i 1 to t i , the log-return of the underlying assets can be written as:
ln ( S t i S t i 1 ) = t i 1 t i ( r t 1 2 v t ) d t + t i 1 t i v t d W ^ t S
Divide both sides of the Equation (27) by N Δ t and add up from 1 to N, then take conditional expectation.
i = 1 N 1 N Δ t E 0 T [ ln 2 ( S t i S t i 1 ) ] = i = 1 N 1 N Δ t E 0 T [ t i 1 t i ( r t 1 2 v t ) d t ] 2 + i = 1 N 1 N Δ t E 0 T [ t i 1 t i v t d W ^ t S ] 2 + i = 1 N 1 N Δ t 2 E 0 T [ t i 1 t i ( r t 1 2 v t ) d t t i 1 t i v t d W ^ t S ]
As the sampling frequency increases, N , then Δ t 0 and Δ u 0 .
For the first item to the right (28), let φ ( t , r t , v t ) = r t 1 2 v t 2 . Since φ ( t , r t , v t ) is a bounded continuous function clearly, hence, according to the integral mean value theorem, there is a ξ [ t i 1 , t i ] satisfied
t i 1 t i ( r t 1 2 v t ) d t = φ ( ξ , r ξ , v ξ ) Δ t
According to (29), the following results can be derived:
lim Δ t 0 i = 1 N 1 N Δ t E 0 T [ t i 1 t i ( r t 1 2 v t 2 ) d t ) 2 ] = lim Δ t 0 i = 1 N 1 N Δ t E 0 T [ φ 2 ( ξ , r ξ , v ξ ) Δ t 2 ] = Δ t · E 0 T [ φ 2 ( ξ , r ξ , v ξ ) ] = 0
Using the properties of the It o ^ ’s integral again, the second item to the right (28) is simply computed:
i = 1 N 1 N Δ t E 0 T [ t i 1 t i v t d W ^ t S ] 2 = i = 1 N 1 N Δ t E 0 T [ t i 1 t i v t d t ] = 1 T E 0 [ 0 T v t d t ]
For the third item on the right (28), by substituting the (A1) and (A2), we can obtain:
i = 1 N 1 N Δ t 2 E 0 T [ t i 1 t i ( r t 1 2 v t ) d t t i 1 t i v t d W ^ t S ] = i = 1 N 1 N Δ t { E 0 T [ 2 t i 1 t i r t d t · t i 1 t i v t d W ^ t S ] E 0 T [ t i 1 t i v t d t · t i 1 t i v t d W ^ t S ] } = i = 1 N 1 N Δ t { E 0 T 2 t i 1 t i η * + ( r t i 1 η * ) e h * ( t t i 1 ) d t · t i 1 t i v t d W ^ t S + E 0 T 2 t i 1 t i ξ t i 1 t e h * ( t u ) r u d W ^ u r d t · t i 1 t i v t d W ^ t S E 0 T t i 1 t i θ + ( v t i 1 θ ) e k ( t t i 1 ) d t · t i 1 t i v t d W ^ t S E 0 T σ t i 1 t i e k ( t u ) v u ( ρ d W ^ u S + 1 ρ 2 d W ^ u v ) d t · t i 1 t i v t d W ^ t S }
We divide (32) into the following four parts: (33), (34), (35), (36) and compute the value of each part:
E 0 T [ 2 t i 1 t i ( η * + ( r s η * ) e h * ( t s ) ) d t · t i 1 t i v t d W ^ t S ] = E 0 T [ ( η * Δ t + r s η * h * ( 1 e h * Δ t ) ) · t i 1 t i v t d W ^ t S ] = 0
E 0 T [ 2 t i 1 t i ξ s t e h * ( t u ) r u d W ^ u r d t · t i 1 t i v t d W ^ t S ] = E 0 T [ 2 t i 1 t i ξ u t i e h * ( t u ) d t r u d W ^ u r · t i 1 t i v t d W ^ t S ] = 0
E 0 T [ t i 1 t i ( θ + ( v s θ ) e k ( t s ) ) d t · t i 1 t i v t d W ^ t S ] = E 0 T [ ( θ Δ t + v s θ k ( 1 e k Δ t ) ) · t i 1 t i v t d W ^ t S ] = 0
E 0 T [ σ s t e k ( t u ) v u d W ^ u v d t · t i 1 t i v t d W ^ t S ] = E 0 T [ t i 1 t i σ u t i e k ( t u ) d t v u ( ρ d W u S ^ + 1 ρ 2 d W ^ u v ) · t i 1 t i v t d W ^ t S ] = E 0 T [ t i 1 t i σ u t i e k ( t u ) d t · v u d u ] = E 0 T [ σ e k ( λ δ ) · v δ · Δ t · Δ u ]
where λ , δ [ t i 1 , t i ] . This is the result of using the Integral Median Theorem again.
lim Δ t 0 , Δ u 0 i = 1 N 1 N Δ t E 0 T [ σ e k ( λ δ ) · v δ · Δ t · Δ u ] = Δ u · E 0 T [ σ e k ( λ δ ) · v δ ] = 0
Finally, when N and Δ t 0 , the price of log-return variance swaps at 0 moment can be simplified as
lim Δ t 0 E 0 T [ σ R 2 ( 0 , N , T ) ] = 1 T E 0 T [ 0 T v t d t ]
Thus, we can write (38) in discrete form (26). ☐
In particular, (38) is the price obtained under the continuous model.
Obviously, compared with the closed-form solution (25), the approximate Formula (26) has a simpler form, and numerical simulation in Section 4 also verifies the reliability of (26). What’s more, the interest rate is the volatility of r t does not appear in the approximate formula. This means the influence of stochastic interest rate becomes smaller and smaller with the increase of sampling frequency.

4. Numerical Analysis

In this section, we show some numerical examples for illustration purpose. We present the results of the closed-form formula, Monte Carlo simulations, the first-order approximation and the continuous model that can help readers understand of our pricing formula intuitively. Compared with Monte Carlo simulation, we verify the validity of our closed-solution and the approximation formula. In addition, according to the approximate formula, we infer that the impact of stochastic interest rate on the results decreases with the increase of sampling frequency. We will also verify this point.
To achieve these purposes, we make the following assumptions:
(1) The values of the parameters involved in the pricing formula are
r 0 = 5 % , v 0 = 0.04 , k = 11.35 , θ = 0.022 , σ = 0.618 , ρ = 0.64 , ξ = 0.01 , h = 1.2 , η = 0.05 , T = 1 .
In particular, the above parameters are obtained from Zhu and Lian [18], and Cao [24].
(2) The asset price S 0 = 10 .
(3) The number of the paths N = 100 , 000 for all the MC simulations presented here.

4.1. Monte Carlo Simulations

We firstly adopt the simple Euler-Maruyama discretization for the CIR–Heston hybrid model in our MC simulations:
d S t = S t 1 + r t 1 S t 1 Δ t + v t 1 S t 1 Δ t W 1 , t d v t = v t 1 + k ( θ v t 1 ) Δ t + σ v t 1 Δ t ( ρ W 1 , t + 1 ρ 2 W 2 , t ) d r t = r t 1 + h ( η r t 1 ) Δ t + ξ r t 1 Δ t W 3 , t
where W 1 , t , W 2 , t , W 3 , t , are three independent standard normal random variables.
Some Monte Carlo simulations, the closed-form solution (25), the approximate solution (26) and the continuous model (38) price curve are carried out with MATLAB software. We show the results in Figure 1. We can clearly observed from Figure 1 that the closed-form solution matches well with the results of some Monte Carlo simulation. However, Monte Carlo simulation is expected to takes 6705.83 s, and our closed-form solution takes only 0.0223 s.
In particular, compared with the semi-solution in Cao [24,26], there is no need to adjust the parameters for our closed-form solution. In the mean time, with the increase of sampling frequency, the price curve of the closed-form solution and the approximate solution is closer and closer, eventually converge to the continuous case. Specifically, we select the prices corresponding to the above curve under different sampling frequencies, and put them together with some Monte Carlo simulation results in Table 1. According to the data in Table 1, when the sampling frequency N = 4 , the error between the closed-form solution and the approximate solution is 21.31. When the sampling frequency N reaches 200, the error is only 0.67.

4.2. Reliability of Approximate Formula

By observing the form of approximate formula, when the sampling frequency is large enough, the influence of stochastic interest rate is very small. This means that there is little difference between the results of different stochastic interest rates and the approximate formula in the case of large sampling. In order to verify this point, we need to compare the results of stochastic interest rates under different parameters. According to Cao [24] and Zhao [27], the long-term mean of stochastic interest rate η has a great influence on the formula, but the mean-reverting speed h and volatility ξ have little impact. Hence, to improve efficiency, we only select different long-term mean values, then take h = 1.2 and ξ = 0.01 . Figure 2 shows that when the sampling frequency is small, there is a large difference between the variance swaps prices calculated by the long-term mean of different stochastic interest rates; However, with the increase of sampling frequency, the difference between these prices and the price of the approximate formula becomes smaller and smaller. This fact shows that the influence of stochastic interest rate is indeed very small when the sampling frequency is large enough. What’s more, when the sampling frequency reaches 150, the maximum error of the approximate formula does not exceed 0.8, which means that the approximate formula proposed in this paper is reliable.

5. Conclusions

In this paper, we studied the problem of pricing log-return variance swaps under the CIR–Heston hybrid model. Compared with Cao’s work, the contribution of our work is to propose a more concise approach from the perspective of stochastic analysis. There is no need to solve complex PDEs, and we obtained a closed-form solution instead of a semi-closed solution. In the mean time, considering the influence of sampling frequency on the pricing formula, we further extend the closed-form solution to an approximate formula. The advantage of the approximate solution is a simpler form and reliable in the case of large sampling. Some numerical simulation shows that our closed-form solution matches well with the results of MC simulation. By comparing the price curves of the closed-form solution and the approximate solution, we conclude that the error become smaller and smaller with the increase of sampling frequency. Moreover, in the case of large samples, the error between the closed-form solution determined from different stochastic interest rates and the approximate solution is quite little. Therefore, we have reason to confirm that the influence of stochastic interest rate on the pricing formula does decrease with the increase of sampling frequency and our approximate formula is reliable.
Finally, it should be noted that our pricing approach is general. It also can be extended to other stochastic interest rate and stochastic volatility models, such as the Heston–Hull–White hybrid model. However, when the model becomes a full coefficient correlation, how to obtain the pricing formula by the method proposed in this paper is a problem worthy of study in the future.

Author Contributions

Conceptualization, Y.W.; Data curation, G.L.; Formal analysis, C.M.; Methodology, C.M. and G.L.; Project administration, G.L. and Y.W.; Software, C.M.; Writing—original draft, C.M.; C.M., Y.W. and G.L.: conceived of the presented idea, developed the theory and performed the computations, discussed the results, wrote the paper, and approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 11471091 and 12101163).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not Applicable.

Acknowledgments

The authors are especially grateful to the anonymous referees for their constructive suggestions and helpful comments to improve our paper.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Proof of Theorem 1. 
For any 0 t i 1 < t , upon applying It o ^ ’ s lemma on d ( e h * ( t u ) r u ) and d ( e k ( t u ) v u ) with respect to u and integrate both sides from t i 1 to t, we obtain
r t = η * + ( r t i 1 η * ) e h * ( t t i 1 ) + ξ t i 1 t e h * ( t u ) r u d W ^ u r
v t = θ + ( v t i 1 θ ) e k ( t t i 1 ) + σ t i 1 t e k ( t u ) v u ( ρ d W ^ t S + 1 ρ 2 d W ^ t v )
η * + ( r t i 1 η * ) e h * ( t t i 1 ) , θ + ( v t i 1 θ ) e k ( t t i 1 ) can be considered as a constant on taking conditional moments under information set F t i 1 , and ξ t i 1 t e h * ( t u ) d W u r and σ t i 1 t e k ( t u ) v u ( ρ d W ^ t S + 1 ρ 2 d W ^ t v ) can be regard as a martingale. Hence, the conditional expectation can be calculated as
E t i 1 T ( r t ) = η * + ( r t i 1 η * ) e h * ( t t i 1 ) E t i 1 T ( v t ) = θ + ( v t i 1 θ ) e k ( t t i 1 )
According to the property of the It o ^ ’ s integral, the variance of r t and v t can be written as
V a r t i 1 ( r t ) = V a r t i 1 ( ξ t i 1 t e h * ( t u ) r u d W ^ u r ) = ξ 2 t i 1 t e 2 h * ( t u ) E t i 1 T ( r u ) du = ξ 2 [ η * 2 h * ( 1 e 2 h * ( t t i 1 ) ) + r t i 1 η * h * ( e h * ( t t i 1 ) e 2 h * ( t t i 1 ) ) ]
V a r t i 1 ( v t ) = V a r t i 1 ( σ t i 1 t e k ( t u ) v u ( ρ d W ^ u S + 1 ρ 2 d W ^ u v ) ) = σ 2 t i 1 t e 2 k ( t u ) E t i 1 T ( v u ) d u = σ 2 [ θ 2 k ( 1 e 2 k ( t t i 1 ) ) + v t i 1 θ k ( e k ( t t i 1 ) e 2 k ( t t i 1 ) ) ]
 ☐

Appendix B

Proof of Theorem 2. 
During the period [ t i 1 , t i ] , integrate (12) on both sides from s to t, the result can be expressed as
ln ( S t i S t i 1 ) = t i 1 t i ( r t 1 2 v t ) d t + t i 1 t i v t d W ^ t s
Plug (A1) and (A2) into (A6), we can obtain
ln ( S t i S t i 1 ) = t i 1 t i ( η * + ( r t i 1 η * ) e h * ( t t i 1 ) + ξ t i 1 t i e h * ( t u ) r u d W ^ u r ) d t 1 2 t i 1 t i ( θ + ( v t i 1 θ ) e k ( t t i 1 ) + σ t i 1 t e k ( t u ) v u ( ρ d W ^ u S + 1 ρ 2 d W ^ u v ) ) d t + t i 1 t i v t d W ^ t S
Noting that t i 1 u < t t i , then exchange the order of integral that is produced by t i 1 t i r t d t and t i 1 t i v t d t . (A7) can be rewritten as
ln ( S t i S t i 1 ) = L 1 + L 2 + L 3 + L 4
where
L 1 = t i 1 t i ( η * + ( r t i 1 η * ) e h * ( t t i 1 ) 1 2 ( θ + ( v t i 1 θ ) e k ( t t i 1 ) ) ) d t = ( η * 1 2 θ ) Δ t + 1 h * ( 1 e h * Δ t ) ( r t i 1 η * ) 1 2 k ( 1 e k Δ t ) ( v t i 1 θ )
L 2 = t i 1 t i ξ t i 1 t i e h * ( t u ) r u d W ^ u r d t = t i 1 t i ξ u t i e h * ( t u ) d t r u d W ^ u r = ξ h * t i 1 t i ( 1 e h * ( t i u ) ) r u d W ^ u r
L 3 = 1 2 t i 1 t i σ 1 ρ 2 t i 1 t i e k ( t u ) v u d W ^ u v d t = 1 2 t i 1 t i σ 1 ρ 2 u t i e k ( t u ) d t v u d W ^ u v = σ 1 ρ 2 2 k t i 1 t i ( 1 e k ( t i u ) ) v u d W ^ u v
L 4 = t i 1 t i [ ρ σ 2 k ( 1 e k ( t i u ) ) + 1 ] v t d W ^ t S
Hence, the conditional expectations can be written in the following form:
E t i 1 [ ln ( S t i S t i 1 ) ] = E t i 1 [ L 1 ] + E t i 1 [ L 2 ] + E t i 1 [ L 3 ] + E t i 1 [ L 4 ]
where
E t i 1 T [ L 1 ] = ( η * 1 2 θ ) Δ t + 1 h * ( 1 e h * Δ t ) ( r t i 1 η * ) 1 2 k ( 1 e k Δ t ) ( v t i 1 θ )
E t i 1 [ L 2 ] = E t i 1 [ L 3 ] = E t i 1 [ L 4 ] = 0
Similarly, the conditional variance can be represented as:
V a r t i 1 [ ln ( S t i S t i 1 ) ] = V a r t i 1 [ L 1 ] + V a r t i 1 [ L 2 ] + V a r t i 1 [ L 3 ] + V a r t i 1 [ L 4 ]
where
V a r t i 1 [ L 1 ] = 0
V a r t i 1 [ L 2 ] = E t i 1 T [ ( ξ h * t i 1 t i ( 1 e h * ( t i u ) ) r u d W ^ u r ) 2 ] E t i 1 T [ ( ξ h * t i 1 t i ( 1 e h * ( t i u ) ) r u d W ^ u r ] 2 = ξ 2 ( h * ) 2 t i 1 t i ( 1 e h * ( t i u ) ) 2 E t i 1 T ( r u ) d u = ξ 2 ( h * ) 2 η * Δ t 2 ξ 2 ( h * ) 3 η * ( 1 e h * Δ t ) + ξ 2 2 ( h * ) 3 η * ( 1 e 2 h * Δ t ) + 2 ξ 2 ( h * ) 2 Δ t e h * Δ t + ξ 2 ( h * ) 3 ( 1 e 2 h * Δ t ) ( r t i 1 η * )
V a r t i 1 [ L 3 ] = E t i 1 T ( σ 1 ρ 2 2 k t i 1 t i ( 1 e k ( t i u ) ) v u d W ^ u v ) 2 E t i 1 T [ 1 2 σ k t i 1 t i ( 1 e k ( t i u ) ) v u d W ^ u v ] 2 = σ 2 ( 1 ρ 2 ) 4 k 2 t i 1 t i ( 1 e k ( t i u ) ) 2 E t i 1 T ( v u ) d u = σ 2 ( 1 ρ 2 ) { 1 4 k 2 θ Δ t 1 2 k 3 θ ( 1 e k Δ t ) + 1 8 k 3 θ ( 1 e 2 k Δ t ) + 1 2 k 2 Δ t e k Δ t + 1 4 k 3 ( 1 e 2 k Δ t ) ( v t i 1 θ ) }
V a r t i 1 [ L 4 ] = E t i 1 T ( t i 1 t i [ ρ σ 2 k ( 1 e k ( t i u ) ) + 1 v t d W ^ t S ) 2 ] ( E t i 1 T [ t i 1 t i [ ρ σ 2 k ( 1 e k ( t i u ) ) + 1 ] v t d W ^ t S ] ) 2 = t i 1 t i [ ρ σ 2 k ( 1 e k ( t i u ) ) + 1 ] 2 E t i 1 T ( v u ) d u = ( θ θ ρ σ k + σ 2 ρ 2 θ 4 k 2 ) Δ t + ( θ ρ σ k 2 σ 2 ρ 2 θ 2 k 3 ) ( 1 e k Δ t ) + θ σ 2 ρ 2 8 k 3 ( 1 e 2 k Δ t ) σ 2 ρ 2 2 k 2 e k Δ t Δ t + ( σ 2 ρ 2 4 k 3 + 1 k σ ρ k 2 ) ( 1 e k Δ t ) + σ 2 ρ 2 4 k 3 ( e k Δ t e 2 k Δ t ) + ρ σ k Δ t e k Δ t ( v t i 1 θ )
The value of C o v ( L 2 , L 3 ) , C o v ( L 3 , L 4 ) , C o v ( L 2 , L 3 ) are equal to 0 since W ^ t S , W ^ t v , W ^ t r are three independent Brownian motions
Finally, organizing the above results as (20) and (21), Theorem 2 is proved. ☐

Appendix C

Proof of Theorem 3. 
Taking the conditional expectation under the information set F 0 on both sides of (23) simultaneously
E 0 T [ ( r t i 1 η * ) ] = E 0 [ ( r 0 η * ) e h * ( t i 1 0 ) + ξ 0 t i 1 e h * ( t i 1 u ) r u d W ^ u r ] = ( r 0 η * ) e h * t i 1
E 0 T [ ( v t i 1 θ ) ] = E 0 [ ( v 0 θ ) e k ( t i 1 0 ) + σ 0 t i 1 e k ( t i 1 u ) v u ( ρ d W ^ u S + 1 ρ 2 d W ^ u v ) ] = ( v 0 θ ) e k t i 1
E 0 T [ ( r t i 1 η * ) 2 ] = E 0 [ ( r 0 η * ) 2 e 2 h * ( t i 1 0 ) + ξ 2 0 t i 1 e 2 h * ( t i 1 u ) ( η * + ( r 0 η * ) e h * u ) d u ] = e 2 h * t i 1 ( r 0 η * ) 2 + η * ξ 2 2 h * ( 1 e 2 h * t i 1 ) + ξ 2 h * ( e h * t i 1 e 2 h * t i 1 ) ( r 0 η * )
E 0 T [ ( v t i 1 θ ) 2 ] = E 0 T [ ( v 0 θ ) 2 e 2 k ( t i 1 0 ) ] + E 0 T [ σ 2 0 t i 1 e 2 k ( t i 1 u ) v u ( ρ d W ^ u S + 1 ρ 2 d W ^ u v ) 2 ] = e 2 k t i 1 ( v 0 θ ) 2 + θ σ 2 2 k ( 1 e 2 k t i 1 ) + σ 2 k ( e k t i 1 e 2 k t i 1 ) ( v 0 θ )
E 0 T [ ( r t i 1 η * ) ( v t i 1 θ ) ] = e h * t i 1 e k t i 1 ( r 0 η * ) ( v 0 θ )
Finally, plug the above results into (23), then (24) can be obtained. ☐

References

  1. Carr, P.; Lee, R. Volatility derivatives. Annu. Rev. Financ. Econ. 2009, 1, 319–339. [Google Scholar] [CrossRef] [Green Version]
  2. Grünbichler, A.; Longstaff, F.A. Valuing futures and options on volatility. J. Bank. Financ. 1996, 20, 985–1001. [Google Scholar] [CrossRef]
  3. Howison, S.; Rafailidis, A.; Rasmussen, H. On the pricing and hedging of volatility derivatives. Appl. Math. Financ. 2004, 11, 317–346. [Google Scholar] [CrossRef] [Green Version]
  4. Heston, S.; Nandi, S. Derivatives on Volatility: Some Simple Solutions Based on Observables. SSRN Electron. J. 2000. [Google Scholar] [CrossRef] [Green Version]
  5. 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]
  6. Swishchuk, A. Modeling and pricing of variance swaps for stochastic volatilities with delay. WILMOTT Mag. 2005, 19, 63–73. [Google Scholar]
  7. Swishchuk, A. Modeling and pricing of variance swaps for multi-factor stochastic volatilities with delay. Can. Appl. Math. Q. 2006, 14, 439–467. [Google Scholar]
  8. Habtemicael, S.; SenGupta, I. Pricing variance and volatility swaps for Barndorff-Nielsen and Shephard process driven financial markets. Int. J. Financ. Eng. 2016, 3, 1650027. [Google Scholar] [CrossRef]
  9. Issaka, A.; SenGupta, I. Analysis of variance based instruments for Ornstein-Uhlenbeck type models: Swap and price index. Ann. Financ. 2017, 13, 401–434. [Google Scholar] [CrossRef]
  10. SenGupta, I.; Wilson, W.; Nganje, W. Barndorff-Nielsen and Shephard model: Oil hedging with variance swap and option. Math. Financ. Econ. 2019, 13, 209–226. [Google Scholar] [CrossRef]
  11. Jarrow, R.; Kchia, Y.; Larsson, M.; Protter, P. Discretely sampled variance and volatility swaps versus their continuous approximations. Financ. Stoch. 2013, 17, 305–324. [Google Scholar] [CrossRef] [Green Version]
  12. Farnoosh, R.; Sobhani, A.; Rezazadeh, H.; Beheshti, M.H. Numerical method for discrete double barrier option pricing with time-dependent parameters. Comput. Math. Appl. 2015, 70, 2006–2013. [Google Scholar] [CrossRef]
  13. Farnoosh, R.; Sobhani, A.; Beheshti, M.H. Efficient and fast numerical method for pricing discrete double barrier option by projection method. Comput. Math. Appl. 2017, 73, 1539–1545. [Google Scholar] [CrossRef]
  14. Lian, G.; Zhu, S.P.; Elliott, R.J.; Cui, Z. Semi-analytical valuation for discrete barrier options under time-dependent lévy processes. J. Bank. Financ. 2017, 75, 167–183. [Google Scholar] [CrossRef]
  15. Little, T.; Pant, V. A finite-difference method for the valuation of variance swaps. J. Comput. Financ. 2001, 5, 81–101. [Google Scholar] [CrossRef]
  16. Heston, S.L. A closed-form solution for options with stochastic volatility with applications to bond and currency options. Rev. Financ. Stud. 1993, 6, 327–343. [Google Scholar] [CrossRef] [Green Version]
  17. Michael, A. Kouritzin, Explicit Heston solutions and stochastic approximation for path-dependent option pricing. Int. J. Theor. Appl. Financ. 2018, 21, 1850006. [Google Scholar]
  18. Zhu, S.P.; Lian, G.H. A closed-form exact solution for pricing variance swaps with stochastic volatility. Math. Financ. 2011, 21, 233–256. [Google Scholar]
  19. Zhu, S.P.; Lian, G.H. Analytically pricing volatility swaps under stochastic volatility. J. Comput. Appl. Math. 2015, 288, 332–340. [Google Scholar] [CrossRef]
  20. Zhang, L.W. A closed-form pricing formula for variance swaps with mean-reverting Gaussian volatility. ANZIAM J. 2014, 55, 362–382. [Google Scholar]
  21. 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]
  22. Liu, W.Y.; Zhu, S.P. Pricing variance swaps under the Hawkes jump-diffusion process. J. Futur. Mark. 2019, 39, 635–655. [Google Scholar] [CrossRef]
  23. Kim, J.H.; Yoon, J.H.; Yu, S.H. Multiscale stochastic volatility with the Hull-CWhite rate of interest. J. Futur. Mark. 2014, 34, 819–837. [Google Scholar] [CrossRef]
  24. Cao, J.; Lian, G.; Roslan, T.R.N. Pricing variance swaps under stochastic volatility and stochastic interest rate. Appl. Math. Comput. 2016, 277, 72–81. [Google Scholar] [CrossRef] [Green Version]
  25. Cox, J.C.; Ingersoll, J.E., Jr.; Ross, S.A. A theory of the term structure of interest rates. Econometrica 1985, 53, 385–407. [Google Scholar] [CrossRef]
  26. Cao, J.; Roslan, T.R.N.; Zhang, W. The valuation of variance swaps under stochastic volatility, stodwastic interest rate and full correlation structure. Korean Math. Soc. 2020, 57, 1167–1186. [Google Scholar]
  27. Han, Y.; Zhao, L. A closed-form pricing formula for variance swaps under MRG-Vasicek model. Comput. Appl. Math. 2019, 38, 142. [Google Scholar] [CrossRef]
  28. Xu, D.X.; Yang, B.Z.; Kang, J.H.; Huang, N.J. Variance and volatility swaps valuations with the stochastic liquidity risk. Phys. A Stat. Mech. Its Appl. 2021, 566, 125679. [Google Scholar] [CrossRef]
  29. He, X.J.; Zhu, S.P. How should a local regime-switching model be calibrated? J. Econom. Dyn. Control 2017, 78, 149–163. [Google Scholar] [CrossRef] [Green Version]
  30. He, X.J.; Zhu, S.P. On full calibration of hybrid local volatility and regime-switching models. J. Futur. Mark. 2018, 38, 586–606. [Google Scholar] [CrossRef]
  31. Zhu, S.P.; Lian, G.H. On the valuation of variance swaps with stochastic volatility. Appl. Math. Comput. 2012, 219, 1654–1669. [Google Scholar] [CrossRef]
  32. Bernard, C.; Cui, Z. Prices and asymptotics for discrete variance swaps. Appl. Math. Financ. 2014, 21, 140–173. [Google Scholar] [CrossRef] [Green Version]
  33. Zhu, S.P.; Lian, G.H. Pricing forward-start variance swaps with stochastic volatility. Appl. Math. Comput. 2015, 250, 920–933. [Google Scholar] [CrossRef]
  34. Cao, J.P.; Fang, Y.B. An analytical approach for variance swaps with an Ornstein-Uhlenbeck process. Aust. N. Z. Ind. Appl. Math. J. 2017, 59, 83–102. [Google Scholar]
  35. Sanae, R. Analytically pricing variance swaps in commodity derivative markets under stochastic convenience yields. Commun. Math. Sci. 2021, 19, 111–146. [Google Scholar]
  36. Larissa, B.; Maran, R.M.; Ioan, B.; Anca, N.; Mircea-Iosif, R.; Horia, T.; Gheorghe, F.; Ema Speranta, M.; Dan, M.I. Adjusted Net Savings of CEE and Baltic Nations in the Context of Sustainable Economic Growth: A Panel Data Analysis. Risk Financ. Manag. 2020, 13, 234. [Google Scholar] [CrossRef]
  37. Brigo, D.; Mercurio, F. Interest Rate Models-Theory and Practice: With Smile, Inflation and Credit; Springer Science and Business Media: Berlin/Heidelberg, Germany, 2007. [Google Scholar]
Figure 1. A comparison of fair strike prices based on the closed-solution, Monte Carlo simulations, the first-order approximation and the continuous model.
Figure 1. A comparison of fair strike prices based on the closed-solution, Monte Carlo simulations, the first-order approximation and the continuous model.
Mathematics 10 00005 g001
Figure 2. The values of variance swaps with different η in our hybrid model.
Figure 2. The values of variance swaps with different η in our hybrid model.
Mathematics 10 00005 g002
Table 1. The numerical results of Monte Carlo simulations, discrete sampling, first-order approximation and continuous sampling.
Table 1. The numerical results of Monte Carlo simulations, discrete sampling, first-order approximation and continuous sampling.
Sampling FrequencyMonte CarloDiscreteFirst OrderContinuous
Quarterly (N = 4)318.11318.13296.82236.31
Monthly (N = 12)274.24274.02258.91236.31
Fortnightly(N = 16)265.33265.55252.53236.31
Weekly (N = 52)247.43247.48241.80236.31
Daily (N = 252)240.14240.12239.45236.31
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mao, C.; Liu, G.; Wang, Y. A Closed-Form Pricing Formula for Log-Return Variance Swaps under Stochastic Volatility and Stochastic Interest Rate. Mathematics 2022, 10, 5. https://doi.org/10.3390/math10010005

AMA Style

Mao C, Liu G, Wang Y. A Closed-Form Pricing Formula for Log-Return Variance Swaps under Stochastic Volatility and Stochastic Interest Rate. Mathematics. 2022; 10(1):5. https://doi.org/10.3390/math10010005

Chicago/Turabian Style

Mao, Chen, Guanqi Liu, and Yuwen Wang. 2022. "A Closed-Form Pricing Formula for Log-Return Variance Swaps under Stochastic Volatility and Stochastic Interest Rate" Mathematics 10, no. 1: 5. https://doi.org/10.3390/math10010005

APA Style

Mao, C., Liu, G., & Wang, Y. (2022). A Closed-Form Pricing Formula for Log-Return Variance Swaps under Stochastic Volatility and Stochastic Interest Rate. Mathematics, 10(1), 5. https://doi.org/10.3390/math10010005

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