# Pricing Arithmetic Asian Options under Hybrid Stochastic and Local Volatility

Authors: Min-Ku Lee; Jeong-Hoon Kim; Kyu-Hwan Jang
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

Recently, hybrid stochastic and local volatility models have become an industry standard for the pricing of derivatives and other problems in finance. In this study, we use a multiscale stochastic volatility model incorporated by the constant elasticity of variance to understand the price structure of continuous arithmetic average Asian options. The multiscale partial differential equation for the option price is approximated by a couple of single scale partial differential equations. In terms of the elasticity parameter governing the leverage effect, a correction to the stochastic volatility model is made for more efficient pricing and hedging of Asian options.


---

## Body

## 1. Introduction

Since the well-known work of Black and Scholes [1] on the classical vanilla European option, there has been concern about the pricing of more complicated exotic options. An exotic option is a derivative which has a payoff structure more complex than commonly traded vanilla options. They are usually traded in over-the-counter market or embedded in structured products. Also, the pricing of them tends to require more complex methods than the classical Black-Scholes approach.This paper is concerned with one of the exotic options called an Asian option. This option is a path dependent option whose final payoff depends on the paths of its underlying asset. More precisely, the payoff is determined by the average value of underlying prices over some prescribed period of time. The name of “Asian” options is known to come from the fact that they were first priced in 1987 by David Spaughton and Mark Standish of Bankers Trust when they were working in Tokyo, Japan (cf. [2]). The main motivation of creating these options is that their averaging feature could reduce the risk of market manipulation of the underlying risky asset at maturity. Since Asian options reduce the volatility inherent in the option, the price of these options is usually lower than the price of classical European vanilla options. Note that there are two types of Asian options depending on the style of averaging: one is the arithmetic average Asian option and the other is the geometric average Asian option.Since there is no general analytical formula for the price of Asian option, a variety of techniques have been developed to approximate the price of this option. Subsequently, there has been quite an amount of literature devoted to studying this option. For instance, Geman and Yor [3] computed the Laplace transform of the price of continuously sampled Asian options. However, there is a problem of slow convergence for low volatility or short time-to-maturity cases as indicated by Fu et al. [4]. Linetsky [5] derived a new integral formula for the price of Asian options but with the same convergence problem. Apart from this approximation technique, there are Monte Carlo simulation approach and partial differential equation (PDE) approach. Monte Carlo simulation methods typically require variance reduction techniques and also have to take into account the discretization error caused by discrete sampling. Refer to Kemna and Vorst [6] for a discussion of the pricing of Asian options with Monte Carlo methods. On the other hand, the PDE methods must deal with an extra state variable representing the running sum of the underlying process, which leads to an issue for reducing the dimension of the PDE. Refer to Ingersoll [7] and Rogers and Shi [8] and Vecer [9] for the PDE methods.It is well known that the constant volatility assumption, on which a review of literature quoted above is based, for the underlying asset price is severely in contrast with many empirical studies which demonstrate the skew or smile effect of implied volatility, fat-tailed and asymmetric returns distributions, and the mean-reversion of volatility. Thus a number of alternative underlying models have been proposed. The constant elasticity of variance model by Cox [10], a stochastic volatility model by Heston [11] or Fouque et al. [12], and a Levy model by Carr et al. [13] are among those representative ones that can reproduce an empirically reasonable outcome.So, it is desirable to study Asian options based on these alternative models. In fact, there are a number of recent studies along the lines of this type of extension. For instance, B. Peng and F. Peng [14] for the CEV model, Fouque and Han [15] for a stochastic volatility model, and Lemmens et al. [16] for a Levy type model are among those works extending the price for the Black-Scholes framework with constant volatility. As long as we understand, up to now, however, there has been no work for Asian options based on a hybrid stochastic and local volatility model. Recently, hybrid stochastic and local volatility models have become an industry standard for the pricing of derivatives and several financial institutions have incorporated those models into their systems [17]. Therefore, it is worth studying Asian options on the hybrid models.This paper studies the pricing of an arithmetic Asian option under a hybrid stochastic and local volatility model which was introduced by Choi et al. [18], where the volatility is given by the product of a multiscale stochastic process and a power (the elasticity of variance) of the underlying's price. The hybrid nature of this volatility enables us to capture the leverage effect produced by the constant elasticity of variance (CEV) model as well as the smile effect of implied volatility, fat-tailed and asymmetric returns distributions, the tendency of the volatility process to revert towards a long-term mean at a certain rate, and a degree of correlation between the randomness of volatility and the randomness of underlying's price produced by a “pure” stochastic volatility (SV) model. This hybrid model is called the SVCEV model. So, this paper will generalize both [14, 15] into an approximation problem for the price of arithmetic average Asian options based on the SVCEV model.This paper is structured as follows. In Section2, we review the SVCEV model introduced by [18]. A multiscale partial differential equation for the price of an arithmetic average Asian option is obtained in Section 3. Section 4 is devoted to obtain an approximated option price under the condition of fast mean-reverting volatility. In Section 5, the effect of hybrid structure on the pricing of the Asian option is illustrated by a numerical study.


## 2. The SVCEV Model

In this section, we establish a partial differential equation for the price of Asian floating strike call option based on the SVCEV model.As introduced by [18], the SVCEV model for the underlying's price is given by the stochastic differential equations (SDEs):
(1)
d
S
t
=
μ
S
t
d
t
+
σ
t
S
t
1
+
γ
d
W
t
,
σ
t
∶
=
f
(
Y
t
)
,
d
Y
t
=
α
(
m
-
Y
t
)
d
t
+
β
d
B
t
,
where μ and m are some constants, γ is greater than -
1, α and β are positive constants, and W
t and B
t are correlated Brownian motions such that d
〈
W
,
B
〉
t
=
ρ
d
t for some ρ. If ρ
>
0, then S
t may fail to be a true martingale since Y
t may go infinite. Refer to [19]. Therefore, ρ
<
0 is assumed here. It is observed in most financial markets that there is a minus correlation or leverage effect between stock price and volatility impact. However, some commodity markets show the opposite effect (the inverse leverage effect). So, the model under the negative correlation condition can apply to many financial markets but has limits for some commodity markets. Also, generally speaking, the correlation ρ may rely on time but it is assumed to be a constant for simplicity. In fact, in most real situations, it is taken to be such. We do not specify the function f but it has to satisfy a growth condition to avoid the nonexistence of moments of S
t. It is assumed in this paper that 0
<
c
1
≤
f
≤
c
2
<
∞ for some constants c
1 and c
2.From theI to formula, the solution Y of the second equation in (1) is an ergodic process given by the Ornstein-Uhlenbeck process:
(2)
Y
t
=
m
+
(
Y
0
-
m
)
e
-
α
t
+
β
∫
0
t
‍
e
-
α
(
t
-
s
)
d
B
s
and so Y
t
~
N
(
m
+
(
Y
0
-
m
)
e
-
α
t
,
(
β
2
/
2
α
)
(
1
-
e
-
2
α
t
)
), which has an invariant distribution given by N
(
m
,
β
2
/
2
α
). Later, we will use notation 〈
·
〉 for the average with respect to this invariant distribution; that is,
(3)
〈
g
〉
=
1
2
π
ν
2
∫
-
∞
+
∞
‍
g
(
y
)
e
-
(
y
-
m
)
2
/
2
ν
2
d
y
,
ν
2
≡
β
2
2
α
for arbitrary function g.Now, we take the processY
t to be a fast mean-reverting process as in Fouque et al. [20]. This means that we take into account a fast time scale volatility factor as a major driving force for the volatility movement of the underlying asset. This assumption can be justified by an empirical analysis of high-frequency S&P 500 index data. The analysis confirms that volatility is fast mean-reverting when looked at over the time scale of a derivative contract although it reverts slowly to its mean in comparison to the tick-by-tick fluctuations of the index value. See [20] for details. Also, this assumption will provide us with analytic tractability for an approximation for the price of the Asian option as also shown in [15]. So, in terms of a parameter, say ϵ, which is taken as positive and small, we assume that
(4)
α
=
1
ϵ
,
β
=
ν
2
ϵ
,
where ν
~
𝒪
(
1
) which implies that the long run magnitude of volatility fluctuations remains fixed as a constant. Then, under a risk-neutral probability measure P
*, we have
(5)
d
S
t
=
r
S
t
d
t
+
f
(
Y
t
)
S
t
1
+
γ
d
W
t
*
,
d
Y
t
=
[
1
ϵ
(
m
-
Y
t
)
-
ν
2
ϵ
Λ
(
Y
t
)
]
d
t
+
ν
2
ϵ
d
B
t
*
,
where W
t
* and B
t
* are two Brownian motions under the measure P
* whose correlation is given by d
〈
W
t
*
,
B
t
*
〉
=
ρ
d
t and Λ
(
y
) denotes the combined market price of risk given by
(6)
Λ
(
y
)
=
ρ
μ
-
r
f
(
y
)
+
λ
(
y
)
1
-
ρ
2
.
Here, λ is the market price of volatility risk, which is assumed to be a bounded function depending on y.


## 3. A PDE for Option Price

In this section, we utilize the generalization of Vecer's dimension reduction technique given by Fouque and Han [15] to derive a two-space dimensional PDE representation for the option price.In this paper, a payoff function for arithmetic average Asian options is given by(7)
h
(
1
T
∫
0
T
‍
S
t
d
t
-
K
1
S
T
-
K
2
)
,
where h is a function satisfying
(8)
h
(
α
x
)
=
α
h
(
x
)
,
for any real number α. Note that, when K
1
=
0, it becomes a payoff for fixed strike Asian options, whereas, when K
2
=
0, it becomes a payoff for floating strike Asian options. So, the option price P
(
t
,
x
,
y
;
T
,
K
1
,
K
2
) is defined by
(9)
P
(
t
,
s
,
y
;
T
,
K
1
,
K
2
)
=
𝔼
*
[
e
-
r
(
T
-
t
)
h
(
1
T
∫
0
T
‍
S
u
d
u
-
K
1
S
T
-
K
2
)
∣
S
t
=
s
,
Y
t
=
y
]
under a risk-neutral measure P
*.First, we would like to replicate the averaged process(
1
/
t
)
∫
0
t
‍
S
u
d
u with a portfolio
(10)
X
t
=
α
t
S
t
+
β
t
e
r
t
.
Here, α
t is assumed to be a nonrandom function. By the self-financing strategy, we note that
(11)
d
X
t
=
α
t
d
S
t
+
β
t
d
(
e
r
t
)
=
α
t
d
S
t
+
r
(
X
t
-
α
t
S
t
)
d
t
=
r
X
t
d
t
+
α
t
(
d
S
t
-
r
S
t
d
t
)
=
r
X
t
d
t
+
α
t
f
(
Y
t
)
S
t
1
+
γ
d
W
t
*
.
The assumption of the nonrandomness of α
t yields
(12)
d
(
e
r
(
T
-
t
)
α
t
S
t
)
=
e
r
(
T
-
t
)
α
t
(
d
S
t
-
r
S
t
d
t
)
+
e
r
(
T
-
t
)
S
t
d
α
t
.
Thus, using (11) and (12), one can obtain
(13)
d
(
e
r
(
T
-
t
)
X
t
)
=
-
r
e
r
(
T
-
t
)
X
t
d
t
+
e
r
(
T
-
t
)
d
X
t
=
e
r
(
T
-
t
)
α
t
(
d
S
t
-
r
S
t
d
t
)
=
d
(
e
r
(
T
-
t
)
α
t
S
t
)
-
e
r
(
T
-
t
)
S
t
d
α
t
.
By integration of (13), we obtain
(14)
X
T
-
e
r
T
X
0
=
∫
0
T
‍
d
(
α
t
e
r
(
T
-
t
)
S
t
)
-
∫
0
T
‍
e
r
(
T
-
t
)
S
t
d
α
t
,
which can be written as
(15)
X
T
=
-
∫
0
T
‍
e
r
(
T
-
t
)
S
t
d
α
t
+
e
r
T
X
0
+
α
T
S
T
-
α
0
e
r
T
S
0
.
For technical reason, if we choose the trading strategy α
t as
(16)
α
t
=
1
-
e
-
r
(
T
-
t
)
r
T
and the initial portfolio price X
0 as
(17)
X
0
=
x
=
α
0
S
0
-
e
-
r
T
K
2
,
then the final portfolio price X
T becomes (
1
/
T
)
∫
0
T
‍
S
t
d
t
-
K
2 and so the general payoff function of arithmetic average Asian options given by (8) becomes h
(
X
T
-
K
1
S
T
). Refer to [9, 21] or [15] for details. In terms of the portfolio process X
t, the price of an arithmetic average Asian option at time t
=
0 is given by
(18)
P
(
0
,
s
,
y
;
T
,
K
1
,
K
2
)
=
𝔼
*
[
e
-
r
T
h
(
X
T
-
K
1
S
T
)
∣
S
0
=
s
,
Y
0
=
y
]
.Next, we change the probability measureP
* into the probability measure P
~
* defined by
(19)
d
P
~
*
d
P
*
=
e
-
r
T
S
T
S
0
=
exp
[
∫
0
T
‍
f
(
Y
t
)
S
t
γ
d
W
~
t
*
-
1
2
∫
0
T
‍
f
2
(
Y
t
)
S
t
2
γ
d
t
]
.Then we are ready to obtain the option price as the solution of a partial differential equation.Theorem 1.
Letψ
t
:
=
X
t
/
S
t. Then the Asian option price at time t
=
0 can be expressed by
(20)
P
(
0
,
s
,
y
;
T
,
K
1
,
K
2
)
=
s
𝔼
~
*
[
h
(
ψ
T
-
K
1
)
∣
ψ
0
=
ψ
,
Y
0
=
y
]
under the measure P
~
*, where
(21)
ψ
=
x
s
=
1
-
e
-
r
T
r
T
-
K
2
s
e
-
r
T
.
If u
(
t
,
ψ
,
y
;
T
,
K
1
,
K
2
)
:
=
𝔼
~
*
[
h
(
ψ
T
-
K
1
)
∣
ψ
t
=
ψ
,
Y
t
=
y
], then u is given by the solution of the PDE:
(22)
u
t
+
1
2
(
α
t
-
ψ
)
2
f
2
(
y
)
s
2
γ
u
ψ
ψ
+
ρ
ν
2
ϵ
(
α
t
-
ψ
)
f
(
y
)
s
γ
u
ψ
y
+
(
1
ϵ
(
m
-
y
)
-
ν
2
ϵ
(
Λ
(
y
)
-
ρ
f
(
y
)
)
)
u
y
+
ν
2
ϵ
u
y
y
=
0
with the terminal condition
(23)
u
(
T
,
ψ
,
y
;
T
,
K
1
,
K
2
)
=
h
(
ψ
-
K
1
)
.Proof.
By the Itô formula, we obtain the following two results:(24)
d
(
S
t
-
1
)
=
1
S
t
[
(
f
2
(
Y
t
)
S
t
2
γ
-
r
)
d
t
-
f
(
Y
t
)
S
t
γ
d
W
t
*
]
,
d
X
t
=
S
t
[
r
ψ
t
d
t
+
α
t
f
(
Y
t
)
S
t
γ
d
W
t
*
]
.
Then we have(25)
d
ψ
t
=
X
t
d
(
S
t
-
1
)
+
S
t
-
1
d
X
t
+
d
X
t
d
(
S
t
-
1
)
=
ψ
t
[
f
2
(
Y
t
)
S
t
2
γ
d
t
-
f
(
Y
t
)
S
t
γ
d
W
t
*
]
-
α
t
[
f
2
(
Y
t
)
S
t
2
γ
d
t
-
f
(
Y
t
)
S
t
γ
d
W
t
*
]
=
f
(
Y
t
)
S
γ
(
α
t
-
ψ
t
)
(
d
W
t
*
-
f
(
Y
t
)
S
t
γ
d
t
)
∶
=
f
(
Y
t
)
S
γ
(
α
t
-
ψ
t
)
d
W
~
t
*
,
where
(26)
W
~
t
*
:
=
W
t
*
-
∫
0
t
‍
f
(
Y
u
)
S
u
γ
d
u
,
which is a Brownian motion under the measure P
~
* by the Girsanov theorem (cf. [22]).
From (8) and (18), we have
(27)
P
(
0
,
s
,
y
;
T
,
K
1
,
K
2
)
=
𝔼
*
[
e
-
r
T
h
(
X
T
-
K
1
S
T
)
∣
S
0
=
s
,
Y
0
=
y
]
=
s
𝔼
*
[
e
-
r
T
S
T
S
0
h
(
ψ
T
-
K
1
)
∣
ψ
0
=
ψ
,
Y
0
=
y
]
=
s
𝔼
~
*
[
h
(
ψ
T
-
K
1
)
∣
ψ
0
=
ψ
,
Y
0
=
y
]
,
where
(28)
ψ
=
x
s
=
1
-
e
-
r
T
r
T
-
K
2
s
e
-
r
T
.
Ifu
(
t
,
ψ
,
y
;
T
,
K
1
,
K
2
)
:
=
𝔼
~
*
[
h
(
ψ
T
-
K
1
)
∣
ψ
t
=
ψ
,
Y
t
=
y
], the Asian option price at time t
=
0 is given by
(29)
P
(
0
,
s
,
y
;
T
,
K
1
,
K
2
)
=
s
u
(
0
,
ψ
,
y
;
T
,
K
1
,
K
2
)
.
Then, by the Feynman-Kac formula (cf. [22]) and (25), one can obtain (22).Once the solutionP
(
0
,
s
,
y
;
T
,
K
1
,
K
2
) is determined, from the result of Fouque and Han [15], the price P
(
t
,
s
,
y
;
T
,
K
1
,
K
2
) is given by
(30)
P
(
t
,
s
,
y
;
T
,
K
1
,
K
2
)
=
T
-
t
T
P
(
0
,
s
,
y
;
T
,
K
1
,
K
2
)
=
T
-
t
T
s
u
(
0
,
ψ
,
y
;
T
,
K
1
,
K
2
)
,
where u is the solution of the PDE (22) at time t
=
0.


## 4. Multiscale Analysis

In this section, we are interested in the solution of the multiscale PDE (22) in the asymptotic form u
=
u
0
+
ϵ
u
1
+
⋯
+
ϵ
i
/
2
u
i
+
⋯. Then one can come up with single scale PDEs much easier to solve than the PDE (22) itself as follows.Substituting the asymptotic form ofu into the PDE (22) yields
(31)
1
ϵ
ℒ
0
u
0
+
1
ϵ
(
ℒ
0
u
1
+
ℒ
1
u
0
)
+
(
ℒ
0
u
2
+
ℒ
1
u
1
+
ℒ
2
u
0
)
+
ϵ
(
ℒ
0
u
3
+
ℒ
1
u
2
+
ℒ
2
u
1
)
+
⋯
+
ϵ
i
/
2
(
ℒ
0
u
i
+
2
+
ℒ
1
u
i
+
1
+
ℒ
2
u
i
)
+
⋯
=
0
with the terminal condition
(32)
∑
i
=
0
∞
‍
ϵ
i
/
2
u
i
=
h
(
ψ
-
K
1
)
,
where
(33)
ℒ
0
:
=
(
m
-
y
)
∂
∂
y
+
ν
2
∂
2
∂
y
2
,
ℒ
1
(
f
(
y
)
)
:
=
ρ
ν
2
(
α
t
-
ψ
)
f
(
y
)
s
γ
∂
2
∂
ψ
∂
y
+
2
ν
(
Λ
(
y
)
-
ρ
f
(
y
)
)
∂
∂
y
,
ℒ
2
(
f
(
y
)
)
:
=
∂
∂
t
+
1
2
(
α
t
-
ψ
)
2
f
(
y
)
2
s
2
γ
∂
2
∂
ψ
2
.Note thatℒ
0 is the infinitesimal generator of the Ornstein-Uhlenbeck process Y
t. Then we have a hierarchy of PDEs as follows:
(34)
1
ϵ
-
term
:
ℒ
0
u
0
=
0
,
1
ϵ
-
term
:
ℒ
0
u
1
+
ℒ
1
u
0
=
0
,
ϵ
-
(
i
+
2
)
/
2
-
term
:
ℒ
0
u
i
+
2
+
ℒ
1
u
i
+
1
+
ℒ
2
u
i
=
0
,
∀
i
≥
0
,
with the following terminal condition:
(35)
u
0
(
T
,
ψ
,
y
;
T
,
K
1
,
K
2
)
=
h
(
ψ
-
K
1
)
,
u
i
(
T
,
ψ
,
y
;
T
,
K
1
,
K
2
)
=
0
,
∀
i
≥
1
.In the following two theorems, we obtain PDEs for the leading order termu
0 as well as the first correction term u
1.Theorem 2.
Assume that the termu
i does not grow as much as u
i
~
e
y
2
/
2 with respect to the variable y for i
=
0,1. Then the leading order term u
0 is independent of y and satisfies the PDE problem:
(36)
ℒ
2
(
σ
¯
)
u
0
(
t
,
ψ
;
T
,
K
1
,
K
2
)
=
∂
u
0
∂
t
+
1
2
(
ψ
-
α
t
)
2
σ
¯
2
s
2
γ
∂
2
u
0
∂
ψ
2
=
0
,
u
0
(
T
,
ψ
;
T
,
K
1
,
K
2
)
=
h
(
ψ
-
K
1
)
,
where
(37)
σ
¯
=
〈
f
2
〉
=
1
2
π
ν
2
∫
-
∞
∞
‍
f
2
(
y
)
e
-
(
y
-
m
)
2
/
2
ν
2
d
y
.Proof.
From the𝒪
(
1
/
ϵ
) term in (34) and the fact that ℒ
0 is the generator of Y
t, the growth condition of u
0 with respect to the variable y leads to the fact that u
0 does not depend on y. So, it is represented as
(38)
u
0
=
u
0
(
t
,
ψ
;
T
,
K
1
,
K
2
)
without dependence on the variable y. From this fact, equation ℒ
0
u
1
+
ℒ
1
u
0
=
0 becomes ℒ
0
u
1
=
0 and so, by the same reason as for u
0, u
1 also does not depend on the variable y:
(39)
u
1
=
u
1
(
t
,
ψ
;
T
,
K
1
,
K
2
)
.
Since ℒ
1 is an operator in which every term contains the derivative with respect to the variable y, ℒ
0
u
2
+
ℒ
1
u
1
+
ℒ
2
u
0
=
0 becomes ℒ
0
u
2
+
ℒ
2
u
0
=
0. By applying the Fredholm alternative (cf. Ramm [23]) to this Poisson equation with respect to the operator ℒ
0, we have 〈
ℒ
2
u
0
〉
=
0; that is,
(40)
ℒ
2
(
σ
¯
)
u
0
=
∂
u
0
∂
t
+
1
2
(
ψ
-
α
t
)
2
σ
¯
2
s
2
θ
∂
2
u
0
∂
ψ
2
=
0
.
By adding the boundary condition (35) for u
0 to this PDE, we obtain the PDE problem (36) as desired.Next, we obtain a PDE for the first correction termu
1.Theorem 3.
Assume that the termu
i does not grow as much as u
i
~
e
y
2
/
2 with respect to the variable y for i
=
0,1. Then the first order correction term u
1 is independent of y and satisfies the PDE problem:
(41)
ℒ
2
(
σ
¯
)
u
1
(
t
,
ψ
;
T
,
K
1
,
K
2
)
=
V
2
(
α
t
-
ψ
)
2
s
2
γ
∂
2
u
0
∂
ψ
2
+
V
3
(
α
t
-
ψ
)
3
s
2
γ
∂
3
u
0
∂
ψ
3
,
u
1
(
T
,
ψ
;
T
,
K
1
,
K
2
)
=
0
,
where V
2 and V
3 are constants, respectively, given by
(42)
V
2
=
ν
2
(
〈
Λ
ϕ
′
〉
-
ρ
〈
f
ϕ
′
〉
)
,
V
3
=
-
ρ
ν
2
〈
f
ϕ
〉
.
Here, ϕ is the solution of
(43)
ℒ
0
ϕ
(
y
)
=
f
(
y
)
2
-
σ
¯
2
.Proof.
Since〈
ℒ
2
u
0
〉
=
0 and ℒ
0
u
2
+
ℒ
2
u
0
=
0, we obtain
(44)
ℒ
0
u
2
=
-
ℒ
2
u
0
=
-
(
ℒ
2
u
0
-
〈
ℒ
2
u
0
〉
)
=
-
1
2
(
ψ
-
α
t
)
2
(
f
(
y
)
2
-
σ
¯
2
)
s
2
θ
∂
2
u
0
∂
ψ
2
.
Then, from the definition of ϕ, we obtain
(45)
u
2
(
t
,
ψ
,
y
;
T
,
K
1
,
K
2
)
=
-
1
2
(
ψ
-
α
t
)
2
(
ϕ
(
y
)
+
c
(
t
,
ψ
)
)
s
2
γ
∂
2
u
0
∂
ψ
2
for some function c
(
t
,
ψ
) independent of y.
By applying the Fredholm alternative to the Poisson equationℒ
0
u
3
+
ℒ
1
u
2
+
ℒ
2
u
1
=
0, we have equation 〈
ℒ
1
u
2
+
ℒ
2
u
1
〉
=
0; that is,
(46)
ℒ
2
(
σ
¯
)
u
1
=
-
〈
ℒ
1
(
f
(
·
)
)
u
2
〉
.
Combining (45) and (46), we obtain (41).


## 5. Approximate Option Price

From the results in Section4, we formally obtain the first order approximation:
(47)
u
(
t
,
ψ
,
y
;
T
,
K
1
,
K
2
)
≈
u
0
(
t
,
ψ
;
T
,
K
1
,
K
2
)
+
ϵ
u
1
(
t
,
ψ
;
T
,
K
1
,
K
2
)
,
where u
0 and u
1 are independent of the unobservable variable y and they are given by the PDEs in Theorems 2 and 3, respectively. Note that u
0 solves a homogeneous equation with the nonzero final condition while u
1 solves a nonhomogeneous equation with the zero final condition. In terms of accuracy of the approximation (47), if the payoff function h is smooth enough, it follows straightforwardly from [20] that the approximation is of order ϵ in the pointwise convergent sense. Otherwise, it requires a regularization of the payoff as done for European vanilla options in [24]. Since the relevant argument (with an exact form of the density function for α
T) is omitted in this paper, we limit ourselves to the case of regularized payoffs here.If we defineP
0 and P
1 by
(48)
P
0
(
t
,
s
;
T
,
K
1
,
K
2
)
∶
=
T
-
t
T
s
u
0
(
0
,
ψ
;
T
,
K
1
,
K
2
)
,
P
1
(
t
,
s
;
T
,
K
1
,
K
2
)
∶
=
T
-
t
T
s
u
1
(
0
,
ψ
;
T
,
K
1
,
K
2
)
,
respectively, then, from (30), the option price P
(
t
,
s
,
y
;
T
,
K
1
,
K
2
) has an approximation given by
(49)
P
(
t
,
s
,
y
;
T
,
K
1
,
K
2
)
≈
P
0
(
t
,
s
;
T
,
K
1
,
K
2
)
+
ϵ
P
1
(
t
,
s
;
T
,
K
1
,
K
2
)
.In this section, we compute numerically the leading order priceP
0 and the first correction term P
~
1
:
=
ϵ
P
1 by using the finite difference method (the Crank-Nicolson method). The solution has the truncation error 𝒪
(
(
Δ
t
)
2
)
+
𝒪
(
(
Δ
ψ
)
2
), where Δ
t
=
0.005 and Δ
ψ
=
0.0104.Figure1 shows P
0 and P
~
1 at time t
=
0 for three different values of γ (the elasticity parameter). The parameter values used in this figure are r
=
0.06, σ
¯
=
0.5, T
=
1, K
1
=
0, K
2
=
2, V
~
2
∶
=
ϵ
V
2
=
-
0.01, and V
~
3
∶
=
ϵ
V
3
=
0.004. Note that the parameter ϵ is absorbed by V
~
2 and V
~
3. We call the price P
0
+
P
~
1 the (approximate) SVCEV price for the arithmetic average Asian option. It can reduce to the price for the well-known models. In particular, the price P
0 with γ
=
0 is the Black-Scholes price. The price P
0
+
P
~
1 with γ
=
0 corresponds to the stochastic volatility model studied by Fouque and Han [15]. Let us call it the SV price. When the leverage effect (γ
<
0) takes place, the SVCEV price becomes lower than the SV price. When the inverse leverage effect (γ
>
0) takes place, the SVCEV price becomes higher than the SV price. Regardless of the value of γ, the first order stochastic volatility correction term P
~
1 is positive and has the maximum value near s
=
K
2.The leading order termP
0 and the first order correction term P
~
1 of the price of an arithmetic average Asian option are drawn for three different values of the elasticity of variance parameter γ, where r
=
0.06, σ
¯
=
0.5, T
=
1, K
1
=
0, K
2
=
2, V
~
2
=
-
0.01, and V
~
3
=
0.004.
(a)
P
0
(b)
P
~
1Figure2 shows the correction term P
~
1 with respect to the group parameters V
~
2 and V
~
3. This figure shows that the correction term P
~
1 decreases as the parameter V
~
2 increases whereas it increases as the V
~
3 increases. Figure 2 also shows that as the parameter γ becomes larger, the correction term P
~
1 becomes larger.The correction termP
~
1 with respect to the group parameters V
~
2 and V
~
3; r
=
0.06, σ
¯
=
0.5, T
=
1, S
=
2, K
1
=
0, and K
2
=
2.
(a)
P
~
1
(b)
P
~
1


## 6. Conclusion

A frequent criticism of the stochastic volatility or local volatility models for path dependent options is that they do not produce deltas precise enough for hedging purposes. So, relevant industry experts recommend using a hybrid stochastic local volatility model of their own development for best pricing option products. See, for instance, [25]. In this paper, by transforming the path dependent problem for Asian options in the hybrid SVCEV model into a European vanilla style problem, we approximate the volatility dependent price by the nonvolatility dependent price. The approximate price provides not only a correction to the price of [15] for different values of the elasticity of variance parameter but also a correction to the CEV price based on the assumption of fast mean-reverting stochastic volatility. The elasticity of variance plays an important role in characterizing volatile markets as well as differentiating commodity markets from financial markets. Refer to [26]. The numerically solved Asian option prices show that both the leading order option price (the CEV price) and the corrected price (the approximate SVCEV price) increase as the elasticity of variance goes up. This result can provide more efficient risk hedging by choosing an appropriate elasticity parameter based on the observed market volatility or the chosen commodity.



---

*Source: 784386-2014-01-08.xml*