Next Article in Journal
Realistic Choice of Annual Matrices Contracts the Range of λS Estimates
Next Article in Special Issue
Research on the Factors of Competition in the Green Procurement Processes: A Case Study for the Conditions of Romania Using PLS-SEM Methodology
Previous Article in Journal
The Riemann-Lebesgue Integral of Interval-Valued Multifunctions
Previous Article in Special Issue
Constructing a Novel Early Warning Algorithm for Global Budget Payments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pricing of Arithmetic Asian Options under Stochastic Volatility Dynamics: Overcoming the Risks of High-Frequency Trading

1
Department of Business Administration, Soochow University, Taipei 100006, Taiwan
2
Department of Financial Engineering and Actuarial Mathematics, Soochow University, Taipei 100006, Taiwan
3
Taiwan/Pervasive Artificial Intelligence Research (PAIR) Labs, Taipei 106031, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(12), 2251; https://doi.org/10.3390/math8122251
Submission received: 30 October 2020 / Revised: 16 December 2020 / Accepted: 18 December 2020 / Published: 20 December 2020
(This article belongs to the Special Issue Mathematical Analysis in Economics and Management)

Abstract

:
This research extended the model developed by Hull and White by integrating Taylor-series expansion into the model for deriving approximate analytical solutions for stochastic volatility forward-starting Asian options. Numerical experiments were performed to compare the proposed model with the Monte Carlo model over numerous simulations and demonstrated that the developed model has a pricing accuracy greater than 99%. Furthermore, the computation time was approximately 10−5 s for each simulation. The model’s outstanding computational performance demonstrates its capability to address the challenges of high-frequency trading.

1. Introduction

High-frequency trading (HFT) is a method of automated trading that involves using computer programs to rapidly perform large numbers of transactions. HFT is controversial and has received some harsh criticism [1,2,3]. The qualitative changes in the functioning principles and structure of the derivatives market have a substantial impact on processes associated with investment activities so much so that major financial markets are regulated by policy measures that aim to reduce the threat from HFT. This evolution has already led to changes in how information is disseminated to traders and has given rise to new forms of trading. Therefore, a set of faultless regulations and hedging tools is necessary to restrict the disturbance from HFT, avoid circuit breakers, and maintain global financial order.
G20 leaders have aimed to avoid disturbances in interest rates and the foreign exchange derivatives market by committing to increasing transparency in derivatives markets through the implementation of mandatory reporting of derivatives contracts [4,5,6]. To address regulators’ concerns and the potential threats of HFT, a forward-starting Asian option (FSAO) evaluation formula was developed in this study to cope with the potential instability of interest rates, equity, and commodity and foreign-exchange markets. The FSAO can be regarded as a special type of Asian option in which the arithmetic average price is calculated from the first trading date of the contract.
The advantages of the Asian option are manifold. First, it is popular in practice in the financial market because it has a lower cost than its counterpart, the plain vanilla option. Second, averaging underlying asset prices can prevent manipulation from the financial market [7,8,9]. Third, it is commonly used in interest rate and foreign exchange markets, providing an expiry price or strike price that remains more stable than plain vanilla over contract periods.
In academia, an FSAO is a path-dependent derivative; its payoff function depends on the average value of the underlying asset over a relatively short time period in which the option has reached partial maturity and is near its expiration date. A precise analytical solution does not exist. Thus, extensive research has investigated such options using numerical approaches. This paper briefly catalogues four main methods described in the literature: The Monte Carlo approach with variance reduction, approximation of the unknown distribution on the average price with a more tractable one, derivation of the exact representation of the Asian option price, and extensions of the Heston model [10].
The first approach, Monte Carlo simulation with variance reduction, is a benchmark method for option valuation. It has often been employed directly for option pricing or to assist in solving complex partial differential equations in option valuations. Well-known examples of the direct application of Monte Carlo simulation for option evaluation are in the studies by [11,12]. They numerically solved pricing partial differential equations, which inspired a series of subsequent studies, such as those by [13,14,15,16,17,18]. The main drawback of this evaluation approach is that it lacks computational efficiency.
In the second approach, the unknown distribution of the Asian option is approximated to obtain an approximate analytical option price formula. Well-known examples of researchers who have adopted this approach are [19,20,21,22,23,24,25,26,27,28]. Generally, the method involves applying Edgeworth expansion or a Taylor series approach to a log-normal reference distribution to approximate the distribution of the arithmetic average of the underlying price or strike price. For example, [29] used a moment-matched compound gamma distribution. These researchers developed analytic option price formulas, but the formulas lack reliable error estimates.
Introduced by [30], the third approach involves the use of the reciprocal of the average, for which the Laguerre series converges. Dufresne improved non–closed form solutions by using an advance mathematical technique. [31] adopted this approach to derive an analytical solution for complex option evaluations. However, this approach does not guarantee that the series will converge to the true price and does not consider the stochastic volatility effect.
First developed by [10], the Heston model has attracted considerable attention, with extensive research exploring its explicit (distributional and pathwise) solutions and exact Monte Carlo simulation. For example, [32] studied the exact simulation of stochastic volatility and other affine jump diffusion processes, and [33] explored explicit Heston solutions and stochastic approximation for path-dependent option pricing. Additionally, [34] recently developed exact option pricing methods for the popular 3/2 model. Finally, an exotic unified approach to pricing Asian options in stochastic jump diffusion models, which was extended to include alternative exotic contracts by [35].
Motivated by these empirical results, the present study extended the research of Hull and White [19,20,25,28] to include a Taylor series expansion technique for deriving the approximate analytical solution of an FSAO in the context of stochastic volatility pricing.
This study contributes to the literature by providing an efficient and accurate analytical solution for a stochastic-volatility FSAO (SVFSAO). Empirical evidence strongly supports the use of two well-separated timescales during evaluation.
The remainder of this paper is structured as follows. Section 2 outlines the problem being addressed, describes the derivation of an approximate analytic SVFSAO with Taylor expansion, and presents the solution for the optimal hedge ratio. Section 3 discusses the model’s pricing and hedging performance in comparison with the solutions of the Monte Carlo simulation, which served as the benchmark. Section 4 presents conclusions.

2. Materials and Methods

This paper first provides the approximate analytical solution for the SVFSAO (Because of the complexity of deriving the SVFSAO, the whole derivation process is not included in this paper. Readers who are interested in the aforementioned derivation may receive more details from the authors upon reasonable request.) and subsequently addresses the hedge ratio of the SVFSAO. Details regarding the derivation process of the approximate analytical solution for the hedge ratio are provided in the Appendix A.

2.1. Approximate Analytical Solution for SVFSAO Pricing

Following Hull and White [1], Tsao, Chang, and Lin [14], and Lin and Chang [28], a European-style FSAO with stochastic volatility written on an asset with maturity date T was investigated.
The price of the underlying asset, which was tradable and uncorrelated with aggregate consumption, was assumed to follow geometric Brownian motion in a risk-neutral world with stochastic volatility and was calculated as follows:
d S t S t = ( r d ) d t + σ t d Z t ,
d σ t 2 σ t 2 = μ d t + ξ d W t ,
where S t is the underlying asset price at time t, r is the annualized risk-free interest rate, d is the continuous yield, σ t is the instantaneous volatility of the underlying asset at time t, μ is the drift term, and ξ is the volatility term for the variance of the underlying asset. Both μ and ξ remain constant over time. In addition, d Z t and d W t are independent Wiener processes.
Following the settings of Tsao, Chang, and Lin [14], this paper assumes that the strike price K of option is the arithmetic average of underlying asset prices over the period [T − A, T]. It follows that the strike price K is
K = 1 A T A T S u d u
where T − A is the time instant after time zero of the option’s issuance. If T − A = 0, then this option becoming a plain vanilla Asian option. It is important to distinguish the following time windows: [0, T − A] and [T − A, T]. Before T − A, investors know nothing about the exercise price K. After time T − A information about K gradually accumulates.
  • Time Window before T − A
The price of a forward starting Asian call option at time t ( 0 < t T A ) is denoted as
C = e r ( T t ) E [ ( S T K ) + | F t ]
where E is the conditional expectation operator with respect to the risk-neutral probability measure and F t is a filtration.
From Tsao, Chang, and Lin [14], it is straightforward to have the analytic approximation formula for forward starting Asian call option as
C = S t e ( r d ) A d ( T t ) [ m 1 N ( m 1 ν 1 ) + ν 1 2 π e m 1 2 2 ν 1 ]
where
m 1 = ( r d σ 2 2 ) A 2 + A σ 2 4 + ( r d σ 2 2 ) 2 A 2 3
ν 1 = A σ 2 3 + A 2 σ 2 4 + 7 ( r d σ 2 2 ) 2 A 3 σ 2 15 + 3 ( r d σ 2 2 ) A 2 σ 2 4
  • Time Window After T − A
For the case of time t ( T A < t T ) , at time t, the investors knew the information of the underlying asset price between [ S T A , S t ] . The strike price can be split into two parts as follows,
K = M t + 1 A t T S u d u
where
M t = 1 A T A t S u d u
The analytic approximation formula for forward starting Asian call option is straightforward as
C = S t e r ( T t ) [ m 2 N ( m 2 ν 2 ) + ν 2 2 π e m 2 2 2 ν 2 ]
where
m 2 = 1 + ( r d σ 2 2 ) ( T t ) ( T t ) A ( r d σ 2 2 ) ( T t ) 2 2 A M t S t + ( T t ) σ 2 2 + ( T t ) 2 ( σ ( r d σ 2 2 ) 2 2 σ 2 4 A ) ( T t ) 3 σ ( r d σ 2 2 ) 2 6 A
ν 2 = ( T t ) σ 2 [ 1 ( T t ) A + ( T t ) 2 3 A 2 ] + σ 2 ( T t ) 2 [ 2 ( r d σ 2 2 ) + σ 2 2 ] + σ 2 ( T t ) 3 [ ( r d σ 2 2 ) 2 σ 2 3 A 5 ( r d σ 2 2 ) 3 A ] + σ 2 ( T t ) 4 [ σ 2 12 A 2 2 ( r d σ 2 2 ) 2 3 A + 5 ( r d σ 2 2 ) 12 A 2 ] + σ 2 ( T t ) 5 2 ( r d σ 2 2 ) 2 15 A 2
Let V ¯ be the mean variance over the time interval [0, T], from Hull and White [1]:
V ¯ = 1 T 0 T σ u 2 d u
The first three moments of V ¯ at μ = 0 can be derived as (Following the settings of Hull and White [19], Lin, and Chang [28], and to simplify the derivation of the analytic approximate solution for options, this paper only discusses the case in which μ = 0 . However, the case in which μ 0 and the moments of V ¯ with nonzero μ can be derived without considerable effort; see Hull and White [19] for further details.)
E ( V ¯ ) = V 0 ,
E ( V ¯ 2 ) = 2 ( e ξ 2 T ξ 2 T 1 ) ξ 4 T 2 V 0 2 ,   and
E ( V ¯ 3 ) = ( e 3 ξ 2 T 9 e ξ 2 T 6 ξ 2 T + 8 ) 3 ξ 6 T 3 V 0 3 .
The derivations of variance and skewness of V ¯ are straightforward as,
V ( V ¯ ) = 2 σ 0 4 ( e k k 1 ) k 2 σ 0 4   and
S k e w ( V ¯ ) = σ 0 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 + 6 k 3 ) 3 k 3 ] ,
where k = ξ 2 T .
Proposition 1.
The approximate analytical solution f ( S 0 , σ 0 2 ) for pricing the SVFSAO when t [ 0 ,     T A ,     T ] can be approximated using the following formula:
f ( S 0 , σ 0 2 ) = S t e ( r d ) A d ( T t ) [ m 1 N ( m 1 v 1 ) + v 1 2 π e m 1 2 2 v 1 ] + 1 2 S t e ( r d ) A d ( T t ) [ M 12 N ( m 1 v 1 ) + 1 2 π e m 1 2 2 v 1 M 11 A ¯ 1 1 32 π e m 1 2 2 v 1 v 1 3 2 V 11 2 + 1 8 π e m 1 2 2 v 1 v 1 1 2 V 11 B ¯ 1 + 1 8 π e m 1 2 2 v 1 v 1 1 2 V 12 ] × [ 2 σ 4 ( e k k 1 ) k 2 σ 4 ] + 1 6 S t e ( r d ) A d ( T t ) [ M 13 N ( m 1 v 1 ) + 1 2 π e m 1 2 2 v 1 α 11 + 1 32 π e m 1 2 2 v 1 α 12 + 1 8 π e m 1 2 2 v 1 ( α 13 + α 14 ) ] × σ 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 6 k 3 ) 3 k 3 ]
where
A ¯ 1 = V ¯ ( m 1 v 1 ) = M 11 v 1 1 2 1 2 m 1 v 1 3 2 V 11
B ¯ 1 = V ¯ ( m 1 2 2 v 1 ) = m 1 v 1 1 M 11 + 1 2 m 1 2 v 1 2 V 11
α 11 = 2 M 12 A ¯ 1 + M 11 A ¯ 1 B ¯ 1 + M 11 ( v 1 1 2 M 12 v 1 3 2 M 11 V 11 1 2 m 1 v 1 3 2 V 12 + 3 4 m 1 v 1 5 2 V 11 2 )
α 12 = 3 2 v 1 5 2 V 11 3 v 1 3 2 V 11 2 B ¯ 1 2 v 1 3 2 V 11 V 12
α 13 = 1 2 v 1 3 2 V 11 V 12 + v 1 1 2 V 12 B ¯ 1 + v 1 1 2 V 13
α 14 = 1 2 v 1 3 2 V 11 2 B ¯ 1 + v 1 1 2 V 11 B ¯ 1 2 + v 1 1 2 V 12 B ¯ 1 + v 1 1 2 V 11 ( v 1 1 M 11 2 + 2 m 1 v 1 2 M 11 V 11 m 1 v 1 1 M 12 m 1 2 v 1 3 V 11 2 + 1 2 m 1 2 v 1 2 V 12 )
M 11 = m 1 V ¯ = 1 3 A 2 ( r d ) + 1 6 A 2 V ¯
M 12 = 2 m 1 V ¯ 2 = 1 6 A 2
M 13 = 3 m 1 V ¯ 3 = 0
V 11 = v 1 V ¯ = [ 1 3 A + 3 4 A 2 ( r d ) + 7 15 A 3 ( r d ) 2 ] + [ 1 4 A 2 14 15 A 3 ( r d ) ] V ¯ + 7 20 A 3 V ¯ 2
V 12 = 2 v 1 V ¯ 2 = [ 1 4 A 2 14 15 A 3 ( r d ) ] + 7 10 A 3 V ¯
V 13 = 3 v 1 V ¯ 3 = 7 10 A 3
When t [ 0 , T A , T ] , then
f ( S 0 , σ 0 2 ) = S t e r ( T t ) [ m 2 N ( m 2 v 2 ) + v 2 2 π e m 2 2 2 v 2 ] + 1 2 S t e r ( T t ) [ M 22 N ( m 2 v 2 ) + 1 2 π e m 2 2 2 v 2 M 21 A ¯ 2 1 32 π e m 2 2 2 v 2 v 2 3 2 V 21 2 + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 21 B ¯ 2 + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 22 ] × [ 2 σ 4 ( e k k 1 ) k 2 σ 4 ] + 1 6 S t e r ( T t ) [ M 23 N ( m 2 v 2 ) + 1 2 π e m 2 2 2 v 2 α 21 + 1 32 π e m 2 2 2 v 2 α 22 + 1 8 π e m 2 2 2 v 2 ( α 23 + α 24 ) ] × σ 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 6 k 3 ) 3 k 3 ]
where
A ¯ 2 = V ¯ ( m 2 v 2 ) = M 21 v 2 1 2 1 2 m 2 v 2 3 2 V 21
B ¯ 2 = V ¯ ( m 2 2 2 v 2 ) = m 2 v 2 1 M 21 + 1 2 m 2 2 v 2 2 V 21
α 21 = 2 M 22 A ¯ 2 + M 21 A ¯ 2 B ¯ 2 + M 21 ( v 2 1 2 M 22 v 2 3 2 M 21 V 21 1 2 m 2 v 2 3 2 V 22 + 3 4 m 2 v 2 5 2 V 21 2 )
α 22 = 3 2 v 2 5 2 V 21 3 v 2 3 2 V 21 2 B ¯ 2 2 v 2 3 2 V 21 V 22
α 23 = 1 2 v 2 3 2 V 21 V 22 + v 2 1 2 V 22 B ¯ 2 + v 2 1 2 V 23
α 24 = 1 2 v 2 3 2 V 21 2 B ¯ 2 + v 2 1 2 V 21 B ¯ 2 2 + v 2 1 2 V 22 B ¯ 2 + v 2 1 2 V 21 ( v 2 1 M 21 2 + 2 m 2 v 2 2 M 21 V 21 m 2 v 2 1 M 22 m 2 2 v 2 3 V 21 2 + 1 2 m 2 2 v 2 2 V 22 )
M 21 = m 2 V ¯ = [ 1 2 ( T t ) 2 ( r d ) + 1 6 A ( T t ) 3 ( r d ) ] + [ 1 4 ( T t ) 2 1 12 A ( T t ) 3 ] V ¯
M 22 = 2 m 2 V ¯ 2 = [ 1 4 ( T t ) 2 1 12 A ( T t ) 3 ]
M 23 = 3 m 2 V ¯ 3 = 0
V 21 = v 2 V ¯ = [ 1 1 A ( T t ) + 2 ( T t ) ( r d ) + 1 3 A 2 ( T t ) 2 5 3 A ( T t ) 2 ( r d )   + ( T t ) 2 ( r d ) 2 + 5 12 A 2 ( T t ) 3 ( r d ) 2 3 A ( T t ) 3 ( r d ) 2   + 2 15 A 2 ( T t ) 4 ( r d ) 2 ] ( T t ) + [ 1 + 1 A ( T t ) 2 ( T t ) ( r d )   1 4 A 2 ( T t ) 2 + 4 3 A ( T t ) 2 ( r d ) 4 15 A 2 ( T t ) 3 ( r d ) ] ( T t ) 2 V ¯   + [ 3 4 1 2 A ( T t ) + 1 10 A 2 ( T t ) 2 ] ( T t ) 3 V ¯ 2
V 22 = 2 v 2 V ¯ 2 = [ 1 + 1 A ( T t ) 2 ( T t ) ( r d ) 1 4 A 2 ( T t ) 2 + 4 3 A ( T t ) 2 ( r d ) 4 15 A 2 ( T t ) 3 ( r d ) ] ( T t ) 2 + [ 3 2 1 A ( T t ) + 1 5 A 2 ( T t ) 2 ] ( T t ) 3 V ¯
V 23 = 3 v 2 V ¯ 3 = [ 3 2 1 A ( T t ) + 1 5 A 2 ( T t ) 2 ] ( T t ) 3

2.2. Approximate Analytical Solution for Hedging the SVFSAO

After the analytical solution for pricing the SVFSAO is derived, the hedge ratio ( Δ = f ( S 0 , σ 0 2 ) S 0 ) can be derived directly through differentiation of the option price f ( S 0 , σ 0 2 ) with respect to the asset price S 0 at the present time:
D e l t a = Δ = f ( S 0 , σ 0 2 ) S 0
Proposition 2.
The approximate analytical solution for hedging the SVFSAO at t [ 0 , T A , T ] can be approximated using the following formula:
D e l t a = Δ = f ( S 0 , σ 0 2 ) S 0 = e ( r d ) A d ( T t ) [ m 1 N ( m 1 v 1 ) + v 1 2 π e m 1 2 2 v 1 ] + 1 2 e ( r d ) A d ( T t ) [ M 12 N ( m 1 v 1 ) + 1 2 π e m 1 2 2 v 1 M 11 A ¯ 1 1 32 π e m 1 2 2 v 1 v 1 3 2 V 11 2 + 1 8 π e m 1 2 2 v 1 v 1 1 2 V 11 B ¯ 1 + 1 8 π e m 1 2 2 v 1 v 1 1 2 V 12 ] × [ 2 σ 4 ( e k k 1 ) k 2 σ 4 ] + 1 6 e ( r d ) A d ( T t ) [ M 13 N ( m 1 v 1 ) + 1 2 π e m 1 2 2 v 1 α 11 + 1 32 π e m 1 2 2 v 1 α 12 + 1 8 π e m 1 2 2 v 1 ( α 13 + α 14 ) ] × σ 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 6 k 3 ) 3 k 3 ]
When t [ 0 , T A , T ] , then
D e l t a = Δ = f ( S 0 , σ 0 2 ) S 0 = e r ( T t ) [ m 2 N ( m 2 v 2 ) + v 2 2 π e m 2 2 2 v 2 + S t N ( m 2 v 2 ) m 2 S t ] + 1 2 { e r ( T t ) [ M 22 N ( m 2 v 2 ) + 1 2 π e m 2 2 2 v 2 M 21 A ¯ 2 1 32 π e m 2 2 2 v 2 v 2 3 2 V 21 2 + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 22 + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 21 B ¯ 2 ] + S t e r ( T t ) [ 1 2 π e m 2 2 2 v 2 M 22 v 2 1 2 m 2 S t 1 2 π e m 2 2 2 v 2 m 2 v 2 1 M 21 A ¯ 2 m 2 S t + 1 32 π e m 2 2 2 v 2 m 2 v 2 5 2 V 21 2 m 2 S t 1 8 π e m 2 2 2 v 2 m 2 v 2 3 2 V 22 m 2 S t 1 8 π e m 2 2 2 v 2 m 2 v 2 3 2 V 21 B ¯ 2 m 2 S t + 1 2 π e m 2 2 2 v 2 M 21 A ¯ 2 S t + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 21 B ¯ 2 S t ] } × [ 2 σ 4 ( e k k 1 ) k 2 σ 4 ] + 1 6 { e r ( T t ) [ M 23 N ( m 2 v 2 ) + 1 2 π e m 2 2 2 v 2 α 21 + 1 32 π e m 2 2 2 v 2 α 22 + 1 8 π e m 2 2 2 v 2 ( α 23 + α 24 ) ] + S t e r ( T t ) [ 1 2 π e m 2 2 2 v 2 M 23 v 2 1 2 m 2 S t 1 2 π e m 2 2 2 v 2 m 2 v 2 1 α 21 m 2 S t + 1 2 π e m 2 2 2 v 2 α 21 S t 1 32 π e m 2 2 2 v 2 m 2 v 2 1 α 22 m 2 S t + 1 32 π e m 2 2 2 v 2 α 22 S t 1 8 π e m 2 2 2 v 2 m 2 v 2 1 ( α 23 + α 24 ) m 2 S t + 1 8 π e m 2 2 2 v 2 ( α 23 S t + α 24 S t ) ] } × σ 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 6 k 3 ) 3 k 3 ]
Please refer to the Appendix A for notations and proofs of the propositions.

3. Results

Monte Carlo simulation was performed to clarify the differences in pricing between a constant-volatility FSAO and SVFSAO according to the pricing results and computer operation time. Next, the accuracy of the approximate analytical solution for the SVFSAO was compared with the solution of Monte Carlo simulation, which served as the benchmark. Finally, the hedge ratio of the SVFSAO was determined and compared with that of a constant-volatility FSAO. The results of the aforementioned three numerical analyses demonstrate the high value of the developed model and its ability to address the potential risk incurred by HFT.

3.1. Differences in Pricing between Constant-Volatility FSAOs and SVFSAOs

To demonstrate the value of the developed SVFSAO, Monte Carlo simulation was performed to price FSAOs under the conditions of constant volatility and stochastic volatility in different scenarios (To ensure consistency with reality, for numerical investigations, it was assumed that there are 365 days in a year.).
As indicated in Table 1, pricing differences between FSAOs under constant volatility and stochastic volatility conditions were evident, exceeding 4% for several conditions. This result demonstrates that the constant-volatility FSAO pricing model could be prone to mispricing and have detrimental effects on users’ justifications for hedging.
Computation time is a crucial factor for managing computer-based trading, especially HFT. The time required for computing a model price for options is at least 11 s for the case of a constant-volatility FSAO and at least 28 s for the case of an SVFSAO. A pricing technique that relies on simulation is incapable of handling HFT. In the present study, this numerical analysis supported the development of an approximate analytical solution for FSAOs.
Table 2 details the pricing performance of the present study’s approximate analytical path-dependent FSAO solutions. Compared with Monte Carlo simulation, for a large sample path, the developed model had a superior pricing speed and accuracy, proving that the FSAO can be applied to manage HFT.
The values in the “AAS” column of Table 2 (The parameter settings were selected with reference to [25]) are the approximate analytical solutions for the SVFSAO. Compared with the solutions of Monte Carlo simulation, the present study’s approximate analytical solutions were closer to the theoretical value, and all pricing errors were less than 1%. Furthermore, the computing time for the approximate analytical solutions was more than 2 million times faster than that of the Monte Carlo solutions; only approximately 20 μs was required for computing a SVFSAO pricing solution. This computing performance is sufficient for dealing with the fast liquidity of HFT.

3.2. Comparison of the Delta Value of FSAO between Constant and Stochastic Volatility Conditions

Hedging is vital for option market participants. Without rebalancing the underlying asset position according to the hedging ratio, delta, option market participants cannot maintain a risk-free portfolio of options and underlying assets and may incur investment loss from the pricing changes of underlying assets.
Table 3 indicates the hedging performance in terms of the delta value under the developed SVFSAO and constant-volatility FSAO. The employed constant-volatility FSAO hedge model refers to the model developed by Tsao, Chang, and Lin [14].
The values in the “delta (SV)” column represent the developed model, which exhibited greater volatility than the values in the “delta (CV)” column for the TCL model. The differences between SV and CV were more obvious at higher ξ values. The results in this table support that the constant-volatility FSAO model may underprice the hedge ratio delta. The market participants with FSAOs should be careful when using the TCL constant-volatility model because it could lead them to take the wrong hedging position.
The accuracy of the present study’s pricing formula is based on two implications. The first is that the analytical solutions of hedge ratios for these options are reliable. The second is that the present study’s pricing formula can be used as a control variable to enhance computational accuracy and efficiency when numerical methods (e.g., Monte Carlo simulation or binomial tree approaches) are being used for pricing American SVFSAOs. Pricing and hedging performance with respect to computation time and pricing accuracy are critical for handling the added risk introduced by the increase in computer-based trading, particularly HFT.

4. Conclusions

An approximate analytical solution for a stochastic-volatility Asian option valuation model based on the model introduced by Hull and White [1] that integrates Taylor series expansion was developed in this study. Greek hedging was also developed to verify the natural hedging ability of derivatives.
Numerical experiments demonstrated that the present study’s approximate analytical solution outperforms performance and a higher computational efficiency compared with a benchmark approach for large-sample simulations. The approximate analytical solution in this study can serve as a practical tool for pricing and hedging SVFSAOs; it is straightforward to implement and is suitable for other complex derivative pricing because of the importance of computational speed. Thus, this research provides evidence that the developed SVFSAO model is capable of overcoming the potential challenges, such as liquidity and volatility toxicity, of HFT.
Other advantages of the developed model are summarized as follows: First, it can easily calculate the hedge ratios of SVFSAOs for use in hedging SVFSAO valuations. Second, the developed FSAO valuation model can be applied as a generalized Asian option valuation model under the condition of stochastic volatility pricing. The developed SVFSAO can be applied for the valuation of any type of Asian option simply through adjustment of the time period variable. Third, despite the fact that a normal probability distribution was assumed in the approximate arithmetic normal Distribution, It Did Not Generate Any Inefficient And Inaccurate Questions On Model Pricing And hedging performance. Hence, this assumption can be relaxed in future research.
However, some limitations were inevitable because, following Hull and White [1], the correlation between the instantaneous underlying price process and the stochastic volatility process was ignored. Moreover, the developed approximate analytical solution is inapplicable for derivation under the condition of mean-reversion pricing.

Author Contributions

C.-C.H. and C.-G.L. conceived of the presented idea. C.-C.H., C.-G.L. and T.-J.K. developed the theory and performed the computations. C.-C.H., C.-G.L. and T.-J.K. verified the analytical methods. C.-C.H. and C.-G.L. encouraged T.-J.K. to investigate “Pricing of Arithmetic Asian Option under Stochastic Volatility Dynamics- Overcome the Risk of High Frequent Trading” and supervised the findings of this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We acknowledge any support given from Department of Business Administration in Soochow University, Taiwan, Department of Financial Engineering and Actuarial Mathematics in Soochow University, Taiwan, and Pervasive Artificial Intelligence Research (PAIR) Labs, Taiwan.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Derivation of Hedge Ratio—Delta

When t [ 0 , T A ] , then
D e l t a = Δ = f S t = C S t + 1 2 C 2 S t V a r ( V ¯ ) + 1 6 C 3 S t S k e w ( V ¯ )
where
C S t = e ( r d ) A d ( T t ) [ m 1 N ( m 1 v 1 ) + v 1 2 π e m 1 2 2 v 1 ]
C 2 S t = e ( r d ) A d ( T t ) [ M 12 N ( m 1 v 1 ) + 1 2 π e m 1 2 2 v 1 M 11 A ¯ 1 1 32 π e m 1 2 2 v 1 v 1 3 2 V 11 2 + 1 8 π e m 1 2 2 v 1 v 1 1 2 V 11 B ¯ 1 + 1 8 π e m 1 2 2 v 1 v 1 1 2 V 12 ]
V a r ( V ¯ ) = 2 σ 4 ( e k k 1 ) k 2 σ 4
C 3 S t = e ( r d ) A d ( T t ) [ M 13 N ( m 1 v 1 ) + 1 2 π e m 1 2 2 v 1 α 11 + 1 32 π e m 1 2 2 v 1 α 12 + 1 8 π e m 1 2 2 v 1 ( α 13 + α 14 ) ]
S k e w ( V ¯ ) = σ 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 6 k 3 ) 3 k 3 ]
in which
m 1 = 1 2 A ( r d ) + 1 3 A 2 ( r d ) 2 1 3 A 2 ( r d ) σ 2 + 1 12 A 2 σ 4
v 1 = [ 1 3 A + 3 4 A 2 ( r d ) + 7 15 A 3 ( r d ) 2 ] σ 2 + [ 1 8 A 2 7 15 A 3 ( r d ) ] σ 4 + 7 60 A 3 σ 6
M 11 = m 1 V ¯ = 1 3 A 2 ( r d ) + 1 6 A 2 V ¯
M 12 = 2 m 1 V ¯ 2 = 1 6 A 2
M 13 = 3 m 1 V ¯ 3 = 0
V 11 = v 1 V ¯ = [ 1 3 A + 3 4 A 2 ( r d ) + 7 15 A 3 ( r d ) 2 ] + [ 1 4 A 2 14 15 A 3 ( r d ) ] V ¯ + 7 20 A 3 V ¯ 2
V 12 = 2 v 1 V ¯ 2 = [ 1 4 A 2 14 15 A 3 ( r d ) ] + 7 10 A 3 V ¯
V 13 = 3 v 1 V ¯ 3 = 7 10 A 3
A ¯ 1 = M 11 v 1 1 2 1 2 m 1 v 1 3 2 V 11
B ¯ 1 = m 1 v 1 1 M 11 + 1 2 m 1 2 v 1 2 V 11
α 11 = 2 M 12 A ¯ 1 + M 11 A ¯ 1 B ¯ 1 + M 11 ( v 1 1 2 M 12 v 1 3 2 M 11 V 11 1 2 m 1 v 1 3 2 V 12 + 3 4 m 1 v 1 5 2 V 11 2 )
α 12 = 3 2 v 1 5 2 V 11 3 v 1 3 2 V 11 2 B ¯ 1 2 v 1 3 2 V 11 V 12
α 13 = 1 2 v 1 3 2 V 11 V 12 + v 1 1 2 V 12 B ¯ 1 + v 1 1 2 V 13
α 14 = 1 2 v 1 3 2 V 11 2 B ¯ 1 + v 1 1 2 V 11 B ¯ 1 2 + v 1 1 2 V 12 B ¯ 1 + v 1 1 2 V 11 ( v 1 1 M 11 2 + 2 m 1 v 1 2 M 11 V 11 m 1 v 1 1 M 12 m 1 2 v 1 3 V 11 2 + 1 2 m 1 2 v 1 2 V 12 )
When t [ T A , T ] ,
D e l t a = Δ = f S t = C S t + 1 2 C 2 S t V a r ( V ¯ ) + 1 6 C 3 S t S k e w ( V ¯ )
where
C S t = e r ( T t ) [ m 2 N ( m 2 v 2 ) + v 2 2 π e m 2 2 2 v 2 + S t N ( m 2 v 2 ) m 2 S t ]
C 2 S t = e r ( T t ) [ M 22 N ( m 2 v 2 ) + 1 2 π e m 2 2 2 v 2 M 21 A ¯ 2 1 32 π e m 2 2 2 v 2 v 2 3 2 V 21 2 + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 22 + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 21 B ¯ 2 ] + S t e r ( T t ) [ 1 2 π e m 2 2 2 v 2 M 22 v 2 1 2 m 2 S t 1 2 π e m 2 2 2 v 2 m 2 v 2 1 M 21 A ¯ 2 m 2 S t + 1 32 π e m 2 2 2 v 2 m 2 v 2 5 2 V 21 2 m 2 S t 1 8 π e m 2 2 2 v 2 m 2 v 2 3 2 V 22 m 2 S t 1 8 π e m 2 2 2 v 2 m 2 v 2 3 2 V 21 B ¯ 2 m 2 S t + 1 2 π e m 2 2 2 v 2 M 21 A ¯ 2 S t + 1 8 π e m 2 2 2 v 2 v 2 1 2 V 21 B ¯ 2 S t ]
V a r ( V ¯ ) = 2 σ 4 ( e k k 1 ) k 2 σ 4
C 3 S t = e r ( T t ) [ M 23 N ( m 2 v 2 ) + 1 2 π e m 2 2 2 v 2 α 21 + 1 32 π e m 2 2 2 v 2 α 22 + 1 8 π e m 2 2 2 v 2 ( α 23 + α 24 ) ] + S t e r ( T t ) [ 1 2 π e m 2 2 2 v 2 M 23 v 2 1 2 m 2 S t 1 2 π e m 2 2 2 v 2 m 2 v 2 1 α 21 m 2 S t + 1 2 π e m 2 2 2 v 2 α 21 S t 1 32 π e m 2 2 2 v 2 m 2 v 2 1 α 22 m 2 S t + 1 32 π e m 2 2 2 v 2 α 22 S t 1 8 π e m 2 2 2 v 2 m 2 v 2 1 ( α 23 + α 24 ) m 2 S t + 1 8 π e m 2 2 2 v 2 ( α 23 S t + α 24 S t ) ]
S k e w ( V ¯ ) = σ 6 [ e 3 k ( 9 + 18 k ) e k + ( 8 + 24 k + 18 k 2 6 k 3 ) 3 k 3 ]
in which
m 2 = [ 1 + ( T t ) ( r d ) 1 A ( T t ) 1 2 A ( T t ) 2 ( r d ) M t S t + 1 2 ( T t ) 2 ( r d ) 2 1 6 A ( T t ) 3 ( r d ) 2 ] + [ 1 2 ( T t ) 2 ( r d ) + 1 6 A ( T t ) 3 ( r d ) ] σ 2 + [ 1 8 ( T t ) 2 1 24 A ( T t ) 3 ] σ 4
v 2 = σ 2 ( T t ) [ 1 1 A ( T t ) + 2 ( T t ) ( r d ) + 1 3 A 2 ( T t ) 2 5 3 A ( T t ) 2 ( r d ) + ( T t ) 2 ( r d ) 2 + 5 12 A 2 ( T t ) 3 ( r d ) 2 3 A ( T t ) 3 ( r d ) 2 + 2 15 A 2 ( T t ) 4 ( r d ) 2 ] + σ 4 ( T t ) 2 [ 1 2 + 1 2 A ( T t ) ( T t ) ( r d ) 1 8 A 2 ( T t ) 2 + 2 3 A ( T t ) 2 ( r d ) 2 15 A 2 ( T t ) 3 ( r d ) ] + σ 6 ( T t ) 3 [ 1 4 1 6 A ( T t ) + 1 30 A 2 ( T t ) 2 ]
M 21 = m 2 V ¯ = [ 1 2 ( T t ) 2 ( r d ) + 1 6 A ( T t ) 3 ( r d ) ] + [ 1 4 ( T t ) 2 1 12 A ( T t ) 3 ] V ¯
M 22 = 2 m 2 V ¯ 2 = [ 1 4 ( T t ) 2 1 12 A ( T t ) 3 ]
M 23 = 3 m 2 V ¯ 3 = 0
V 21 = v 2 V ¯ = [ 1 1 A ( T t ) + 2 ( T t ) ( r d ) + 1 3 A 2 ( T t ) 2 5 3 A ( T t ) 2 ( r d ) + ( T t ) 2 ( r d ) 2 + 5 12 A 2 ( T t ) 3 ( r d ) 2 3 A ( T t ) 3 ( r d ) 2 + 2 15 A 2 ( T t ) 4 ( r d ) 2 ] ( T t ) + [ 1 + 1 A ( T t ) 2 ( T t ) ( r d ) 1 4 A 2 ( T t ) 2 + 4 3 A ( T t ) 2 ( r d ) 4 15 A 2 ( T t ) 3 ( r d ) ] ( T t ) 2 V ¯ + [ 3 4 1 2 A ( T t ) + 1 10 A 2 ( T t ) 2 ] ( T t ) 3 V ¯ 2
V 22 = 2 v 2 V ¯ 2 = [ 1 + 1 A ( T t ) 2 ( T t ) ( r d ) 1 4 A 2 ( T t ) 2 + 4 3 A ( T t ) 2 ( r d ) 4 15 A 2 ( T t ) 3 ( r d ) ] ( T t ) 2 + [ 3 2 1 A ( T t ) + 1 5 A 2 ( T t ) 2 ] ( T t ) 3 V ¯
V 23 = 3 v 2 V ¯ 3 = [ 3 2 1 A ( T t ) + 1 5 A 2 ( T t ) 2 ] ( T t ) 3
A ¯ 2 = V ¯ ( m 2 v 2 ) = M 21 v 2 1 2 1 2 m 2 v 2 3 2 V 21
B ¯ 2 = V ¯ ( m 2 2 2 v 2 ) = m 2 v 2 1 M 21 + 1 2 m 2 2 v 2 2 V 21
α 21 = 2 M 22 A ¯ 2 + M 21 A ¯ 2 B ¯ 2 + M 21 ( v 2 1 2 M 22 v 2 3 2 M 21 V 21 1 2 m 2 v 2 3 2 V 22 + 3 4 m 2 v 2 5 2 V 21 2 )
α 22 = 3 2 v 2 5 2 V 21 3 v 2 3 2 V 21 2 B ¯ 2 2 v 2 3 2 V 21 V 22
α 23 = 1 2 v 2 3 2 V 21 V 22 + v 2 1 2 V 22 B ¯ 2 + v 2 1 2 V 23
α 24 = 1 2 v 2 3 2 V 21 2 B ¯ 2 + v 2 1 2 V 21 B ¯ 2 2 + v 2 1 2 V 22 B ¯ 2 + v 2 1 2 V 21 ( v 2 1 M 21 2 + 2 m 2 v 2 2 M 21 V 21 m 2 v 2 1 M 22 m 2 2 v 2 3 V 21 2 + 1 2 m 2 2 v 2 2 V 22 )
m 2 S t = M t S t 2
A ¯ 2 S t = 1 2 v 2 3 2 V 21 m 2 S t
B ¯ 2 S t = v 2 1 M 21 m 2 S t + m 2 v 2 2 V 21 m 2 S t
α 21 S t = 2 M 22 A ¯ 2 S t + M 21 B ¯ 2 A ¯ 2 S t + M 21 A ¯ 2 B ¯ 2 S t + M 21 ( 3 4 v 2 5 2 V 21 2 m 2 S t 1 2 v 2 3 2 V 22 m 2 S t )
α 22 S t = v 2 3 2 V 21 2 B ¯ 2 S t
α 23 S t = v 2 1 2 V 22 B ¯ 2 S t
α 24 S t = 1 2 v 2 3 2 V 21 2 B ¯ 2 S t + 2 v 2 1 2 V 21 B ¯ 2 B ¯ 2 S t + v 2 1 2 V 22 B ¯ 2 S t + v 2 1 2 V 21 ( 2 v 2 2 M 21 V 21 m 2 S t v 2 1 M 22 m 2 S t 2 m 2 v 2 3 V 21 2 m 2 S t + m 2 v 2 2 V 22 m 2 S t )

References

  1. Bernales, A. Make-take Decisions under High-frequency Trading Competition. J. Finance Mark. 2019, 45, 1–18. [Google Scholar] [CrossRef]
  2. Shafi, K.; Latif, N.; Shad, S.A.; Idrees, Z. High-frequency trading: Inverse relationship of the financial markets. Phys. A Stat. Mech. Appl. 2019, 527, 121067. [Google Scholar] [CrossRef]
  3. Bibinger, M.; Hautsch, N.; Malec, P.; Reiss, M. Estimating the Spot Covariation of Asset Prices—Statistical Theory and Empirical Evidence. J. Bus. Econ. Stat. 2019, 37, 419–435. [Google Scholar] [CrossRef]
  4. Said, E.; Ayed, A.B.H.; Thillou, D.; Rabeyrin, J.-J.; Abergel, F. Market Impact: A Systematic Study of the High Frequency Options Market. Quant. Finance 2020, 21, 1–16. [Google Scholar] [CrossRef]
  5. Chaboud, A.P.; Chiquoine, B.; Hjalmarson, E.; Vega, C. Rise of the Machines: Algorithmic Trading in the Foreign Exchange Market. J. Finance 2014, 69, 2045–2084. [Google Scholar] [CrossRef] [Green Version]
  6. Baldauf, M.; Mollner, J. High-Frequency Trading and Market Performance. J. Finance 2020, 75, 1495–1523. [Google Scholar] [CrossRef] [Green Version]
  7. Corsaro, S.; Kyriakou, I.; Marazzina, D.; Marino, Z. A general framework for pricing Asian options under stochastic volatility on parallel architectures. Eur. J. Oper. Res. 2019, 272, 1082–1095. [Google Scholar] [CrossRef] [Green Version]
  8. Wang, X. Analytical valuation of Asian options with counterparty risk under stochastic volatility models. J. Futures Mark. 2020, 40, 410–429. [Google Scholar] [CrossRef]
  9. Yao, N.; Ling, Z.; Zhang, J.; Xizo, M. Short maturity conditional Asian options in local volatility models. Math. Finance Econ. 2020, 14, 307–328. [Google Scholar] [CrossRef]
  10. Heston, S.L. A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options. Rev. Finance Stud. 1993, 6, 327–343. [Google Scholar] [CrossRef] [Green Version]
  11. Kemna, A.; Vorst, A. A pricing method for options based on average asset prices. J. Bank. Finance 1990, 2, 52–66. [Google Scholar]
  12. Rogers, L.C.G.; Shi, Z. The value of an Asian option. J. Appl. Probab. 1995, 32, 1077–1088. [Google Scholar] [CrossRef]
  13. Hoogland, J.K.; Neumann, C.D.D. Local scale invariance and contingent claim pricing. Int. J. Theor. Appl. Finance 2001, 4, 1–21. [Google Scholar] [CrossRef] [Green Version]
  14. Hoogland, J.K.; Neumann, C.D.D. ; Neumann, C.D.D. Local scale invariance and contingent claim pricing II. Int. J. Theor. Appl. Finance 2001, 4, 23–43. [Google Scholar] [CrossRef] [Green Version]
  15. Vecer, J. A new PDE approach for pricing arithmetic average Asian options. J. Comput. Finance 2001, 4, 105–113. [Google Scholar] [CrossRef] [Green Version]
  16. Vecer, J. Unified pricing of Asian options risk. J. Comput. Finance 2001, 4, 113–116. [Google Scholar]
  17. Linetsky, V. Pricing Equity Derivatives Subject to Bankruptcy. Math. Finance 2006, 16, 255–282. [Google Scholar] [CrossRef] [Green Version]
  18. Vanmaele, M.; Deelstra, G.; Liinev, J.; Dhaene, J.; Goovaerts, M.J. Bounds for the price of discrete arithmetic Asian options. J. Comput. Finance 2006, 185, 51–90. [Google Scholar] [CrossRef]
  19. Hull, J.; White, A. The Pricing of Options on Assets with Stochastic Volatilities. J. Finance 1987, 42, 281–300. [Google Scholar] [CrossRef]
  20. Stein, E.M.; Stein, J.C. Stock Price Distributions with Stochastic Volatility: An Analytic Approach. Rev. Financ. Stud. 1991, 4, 727–752. [Google Scholar] [CrossRef]
  21. Turnbull, S.; Wakeman, L. A quick algorithm for pricing European average options. J. Financ. Quant. Anal. 1991, 26, 377–389. [Google Scholar] [CrossRef]
  22. Ritchken, P.; Sankarasubramanian, L.; Vijh, A.M. The valuation of path dependent contracts on the average. Manag. Sci. 1993, 39, 1202–1213. [Google Scholar] [CrossRef]
  23. Milevsky, M.A.; Posner, S.E. Asian options, the sum of lognormals, and the reciprocal gamma distribution. J. Financ. Quant. Anal. 1998, 33, 409–422. [Google Scholar] [CrossRef]
  24. Fouque, J.P.; Han, C.H. Pricing Asian options with stochastic volatility. Quant. Finance 2003, 3, 353–362. [Google Scholar] [CrossRef]
  25. Tsao, C.Y.; Chang, C.C.; Lin, C.G. Analytic Approximation Formulae for Pricing Forward Starting Asian option. J. Futures Mark. 2003, 23, 487–516. [Google Scholar] [CrossRef]
  26. Sun, J.; Chen, L.; Li, S. A quasi-analytical pricing model for arithmetic Asian options. J. Futures Mark. 2013, 33, 1143–1166. [Google Scholar] [CrossRef]
  27. Li, W.; Chen, S. Pricing and hedging of arithmetic Asian options via the Edgeworth series expansion approach. J. Financ. Data Sci. 2016, 2, 1–25. [Google Scholar] [CrossRef] [Green Version]
  28. Lin, C.G.; Chang, C.C. Approximate Analytic Solution for Asian Options with Stochastic Volatility. North Am. J. Econ. Finance 2019, 54, 100949. [Google Scholar] [CrossRef]
  29. Aprahamian, H.; Maddah, B. Pricing Asian options via compound gamma and orthogonal polynomials. Appl. Math. Comput. 2015, 264, 21–43. [Google Scholar] [CrossRef]
  30. Dufresne, D. Laguerre series for Asian and other options. Math. Finance 1999, 10, 407–428. [Google Scholar] [CrossRef]
  31. Willems, S. Asian Option Pricing with Orthogonal Polynomials. Quant. Finance 2019, 19, 605–618. [Google Scholar] [CrossRef] [Green Version]
  32. Broadie, M.; Kaya, O. Exact Simulation of Stochastic Volatility and Other Affine Jump Diffusion Processes. Oper. Res. 2006, 54, 217–231. [Google Scholar] [CrossRef]
  33. Kouritzin, M.A. Explicit Heston solutions and stochastic approximation for path-dependent option pricing. Int. J. Theor. Appl. Finance 2018, 21, 217–231. [Google Scholar] [CrossRef]
  34. Kouarfate, R.; Kouritzin, A.; Mackay, A. Explicit Solution Method for the 3/2 Model; Cornell University: Ithaca, NY, USA, 2020. [Google Scholar]
  35. Kirkby, J.L.; Nguyen, D. Efficient Asian option pricing under regime switching jump diffusions and stochastic volatility models. Ann. Finance 2020, 16, 307–351. [Google Scholar] [CrossRef]
Table 1. Comparison of FSAOs with stochastic and constant volatilities.
Table 1. Comparison of FSAOs with stochastic and constant volatilities.
(1)(2)(3)(4)(5)
μ ξ MC(SV)MC(CV)Diff Ratio
(2) − (1)/(2)
Time
(SV)
Time
(CV)
00.34.962404.989230.537758%28.8371811.98062
(0.02769)(0.02047)
0.64.886334.989232.062443%28.6793611.97470
(0.02786)(0.02047)
0.94.763764.989234.519134%28.4134411.96720
(0.028)(0.02047)
0.10.35.088544.98923−1.990488%28.7806211.96250
(0.02862)(0.02047)
0.65.009714.98923−0.410484%28.6831211.96968
(0.02881)(0.02047)
0.94.882684.989232.135600%28.4134411.96500
(0.02895)(0.02047)
0.20.35.218884.98923−4.602915%28.7697011.96500
(0.0296)(0.02047)
0.65.137214.98923−2.965989%28.6881211.96812
(0.02978)(0.02047)
0.95.005584.98923−0.327706%28.4103211.95906
(0.02994)(0.02047)
Note: The parameters were set as follows: S 0 = 100 , r = 10 % , d = 0 % , σ 0 = 30 % , T = 180 days, A = 90 days. Option prices in columns (1) and (2) reflect stochastic and constant volatility conditions, respectively, of Monte Carlo simulation with 100,000 paths and 50 repetitions. The standard error is presented in parentheses. Column (3) displays the difference ratio, which denotes the percentage difference in pricing. Columns (4) and (5) respectively indicate the time required for simulation under stochastic and constant volatility conditions. Monte Carlo simulation was performed using a personal computer with the following specifications: P4-3.4 GHz with 1 GB RAM.
Table 2. Comparison of the analytical solutions of the proposed FSAO and the Monte Carlo method.
Table 2. Comparison of the analytical solutions of the proposed FSAO and the Monte Carlo method.
(1)(2)(3)(4)(5)(6)
T
(Days)
MCASCDiff Ratio
(2) − (1)/(1)
Time
(MC)
Time
(AAS)
Efficiency = (4)/(5)
301.834101.8374160.180797%57.288142.343 × 10−52.445 × 106
(0.0097)
903.364403.3772930.383218%57.516262.281 × 10−52.522 × 106
(0.018)
1203.963573.9816090.455120%57.486562.312 × 10−52.486 × 106
(0.02133)
2405.945185.9808350.599730%57.197822.500 × 10−52.288 × 106
(0.03275)
Note: The parameters were set as follows: S 0 = 100 , r = 10 % , d = 0 % , μ = 0 , ξ = 0.15 , σ 0 = 30 % , A = T / 2 trading days. Column (1) is the results of the Monte Carlo (MC) method, with 100,000 paths and 50 repetitions. The standard error is presented in parenthesis. Column (2) displays the option value obtained using the analytical solution proposed in this study. Column (3) presents the difference ratio, which denotes the pricing difference as a percentage. Columns (4) and (5) indicate the time required for simulation using the MC method and the approximate analytical solution (AAS), respectively. Column (6) presents the efficiency (the required simulation time for the benchmark method divided by that required for the present study’s formula). MC simulation was performed on a personal computer with the following specifications: P4-3.4 GHz with 1 GB RAM.
Table 3. Comparison of the delta value of the FSAO between constant and stochastic volatility conditions.
Table 3. Comparison of the delta value of the FSAO between constant and stochastic volatility conditions.
(1)(2)(3)
σ ξ Delta
(CV)
Delta
(SV)
Diff Ratio
(2) − (1)/(1)
0.20.21.4381 × 10−193.4619 × 10−19140.7296%
0.41.4381 × 10−191.1246 × 10−18682.0450%
0.61.4381 × 10−193.0113 × 10−181993.9933%
0.81.4381 × 10−196.9566 × 10−184737.4300%
0.30.21.0431 × 10−91.3115 × 10−925.7356%
0.41.0431 × 10−92.2116 × 10−9112.0294%
0.61.0431 × 10−94.0376 × 10−9287.0919%
0.81.0431 × 10−97.3144 × 10−9601.2435%
0.40.23.3847 × 10−63.6421 × 10−67.6060%
0.43.3847 × 10−64.4579 × 10−631.7077%
0.63.3847 × 10−65.9667 × 10−676.2859%
0.83.3847 × 10−68.4084 × 10−6148.4258%
0.50.21.5510 × 10−41.5954 × 10−42.8660%
0.41.5510 × 10−41.7323 × 10−411.6924%
0.61.5510 × 10−41.9726 × 10−427.1859%
0.81.5510 × 10−42.3357 × 10−450.5983%
Note: The parameters were set as follows: S 0 = 100 , r = 10 % , d = 0 % , μ = 0 , T = 0.5 ,   A = T / 2 , t = 0.75 T . Column (1) presents the delta value of the call option when the analytical solution of the TCL model is used. Column (2) presents the delta value of the call option when the analytical solution of the proposed model is used. Column (3) provides the difference ratio, which denotes the delta difference as a percentage.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hsu, C.-C.; Lin, C.-G.; Kuo, T.-J. Pricing of Arithmetic Asian Options under Stochastic Volatility Dynamics: Overcoming the Risks of High-Frequency Trading. Mathematics 2020, 8, 2251. https://doi.org/10.3390/math8122251

AMA Style

Hsu C-C, Lin C-G, Kuo T-J. Pricing of Arithmetic Asian Options under Stochastic Volatility Dynamics: Overcoming the Risks of High-Frequency Trading. Mathematics. 2020; 8(12):2251. https://doi.org/10.3390/math8122251

Chicago/Turabian Style

Hsu, Chih-Chen, Chung-Gee Lin, and Tsung-Jung Kuo. 2020. "Pricing of Arithmetic Asian Options under Stochastic Volatility Dynamics: Overcoming the Risks of High-Frequency Trading" Mathematics 8, no. 12: 2251. https://doi.org/10.3390/math8122251

APA Style

Hsu, C. -C., Lin, C. -G., & Kuo, T. -J. (2020). Pricing of Arithmetic Asian Options under Stochastic Volatility Dynamics: Overcoming the Risks of High-Frequency Trading. Mathematics, 8(12), 2251. https://doi.org/10.3390/math8122251

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