Discrete Dynamics in Nature and Society

Volume 2014 (2014), Article ID 621895, 12 pages

http://dx.doi.org/10.1155/2014/621895

Research Article

Fractional Order Stochastic Differential Equation with Application in European Option Pricing

Qing Li,¹ Yanli Zhou,2,3 Xinquan Zhao,¹ and Xiangyu Ge¹

¹School of Statistics and Mathematics, Zhongnan University of Economics and Law, Wuhan 430073, China

²School of Finance, Zhongnan University of Economics and Law, Wuhan 430073, China

³Department of Mathematics and Statistics, Curtin University, Perth, WA 6845, Australia

Received 20 January 2014; Revised 3 July 2014; Accepted 17 July 2014; Published 13 August 2014

Academic Editor: Jianping Li

Copyright © 2014 Qing Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Memory effect is an important phenomenon in financial systems, and a number of research works have been carried out to study the long memory in the financial markets. In recent years, fractional order ordinary differential equation is used as an effective instrument for describing the memory effect in complex systems. In this paper, we establish a fractional order stochastic differential equation (FSDE) model to describe the effect of trend memory in financial pricing. We, then, derive a European option pricing formula based on the FSDE model and prove the existence of the trend memory (i.e., the mean value function) in the option pricing formula when the Hurst index is between 0.5 and 1. In addition, we make a comparison analysis between our proposed model, the classic Black-Scholes model, and the stochastic model with fractional Brownian motion. Numerical results suggest that our model leads to more accurate and lower standard deviation in the empirical study.

1. Introduction

Time series incorporating memory structure has been widely used in biological, chemical, and physical system. Memory effects also exist in financial systems. For example, the decision will be effected spontaneously by the past experience of decision makers. Plenty of financial variables with long memory effects have been found [1–4], such as the gross domestic product (GDP), interest rate, foreign exchange rates, stock price, and futures price. Garzareli et al. have proved the existence of memory effects in the stock price series by the conditional probability approach and measured the extent of long memory (autocorrelation) [5].

Memory effect is often measured by the autocorrelation function, and, recently, the Hurst index as an effective tool was introduced to measure the memory effect [6]. The Hurst index is often denoted by   . In the case of , time series has negative correlation and antipersistent behavior, which is called short-dependence memory. When , the time series has no dependence. However, in the case of , time series has positive correlation and persistent behavior, which is long-dependence memory. The persistent behavior was also called “Joseph Effect” by Mandelbrot and Wallis [7]. Cajueiro and Tabak [8, 9] have also found that memory effect exists in financial markets.

A number of researchers used fractional Brownian motion to depict the characteristic of memory. Mandelbrot and Van Ness first found that long memory effects exist in stock returns and gave the definition of fractional Brownian motion [10]. Since then, describing the memory by the fractional Brownian motion in financial market becomes more and more popular. For instance, Bȩben and Orłowski [11], Huang and Yang [12], Evertsz [13], Lo [14], and Wen et al. [15, 16] have shown that the returns are of long-term (or short-term) dependence in the markets. After Black and Scholes [17] developed the option pricing theory based on the classical stochastic differential equation, a large number of literatures studied the option price based on the fractional Brownian motion. For example, Necula [18], Rostek [19], and Hu and Øksendal [20] obtained the Black-Scholes option pricing formula under fractional Brownian motion. Ren et al. [21] have considered the option pricing model for . In the case of , the option pricing formula was studied by Wang et al. [22]. Chen et al. [23] established the mixed fractional-fractional version of Black-Scholes model with and gave the Ito’s formula correspondingly.

However, the memory effects contain not only the noise memory effect but also the trend memory effect. Stochastic differential equation with fractional Brownian motion only describes the noise memory but cannot be used to study the trend memory effect of stock price. So we will describe the trend memory process by using the fractional derivative, which is another effective instrument to describe the memory effect. In particular, fractional calculus has been successfully applied in biology, physics, chemistry, and hydrology. Recently, the concept of fractal has been extended in financial mathematics [24]. This is due to the fact that fractional integral and derivatives can depict the memory and inherent process [25]. It has been realized that fractional derivative provides an excellent mathematical instrument for the description of complex process, irregular increment, memory properties, and intermediate process [25–28].

The fractional derivative is given as below:

where is a fraction. This fractional differential equation is an appropriate mathematical approach to depict memory process of the increment. However, the fractional order derivative above only denotes the memory effect of a fixed process. Since the process in financial market has stochastic effect, we add stochastic process into fractional order ordinary differential equation. In this work, we propose a new model constructed by stochastic differential equation with fractional order. We denote the stochastic process of the asset price by fractional order stochastic differential equation as follows:

In (2), is the Hurst index, which is an exponent describing the memory of the time series, and can be calculated by the analysis approach [6]. In the special case of (i.e., ), the equation is reduced to the classic stochastic differential equation. Jumarie gave the Taylor’s series of fractional order, expressed in terms of fractional differential by using Taylor’s series of fractional order, and, hence, obtained the expression of , which involves the so-called Mittag-Leffler function [29, 30]. Momani and Odibat presented the numerical approach of differential equation of fractional order [31]. Odibat proposed algorithms to compute the functions of fractional derivative [32].

The rest of this paper is organized as follows. Section 2 gives some basic concepts and theories on the fractional order ordinary differential equations and Hurst index and then establishes the fractional order stochastic differential equation in the financial market. In Section 3, based on the proposed stochastic differential equation with fractional order derivative, we give the corresponding Ito formula under the FSDE and then derive the fractional European option pricing formula. In Section 4, we conduct the empirical analysis of fractional order formula of stock price process by using the Monte Carlo simulation method, and we also make comparison analysis of option pricing formula under FSDE with the classic option pricing formula and option pricing formula based on fractional Brownian motion. The conclusions drawn from this study are presented in Section 5.

2. Fractional Order Stochastic Differential Equation

In this section, we first give some preliminaries about the fractional order ordinary differential equation and then expand them to the field of the stochastic differential equations. Thus, based on these previous research results, we can construct the generalized fractional order stochastic differential equation.

2.1. Fractional Order Ordinary Differential Equations (FODE)

Now we introduce the definitions of fractional order integration and fractional order derivative. There exist several definitions of fractional derivatives, which are related to different applications. In our paper, we consider these two definitions, which are Riemann-Liouville integral and Caputo derivative [30].

Definition 1. is a continuous function. Its Riemann-Liouville fractional integral of order of function is defined as follows:

where is a fraction and is the Gamma function with .

Definition 2. Consider the function of Definition 1, and Caputo fractional derivative of order of function is defined as

where is a fraction, is an integer and is the value of rounded up to the nearest integer, and is the ordinary derivative of .

Based on the definitions above, the following equality holds [33, 34]:

In order to get the relations between the fractional derivative and ordinary derivative, we introduce the Taylor expansion of fractional order.

Proposition 3. Assume that the continuous function has fractional derivative of fractional order , for any positive integer at any , ; then the following equality holds:

where is the derivative of order of , which can be denoted by .

Lemma 4. Assume that , ; then,

Let be equal to 1 in (7), and take integration with respect to ; we then have the following result:

The proof of the lemma above can be found in [29].

By employing the fractional order Taylor formula and (5), we get the applications below. Given that is an integer with , the following results hold:

We then compare the two equations, (8) and (9), when ; thus, the relationship between fractional difference and finite difference is obtained as follows:

For the purpose of constructing the fractional order stochastic differential equations in this section, now we give some results of the integral with respect to in Lemma 5 presented below. Its detailed proof can be obtained in [29, 30].

Lemma 5. Let denote a continuous function; then its integral with respect to is defined by the following equalities:

where ; on making , we can have the result: .

2.2. Memory Effect and the Hurst Index

Time series is a stochastic process with recorded at the discrete times . A time series has the memory structure, if the lag period information affects the future changes. Time series displays long memory when the correlation between current and lag observations does not weaken to zero quickly over time.

Let be a stationary stochastic process with autocorrelation function , , where denotes the time lag. If , is called a long memory process; if , is called a short memory process, and, otherwise, if , for , has no memory effect. The classical approach to measure the stochastic memory process is the autocorrelation function. Now, the Hurst index is widely used as an effective substitute of the autocorrelation function to determine long-range or short-range dependence.

The memory effect can be described by the memory parameter, namely, the Hurst index. Hurst index measures the smoothness of time series based on the asymptotic behavior of the rescale of the stochastic process. A key property of memory process is self-similarity, which is denoted by the Hurst index.

Definition 6. Stochastic process is self-similar with Hurst index for any and at any time ; then we denote it by , where Hurst index describes the self-similarity of stochastic process, and represents equality of the distribution.

In the following lemma, some basic properties are given and the corresponding proofs can be obtained in [6, 8].

Lemma 7. Suppose a time series is self-similar with strictly stationary increment; then this time series has the following properties.

(1)The expectation of is and, thus, for all .(2)The covariance function , which has the following result:

(3)The autocovariance function of is given by , , where is the lag period:

(4)If , then we get the relationship between autocovariance function and Hurst index:

which means in the case of ; similarly, in the case , , and in the case , . According to the autocovariance function, we have that, in the case of , the times series exhibit short-range dependence; in the case of , the times series has no dependence, which is a perfect random walk; and in the case of , time series has long-range dependence.

The Hurst indexis usually estimated by the statistic approach. Given a stochastic process of length , we divide the time interval into contiguous subintervals of length such that . For each subinterval, the average value is .

The running sum of the accumulated deviations from the mean is given as

The range over the time period is

The standard deviation of , is

The rescaled range is , and the relationship between statistic and is

Thus, we can get the result:

where is a constant and is the Hurst index.

As a consequence, we can get the Hurst index of the observed time by linear regression:

2.3. Fractional Order Stochastic Differential Equation (FSDE)

Here, we generalize the classic stochastic differential equation to establish the fractional order stochastic differential equation based on the results presented before and then apply it to the option pricing in the next section.

Definition 8. Assuming that a financial asset price is , according to the fractional ordinary differential equation, and considering the stochastic process, we can get the FSDE as follows:

where is the drift parameter, is the diffusion parameter, is the Wiener process, , (normal distribution), and and are uncorrelated, , , .

In a special case, suppose , , and then we have the linear stochastic differential equation:

By using the results of (11), we can rewrite (22) into the following form of with respect to :

where is the first order derivative of about time .

3. European Call Option Pricing Based on FOSDE

In this section, the corresponding Ito’s formula and European call option pricing formula are derived based on the fractional order stochastic differential equation.

3.1. Ito’s Lemma Based on FSDE

Lemma 9. Assume that the stock price follows the fractional order stochastic differential equation as below:

then, the function is still an Ito stochastic process, and the following expressions hold.

When

When

Proof. According to the Ito formula, we notice that

and the discrete form of is .

In this paper, we only consider the case that . There are two reasons for this consideration: first, the Hurst index is much larger than 0 generally; second, when , , and are infinitesimal. Hence, we do not need to consider the case of .

In the case of , since , , we have

According to the Ito formula presented above, we can get

Thus, the differential form is given below:

In the case of

Similarly, we obtain the differential form as follows:

To price a European option, we first introduce Lemma 10, which connects the fractional order stochastic differential equations to the partial differential equations.

Lemma 10. is the solution of the partial differential equations:

Proof. First, make portfolios and .

In the case of ,

When , we can get the riskless asset portfolio

And because the portfolio is riskless, according to the Bellman Equation, we have , where is the riskless rate. Thus, we get the equation . Consequently, we obtain the first partial differential equation

In the case of ,

When , we can also get the riskless asset portfolio

And again because is riskless, we can get the equation

Similarly, the second partial differential equation can be obtained as below:

3.2. European Call Option Based on FSDE

Before we proceed to price the European call option, we make the assumptions as below:

(1) is the riskless rate and is a constant;(2)the exchange of the stock is continuous and the stock can be divided;(3)the tax of the stock exchange is free;(4)the bonus of the stock cannot be paid within the duration of derivatives;(5)no arbitrage exists in the market;(6)the price of stock follows a fractional order stochastic differential equation

(7)the strike price is ;(8)the maturity is ,

where is the price of the stock and is the riskless interest rate; is the volatility of the price of stock; is Hurst parameter of the stock.

In the following work, we will derive the fractional option pricing formula based on the risk-neutral assumption. If the price of underlying asset is subject to the geometric Brownian motion and the return is equal to the riskless interest rate (i.e., ), we have

In the case of , according to Ito’s Lemma 9, we can get the price of the stock as

Integrate (44) and use Lemma 5; then, we can get the solution of :

Therefore, the European call option pricing formula follows:

where

Proof. The price of the European option is given by , where is the expectation of the option price based on risk-neutral, and the price of the asset obeys the lognormal distribution:

Let ; obviously, , and the probability density function , where , and . Hence,

where

So we get the European option pricing formula as follows:

In the case of , in a similar way, according to Ito’s Lemma 9, the price of the stock is

Notice that , represents the daily logarithm returns of stock , and , is the returns of one year; thus, , and (52) can be written as below:

By integrating (52) and employing Lemma 5, we get the solution of :

Consequently, the European call option pricing formula is obtained:

where

Proof. Let , and it is obvious that , so the probability density function , where

Hence,

where

Finally, the European option pricing formula is given as below:

From the result we derived, the option price formula contains mean value function of the logarithmic returns of stock price, which is the effect of trend memory. Therefore, we proved that trend memory exists in the financial systems.

Now, we give the European call option pricing formula under the risk-neutral measure. Let the mean returns of stock be equal to the riskless rate ; by taking the expectation of the returns in case , we have , where is the riskless returns. Then, we simplify the mean value function and have ; thus, we get the option pricing formula

where

4. Results and Discussion

To explain the memory effects in financial market, we make some comparisons in this section between our proposed European option pricing model and its underlying stock price equation and the well-known classic models, such as the Black-Scholes model (Black and Scholes (1973) [17]) and Black-Scholes model under fractional Brownian motion (Necula (2002) [18], Hu and Øksendal (2003) [20]).

4.1. Comparing European Pricing Formula with Other Models

Classical Black-Scholes model: [17].

The European option pricing formula is , where

The classical Black-Scholes model was established under the assumption that the price process is Markov process and that the price process is independent and has no memory effect; however, the memory effects exist in price process.

SDE with fractional Brownian motion: [18, 20].

The European option pricing formula is , where

The fractional Brownian motion model has improved the Black-Scholes model by considering the memory effect of the asset price but only considered the memory effect of the noise.

Fractional order SDE (FSDE model): , .

In the case of , the European call option pricing formula is

where

When , the option formula is reduced to the classic option formula.

In the case of , the European call option pricing formula is

where

In the paper, our model takes a new memory effect into consideration, which is called the trend memory effect of the asset price.

4.2. Comparing the Asset Price Equation with Other Models

The underlying asset price equations to be used to make comparison are given as follows:

(1)the SDE model: [17]. The stock price equation is ;(2)the SDE model with fractional Brownian motion: [18, 20]. The stock price equation is ;(3)the fractional order SDE model: , ; when , the stock price equation is

However, when , the stock price equation is

To illustrate the proposed FSDE model, we simulate the three types of stochastic differential equations presented above by using the Monte Carlo simulation method and then make comparison of these three different models. The data used in the empirical analysis is the daily closing price index series of CSI300 index of China. The time range is from January 4, 2012, to October 27, 2012, with the initial value (the price index on January 4, 2012) and the final value . We choose the one-year bonds interest rate in China as the riskless rate and the mean yield . The Hurst parameter is , which is estimated by analysis approach. Given a 95 percent confidence interval, the simulation results are shown in Tables 1, 2, and 3.

Table 1: Simulation results of SDE model.

Simulation times
Real value
Mean value
Standard deviation
Error rate
Confidence interval

100
2445
2351.9
198.6958
3.81%
[2313, 2390.8]

1000
2445
2358.7
195.7074
3.53%
[2346.5, 2370.8]

10000
2445
2357
200.5667
3.60%
[2353.1, 2360.9]

100000
2445
2360
202.2802
3.48%
[2358.7, 2361.2]

Table 2: Simulation results of SDE with FBM.

Simulation times
Real value
Mean value
Standard deviation
Error rate
Confidence interval

100
2445
2358.2
77.4242
3.55%
[2343.0, 2373.3]

1000
2445
2360.8
76.3254
3.44%
[2356.1, 2365.6]

10000
2445
2360
78.0423
3.48%
[2358.6, 2361.6]

100000
2445
2361.2
78.5882
3.43%
[2360.8, 2361.7]

Table 3: Simulation results of FSDE.

Simulation times
Real value
Mean value
Standard deviation
Error rate
Confidence interval

100
2445
2362.7
19.0566
3.37%
[2359.1, 2366.5]

1000
2445
2368.4
21.1279
3.13%
[2367.2, 2369.8]

10000
2445
2367.0
19.7396
3.19%
[2366.7, 2367.5]

100000
2445
2369.5
20.4312
3.09%
[2369.5, 2369.7]

From the results in Tables 1–3, by using Monte Carlo simulation, we conclude that the error of our proposed FSDE model is smaller than the conventional SDE model and SDE with FBM model. If we take the simulating process 100000 times, which is large enough for the error analysis, we obtain that the error rate of SDE model is 3.48%, the SDE driven by fractional Brownian motion model is 3.43%, and the FSDE model is 3.09%, respectively.

In addition, the standard derivation of simulation is also much lower than the SDE model and SDE with FBM model, and the confidence interval is smaller than those two classic models. In the same way, when we simulate 100000 times, the standard deviation of SDE model is 202.2802, the SDE driven by fractional Brownian motion model is 78.5882, and FSDE model is 20.4312. Thus, we get the conclusion that the FSDE has about 10 times lower standard derivation than the SDE model and about 2.5 times lower standard derivation than the SDE driven by fractional Brownian motion model.

5. Conclusions and Future Research

Because the fractional order ordinary differential equations can capture the memory effect in the financial system, we established the fractional order stochastic differential equation by adding the stochastic process into the fractional ordinary differential equation. Based on this stochastic differential equation with fractional order, we apply the fractional order stochastic differential equation to the financial market. We constructed the stock price , where , is Hurst index, and derived the stock price process in the cases of and , respectively, and the European call option pricing formula under the fractional order stochastic differential equation. From the European option pricing formula, we find the trend memory in stock price process when Hurst index is between 0.5 and 1.

In addition, we made some comparisons in terms of the pricing option formula and its underlying stock price process between our proposed approach and the other two classical models. We find that the new approach leads to a better result than the classic approach and fractional Brownian motion approach when we simulate the stock prices by Monte Carlo simulation.

It would be an interesting work if we improve our model by connecting fractional ordinary differential equation with fractional Brownian motion, which can describe both the trend memory and the noise memory.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The first author acknowledges the scholarship funds for the funds for fostering one hundred best Ph.D. degree papers (2013YB1901) provided by Zhongnan University of Economics and Law and also acknowledges the support of the School of Mathematics Science in Monash University for International Joint Training for Ph.D. This research was supported by the Chinese National Social Science Foundation (10BJY104).

References

E. Panas, “Long memory and chaotic models of prices on the London metal exchange,” Resources Policy, vol. 27, no. 4, pp. 485–490, 2001.

B. J. West and S. Picozzi, “Fractional Langevin model of memory in financial time series,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, vol. 65, no. 3, Article ID 037106, 2002.

F. Wen and Z. Liu, “A copula-based correlation measure and its application in chinese stock market,” International Journal of Information Technology and Decision Making, vol. 8, no. 4, pp. 787–801, 2009.

F. Wen, X. Gong, Y. Chao, and X. Chen, “The effects of prior outcomes on risky choice: evidence from the stock market,” Mathematical Problems in Engineering, vol. 2014, Article ID 272518, 8 pages, 2014.

F. Garzareli, M. Cristelli, A. Zaccaria, and L. Pietronero, “Memory effects in stock price dynamics: evidence of technical trading,” Scientific Reports, vol. 4, article 4487, 2014.

H. E. Hurst, “Long-term storage capacity of reservoirs,” Transactions of the American Society of Civil Engineers, vol. 116, pp. 770–799, 1951.

B. B. Mandelbrot and J. R. Wallis, “Noah, Joseph and operational hydrology,” Water Resources Research, vol. 4, no. 5, pp. 909–918, 1968.

D. O. Cajueiro and B. M. Tabak, “The Hurst exponent over time: testing the assertion that emerging markets are becoming more efficient,” Physica A: Statistical Mechanics and its Applications, vol. 336, no. 3-4, pp. 521–537, 2004.

D. O. Cajueiro and B. M. Tabak, “Ranking efficiency for emerging markets,” Chaos, Solitons and Fractals, vol. 22, no. 2, pp. 349–352, 2004.

B. B. Mandelbrot and J. W. Van Ness, “Fractional brownian motions, fractional noises and applications,” Society for Industrial and Applied Mathematics, vol. 10, no. 4, pp. 422–437, 1968.

M. Bȩben and A. Orłowski, “Correlations in financial time series: established versus emerging markets,” The European Physical Journal B, vol. 20, no. 4, pp. 527–530, 2001.

B. Huang and C. W. Yang, “The fractional structure in multinational stock returns,” Applied Economics Letters, vol. 2, pp. 67–71, 1995.

C. J. G. Evertsz, “Fractional geometry of financial time series,” Fractals, vol. 3, pp. 609–616, 1995.

A. W. Lo, “Long term memory in stock market prices,” Econometrica, vol. 59, pp. 1279–1313, 1991.

F. Wen and X. Yang, “Skewness of return distribution and coefficient of risk premium,” Journal of Systems Science and Complexity, vol. 22, no. 3, pp. 360–371, 2009.

F. Wen, Z. He, and X. Chen, “Investors’ risk preference characteristics and conditional skewness,” Mathematical Problems in Engineering, vol. 2014, Article ID 814965, 14 pages, 2014.

F. Black and M. Scholes, “The pricing of option and corporate liabilities,” Journal of Political Economy, vol. 81, pp. 637–659, 1973.

C. Necula, “Option pricing in a fractional Brownian motion environment,” Social Science Research Network, 2002.

S. Rostek, Option Pricing in Fractional Brownian Markets, vol. 622 of Lecture Notes in Economics and Mathematical Systems, Academic Press, Springer, Berlin, Germany, 2009.

Y. Hu and B. Øksendal, “Fractional white noise calculus and applications to finance,” Infinite Dimensional Analysis, Quantum Probability and Related Topics, vol. 6, no. 1, pp. 1–32, 2003.

F. Ren, X. Wang, and J. Liang, “A proof for French's empirical formula on option pricing,” Chaos, Solitons and Fractals, vol. 12, no. 13, pp. 2441–2453, 2001.

X. Wang, W. Qiu, and F. Ren, “Option pricing of fractional version of the Black-Scholes model with Hurst exponent H being in (1/3, 1/2 ),” Chaos, Solitons and Fractals, vol. 12, no. 3, pp. 599–608, 2001.

J. Chen, F. Ren, and W. Qiu, “Option pricing of a mixed fractional-fractional version of the Black-Scholes model,” Chaos, Solitons and Fractals, vol. 21, no. 5, pp. 1163–1174, 2004.

W. Chen, “Nonlinear dynamics and chaos in a fractional-order financial system,” Chaos, Solitons and Fractals, vol. 36, no. 5, pp. 1305–1314, 2008.

I. Pollubny, Fractional Differential Equations, Academic Press, San Diego, Calif, USA, 1999.

S. Das, Functional Fractional Calculus for System Identification and Controls, Springer, Berlin, Germany, 2008.

R. Caponetto, G. Dongola, L. Fortuna, and I. Petras, Fractional Order Systems: Modeling and Control Applications, World Scientific, River Edge, NJ, USA, 2010.

C. A. Monje, Y. Chen, B. Vinagre, D. Xue, and V. Feliu, Fractional-Order Systems and Controls: Fundamentals and Applications, Springer, London, UK, 2010.

G. Jumarie, “On the representation of fractional Brownian motion as an integral with respect to (dt)a,” Applied Mathematics Letters, vol. 18, no. 7, pp. 739–748, 2005.

G. Jumarie, “On the solution of the stochastic differential equation of exponential growth driven by fractional Brownian motion,” Applied Mathematics Letters: An International Journal of Rapid Publication, vol. 18, no. 7, pp. 817–826, 2005.

S. Momani and Z. Odibat, “Numerical comparison of methods for solving linear differential equations of fractional order,” Chaos, Solitons and Fractals, vol. 31, no. 5, pp. 1248–1255, 2007.

Z. M. Odibat, “Computational algorithms for computing the fractional derivatives of functions,” Mathematics and Computers in Simulation, vol. 79, no. 7, pp. 2013–2020, 2009.

G. Jumarie, “Stochastic differential equations with fractional Brownian motion input,” International Journal of Systems Science: Principles and Applications of Systems and Integration, vol. 24, no. 6, pp. 1113–1131, 1993.

G. Jumarie, “Fractional Brownian motions via random walk in the complex plane and via fractional derivative. Comparison and further results on their Fokker-Planck equations,” Chaos, Solitons and Fractals, vol. 22, no. 4, pp. 907–925, 2004.