# Pricing Currency Option in a Mixed Fractional Brownian Motion with Jumps Environment

Authors: Foad Shokrollahi; Adem Kılıçman
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

A new framework for pricing the European currency option is developed in the case where the spot exchange rate fellows a mixed fractional Brownian motion with jumps. The jump mixed fractional partial differential equation is obtained. Some Greeks and properties volatility are discussed. Finally the numerical simulations illustrate that our model is flexible and easy to implement.


---

## Body

## 1. Introduction

A currency option is a contract, which gives the owner the right but not the obligation to purchase or sell the indicated amount of foreign currency at a specified price within a specified period of time (American Option) or on a fixed date (European Option). Since the currency option can be used as a tool for investment and hedging, it is one of the best tools for corporations or individuals to hedge against adverse movements in exchange rates. Thus in the present work the theoretical models for pricing currency options have been carried out.Option pricing was introduced by Black-Scholes in 1973. However, Dravid et al. [1], Ho et al. [2], Toft and Reiner [3], and Kwok and Wong [4] also worked on that direction. Duan and Wei [5] indicated that option pricing by Black-Scholes model which is based on Brownian motion cannot illustrate clearly two phenomena from stock markets: first asymmetric leptokurtic features and second the volatility smile. In a work by Garman and Kohlhagen (Hereafter G-K) [6] was extended the Black-Scholes model in order to make valuation European currency option, having two fundamental features: (1) estimating the market volatility of an underlying asset generally as a function of price and time without direct reference to the specific investor characteristics such as expected yield, risk aversion measures, or utility functions; (2) self-replicating strategy or hedging. However, some researchers (see [7]) presented some evidence of the mispriced currency options by the G-K model. The significant causes of why this model is not suitable for stock markets are that the currencies are different from the stocks in main respects and geometric Brownian motion cannot resolve the conduct of currency return; see [8]. Since then, in order to overcome these problems, many systems for pricing currency options were proposed by using amendments of the G-K model, such as, Rosenberg [9], Sarwar and Krehbiel [10],and Bollen and Rasiel [11].However all these researches consider that the logarithmic returns of the exchange rate were independent and identically distributed normal random variables. Nevertheless, the experimental investigation on asset return shows that discontinuities or jumps are believed to be an indispensable element of financial stock pricing; see [12–17]. Merton, in [18], proposed a jump diffusion process with Poisson jump to match the abnormal fluctuations of stock price. Based on this theory, Kou [19], Cont and Tankov [20] also considered the problems of pricing options under a jump diffusion environment in a larger setting. Moreover, the empirical studies also demonstrated that the distributions of the logarithmic returns in the asset market generally reveal excess kurtosis with more probability mass around the origin and in the tails and less in the flanks than what would occur for normally distributed data [20]. It can be said that the properties of financial return series are nonnormal, nonindependent, and nonlinear, self-similar, with heavy tails, in both autocorrelations and cross-correlations, and volatility clustering [21–32]. Since fractional Brownian motion has two substantial features such as self-similarity and long-range dependence, thus using it is more applicable to capture behavior from financial asset [32–37]. Although, the fractional Brownian motion is neither a Markov process nor a semimartingale, we cannot apply the common stochastic calculus to analyze it. Fortunately, Xiao et al. [38] employed the Wick product rather than the pathwise product to define a fractional stochastic integral whose mean is indeed zero. This property was very convenient for both theoretical developments and practical applications. Further, in [38], it was stated that if one uses the Wick-Ito-Skorohod integral, then can obtain an arbitrage-free model. Björk and Hult [39] showed recently that utilizing fractional Brownian motion in finance does not make much economic sense because, while Wick integration leads to no arbitrage, the definition of the corresponding self-financing trading strategies is quite restrictive, see for example [40], that is the simple buying-and-hold strategy is not self-financing. Therefore, the fractional market based on Wick integrals is considered which is a beautiful mathematical construction but with restricted applicability in finance.From the above analysis one can conclude that the classical Ito calculus could not help in fractional Brownian motion and to describe an appropriate stochastic integral with respect to fractional Brownian motion is very strict. Indeed, the principal obstacle of applying fractional Brownian motion in finance is that it is not semimartingale. To capture this problem, to take into account the long memory property, and to get fluctuations form financial markets, it is suitable to apply the mixed fractional Brownian motion to take fluctuations from financial asset (see [41–44]). The mixed fractional Brownian motion is a family of Gaussian processes that is a linear combination of Brownian motion and fractional Brownian motion. It is considered as one of special class of long memory processes with Hurst parameter H
>
1
/
2. In economics the first work of using mixed fractional Brownian motion could be found in [45]. Moreover, Cheridito [45] had proved that, for H
∈
(
3
/
4
,
1
), the mixed model with dependent Brownian motion and fractional Brownian motion is equivalent to one with Brownian motion and hence it was arbitrage-free. Recent further applications have been accepted in [46]. Throughout the present work, we follow the idea in [45]; that, is we will assume that H
∈
(
3
/
4
,
1
). Moreover, some empirical studies have proved that this hypothesis is useful and applicable; see [46–51].Further, to capture jumps or discontinuities, fluctuations and to take into account the long memory property, combination of Poisson jumps and mixed fractional Brownian motion is introduced in this paper. The jump mixed fractional Brownian motion is based on the assumption that exchange rate returns are generated by a two-part stochastic process: (1) small, continuous price movements are generated by a mixed fractional Brownian motion and (2) large, infrequent price jumps are generated by a Poisson process (see [52]). This two-part process is intuitively appealing, as it is consistent with an impressive market in which major information arrives infrequently and randomly. In addition, this process may provide a description for empirically observed distributions of exchange rate changes that are skewed, leptokurtic, long memory and have fatter tails than comparable normal distributions and for the apparent nonstationarity of variance. Although, various models have been applied for pricing currency options, utilizing jump fractional Brownian motion to pricing currency options has not been investigated. Then we illustrate how to price European currency options using the G-K type model derived in a jump mixed fractional Brownian environment. The comparative results of our model and other available valuation models indicate that our model is easy to implement.The rest of this paper is organized as follows. In Section2, we briefly state the definitions and properties related to mfBm that will be used in forthcoming sections and we prove some results regarding the quasi-conditional expectation. We present an analytical pricing formula for European currency option in a mixed fractional Brownian motion with jumps environment, in Section 3. Section 4 deals with the jump mixed fractional partial differential equation and discussed some Greeks of this jump mixed model. In Section 5, we indicate how to use our model to price currency options by numerical simulations. The comparison of our jump mixed fractional Brownian motion model and traditional models is presented. Section 6 draws the concluding remarks.


## 2. Preliminaries

In this section we recall some definitions and results which we need for the rest of paper. This outcome can be found in [45, 46, 53, 54].Definition 1.
A mixed fractional Brownian motion of parametersα
,
β, and H is a linear combination of fractional Brownian motion with Hurst parameter H and Brownian motion, defined on the probability (
Ω
,
F
,
P
) for any t
∈
R
+ by
(1)
M
t
H
=
α
B
t
+
β
B
t
H
,
where B
t is a Brownian motion, B
t
H is an independent fractional Brownian motion with Hurst parameter H
∈
(
0,1
), and α and β are two real constants such that (
α
,
β
)
≠
(
0,1
).Now we list some properties in [52, 54] by the following propositions.Proposition 2.
The mixed fractional Brownian motionM
t
H satisfies the following properties: (i)
M
t
H
is a centered Gaussian process and not a Markovian one, for
H
∈
(
0,1
)
∖
1
/
2
;
(ii)
M
0
H
=
0 P-almost surely;
(iii)
the covariation function of
M
t
H
(
α
,
β
)
and
M
t
H
(
a
,
b
)
for any
t
,
s
R
+
is given by
(2)
C
o
v
(
M
t
H
,
M
s
H
)
=
α
2
(
t
∧
s
)
+
β
2
(
t
2
H
+
s
2
H
-
|
t
-
s
|
2
H
)
,
where
∧
denote the minimum of two numbers;
(iv)
the increments of
M
t
H
(
α
,
β
)
are stationary and mixed-self similar for any
h
>
0
,
(3)
M
h
t
H
(
α
,
β
)
≜
M
t
H
(
α
h
1
/
2
,
β
h
H
)
,
where
≜
means “to some law”;
(v)
the increments of
M
t
H
are positively correlated if
1
/
2
<
H
<
1,uncorrelated if
H
=
1
/
2
, and negative correlated if
0
<
H
<
1
/
2
;
(vi)
the increments of
M
t
H
are long range dependent if and only if
H
>
1
/
2
;
(vii)
for all
t
∈
R
+,one has
(4)
E
[
(
M
t
H
(
α
,
β
)
)
n
]
=
{
0
,
n
=
2
l
+
1
,
(
2
l
)
!
2
l
l
!
(
α
2
t
+
β
2
t
2
H
)
l
,
n
=
2
l
.Now, let(
Ω
,
F
,
P
) be a probability field such that B
t is Brownian motion with respect to P and B
t
H is an independent fractional Brownian motion with respect to P. Now we present some results regarding the quasi-conditional expectation that we will need for the rest of the paper (see [16]).Lemma 3.
For every0
<
t
<
T and σ
∈
C one has
(5)
E
~
t
[
e
σ
(
B
T
+
B
T
H
)
]
=
e
σ
(
B
T
+
B
T
H
)
+
(
σ
2
/
2
)
(
T
-
t
)
+
(
σ
2
/
2
)
(
T
2
H
-
t
2
H
)
,
where E
~
t denotes the quasi-conditional expectation with respect to risk-neutral measure.Proof.
See [44].Lemma 4.
Letf be a function such that E
~
t
[
f
(
B
T
,
B
T
H
)
]
<
∞. Then, for every 0
<
t
≤
T and σ
∈
C, one has
(6)
E
~
t
[
f
(
σ
B
T
+
σ
B
T
H
)
]
=
∫
R
‍
1
2
π
[
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
]
×
exp
⁡
[
-
(
x
-
σ
B
t
-
σ
B
t
H
)
2
2
[
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
]
]
f
(
x
)
d
x
.Let f
(
x
)
=
1
A; one can easily obtain the following corollary.Proof.
See [44].Corollary 5.
LetA
∈
B
(
R
). Then
(7)
E
~
t
[
1
A
(
σ
B
T
+
σ
B
T
H
)
]
=
∫
R
(
(
exp
⁡
⁡
[
-
(
x
-
σ
B
t
-
σ
B
t
H
)
2
2
[
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
]
]
)
×
(
2
π
[
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
]
)
-
1
-
(
x
-
σ
B
t
-
σ
B
t
H
)
2
σ
2
(
T
2
H
-
t
2
H
)
)
h
h
h
h
h
h
h
×
1
A
(
x
)
d
x
.
Let θ
,
w
∈
R. Now, one considers the process
(8)
Z
t
*
=
θ
B
t
*
+
w
(
B
t
H
)
*
=
θ
B
t
+
θ
2
t
+
w
B
t
H
+
w
2
t
2
H
,
0
≤
t
≤
T
.
From the Girsanov theorem, there exists a measure P
* such that Z
t
* is a new mixed fractional Brownian motion. One will denote E
t
*
[
·
] the quasi-conditional expectation with respect to P
*. Consider
(9)
X
t
=
exp
⁡
(
-
θ
B
t
-
θ
2
2
t
-
w
B
t
H
-
w
2
2
t
2
H
)
.Proof.
See [44].Lemma 6.
Letf be a function such that E
~
t
[
f
(
θ
B
t
+
w
B
t
H
)
]
≤
∞. Then for every t
≤
T, one has
(10)
E
~
t
*
[
f
(
θ
B
T
+
w
B
T
H
)
]
=
1
X
t
E
~
t
[
f
(
θ
B
T
+
w
B
T
H
)
X
T
]
.Proof.
See [44].Theorem 7.
The price at every timet
∈
[
0
,
T
] of a bounded F
T
H-measurable claim F
∈
L
2 is given by
(11)
F
t
=
e
-
r
(
T
-
t
)
E
~
t
[
F
]
,
where r represents the constants risk-free interest rate.Proof.
See [44].


## 3. Pricing Currency Options in Mixed Fractional Brownian Motion with Jumps

The aim of this section is to derive the pricing formula for European currency options. The mixed fractional Brownian motion with jump model by combining the mixed fractional Brownian motion and jump process (see [43, 52]) is obtained and some properties are analyzed.To derive the new currency option pricing formula in a jump mixed fractional market. The following hypothesis will be provided:(i)
there is no transaction costs or taxes and all securities are perfectly divisible;
(ii)
security trading is continuous;
(iii)
the short-term domestic interest rater
d and foreign interest rate r
f are known and constant thought time;
(iv)
there are no risk-free arbitrage opportunities.Now consider a jump fractional Brownian motion Black-Scholes currency market that has two investments:(a)
a money market account:(12)
d
Π
t
=
r
d
Π
t
d
t
,
wherer
d represent the domestic interest rate;
(b)
a stock whose price satisfies the following equation:(13)
d
S
t
=
S
t
(
μ
-
λ
μ
J
(
t
)
)
d
t
+
σ
S
t
d
B
t
+
σ
S
t
d
B
t
H
+
S
t
(
e
J
(
t
)
-
1
)
d
N
t
,
hhhhh
0
<
t
≤
T
,
S
0
=
S
>
0
,
whereS
t denote the spot exchange rate at time t of one unit of the foreign currency measured in the domestic currency; the drift μ and volatility σ are supposed to be constants; B
t
H is a fractional Brownian motion with Hurst parameter H
>
3
/
4; N
t is a poisson process with rate λ; J
(
t
) is jump size percent at time t  which is a sequence of independent identically distributed; and (
e
J
(
t
)
-
1
)
~
N
(
μ
J
(
t
)
,
δ
2
(
t
)
). Moreover, all three sources of randomness, the fractional Brownian motion B
t
H, the poisson process N
t, and the jump size e
J
(
t
)
-
1, are supposed to be independent.By using the Ito formula, the solution for stochastic differential equation (13) is
(14)
S
t
=
S
∏
i
=
1
N
t
‍
e
J
(
t
i
)
exp
⁡
[
1
2
σ
2
t
2
H
(
μ
-
λ
μ
J
(
t
)
)
t
+
σ
B
t
H
+
σ
B
t
h
h
h
h
h
h
h
h
h
h
h
h
h
-
1
2
σ
2
t
-
1
2
σ
2
t
2
H
]
.
Let C
(
t
,
S
t
) be the price of a European call currency option at time t with strike price K that matures at time T. Then we obtain the pricing formula for currency option by the following theorem.Theorem 8.
The price at everyt
∈
[
0
,
T
] of a European call currency option with strike price K that matures at time T is given by
(15)
C
(
t
,
S
t
)
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
ε
n
×
[
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
(
d
1
)
h
h
h
h
-
K
e
-
r
d
(
T
-
t
)
N
(
d
2
)
∏
i
=
1
n
‍
]
,
where
(16)
d
1
=
(
ln
⁡
(
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
K
)
(
r
d
-
r
f
-
λ
μ
J
(
t
)
)
(
T
-
t
)
+
1
2
σ
2
(
T
-
t
)
+
1
2
σ
2
(
T
2
H
-
t
2
H
)
e
J
(
t
i
)
K
)
×
(
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
)
-
1
,
d
2
=
d
1
-
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
,
ε
n denote the expectation operator over the distribution of ∏
i
=
1
n
‍
e
J
(
t
i
), and N
(
·
) is the cumulative normal distribution. In addition, we can calculate price of a put currency option which is made by the following corollary.Corollary 9.
The price at everyt
∈
[
0
,
T
] of a European put currency option with strike price K and maturity T is given by
(17)
P
(
t
,
S
t
)
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
ε
n
×
[
∏
i
=
1
n
‍
K
e
-
r
d
(
T
-
t
)
N
(
-
d
2
)
h
h
h
-
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
(
-
d
1
)
]
,
where
(18)
d
1
=
(
ln
⁡
(
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
K
)
(
r
d
-
r
f
-
λ
μ
J
(
t
)
)
(
T
-
t
)
&
h
h
h
+
1
2
σ
2
(
T
-
t
)
+
1
2
σ
2
(
T
2
H
-
t
2
H
)
∏
i
=
1
n
‍
e
J
(
t
i
)
K
)
h
h
h
h
h
×
(
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
)
-
1
,
d
2
=
d
1
-
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
.


## 4. Property of Pricing Formula

Assume thatV is the value of a whole portfolio of different option. The value of whole portfolio satisfies in the jump mixed fractional Black-Scholes partial differential equation that present in this theorem.Theorem 10.
The price of a currency option with a bounded payofff
(
S
t
) is given by V
(
S
t
,
t
), where V
(
S
t
,
t
) is the solution of the PDE:
(19)
∂
V
∂
t
+
1
2
σ
2
S
2
∂
2
V
∂
S
2
+
H
σ
2
t
2
H
-
1
S
2
∂
2
V
∂
S
2
hhhhh
+
(
r
d
-
r
f
-
λ
μ
J
(
t
)
)
S
∂
V
∂
S
-
r
d
V
hhhhh
+
λ
E
[
V
(
e
J
(
t
)
S
t
,
t
)
-
V
(
S
t
,
t
)
]
=
0
.
Now, we discuss the properties of jump mixed fractional Brownian motion such as Greeks which are basic tools in risk management and trading options without the knowledge of Greeks can result in high loses.Theorem 11.
SupposeC
=
C
(
t
,
S
t
) be the price of European call currency option the Greeks are given by
(20)
Δ
=
∂
C
∂
S
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
ε
n
h
h
h
h
×
[
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
(
d
1
)
]
.
∇
=
∂
C
∂
K
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
h
h
h
h
×
ε
n
[
-
e
-
r
d
(
T
-
t
)
N
(
d
2
)
]
.
ρ
r
d
=
∂
C
∂
r
d
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
h
h
h
h
×
ε
n
[
K
(
T
-
t
)
e
-
r
d
(
T
-
t
)
]
N
(
d
2
)
.
ρ
r
f
=
∂
C
∂
r
f
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
h
h
h
h
×
ε
n
[
-
S
t
(
T
-
t
)
∏
i
=
1
n
‍
e
J
(
t
i
)
h
h
h
h
h
h
h
h
h
×
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
(
d
1
)
∏
i
=
1
n
‍
]
.
Θ
=
∂
C
∂
t
=
∑
n
=
0
∞
(
e
-
λ
(
T
-
t
)
λ
n
+
1
(
T
-
t
)
n
h
h
h
h
h
h
h
h
h
h
h
h
h
-
n
e
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
-
1
)
(
n
!
)
-
1
ε
n
h
h
h
h
h
×
[
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
(
d
1
)
]
h
h
h
h
h
-
K
e
-
r
d
(
T
-
t
)
N
(
d
2
)
+
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
ε
n
h
h
h
h
h
×
[
(
r
f
+
λ
μ
J
(
t
)
)
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
(
d
1
)
]
h
h
h
h
h
-
r
d
K
e
-
r
d
(
T
-
t
)
N
(
d
2
)
h
h
h
h
h
-
1
2
σ
2
+
2
H
σ
2
t
2
H
-
1
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
h
h
h
h
h
×
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
′
(
d
1
)
.
Γ
=
∂
2
C
∂
S
2
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
h
h
h
h
h
×
ε
n
[
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
S
t
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
N
′
(
d
1
)
]
,
ν
σ
=
∂
C
∂
σ
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
ε
n
h
h
h
h
h
×
[
σ
(
T
-
t
)
+
σ
(
T
2
H
-
t
2
H
)
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
S
t
∏
i
=
1
n
‍
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
h
h
h
h
h
h
h
h
h
×
σ
(
T
-
t
)
+
σ
(
T
2
H
-
t
2
H
)
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
N
′
(
d
1
)
]
.
The following theorem represents the influence of Hurst parameter H in jump mixed fractional Brownian motion.Theorem 12.
The influence of the Hurst parameter can be written as(21)
∂
C
∂
H
=
∑
n
=
0
∞
‍
e
-
λ
(
T
-
t
)
λ
n
(
T
-
t
)
n
n
!
ε
n
×
[
S
t
σ
2
(
T
2
H
ln
⁡
T
-
t
2
H
ln
⁡
t
)
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
h
h
h
×
∏
i
=
1
n
e
J
(
t
i
)
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
N
′
(
d
1
)
σ
2
(
T
2
H
ln
⁡
T
-
t
2
H
ln
⁡
t
)
σ
2
(
T
-
t
)
+
σ
2
(
T
2
H
-
t
2
H
)
]
.Suppose the forward price of an exchange rate is as(22)
f
≜
S
(
t
)
∏
i
=
1
n
‍
e
J
(
t
i
)
e
(
r
d
-
r
f
-
λ
μ
J
(
t
)
)
(
T
-
t
)
,
d
±
=
ln
⁡
(
f
/
K
)
±
(
σ
2
/
2
)
(
T
2
H
-
t
2
H
)
±
(
σ
2
/
2
)
(
T
-
t
)
σ
T
2
H
-
t
2
H
+
T
-
t
.Remark 13.
The relationship of call-put parity satisfies as(23)
C
(
t
,
S
t
)
-
P
(
t
,
S
t
)
=
S
t
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
-
K
e
-
r
d
(
T
-
t
)
.This is just more complicated way to write the trivial equationS
t
=
S
t
+
-
S
t
-.Remark 14.
The relationship of put-call parity can be written as(24)
∂
C
(
t
,
S
t
)
∂
S
t
-
∂
P
(
t
,
S
t
)
∂
S
t
=
e
-
(
r
f
+
λ
μ
J
(
t
)
)
(
T
-
t
)
.


## 5. Simulation Studies

This section indicate how to performance our jump mixed fractional Brownian motion model and to illustrate the effects of jump parameters of our pricing model. For these aims, we report on two sets of numerical experiments. In the first set, we compare theoretical prices of some assumptive options among of the following models: the G-K, the pure fractional Brownian motion (hereafter PFBM), the pure mixed fractional Brownian motion (hereafter PMFBM), and our jump mixed fractional Brownian motion (hereafter JMFBM). These tests will not be based on empirical data, but they will consist of some simulations of different pricing models with some chosen parameters of some simulations of different pricing models with some chosen parameters that will not be based on empirical data. The second set shows the influence of different parameters of JMFBM. The code line is written in MATLAB.


### 5.1. Comparison of Option Prices

Now, for description of discrepancies among these methods: the G-K, the PFBM, the PMFBM, and our JMFBM methods, we present the theoretical prices of some assumptive options applying several models. All parameters which are needed for computing the assumptive currency call options present by the G-K, PFBM, PMFBM, and JMFBM methods are present in Table1, respectively. The fourth and fifth line, which has low and high jump parameters, respectively, provides the parameters for calculating the prices by the JMFBM method. Tables 2 and 3 present the prices computed by different models, where S denotes the price of exchange rate; P
G-K denotes the prices calculated by the G-K method; P
P-F shows the prices calculated by the PFBM; P
P-MF indicates the prices simulated corresponding to PMFBM; and P
J-MF shows the prices calculated by the JMFBM.Table 1
The valuation of chosen parameters used in these models.
Model type
r
d
r
f
σ
K
H
J
λ
μ
J
δ
G-K
0.0210
0.0320
0.1050
100
—
—
—
—
—
PFBM
0.0210
0.0320
0.1050
100
0.76
—
—
—
—
PMFBM
0.0210
0.0320
0.1050
100
0.76
—
—
—
—
JMFBM
0.0210
0.0320
0.1050
100
0.76
0.008
2.3
0.00065
0.0015
JMFBM
0.0210
0.0320
0.1050
100
0.76
0.0031
7.3
−0.00065
0.0018Table 2
Results by different pricing models.
S
Low time to maturity,T
=
0.5, t
=
0
P
G-K
P
P
-
F
P
P
-
MF
P
J
-
MF
*
,
a
P
J
-
MF
*
,
b
80
0.0019
0.0001
0.0259
0.0348
0.0399
85
0.0275
0.0058
0.1469
0.1865
0.2094
90
0.2068
0.0854
0.5664
0.6857
0.7550
95
0.9124
0.5814
1.6014
1.8636
2.0158
98
1.8093
1.3654
2.6454
3.0172
3.2325
100
2.6590
2.1788
3.5415
3.9912
4.2510
102
3.7188
3.2429
4.5993
5.1273
5.4313
110
9.7210
9.5097
10.2519
11.0509
11.5015
120
19.1629
19.1441
19.2787
20.2763
20.8171
130
28.9816
28.9811
28.9944
30.1055
30.6812
140
38.8223
38.8223
38.8232
40.0284
40.6266
*The maximum number of iteration is 100.
aThe parameters for calculations are in the third line of Table 1.
bThe parameters for calculations are in the forth line of Table 1.Table 3
Results by different pricing models.
S
Low time to maturity,T
=
2, t
=
0
P
G-K
P
P
-
F
P
P
-
MF
P
J
-
MF
*
,
a
P
J
-
MF
*
,
b
80
0.2589
0.5780
1.4574
1.9778
2.3126
85
0.6657
1.2020
2.4403
3.2183
3.7063
90
1.4466
2.2186
3.7983
4.8833
5.5476
95
2.7384
3.7120
5.5627
6.9919
7.8466
98
3.7991
4.8557
6.8209
8.4684
9.4401
100
4.6315
5.7272
7.7420
9.5380
10.5878
102
5.5643
6.6810
8.7277
10.6735
11.8008
110
10.2397
11.2873
13.2787
15.8219
17.2446
120
17.7735
18.5190
20.1458
23.3889
25.1289
130
26.4100
26.8349
27.9874
31.8470
33.8372
140
35.5361
35.7410
36.4750
40.8680
43.9468
*The maximum number of iteration is 100.
aThe parameters for calculations are in the third line of Table 1.
bThe parameters for calculations are in the forth line of Table 1.By comparing columnsP
G-K, P
P-F, P
P-MF, and P
J-MF
*
,
a in Tables 2 and 3 for the small maturity cases, we can conclude that the option prices given by four assessment methods almost are the same. The main reason in this case is that jump parameters are very small. Furthermore, the prices obtained by MFBM and JMFBM with small parameters in high and low expiration are close, since the jump parameters are very small. Also, Tables 2 and 3 show that the parameters in column P
J-MF
*
,
a are greater than prices which is obtained by the G-K, PFBM, and MFBM. Next, the prices of assumptive options are investigated in large jump parameters case. Whenever the expiration time increases, the discrepancy among values calculated by JMFBM model and other models increases in large jump parameters case. Tables 2 and 3 show this result by comparing columns P
G-K, P
P-F, P
P-MF, and P
J-MF
*
,
b for the small and large time to maturity cases. The findings from the columns P
G-K, P
P-F, P
J-MF
*
,
b, and P
P-MF with S
∈
[
80,140
] indicate that the amount of discrepancy proportion in prices is larger for out-of-the-money options in the large expiration time case and this discrepancy proportion declines with the increase the exchange rate.


### 5.2. The Influence of Parameters

In what follows, the value of call currency options is displayed by using JMFBM for different parameters. We just investigate the out-of-the-money case. In fact, the other cases in-the-money and at-the-money can also be consideres by applying the same method. Now, we consider the values of our jump mixed fractional Brownian motion currency options for various Hurst parametersH and then investigate the values for different jump parameters. Figure 1 shows the values of call currency option against its parameters, H, λ, and μ
J. The defaulting parameters are S
=
90, K
=
100, r
d
=
0.0210, r
f
=
0.320, σ
=
0.1050, t
=
0, T
=
2, H
=
0.76, J
=
0.0031, λ
=
7.3, μ
J
=
-
0.00065, and δ
=
0.0018. No wonder that Figure 1 shows that (1) the option value is a decreasing function of μ
J, δ and (2) increasing parameter of H and λ comes along with an increase of the option value.Figure 1
Value of the currency call option.Now, we compare the three prices derived from the G-K, PFBM, and JMFBM models for out-of-the-money. These simulation parameters are selected:S
=
89.8, r
d
=
0.0210, r
f
=
0.0320, σ
=
0.1050, H
=
0.76, J
=
0.008, λ
=
2.3, μ
j
=
0.00065, and t
=
0.1 with different exercise price K
∈
[
90,120
] and expiration T
∈
[
0.11,10
]. For out-of-the-money Figure 2 indicates the option price discrepancies by the G-K, PMFBM, and JMFBM methods, respectively. From the figures it is concluded that the prices by JMFBM are better fitted to the G-K model. Therefore, we can see that our JMFBM seems reasonable.Figure 2
Relative difference among the G-K model, PMFBM, and JMFBM for out-of-the-money.


## 5.1. Comparison of Option Prices

Now, for description of discrepancies among these methods: the G-K, the PFBM, the PMFBM, and our JMFBM methods, we present the theoretical prices of some assumptive options applying several models. All parameters which are needed for computing the assumptive currency call options present by the G-K, PFBM, PMFBM, and JMFBM methods are present in Table1, respectively. The fourth and fifth line, which has low and high jump parameters, respectively, provides the parameters for calculating the prices by the JMFBM method. Tables 2 and 3 present the prices computed by different models, where S denotes the price of exchange rate; P
G-K denotes the prices calculated by the G-K method; P
P-F shows the prices calculated by the PFBM; P
P-MF indicates the prices simulated corresponding to PMFBM; and P
J-MF shows the prices calculated by the JMFBM.Table 1
The valuation of chosen parameters used in these models.
Model type
r
d
r
f
σ
K
H
J
λ
μ
J
δ
G-K
0.0210
0.0320
0.1050
100
—
—
—
—
—
PFBM
0.0210
0.0320
0.1050
100
0.76
—
—
—
—
PMFBM
0.0210
0.0320
0.1050
100
0.76
—
—
—
—
JMFBM
0.0210
0.0320
0.1050
100
0.76
0.008
2.3
0.00065
0.0015
JMFBM
0.0210
0.0320
0.1050
100
0.76
0.0031
7.3
−0.00065
0.0018Table 2
Results by different pricing models.
S
Low time to maturity,T
=
0.5, t
=
0
P
G-K
P
P
-
F
P
P
-
MF
P
J
-
MF
*
,
a
P
J
-
MF
*
,
b
80
0.0019
0.0001
0.0259
0.0348
0.0399
85
0.0275
0.0058
0.1469
0.1865
0.2094
90
0.2068
0.0854
0.5664
0.6857
0.7550
95
0.9124
0.5814
1.6014
1.8636
2.0158
98
1.8093
1.3654
2.6454
3.0172
3.2325
100
2.6590
2.1788
3.5415
3.9912
4.2510
102
3.7188
3.2429
4.5993
5.1273
5.4313
110
9.7210
9.5097
10.2519
11.0509
11.5015
120
19.1629
19.1441
19.2787
20.2763
20.8171
130
28.9816
28.9811
28.9944
30.1055
30.6812
140
38.8223
38.8223
38.8232
40.0284
40.6266
*The maximum number of iteration is 100.
aThe parameters for calculations are in the third line of Table 1.
bThe parameters for calculations are in the forth line of Table 1.Table 3
Results by different pricing models.
S
Low time to maturity,T
=
2, t
=
0
P
G-K
P
P
-
F
P
P
-
MF
P
J
-
MF
*
,
a
P
J
-
MF
*
,
b
80
0.2589
0.5780
1.4574
1.9778
2.3126
85
0.6657
1.2020
2.4403
3.2183
3.7063
90
1.4466
2.2186
3.7983
4.8833
5.5476
95
2.7384
3.7120
5.5627
6.9919
7.8466
98
3.7991
4.8557
6.8209
8.4684
9.4401
100
4.6315
5.7272
7.7420
9.5380
10.5878
102
5.5643
6.6810
8.7277
10.6735
11.8008
110
10.2397
11.2873
13.2787
15.8219
17.2446
120
17.7735
18.5190
20.1458
23.3889
25.1289
130
26.4100
26.8349
27.9874
31.8470
33.8372
140
35.5361
35.7410
36.4750
40.8680
43.9468
*The maximum number of iteration is 100.
aThe parameters for calculations are in the third line of Table 1.
bThe parameters for calculations are in the forth line of Table 1.By comparing columnsP
G-K, P
P-F, P
P-MF, and P
J-MF
*
,
a in Tables 2 and 3 for the small maturity cases, we can conclude that the option prices given by four assessment methods almost are the same. The main reason in this case is that jump parameters are very small. Furthermore, the prices obtained by MFBM and JMFBM with small parameters in high and low expiration are close, since the jump parameters are very small. Also, Tables 2 and 3 show that the parameters in column P
J-MF
*
,
a are greater than prices which is obtained by the G-K, PFBM, and MFBM. Next, the prices of assumptive options are investigated in large jump parameters case. Whenever the expiration time increases, the discrepancy among values calculated by JMFBM model and other models increases in large jump parameters case. Tables 2 and 3 show this result by comparing columns P
G-K, P
P-F, P
P-MF, and P
J-MF
*
,
b for the small and large time to maturity cases. The findings from the columns P
G-K, P
P-F, P
J-MF
*
,
b, and P
P-MF with S
∈
[
80,140
] indicate that the amount of discrepancy proportion in prices is larger for out-of-the-money options in the large expiration time case and this discrepancy proportion declines with the increase the exchange rate.


## 5.2. The Influence of Parameters

In what follows, the value of call currency options is displayed by using JMFBM for different parameters. We just investigate the out-of-the-money case. In fact, the other cases in-the-money and at-the-money can also be consideres by applying the same method. Now, we consider the values of our jump mixed fractional Brownian motion currency options for various Hurst parametersH and then investigate the values for different jump parameters. Figure 1 shows the values of call currency option against its parameters, H, λ, and μ
J. The defaulting parameters are S
=
90, K
=
100, r
d
=
0.0210, r
f
=
0.320, σ
=
0.1050, t
=
0, T
=
2, H
=
0.76, J
=
0.0031, λ
=
7.3, μ
J
=
-
0.00065, and δ
=
0.0018. No wonder that Figure 1 shows that (1) the option value is a decreasing function of μ
J, δ and (2) increasing parameter of H and λ comes along with an increase of the option value.Figure 1
Value of the currency call option.Now, we compare the three prices derived from the G-K, PFBM, and JMFBM models for out-of-the-money. These simulation parameters are selected:S
=
89.8, r
d
=
0.0210, r
f
=
0.0320, σ
=
0.1050, H
=
0.76, J
=
0.008, λ
=
2.3, μ
j
=
0.00065, and t
=
0.1 with different exercise price K
∈
[
90,120
] and expiration T
∈
[
0.11,10
]. For out-of-the-money Figure 2 indicates the option price discrepancies by the G-K, PMFBM, and JMFBM methods, respectively. From the figures it is concluded that the prices by JMFBM are better fitted to the G-K model. Therefore, we can see that our JMFBM seems reasonable.Figure 2
Relative difference among the G-K model, PMFBM, and JMFBM for out-of-the-money.


## 6. Conclusions

This paper compounds the jump process and mixed fractional Brownian motion. Furthermore, some particular features of currency pricing formula are also indicated. Furthermore, experimental results obtained in Section5 indicate that when the expiration date is great sufficient the JMFBM model is essential to use. However, our jump mixed fractional Brownian motion model is simple to employ and has the especial properties to describe the unusual movements of financial markets. In conclusion, since mixed fractional Brownian motion is a suitable mathematical model of profoundly correlated random processes and jumps are undeniable component in financial markets, jump mixed fractional Brownian model will be a more appropriate method for pricing currency options.



---

*Source: 858210-2014-04-09.xml*