1. Introduction
An option is a contract that gives its owner the right to purchase (call option) or sell (put option) a fixed number of assets of a specific common stock at a fixed price (exercise price) on or before a specified date (expiration date) (see Vazquez [1] ). Options that can only be exercised on the expiration date are called European options, while American options can be exercised at any time before the expiration date. Option pricing is one of the key issues in financial engineering research. It is not only used for investment decision-making, but also for the entire chain of enterprise operations, such as production pricing, marketing, material and accessory supply, after-sales service, and even for macroeconomic research (see Song [2] ).
The European option pricing model is the basis for pricing other financial derivatives. Therefore, providing a reasonable pricing method for European options is of great significance. Davis et al. [3] considered the pricing of European options when the underlying asset needs to be charged transaction fees during the trading process. They proved that the value function of this problem is a unique viscous solution of a fully nonlinear quasi-variational inequality under different boundary conditions. Wu [4] provided a fuzzy pattern of Black-Scholes formula and the corresponding put-call parity relationship when the related input parameters are fuzzy values. Based on the model proposed by them, financial analysts can be allowed to select any European option price for later usage with an acceptable confidence level of belief. The fuzzy price of European options was also proposed in this study of Wu [5] using the fuzzy set theory concept of “Resolution Identity”. Under these assumptions, the European option price will become a fuzzy number. When a continuous-time Markov chain with a limited number of “economic states” modulates the dynamics of the short rate and volatility of the underlying price process, Mamon et al. [6] offered closed-form solutions for European option values. Fang et al. [7] provided a novel European option pricing method (called COS method) based on the Fourier-cosine series. By utilizing the equivalent martingale measure transformation method, Lin et al. [8] considered the pricing issue for European options for two underlying assets with delays Numerical analysis showed that there is a greater chance of a price difference between delayed options and Black Scholes options as the delay lengthens. He et al. [9] proposed a new random volatility model. By assuming that the long-term average of volatility in the Heston model is random, a closed pricing formula for European options was derived. Empirical research showed that the current model can generally obtain more accurate option prices than the Heston model. Kil and Kim [10] used the scale version of the double-mean-reverting model and got the closed-form formula for European option pricing. Nzokem [11] studied a novel approach for pricing European option, and he obtained the corresponding pricing model based on the stochastic volatility method. The above literature did not investigate which distribution of underlying securities’ returns is most realistic.
In order to price options more accurately, we must construct a reasonable underlying asset price model. Black et al. [12] assumed that the asset price follows the standard Geometric Brownian motion model and obtained a partial differential equation, called B-S equation, to describe the option price function. By solving the partial differential equation, they obtained an analytical solution for the option pricing model. Later empirical research showed that return on assets has the characteristics of “peak and fat tail”, and the standard Geometric Brownian motion model of asset price return is inconsistent with the actual characteristics of market prices (see Mandelbrot [13] , Cassidy et al. [14] , Fam [15] , Zhu and Zhang [16] ). To explain the fat tail, Lim et al. [17] studied the pricing of currency options using a generalized student distribution. The generalized t-distribution they adopt is a function of a distribution multiplied by a function of the form. Cassidy et al. [14] believed that stock returns fit the t-distribution, and then used the t-distribution to price European options.
Although the student distribution has the feature of “fat tail” and is more suitable for asset returns than the normal distribution, the difficulty of pricing options with t-distribution is that fat tail will lead to a deviation in an integral required for option pricing. In this paper, we consider a new distribution, called logarithmic truncated t-distribution, to describe stock returns, and use risk neutral valuation method (see Hull [18] ) to provide a stochastic model for stock prices, and then price European options, and then obtain a European call option pricing model based on logarithmic truncated t-distribution. Finally, we use the put-call parity relationship to find the price of European put option.
2. Stock Prices Model
As we all known, there are a close correlation between the value of European option and the underlying stock prices. A reasonable stock prices model is helpful to estimate the option value more precisely. A new probability distribution needs to be defined in order to create the stock price model.
Definition 2.1. A continuous random variable follows a truncated t-distribution, if it has a probability density function as follows:
(2.1)
where
,
are constant,
,
,
is the Gamma function, and
.
is called the probability density function of X with n degrees of freedom.
For simplification, we can rewrite the above probability density function (2.1) as follows:
(2.2)
Definition 2.2. If X follows a truncated t-distribution, then
is called be a logarithmic truncated t-distributed.
We denote the strike prices by the symbol K, and the expiration data by T. The stock price at the time t is denoted by
. Let
be the distribution function of X with n degrees of freedom. We suppose the stock price satisfies
where u is the compound rate of stock return,
is the volatility of stock return and
can be obtained by historical data.
is defined as a random variable, with the probability density function
.
(2.3)
From Equation (2.3), the stock price at time T can be written as
(2.4)
3. European Option Pricing Model
Risk neutral valuation approach can be used to price European call option. In a risk neutral world, the returns of all assets equal the interest rate r. In the risk neutral world, the stock price can be drove by Equation (3.5), which is similar to Equation (2.4).
(3.5)
where
is the compound rate of stock return in the risk neutral world, and
is the modified stock price at the time T. The expected value of
in the risk neutral world is given as follows:
The above equation can be simplified to the following form
(3.6)
where
,
,
.
In the risk neutral world, the following equation is established.
(3.7)
Comparing Equation (3.6) with Equation (3.7), we have
(3.8)
From Equation (3.8), we get
Let
be the option value at the time t. The cost of an European call option, calculated at the time of expiration T, is
. According to the risk-neutral valuation approach,
can be written as
(3.9)
From Equation (3.9), we have
(3.10)
Considering the inequality
, and using the model given by Equation (3.5), we obtain
(3.11)
Solve the inequality (3.11), so that
. Then Equation (3.10) can be rewritten as
(3.12)
If
, from Equation (3.12), the option value at time t is zero, i.e.,
. In fact, when
, we get the inequality
. That is to say, K exceeds
times
. Since the value of the European call option is approximated zero as K is large enough,
accords with the actual financial meaning.
3.1. The Case of
By Equation (2) and Equation (12), we get the European call option value at the time t, which can be calculated by using the following equation:
Simplifying the above formula, we have
(3.13)
where
,
.
3.2. The Case of
If
, substituting Equation (2.2) into Equation (3.12), we have
Let
,
, and
, and we get
(3.14)
3.3. The Case of
If
, the first integral of the Equation (3.12) can be calculated by using the following method
Then the first integral of Equation (3.12) becomes
(3.15)
where
. The second integral of Equation (3.12) becomes the sum of the three integrals as follows:
Then we get
(3.16)
where
,
. Substituting Equation (3.15) and Equation (3.16) into Equation (3.12), we obtain
(3.17)
3.4. The Case of
When
, the first integral of Equation (3.12) can be simplified as follows:
When
, the second integral of Equation (3.12) equals the integral of the probability density function
from
to
. That is
Since the two integrals above can be solved, the option value at the time t can be calculated by using the following formula:
(3.18)
3.5. Pricing for European Call and Put Options
When the parameters
and the volatility of stock prices
are given, the integrals A1, A2, A3, B1, B2, C1, C2, C3, D1, D2 and D3 are easy to be calculated. The expiration date T, the present moment t and the interest rate r always are predictable, and the values of the parameters n,
can be obtain by using data mining technology.
According to Equations (3.13) (3.14) (3.17) and (3.18), the value of the European call option at the time t can be given as follows:
(3.19)
We use put-call parity relation to find the price of the European put option with the same parameters as earlier. The put-call parity relation can be written as
(3.20)
Substitute Equation (3.19) into Equation (3.20), we obtain the value of the European put option at the time t, namely
4. Conclusions
Previous empirical tests have shown that returns do not follow a Gaussian distribution. Through data fitting, it can be found that its true distribution has a fat tail feature. In this paper, we assumed that stock returns follow a logarithmic truncated t-distribution and proposed a stock price model. By discussing the parameters, we got the mathematical model of the value of European call option. Finally, we used the put-call parity relationship to find the value of European put option.
However, we only studied the European option when the underlying stock returns follow a logarithmic truncated t-distribution. Many issues have not been resolved yet. Implied volatility parameters may change, interest rates may fluctuate, and option sensitivity may exhibit unexpected behavior. These risks may not be the “cost” of maintaining positions, but they can affect pricing and play an important role in option trading. In addition, it is also worth studying how to price other financial derivatives such as American options and exotic options under the assumption of logarithmic truncated t-distribution.
Acknowledgements
This work was supported in part by the key scientific and technological research projects in Henan province of China (Grant No. 172102210271), in part by the Fundamental Research Funds for the Universities of Henan Province (Grant No. NSFRF220442) and the Doctoral Fund of Henan Polytechnic University (Grant No. B2020-36).