1. Introduction
Since the famous Black–Scholes [
1] model was put forward in 1973, many scholars have begun to deal with financial investment problems in terms of the Black–Scholes model. Wang [
2] studied the Black–Scholes option pricing model based on the dynamic investment strategy. Glazyrina [
3] showed another derivation of how the normal approximation of the binomial distribution leads to the Black–Scholes formula from the binomial option pricing model. Ulyah et al., [
4] proposed a new pricing model for short-term foreign stock options by Black–Scholes theory.
It is well-known that the Black–Scholes model assumed that stock prices follow a logarithmic normal distribution and have independent increments. However, Lo [
5] and Lux [
6] have confirmed that the returns of financial assets have some fractal properties, such as the sharp peak, fat tail, and long-term memory. Therefore, the hypothesis in Black–Scholes model is inconsistent with the empirical results. In addition, the underestimation of the options based on Black–Scholes model results in the volatility smile curve. Hence, many scholars have begun to correct the Black–Scholes model. Merton [
7] proposed a jump model to describe the intermittent fluctuation of price. Hull [
8] introduced the stochastic volatility model. Necula [
9] established a new price model by using fractional Brownian motion. Since these models can well depict the characteristics of long memory or sharp peak and fat tail of the return distribution, Hubalek et al., [
10], Xiao [
11] and Gu [
12] popularized the application of these models in the financial field.
In 1988, Tsallis [
13] put forward the theory of non-extensive Tsallis entropy. Then Borland [
14] studied the fractal characteristics of Tsallis distribution from the perspective of microdynamics. Liu [
15] systematically discussed the self-similarity, non-Markovian and long dependence of this distribution. Since the non-extensive Tsallis has the fractal characteristics, it can well model the distributions of asset returns of many financial assets. In 2002, Borland [
16] began to apply Tsallis theory to the pricing problem and achieved a quiet better result. Later, some scholars used this theory to asset portfolio, option pricing and other aspects and achieved better results [
17,
18,
19,
20,
21].
In the process of option pricing, the key issue is to solve the integral
. To address this problem, Borland transformed this integral into a function in terms of
and
T by substituting
with
for the reason that
and
have identical distribution. Zhao [
22] extended this method to the average geometric Asian option pricing problem. In 2004, Borland [
23] used the Feynman–Kac equation and
ansatz operator to reconsider the integral approximation and got the price formula of European call option. Then Wang [
24] extended this method to the average geometric Asian option pricing problem. However, the above scholars did not discuss the effectiveness of the two approximations of Borland. In this paper, we used a new least-square method to reconsider the approximation problem and discussed the effectiveness of our method, by which to promote the wide application of Tsallis theory in the financial field.
The rest of the paper is organized as follows. In
Section 2, we establish the asset price model by non-extensive statistics theory. In
Section 3, we use the least-square method to solve the integral approximation problem. In
Section 4, we deduce the closed solution of European option which is similar to the Black–Scholes formula. In
Section 5, the effectiveness of our method is discussed. In
Section 6, we summarize our paper.
2. Price Model Based on the Non-Extensive Statistics
This section considers a price model which is same as Borland’s model in [
16,
23]. The characteristic of this model is that the noise of stock price fluctuation is a non-Gaussian process.
Suppose that there are two kinds of assets in the market. One is the risk-free bond with a price satisfying the equation below:
where
r is risk-free rate,
is the bond price. The other asset in the market is the stock, its price
follows
with
where
is a non-Gaussian noise source,
is a standard Brownian motion defined on the probability space
.
comes from the non-extensive entropy theory in statistical physics, which refers to the probability distribution when Tsallis entropy is maximized under the constraints of regularization and normalization.
satisfy the following nonlinear Fokker-Planck equation:
where
q is the non-extensive parameter. Plastino and Plastino [
25] have proved that the conditional probability in Equation (
4) has a form as
with
where
is the Gamma function. In addition, as
,
becomes a Gaussian process. Let
,
, we can get the following Tsallis distribution [
26]:
From Equation (
6), it is easily to verify that the mean function of
is zero and the variance function is
Since we are only interested in the process with finite variance, it is reasonable to assume .
When there is no arbitrage opportunity in the market, for a risk-neutral investor, the expected return of the asset in the future must be the risk-free rate. In other words, the fair market value should discount accordingly in the risk-neutral framework at the risk-free rate. The following Theorem deduces the formula of stock price based on Tsallis distribution in the fair market, which is the reorganization and restatement of Borland [
16].
Theorem 1. In the risk-neutral market, the stock price satisfying the stochastic differential Equation (2) iswhere is the price at time 0, and Proof. Define the discounted stock price
is
Since there is generally no arbitrage opportunity in a risk-neutral financial market, Equation (
10) is required to be a martingale process. According to Girasanov’s theorem, we can find an equivalent measure
corresponding to another noise term
that transforms Equation (
10) into a martingale. Let
and
where
Then the new measure
is related to original
by the following Radon–Nikodym derivative
Using
formula on Equation (
10), we obtain
Inserting (
11) into (
12) yields
Thus, becomes a martingale process.
Using
formula we get
which implies that
where
is defined by (
9). ☐
3. Approximate Solution Based on Least-Square Method
Since integral has the form of , it cannot be solved directly. In this section, we use the least-square method (LSM) to find a simple regression equation to approximate this integral.
As mentioned in Equation (
7), the stochastic process
has the mean of 0 and the variance of
. Suppose the standardized variable of
is
expressed as
then
follows a standardized Tsallis distribution with the mean of 0 and the variance of 1. Similarly, for any time
s, it is easily to obtain
This means that
and
are identically distributed, i.e.
It can be clearly seen that if we use the
to substitute the
in Equation (
9), then the form of
can be transformed more simple so that we can solve. Following this intuition, inserting Equation (15) into Equation (
9) yields
By observing Equation (16), we find
is a function about
and
. Hence, the basis functions are
,
and
. Suppose the approximate expression of
is
Definition 1. The estimated error E between and δ iswhere represents the mean under measure . Borland [
16] in 2002 mapped
onto the
by Equation (15) and derived an approximate expression for
.
Remark 1. The approximation of derived by Borland in 2002 is By comparing with Equation (17), it can be seen that
in Equation (17) take the following form respectively,
However, the disadvantage of the method (15) is that
is a random process not a simple distribution, the substitution of (15) will cause a large deviation of the approximation of
. Then Borland proposed a new method in 2004 taken in [
23], in which he used the Feynman–Kac equation and
ansatz to reconsider the approximate of
and got another approximation result.
Remark 2. The approximation of derived by Borland in 2004 is Similar comparison with Equation (17), it can be seen that the
in (17) are
4. European Option Price
In this section, the approximate function is used to solve the pricing problem of European options.
We first consider the pricing of European call options. A European call option has a payoff of
where
represents the price of a risky asset at maturity
T and
K is the strike price. Under the measure
, the price of this option is
with
where
represents the mean value under measure
and
is an indicator function.
Theorem 2. The price of the European call option is given bywhere Proof. By Theorem 1, the approximate expression of the stock price is
To calculate
and
in Equation (22), we should first solve the inequality
. Using Equation (24) yields
Taking logarithms of both sides of Equation (25) obtain
where
The quadratic equation
has the discriminant as
Assuming that
, there are two roots of quadratic Equation (27)
Therefore, we get the solution set of
as
. Substituting the roots into
and
, we finally get
where
Hence the value for a European call option is
☐
Corollary 1. According to the parity formula of European call-put option, the price of European put option is 5. Numerical Results
In this section, we use numerical simulation to compare LSM method and Borland’s methods from two aspects, one is the values of integral approximation, the other is the option prices.
To compare the three methods, we first generate the paths of this process by using the following Euler’s iterative formula:
Moreover, the values of
is calculate by rectangle integral principle that are expressed as
where
n in Equations (28) and (29) represents the number of nodes in the path.
refers to the time point corresponding to each node.
The steps of the numerical simulation mainly include three aspects represented as follows.
Step 1. Use the Euler’s iterative formula (28) to generate the path of .
Step 2. Use the rectangle integral principle (29) to generate .
Step 3. Perform the multiple nonlinear regression based on LSM via (17) to get the regression parameters.
Use 1000 paths, and each path has 100 nodes. The following
Table 1 reports the corresponding simulation results of regression parameter values. As can be seen from
Table 1, the value of each parameter
varies with the difference of
q.
Assuming
,
, we get the following
Figure 1 which depicts the quadratic relationship between
and
. Scatter points represent Monte Carlo simulation results. It can be seen that the approximate functions obtained by the three methods are close to each other in the place with dense scatter points, but differ greatly in the place with few scatter points. The approximation based on LSM method can better fit the sparse points distributed at both ends implying that our method is much closer to the results of Monte Carlo.
The following
Table 2 is the results of errors of
calculated by LSM and the other two methods, respectively. From
Table 2, we can see that the values of error vary with the values of
q. It is obviously to find that the
calculated with the LSM has the smallest errors in the three methods.
Before examining the numerical simulation of the call option price by three method, we first give the price formula of European call option deduced by Borland in 2002 and 2004, respectively.
Remark 3. The price of European call option based on Borland’s method in 2002 iswhere is same as Equation (19) and , take the form aswith Remark 4. The price based on Borland’s method in 2004 iswhere and , are defined byand Letting
,
,
,
and
, we use the price Formulas (23), (30) and (32) to calculate the option prices and absolute errors based on three method respectively. The absolute errors represent the absolute values of the difference of prices between each method and the Monte Carlo simulation. In
Table 3, as can be seen, the price based on LSM is closest to the result of Monte Carlo simulation. Therefore, the obtained option price based on LSM is better than the two methods of Borland’s.
Figure 2 below is an implied volatility curve based on LSM model, which is plotted as a function of strike price
K. We substitute the prices calculated by Equation (23) (
,
,
, T = 0.8 and
) into Blake-Scholes model to back out the implicit volatility. Obviously, these implied fluctuations form a smiling shape, very similar to the shape implied by real market data. Moreover, the downward sloping smile curve reproduces well-known systematic features of the volatility smile that appears when using the standard Black–Scholes to price real options.