Abstract

This paper presents an extension of double Heston stochastic volatility model by incorporating stochastic interest rates and derives explicit solutions for the prices of the continuously monitored fixed and floating strike geometric Asian options. The discounted joint characteristic function of the log-asset price and its log-geometric mean value is computed by using the change of numeraire and the Fourier inversion transform technique. We also provide efficient approximated approach and analyze several effects on option prices under the proposed model. Numerical examples show that both stochastic volatility and stochastic interest rate have a significant impact on option values, particularly on the values of longer term options. The proposed model is suitable for modeling the longer time real-market changes and managing the credit risks.

1. Introduction

Asian option is a special type of option contract in which the payoff depends on the average of the underlying asset price over some predetermined time interval. The averaging feature allows Asian options to reduce the volatility inherent in the option. There are some advantages to trading Asian options in a financial market. One is that these decrease the risk of market manipulation of the financial derivative at expiry. Another is that Asian options have lower relative charge than European or American options. In general, the average considered can be a arithmetic or geometric one and it can be calculated either discretely, for which the average is taken over the underlying asset prices at discrete monitoring time points, or continuously, for which the average is calculated via the integration of the underlying asset price over the monitoring time period. Asian options can be differentiated into two main classes according to their payoff: fixed strike price options (sometimes called “average price”) and floating strike price options (sometimes called “average strike”). All these details are specified by the contracts stipulated by two counterparts, as Asian options are traded actively on the OTC market among investors or traders for hedging the average price of a commodity. For a brief introduction to the development of Asian options, see Boyle and Boyle [1].

As the probability distribution of the average prices of the underlying asset generally does not have a simple analytical expression, it is difficult to obtain the analytical pricing formula for Asian option. Since the best-known closed-form pricing formula for the European vanilla option derived by Black and Scholes [2]), many researchers have devoted themselves to developing the Asian options pricing based on the Black-Scholes assumptions; see, e.g., Kemna and Vorst (1990), Turnbull and Wakeman [3], Ritchken et al. [4], Geman and Yor [5], Rogers and Shi [6], Boyle et al. [7], Angus [8], Linetsky [9], Cui et al. [10], and the references therein. For a recent review, one can refer to Fusai and Roncoroni [11] and Sun et al. [12].

In practice, the Black-Scholes assumptions are hardly satisfied, especially the constant volatility and constant interest rate hypothesis. As the empirical behaviors of the implied volatility smile and heavy tailed in the distribution of log-returns are commonly observed in financial markets. For this reason, stochastic volatility (hereafter SV) models have been proposed in finance (see Hull and White [13], Stein and Stein [14], Heston [15], and others). These models have been applied to value the Asian options (see, e.g., Wong and Cheung [16], Hubalek and Sgarra [17], Kim and Wee [18], and Shi and Yang [19]). In addition, interest rates are stochastic and stock returns are negatively correlated with interest rate changes, which have been examined in previous research.

Although these models mentioned above are able to account for the empirical behaviors, they are still based on a single-factor for volatility dynamics that is inconsistent with the long range memory characteristic of the volatility corrections and the stiff volatility skews. See Alizadeh et al. [20], Fiorentini et al. [21], Chernov et al. [22], Gourieroux [23], Christoffersen et al. [24], Romo [25], and Nagashima et al. [26] for the empirical results. To address this issue, multifactor SV models have recently generated attention in the option pricing literature. For instance, Duffie et al. [27] proposed multifactor affine stochastic volatility models. Based upon the Black-Scholes framework, Fouque et al. [28] introduced a multiscale SV model, in which the volatility processes are driven by two mean-reverting diffusion processes. Gourieroux [23] proposed a multivariate model in which the volatility-covolatility matrix follows a Wishart process.

On the basis of the findings of Christoffersen et al. [24], a double Heston (dbH) model, which consists of two independent variance processes, has recently been reported better than the plain Heston [15] model in the performances of hedging (see Sun [29]) and has also been applied to arithmetic Asian option under discrete monitoring (Mehrdoust and Saber [30]) and forward starting option (Zhang and Sun [31]). However, its extension to continuously monitored geometric Asian option is yet to be considered. On the other hand, many of Asian options often have long-dated maturities since they are used as part of the structured notes which has a long maturity. The movement of interest rates becomes an issue in such cases and constant interest rate assumption should be replaced by an appropriate dynamic interest rate model. Several results are available on the Asian option in the stochastic interest rate framework; see, e.g., Nielsen and Sandmann [32, 33], Zhang et al. [34], Donnelly et al. [35], and He and Zhu [36]. In the above stochastic interest rate framework, the short-term interest rate is assumed to follow a specific parametric one-factor model (see, e.g., Cox et al. [37], Hull and White [13], and Vasicek [38]), which tends to oversimplify the true behavior of interest rate movement. However, empirical tests reported in Lonstaff and Schwartz [39] and Pearson and Sun [40] show that the term structure for the interest rate should involve several sources of uncertainty, and introducing additional state variables (such as the rate of inflation, GDP, etc.) significantly improves the fit.

In this paper, we study the pricing of the continuously monitored geometric Asian options under dbH stochastic volatility model with stochastic interest rate framework (hereafter, dbH-SI model). The contribution of the present paper is twofold. Firstly, this paper extends the dbH model by introducing stochastic interest rate, which is assumed to follow two-factor model with two state variables. Secondly, this paper provides a semiexplicit valuation formula for the geometric Asian options with fixed or floating strike price, which is extremely useful also for the arithmetic average option valuation via Monte Carlo methods with control variables.

The rest of the paper is organized as follows. Section 2 develops the underlying pricing model and describes the geometric Asian option. Section 3 derives the joint characteristic function of a log-return of the underlying asset and its geometric average. Section 4 obtains the analytic expressions for the prices of the fixed strike geometric Asian call option and the floating strike Asian call option under continuous monitoring. Section 5 provides some numerical examples for the proposed approach. Section 6 concludes the paper.

2. Model Formulation

We consider an arbitrage-free, frictionless financial market where only riskless asset and risky asset are traded continuously up to a fixed horizon date . Let be a complete probability space equipped with a filtration satisfying the usual conditions, where is a risk-neutral probability measure. Suppose and are all standard Brownian motions defined on the probability space, and the filtration is generated by these Brownian motions. Moreover, , , and any other Brownian motions are pairwisely independent. Assume that the asset price process , without paying any dividend, satisfies the following stochastic differential equation under :where are all nonnegative constants, which represent the mean-reverting rates, long-term mean levels, and volatilities of variance processes , respectively. We suppose that . The instantaneous interest rate, , is assumed to be a linear combination of and , i.e., , which designates the interest rate as an affine function of two-factor economic variables and and offers the analytic tractability (see Duffie et al. [27]).

In financial market, there are four types of European style continuously monitoring geometric Asian options: fixed strike geometric Asian calls, fixed strike geometric Asian puts, floating strike geometric Asian calls, and floating strike geometric Asian puts. The payoffs at the expiration date for these options are as follows: where is a fixed strike price and is the geometric average of the underlying asset price until time ; i.e., . In the following, we consider only the pricing problem of the geometric Asian call options (hereafter, GAC), while the put options can be dealt with similarly.

For the instantaneous interest rate , one can express the price at time of a zero-coupon bond with maturate as follows (see Cox et al. [37]):where

3. The Joint Characteristic Function

Given the dynamic of the underlying asset price, it is possible to obtain the discounted joint characteristic function for the log-asset value and the log-geometric mean value of the asset price over a certain time period.

Let be the discounted joint characteristic function of two-dimensional random variable, , conditioned on under , where and is the conditional expectation under for . Denote .

Proposition 1. Suppose that , and follow the dynamics in (1). If and , then andwhere

Proof. (i) We first prove that the integrability condition guarantees the existence of the cumulant function in . If and , then where is the forward measure given by the Radon-Nikodym derivative: , and is given above in (3). In the case of , it is triviality.
(ii) In order to determine (5), we start from model (1) and develop with the Brownian motions and expressed as and , respectively, where are 4-dimensional Brownian motion defined on the probability space. For the process, , we haveOn the other hand, we have Using the fact for , thenLet be the field generated by and . By (8) and (11), for , we havewhereSincesubstituting (14) into (12) and applying the Markov property of , lead to (5), which completes the proof.

From Proposition 1, it is clear that we need to search for an exact formula for the discounted joint characteristic function of and and the three different integrals of appearing in (5). We use the same approach introduced by Kim and Wee [18] to obtain an explicit formula for the discounted joint characteristic function . Therefore, two series of functions are introduced as follows. Define and aswhere , are functions of , , , and defined as for .

For , denote

Proposition 2. (1) The argument of : for every . In particular, is continuous on , and for , and are all real numbers in for .
(2) For , we haveThe proof of Proposition 2 is similar to that of Kim and Wee [18].
Using Proposition 2 to Proposition 1 leads to the explicit expression of in Proposition 3. To describe the simplicity of the result, we need to introduce new functions. Define and : aswith .

Proposition 3. (1) If , then and is continuous on . In particular, if and are real numbers for .
(2) For with defined as above, thenHere , for .

Proof. Assume that , using the definitions of and for and . Note that for every and . Therefore, (19) remains valid for any . The time homogeneous Markov property of implies that (22) holds when Propositions 1 and 2 are satisfied.
Now Substituting the expressions of into (15) and (16), and are expressed as follows, i.e., for ,where for .

4. Pricing Geometric Asian Option

Once the discounted joint characteristic function is found, the European continuously monitored geometric Asian option can be valued using numeraire change technique and the inverse Fourier transform approach which are applied in many research works (see, e.g., Geman et al. [41] and Deng [42]). This section derives the pricing formulas for the continuously monitored fixed and floating strike geometric Asian call options.

Theorem 4. Suppose that , and follow the dynamics in (1), then the price at time of the continuously monitored fixed strike geometric Asian call option with maturity and the strike price is given bywhere is given above in (3) and where is the imaginary unit ().

Proof. Sincewhere is defined by the following Radon-Nikodym derivative: and is the forward measure given above, it is well known that the probability distribution functions can be calculated by using the Fourier inversion transform, and then the above two probabilities and in (29) are given bywhere From (29), (31), and (32), we can obtain the required Theorem 4.

Theorem 5. Suppose that , and follow the dynamics in (1), then the price at time of the continuously monitored floating strike geometric Asian call option with maturity is given bywhere

Proof. Sincewhere is defined by the following Radon-Nikodym derivative:under the two probability measures and , the conditional characteristic functions of are given byTherefore Plugging (41) and (42) in (37) yields (35).

In addition to option prices, one can compute derivatives to hedge against changes in the underlying asset price and volatilities and . We omit it due to its triviality.

5. Numerical Examples

In this section, we use the dbH-SI model to analyze the valuation of the continuously monitored fixed strike geometric Asian call option using some numerical examples. To implement the analytic formula given in (27) numerically, we first determine the number of terms taken for the computation of the infinite series expansions for and for in (24) and (25). Second, we investigate the accuracy and efficiency of the approximated analytic formula given in (27). In the end, we compare the option prices varying by and under different models including the dbH-SI, H-SI (i.e., Heston stochastic volatility model with stochastic interest rate), dbH, and Heston models. Furthermore, we use this proposed model (1), i.e., the dbH-SI model, to examine the impacts of the model parameters on option prices. Here we implement the integral formulas given in (27) without truncation, modification, and approximation in Mathematica and Matlab software. Numerical integration was performed using the “NIntegrate()” or “quadgk()” commands which can handle infinite ranges, oscillatory integrands, and singularities in their default version, which employs automatic adaptive integration.

Model (1) parameters are set as follows: , , , , , , , , , , , and . The results are shown in Table 1.

Table 1 investigates the influence of the number of terms taken in (24) and (25). It is shown that the numerical values tend to be quite stable as the number of terms increases. In particular, the option prices stay unchanged after taking and terms, at most. In addition, Table 1 provides option prices with various maturities , strike prices , and the required CPU times (in seconds). From Table 1, we can see that some characteristics of the option are similar to the plain European call option. For example, the values of the option are decreasing functions of the strike price. In addition, under the same parameter values setting in dbH-SI model, longer maturating date will result in higher option values.

It would be interesting to see the performance of our approximated analytic formula approach implemented to the fixed strike GAC option under dbH-SI model. Our numerical example uses the number of terms taken in (24) and (25); other parameter values are similar to those settings of Table 1. We compute the 0.5-,1.5-, and 3-year maturity GAC options by our approximated analytic formula approach and compare the results with Monte Carlo simulation with 10000 sample paths and 100 points for time axis. The numerical results are displayed in Table 2.

Table 2 shows that the approximated analytic formula approach is considerably faster than those of the Monte Carlo simulation (MC). For a given set of parameters, the approximated analytic formula approach calculates the option prices for 5 different strikes and 3 different maturates in approximately 5 seconds. The Monte Carlo simulation takes around 110 seconds for each option price. Table 2 also compares their pricing accuracy. It can be seen that the relative error (RE) in prices is less than 0.4% for all cases. If we regard the Monte Carlo price as the benchmark, then this numerical example confirms that the approximated analytic formula approach is accurate and efficient.

After examining accuracy and efficiency. We shall turn to investigate the comparison of the option prices under different models. Figure 1 compares the dbH-SI model with Heston, dbH, and H-SI models in pricing the fixed strike GAC option against the initial asset price and maturating date with the strike price . From Figure 1, it is interesting to notice that the option prices of the fixed strike GAC option calculated from the dbH-SI model are much higher than those of the dbH model, H-SI model, and Heston model. One possible reason for this phenomena is that the volatilities are always bigger for the dbH-SI model than those for other models. This implies that the proposed model, i.e., the dbH-SI model, has a more significant influence than the dbH model, H-SI model, and Heston model on option prices. We can also clearly observe that the bigger the value of (or ) is, the higher the GAC option price is.

In Figure 2, we use the dbH-SI model to examine the effects of the mean-reverting rates and volatilities of variance processes on the GAC option prices (see Figures 2(a)2(d)). We also examine the impacts of the correlation coefficients (see Figures 2(e) and 2(f)), which determine both the correlation between the underlying asset and its volatilities, and the correlation between the underlying asset and the short interest rates.

From Figures 2(a) and 2(b), we can observe that the option prices increase as the mean-reverting rate increases. Particularly, it is quite remarkable that the mean-reverting rates have a significant effect on the longer term option values. The option price decreases as the volatilities of variance processes increase (see Figures 2(c) and 2(d)). The correlation coefficients have several effects depending on the relation between the strike price and the initial asset price. A negative tends to produce higher values for ITM (in-the-money) options and lower values for OTM (out-of-the-money) options (see Figure 2(e)). A similar result holds with respect to the expiration date: a negative tends to produce higher value for the longer term options and lower values for the shorter options (see Figure 2(f)).

6. Conclusions

The proposed model (the dbH-SI model) incorporates several important features of the underlying asset returns variability. We derive the discounted joint characteristic function of the log-asset price and its log-geometric mean value over time period and obtain approximated analytic solutions to the continuously monitored fixed and floating strike geometric Asian call options using the change of numeraire technique and the Fourier inversion transform approach. Some numerical examples are provided to examine the effects of the proposed model, which reveals some additional features having a significant impact on option values, especially long-term options. The proposed model can be tested empirically by using the option price data from the option market.

Data Availability

All data used to support of the findings for this study are included within this article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The authors acknowledge the financial support provided by the Natural Sciences Foundation of China (11461008), the Guangxi Natural Science Foundation (2018JJA110001), and the Guangxi Graduate Education Innovation Program project (XYCSZ2018059).