Bifractional Black-Scholes Model for Pricing European Options and Compound Options

Feng XU

Journal of Systems Science and Information ›› 2020, Vol. 8 ›› Issue (4) : 346-355.

PDF(167 KB)
PDF(167 KB)
Journal of Systems Science and Information ›› 2020, Vol. 8 ›› Issue (4) : 346-355. DOI: 10.21078/JSSI-2020-346-10
 

Bifractional Black-Scholes Model for Pricing European Options and Compound Options

Author information +
History +

Abstract

Recent empirical studies show that an underlying asset price process may have the property of long memory. In this paper, it is introduced the bifractional Brownian motion to capture the underlying asset of European options. Moreover, a bifractional Black-Scholes partial differential equation formulation for valuing European options based on Delta hedging strategy is proposed. Using the final condition and the method of variable substitution, the pricing formulas for the European options are derived. Furthermore, applying to risk-neutral principle, we obtain the pricing formulas for the compound options. Finally, the numerical experiments show that the parameter HK has a significant impact on the option value.

Key words

bifractional Brownian motion / compound options / long memory property / pricing model

Cite this article

Download Citations
Feng XU. Bifractional Black-Scholes Model for Pricing European Options and Compound Options. Journal of Systems Science and Information, 2020, 8(4): 346-355 https://doi.org/10.21078/JSSI-2020-346-10

1 Introduction

The traditional study of option pricing is based on the assumption that the underlying asset obeys the geometric Brownian motion[1]. However, empirical research on the return on financial assets in recent years shows that financial asset prices are not random walks, but rather fractal features such as self-similarity, long memory and non-periodic cycle. In 1989, Peters[2] proposed the fractal market hypothesis, applied R/S to analyze different capital markets, and discovered the existence of fractal structures and non-periodic cycles, that is, the price of financial assets depends not only on the current time, but also on the previous historical moment. In general, we use fractional Brownian motion to characterize this self-similarity and long memory.
However, the literature [3] points out that the fractional Brownian motion is not a semi-martingale, so the classical stochastic analysis theory cannot be directly applied to the fractional Brownian motion, and the direct application of fractional Brownian motion to the financial environment will generate arbitrage opportunities[4, 5]. Therefore, many scholars began to use modified fractional Brownian motion to describe the behavior pattern of financial asset price changes, such as mixed fractional Brownian motion[6-8] and bifractional Brownian motion[9, 10]. Since the bifractional Brownian motion is a Gaussian process that is more extensive than the fractional Brownian motion, it not only has the characteristics of self-similarity and long memory, but also is a semi-martingale under certain restricted conditions, so it can be applied to finance field. The paper [11] gave the stochastic integral of the bifractional Brownian motion and pointed out that the bifractional Brownian motion can be used to characterize the random volatility of financial assets. In the paper [12], the author used the geometric bifractional Brownian motion to capture the underlying asset of equity warrants. Moreover, a partial differential equation formulation for valuing equity warrants based on Wick integration was proposed. Meanwhile, the pricing results of different models illustrate that the financial asset had the property of long memory. Zhao[13] set the asset price followed bifractional Brownian motion, and constructed quasi-martingale method under the risk neutral measure to solve bifractional Black-Scholes model and two kinds of depressed option by the same way. Using the measure transformation method and quasi-martingale technique, the paper [14] obtained the geometric average Asian call and put option pricing formulas with fixed-price driven by bifractional Brownian motion are generalized. For other studies on bifractional Brownian motion and long memory process can see [15-17].
A compound option[18, 19] is simply an option on an option. The exercise payoff of a compound option involves the value of another option. A compound option then has two expiration dates and two strike prices. There are four main types of compound option, namely, a call on a call, a call on a put, a put on a call and a put on a put. This option has widely been applied in many financial practices. For example, the paper [20] introduced the concept of multi-stage compound options to the valuation of convertible bonds, and found that adopting the finite difference method to solve the Black-Scholes equation for each stage actually resulted in a better numerical efficiency. In addition, the article[21] studied the problem of European-style option pricing in time-changed Lévy models in the presence of compound Poisson jumps.
Our paper is organized as follows. Section 2 contains some preliminaries about bifrational Brownian motion. In section 3, a bifractional Black-Scholes equation for European option is derived by the method of Delta hedging strategy and valuation formulas for the European options by variable substitution. It is derived the pricing formula of a call on a call under bifractional Brownian motion by the method of risk neutral pricing when the parameter HK is greater than 0.5 in section 4, and it can derive the pricing formulas of other compound options in the same manner. Finally, we show that parameter HK affect risk characteristics of the European options by numerical analysis.

2 Preliminaries

In this section, we first describe some basic facts on the bifractional Brownian motion.
The so-called bifractional Brownian motion with H(0,1),K(0,1] is a mean zero Gaussian process BH,K=(BtH,K,t0) with B0H,K=0 and the covariance
E(BtH,K,BsH,K)=12K[(t2H+s2H)K|ts|2HK],s,t0,
(1)
if K=1, then BtH,1 is the fractional Brownian motion with Hurst parameter H(0,1).
Now we list some properties of the bifractional Brownian motion BtH,K by the following proposition.
● The bifractional Brownian motion is self-similar with parameter HK, i.e., for all α>0,
Law(αHKBαtH,K)=Law(BtH,K);
● For HK>12, the bifractional Brownian motion has long-range dependence in the case that
ρ(n)=Cov(BsH,KBs1H,K,Bs+nH,KBs+n1H,K)
satisfies
n=1|ρ(n)|=+;
● If 2HK=1,1<K<2 and H(0,1), the bifractional Brownian motion is a semi-martingale;
● The paths of the bifractional Brownian motion are a.s. Hölder continuous with parameter δ for any δ<HK.

3 The Bifractional Black-Scholes Equation for European Options

Since bifractional Brownian motion is a well-developed mathematical model of strongly correlated stochastic processes, it is the most efficient tool for capturing the long memory behaviour for the financial asset. firstly, we state some basic assumptions that will be used in this paper.
Consider a financial market with two primitive securities, namely, a risky asset St and a risk-free bond Mt. We will need the following assumptions. It is important to note that this article always assumes that the parameters HK are between 12 and 1.
● The governing stochastic differential equation for the price of a risky asset (such as a stock) St is given by
dSt=μStdt+σStdBtH,K,
(2)
where μ is the instantaneous expected return of St, σ is its volatility, BtH,K is the process of the bifractional Brownian motion.
● The price of a risk-free bond Mt is given by
dMt=rMtdt,M0=1,  1tT,
(3)
where r is a constant risk-free rate.
● There are no transaction costs or taxed and dividends are not paid during the time of outstanding options.
● The market is complete.
Under the risk-neutral measure, we have that
dSt=rStdt+σStdBtH,K,
and it is easy to derive that
St=S0exp(rt12σ2t2HK+σBtH,K).
(4)
Lemma 1 Considering a binary differentiable function Vt=V(t,St), St is determined by the above stochastic differential equation (2), we have
dVt=(Vt+μStVS+HKσ2t2HK1St22VS2)dt+σStVSdBtH,K.
(5)
Proof According to Taylor's formula, we have
dVt=Vtdt+VSdSt+122VS2(dSt)2+o(dtdSt),
since E(BtH,K)=t2HK, we can approximately assume that (dBtH,K)2=2HKt2HK1dt, then we have that
(dSt)2=(μStdt+σStdBtH,K)2=μ2St2(dt)2+σ2St2(dBtH,K)2+2μσSt2dtdBtH,K=2HKσ2t2HK1St2+o(dt).
(6)
Ignoring higher order infinitesimal o(t), then the proof is completed.
In the following part, the asset portfolio is constructed, and Black-Scholes partial differential equation driven by bifractional Brownian motion is derived through Delta hedging strategy, then the pricing formula of European options is obtained by the method of boundary-value conditions and variable substitution.
Theorem 1 (Bifractional Black-Scholes Equation) Suppose Ct=C(t,St) is the price of the European call option at time t, and the price of the underlying asset satisfies (2), then Ct satisfies the following partial differential equation
Ctt+rStCtSt+HKσ2t2HK1St22CtSt2rCt=0.
(7)
Proof Construct a self-financing portfolio Πt: Buy one unit a European call option and sale Δ units underlying assets, i.e., Πt=CtΔSt. According to lemma 1, we can know that
dΠt=dCtΔdSt=(Ctt+μStCtSt+HKσ2t2HK1St2Ct2St2)dt+σStCtStdBtH,KΔStμdtΔStσdBtH,K=[Ctt+μSt(CtStΔ)+HKσ2t2HK1St2Ct2St2]dt+(CtStΔ)StσdBtH,K.
(8)
Choose Δ so that the portfolio is risk-free at (t,t+dt), that is dΠt=rΠtdt.
Let Δ=CtSt, we have
(Ctt+HKσ2t2HK1St2Ct2St2)dt=r(CtCtStSt)dt.
(9)
Then (7) is obtained, the proof is completed.
Theorem 2 Assuming that the maturity date is T and the strike price is k, the price Ct of the European call option under the bifractional Brownian motion at any time t[0,T] is
C(t,St)=StΦ(d1)er(Tt)kΦ(d2),
(10)
where Φ() is the cumulative standard normal distribution function and
d1=lnStk+r(Tt)+12σ2(T2HKt2HK)σT2HKt2HK,d2=d1σT2HKt2HK.
Proof According to Theorem 1, Ct is defined by Theorem 1, and the boundary-value condition is C(T,S)=(Sk+).
Let S=ex,C=V(t,x), then we get
x=lnS,Ct=Vt,CS=1SVx,2CS2=1S2(2Vx2Vx).
(11)
Substituting (11) into (7), we get
Vt+(rHKσ2t2HK1)Vx+HKσ2t2HK12Vx2rV=0,
(12)
and the boundary-value condition becomes V(T,x)=(exk)+.
Let u(s,z)=V(t,x)er(Tt),s=12σ2(T2HKt2HK) and z=x+r(Tt)12σ2(T2HKt2HK), then we have
Vt=rer(Tt)uer(Tt)HKσ2t2HK1us+er(Tt)(r+HKσ2t2HK1)uz,
(13)
Vx=er(Tt)uz,2Vx2=er(Tt)2uz2.
(14)
Substituting (13) and (14) into (12), the following heat equation is obtained.
us=2uz2,
(15)
where the boundary-value condition is u(0,z)=(ezk)+.
According to the classical solution theory of heat equation, (15) has a unique strong solution, which is described by the following form
u(s,z)=12πs+(eξk)+e(ξz)24sdξ,
(16)
substituting the boundary-value condition into (16), we can get
u(s,z)=es+zΦ(z+2slnk2s)kΦ(zlnk2s).
(17)
By the inverse transformation, we can derive the formula of the European call option in the bifractional Black-Scholes model, and the proof is completed.
At the same way, we can obtain the formula of the European put option.
Theorem 3 Assuming that the maturity date is T, the strike price is k and St is defined by (2), then the price P(t,St) of European put option at any time t[0,T] is
P(t,St)=er(Tt)kΦ(d2)StΦ(d1),
(18)
where d1,d2 and Φ() are the same as the above.
Corollary 1 The put-call parity formula for the bifracitonal Brownian motion can be written as the following form
P(t,St)C(t,St)=er(Tt)kSt.
(19)

4 The Pricing Formulas of Compound Options

In this section we consider a compound option(a call on a call) with strike price k and maturity T1, which terminal payoff is
CCT1=max(C(T1,ST1)k,0),
(20)
where St denotes the stock price at time t(0<t<T), C(t,St) is the price value of a European call option and 0<T1<T.
Based on these assumptions, we know that at time T1 this option gives the right to buy another call option and at time T this call option gives the right to buy the stock. Now we derive the pricing formulas of this compound option.
Lemma 2 (see [22]) Let E~t[] denotes the quasi-conditional expectation with respect to the risk-neutral measure, the price at every time t[0,T] of the claim FL2(μ) is given by
F(t)=er(Tt)E~t[F].
(21)
Theorem 4 (a call on a call) The value of this compound option at time t is
CC(t,St)=StΦ2(d3+m,d2+n;ρ)ker(Tt)Φ2(d3,d2;ρ)ker(T1t)Φ(d3),
where
d3=lnStX+r(T1t)12σ2(T12HKt2HK)σT12HKt2HK,m=σT12HKt2HK,n=σT2HKt2HK,ρ=T12HKt2HKT2HKt2HK,
d2 and Φ() are the same as the above, Φ2(,;ρ) denotes the two dimensional cumulative standard normal distribution function with correlation coefficient ρ.
Proof Let C(T1,ST1) denotes the value of a European call option, according to Theorem 2, we know that
C(T1,ST1)=ST1Φ(y1)+ker(Tt)Φ(y1σT2HKT12HK),
(22)
where y1=lnST1k+r(TT1)+12σ2(T2HKT12HK)σT2HKT12HK.
Let C(T1,X)=k, then X satisfies the following equation
XΦ(y1¯)ker(Tt)Φ(y1¯σT2HKT12HK)=k,
(23)
where
y1¯=lnXk+r(TT1)+12σ2(T2HKT12HK)σT2HKT12HK,
and the exercise condition of this compound option obviously is S(T1)>X.
According to risk-neutral principle and Lemma 2, the value of this compound option is the expected present value as follows:
CC(t,St)=E~t[er(T1t)(C(T1,ST1)k)+]=er(T1t)E~t[(C(T1,ST1)k)1A]=er(T1t)E~t[C(T1,ST1)1A]ker(T1t)E~t[1A]=I1I2.
(24)
where A={ST1|ST1>X} and 1{} denotes the indicator function.
Noted that
ST=Stexp{r(Tt)12σ2(T2HKt2HK)+σ(BTH,KBtH,K)},
(25)
then we obtain that
A={ST1|ST1>X}={BT1H,KBtH,KT12HKt2HK|BT1H,KBtH,KT12HKt2HK<d3}.
(26)
So we can derive that I2=ker(T1t)Φ(d3).
Due to C(T1,ST1)=er(TT1)E~T1[(STk)1{ST>k}], where
ST=ST1exp{r(TT1)12σ2(T2HKT12HK)+σ(BTH,KBT1H,K)},
(27)
so we have that
C(T,ST1)=er(TT1)E~T1[ST1B]ker(TT1)E~T1[1B],
(28)
where
B={(ST1,ST)|ST1>X,ST>k}={BT1H,KBtH,KT12HKt2HK<d3,BTH,KBtH,KT2HKt2HK<d2},
(29)
then
I1=E~t[ST1B]ker(Tt)Φ2(d3,d2;ρ)=StΦ2(d3+m,d2+σT2HKt2HK;ρ)ker(Tt)Φ2(d3,d2;ρ).
and the proof is completed.
By the similar argument, we may obtain the results on the other compound options(a call on a put, a put on a call and a put on a put, respectively). The following theorems are given without proof.
Theorem 5 (a put on a call) The value of this compound option at time t is
PC(t,St)=StΦ2(d3+m,d2+n;ρ)+ker(Tt)Φ2(d3,d2;ρ)+ker(T1t)Φ(d3).
Theorem 6 (a call on a put) The value of this compound option at time t is
CP(t,St)=StΦ2(d3m,d2+n;ρ)ker(Tt)Φ2(d3,d2;ρ)+ker(T1t)Φ(d3).
Theorem 7 (a put on a put) The value of this compound option at time t is
PP(t,St)=StΦ2(d3+σT12HKt2HK,d2n;ρ)ker(Tt)Φ2(d3,d2;ρ)+ker(T1t)Φ(d3).

5 Numerical Simulation Analysis of Risk Characteristics

Among many risk characteristic parameters, Δ=CtSt and Γ=2CtSt2 are the two most important indexes of option pricing sensitivity. Δ is the slope of the option value, also known as the hedging ratio, which accurately defines how many units the option price will change when the underlying asset price changes by one unit. This sensitivity is what option investors are most concerned about. When constructing the risk-free arbitrage asset portfolio, the ratio of the underlying asset held by the investor to the option position is Δ. Γ is the curve of the option value line. The larger Γ is, the harder it is to avoid risk. In this paper, the Black-Scholes option pricing model is taken as an example to analyze the risk characteristics. According to formula (10), the two hedging parameters can be deduced as follows:
Δ=Φ(d1),
(30)
Γ=1StσT2HKt2HKφ(d1),
(31)
where φ(x)=12πex22.
We take a first look at the value of Δ for different parameters HK. Apparently, in the case displayed in Figure 1, an increasing parameter comes along with a increasing of the value of Δ. Noticed that the value of Δ is positive, which is consistent with the characteristics of European call options. It's easy to see that the value of Δ is negative when the option is European put option.
Figure 1 The value of Δ with varying parameter HK (r=0.03,k=95,σ=0.1,T=2,t=0.5)

Full size|PPT slide

Figure 2 indicates that Γ decreases with the increase of HK when the price of the underlying asset St is not much different from the exercise price k. When St differs greatly from k, Γ is proportional to HK. This indicates that the investors are more willing to exercise this option when the difference between St and k is large, otherwise, this option will not be exercised.
Figure 2 The value of Γ with varying parameter HK (r=0.03,k=95,σ=0.1,T=2,t=0.5)

Full size|PPT slide

6 Conclusions

A large number of scholars have shown that the price volatility of financial assets has a long memory, which leads to the improvement of option pricing model based on Brownian motion hypothesis. In this paper, the Black-Scholes option pricing formula driven by bifractional Brownian motion is solved by Delta hedging strategy. Furthermore, we study the European options and the compound options driven by the bifractional Brownian motion respectively. By comparing the fluctuation of risk characteristics under different HK values, it is found that European call options with higher long memory parameter HK are more conducive to investors to implement effective hedging strategies.

References

1
Black F, Scholes M. The pricing of options and corporate liabilities. Journal of Political Economy, 1973, 81 (3): 637- 657.
2
Peters E E. Fractal structure in the capital markets. Financial Analysts Journal, 1989, 7 (32): 434- 453.
3
Lin S J. Stochastic analysis of fractional Brownian motions. Stochastics—An International Journal of Probability and Stochastic Processes, 1995, 55 (1/2): 121- 140.
4
Bender C, Elliott R J. Arbitrage in a discrete version of the Wick-fractional Black-Scholes market. Mathematics of Operations Research, 2004, 29 (4): 935- 945.
5
Björk T, Hult H. A note on Wick products and the fractional Black-Scholes model. Finance and Stochastics, 2005, 9 (2): 197- 209.
6
Xiao W L, Zhang W G, Zhang X L, et al. Pricing model for equity warrants in a mixed fractional Brownian environment and its algorithm. Physica A: Statistical Mechanics and Its Applications, 2012, 391 (24): 6418- 6431.
7
Sun L. Pricing currency options in the mixed fractional Brownian environment. Physica A: Statistical Mechanics and Its Applications, 2013, 392 (16): 3441- 3458.
8
Shokrollahi F, Kılıçman A. Actuarial approach in a mixed fractional Brownian motion with jumps environment for pricing currency option. Advances in Difference Equations, 2015, 2015 (1): 257.
9
Lei P, Nualart D. A decomposition of the bifractional Brownian motion and some applications. Statistics and Probability Letters, 2008, 79 (5): 619- 624.
10
Russo F, Tudor C. On the bifractional Brownian motion. Stochastic Processes and their Applications, 2006, 116 (5): 830- 856.
11
Es-Sebaiy K, Tudor C. Multidimentsional bifractional Brownian motion: Itô and Tanaka's formulas. Stochastics and Dynamics, 2007, 7 (3): 365- 388.
12
Xiao W L, Zhang W G, Xu W D. Pricing equity warrants in a bifractional Brownian motion. Journal of Systems Engineering, 2013, 28 (3): 348- 354.
13
Zhao W. Research on depressed option pricing driven by bifractional Brownian motion. Journal of Nanjing University (Natural Science Edition), 2017, 40 (4): 21- 25.
14
Song R L, Li X, Wang W. Asian option pricing in double fractional Brownian model with Markovian switching. Journal of Chongqing Technology and Business University (Natural Science Edition), 2019, 36 (1): 73- 77.
15
Chen Z L, Sang L H, Hao X Z. Renormalized self-intersection local time of bifractional Brownian motion. Journal of Inequalities and Applications, 2018, 2018 (1): 1- 20.
16
Ait Ouahra M, Ouahhabi H, Sghir A. Continuity in law of some additive functionals of bifractional Brownian motion. Stochastics, 2019, 91 (4): 613- 628.
17
Qin X, Lin X, Shang Q. Fuzzy pricing of binary option based on the long memory property of financial markets. Journal of Intelligent and Fuzzy Systems, 2020, 2020 (3): 1- 12.
18
Geske R. The valuation of compound options. Journal of Financial Economics, 1979, 7 (1): 63- 81.
19
Gukhal C R. The compound option approach to American options on jump-diffusions. Journal of Economic Dynamics and Control, 2004, 28 (10): 2055- 2074.
20
Gong P, He Z W, Zhu S P. Pricing convertible bonds based on a multi-stage compound option model. Physica A: Statistical Mechanics and Its Applications, 2006, 366 (1): 449- 462.
21
Ivanov R V, Ano K. Option pricing in time-changed Lévy models with compound Possion jumps. Modern Stochastics: Theory and Applications, 2019, 6 (1): 81- 107.
22
Bender C, Sottine T, Valkeila E. Arbitrage with fractional Brownian motion. Theory of Stochastic Processes, 2006, 12 (3): 1- 12.

Funding

the Fundamental Research Funds for Suzhou Vocational University of China(SVU2018YY01)
"QINGLAN" project of SVU
PDF(167 KB)

197

Accesses

0

Citation

Detail

Sections
Recommended

/