Numerical Method for a System of PIDEs Arising in American Contingent Claims under FMLS Model with Jump Diffusion and Regime-Switching Process

Fan, Congyin; Chen, Peimin

doi:https://doi.org/10.1155/2021/8818876

Discrete Dynamics in Nature and Society

On this page

Abstract Introduction Discussion Conclusions Appendix Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 8818876 | https://doi.org/10.1155/2021/8818876

Numerical Method for a System of PIDEs Arising in American Contingent Claims under FMLS Model with Jump Diffusion and Regime-Switching Process

Congyin Fan¹and Peimin Chen²

Academic Editor: Giancarlo Consolo

Received25 Aug 2020

Accepted20 Jun 2021

Published09 Jul 2021

Abstract

This paper investigates a numerical method for solving fractional partial integro-differential equations (FPIDEs) arising in American Contingent Claims, which follow finite moment log-stable process (FMLS) with jump diffusion and regime switching. Mathematically, the prices of American Contingent Claims satisfy a system of problems with free-boundary values, where is the number of regimes of the market. In addition, an optimal exercise boundary is needed to setup with each regime. Therefore, a fully implicit scheme based on the penalty term is arranged. In the end, numerical examples are carried out to verify the obtained theoretical results, and the impacts of state variables in our model on the optimal exercise boundary of American Contingent Claims are analyzed.

1. Introduction

More and more researchers pay much attention to the pricing problem of American Contingent Claims since the pricing model is free-boundary problem and it is of great academic value to solve this kind of model. Hence, we consider the stock loan, which is a kind of famous American contingent claim, in this paper. In fact, Xia and Zhou had studied the stock loan with infinite maturity based on the Black–Scholes framework in 2007 [1]. Following this contribution, many scholars also consider the same topic by the Black–Scholes model under other different conditions. For instance, Chen et al. studied the stock loan with finite maturity under the framework that risk-free interest rate follows the Rednleman–Bartter model without drift term [2]. Cai et al. investigated the value and optimal redemption price of stock loan with infinite and finite maturity under the framework of the hyperexponential jump diffusion model, respectively [3]. Liang et al. provided the formulae of stock loan with automatic termination clause, cap, and margin [4]. Prager and Zhang [5] considered the finite horizon loan valuation under various stock models that include classical aerometric Brownian motion, mean-reverting, and two-state regime switching with both mean-reverting and geometric Brownian motion states.

Most models in literatures use a continuous-time Brownian motion as their primary source of stochasticity. These stock models include geometric Brownian motion with stochastic interest rate and diffusion with hyperexponential jump. However, even if considering the actual phenomenon, the continuous-time Brownian motion does not fully reflect the stochastic nature of financial markets. First of all, empirical evidence has showed that the distribution of asset return can appear “leptokurtic distribution” (see, e.g., [6–8]), which has a higher peak and two heavier tails compared with the standard normal distribution. Secondly, from an economic perspective, regime-switching behavior captures the changing preferences and beliefs of investors concerning stock prices as the state of financial market changes. The presence of regime-switching dynamics in financial market has been well acknowledged. A frequently observed phenomenon is that a transition in a real business cycle, expansion, and contraction usually leads to significant changes of stock returns, interest rates, and financial indices. These changes exhibit certain cyclic or periodic patterns. Consequently, there is a need for more realistic models that better capture asymmetric distribution, phenomenon of jumps, and the dynamic changes of risk asset. Naturally, the finite moment log-stable model with jump diffusion and regime switching (FMLSJR) [9] is employed to capture dynamical process of underlying asset price.

Under framework of FMLSJR, the values of stock loan follow a -dimension coupled partial integral-differential equations (PIDEs). And there are -free boundaries in the pricing model. Associated with each regime, there is a function of optimal redemption prices and a space-fractional partial integral-differential equation. And for the fractional model, the problem of how to solve it has attracted the interest of many researchers. For instance, Jumarie [10] presented the solutions of time-fractional and space-fractional BCS models about the European option by Fourier transforms and Mittag-Leffler function. Liang et al. [4] used the Laplace transform technique to obtain a completed pricing formula of option in the case of the bifractional Black–Scholes–Merton model. Following this work, Kumar et al. [11] investigated the time-fractional BS European option pricing model and obtained an analytical solution. Under the framework of the FMLS model, Chen et al. [12] considered the analytical solution of European option by using a Fourier integral transform and Fox functions. According to the literature review above, it appears that the favored methods used to consider the analytical solution of the fractional models of European option are via integral transform methods. The solutions obtained by means of these methods usually take the form of a convolution of some special functions, which make it difficult to compute. Moreover, these integral transform methods may not be effective for the American Contingent Claims. Hence, studying the numerical approximate solutions of these models appears to be a very practical and important research objective.

There are a lot of recent publications investigating the numerical method for the fractional models of European and American Contingent Claims. For example, Cartea and del-Castillo-Negrete [13] employed the shifted Grawld–Letnik scheme to discrete the fractional derivatives and set an implicit scheme to solve the pricing model of exotic option. Marom and Momoniat [14] presented a comparison of numerical solutions of the FMLS, KoBol, and CGMY models. Li [15] considered a difference scheme with first-order accuracy for a space-fractional BS model. By using a wavelet technique, G. Hariharan [16] gave a numerical scheme based on the wavelet technique to solve the time-fractional BS model arising from European option pricing problem. Chen [17] investigated a second-order finite difference method for the one-dimensional FMLS model governing the valuation of European options and a power penalty method for a space-fractional differential linear complementarity problem arising in the valuation of American options on a single asset and then extended the above methods to the corresponding two dimensional models arising in the valuation of European and American options on two assets. The author also gave the detailed numerical analysis of the methods. Zhang et al. [18] considered the European double barrier option in the case of tempered fractional Black–Scholes equation, and they employed a fast biconjugate gradient stabilized method to solve the linear system. Fan et al. [19] set a fully implicit scheme with first-order accuracy based on the penalty function for the stock loan pricing model, and their method can be used to other pricing model of American Contingent Claims. Zhou et al. [9] investigated the iterative Laplace transform methods for system of fractional PDEs and PIDEs arising in option pricing.

However, from the existing literature, it appears that research on the numerical simulation of the coupled fractional differential equations for American Contingent Claims is also relatively limited. So, in this paper, we give a numerical scheme based on penalty term and present a technique to solve the coupled system.

There are two contributions of this paper. First of all, we prove the lower bound of stock loan values at the case of fractional order with jumps and mechanism transfer. Especially, the involvement of mechanism transfer increases the difficulty of proof. Because of the existence of mechanism transfer term, the objective equation in the model is coupled by partial differential equations. A skill is employed to solve the nonlinear equation in this paper to avoid repeated iteration. Secondly, we present a technique to solve the final coupled system, which makes the numerical results more accurate than the results calculated by the iterative algorithm.

The rest of the paper is organized as follows. In Section 2, the completed pricing model is given. And we introduce our numerical method by extending the penalty method and prove the effectiveness of this method in Section 3. The numerical examples and analysis of parameter impact are displayed in Section 4, and we give the conclusions in the final section.

2. The Mathematical Model

Let be a complete filtering probability space, where and is a fixed time expiration. Assume that there exists an equivalent martingale measure , under which the dynamics of the logarithmic price of the risky asset are given by the following stochastic differential equation as in [9]:where is a continuous-time Markov chain with states , . Assume that at each state , the risk-free interest rate , , and the convexity adjustment . The term is the maximally skewed Lévy process, and the tail index . is a Poisson process, which is characterized by a jump intensity parameter . is a sequence of independent and identically distributed hyperexponential random variables with probability density function as follows:where , and , are the probabilities of different kinds of positive and negative jumps, respectively. Moreover, it satisfies . Similarly, the parameters , , and , , are magnitude parameters of different types of upwards and downward random jumps, respectively. The average jump size is given aswhere is the expectation operator under probability measure . Let the following be the generator matrix of the Markov chain process:with constants , for all , and

According to [9], the value of stock loans satisfies the following coupled FPIDEs:where , . And,which is left-sided Riemann–Liouville fractional derivatives.

So, the boundary conditions of the mathematical model can be obtained as follows ([2, 20, 21]):where denotes the optimal redeem price of stock loan contract. For every , the model is a free-boundary problem. Associated with each regime is an optimal exercise and a FPIDE, and it is impossible to obtain the analytic solution; therefore, the numerical method should be considered. Moreover, according to the conditions of stock loan contract, the investor should continue to hold the contract if the redemption value is less than the holding value. Therefore, the function must satisfy the following inequality:for all and .

3. Fully Implicit Scheme Based on Penalty Method

3.1. Model Transformation

To avoid the effect of time-factor in our model on the numerical results, we first introduce a new variable system as follows:

By variable transformation (the details are displayed in Appendix A), the pricing systems (6) and (8) can be written aswhere and . And, inequality (9) should be rewritten asfor all and .

It is straightforward to obtain that the function in system (11)–(15) can be viewed as an American call option with strike price and free-boundary . The analytic solution of is laborious and even impossible to achieve; therefore, we next employ a finite difference scheme based on the penalty function to solve it.

Now, we extend the penalty method to model (11)–(15) and develop an efficient numerical method via the implementation of a fully implicit scheme. So, we let be a small regularization parameter and is a constant. By adding the penalty terms [22],to the corresponding equations in (11), we obtain the following system of nonlinear space-fractional partial integral-differential equations on a fixed domain:where denotes the maximum stock price, , and . According to the numerical experiments conducted by Kandilarov and Valkov [23], the maximum value of stock is equal to 3 or 4 times of strike price, and we take in this paper. In addition, we will omit the subscript of function for convenience.

3.2. Implicit Scheme

We first take as spatial step, such that , where is a positive integer and then place uniform grids in the direction, namely, . That is,where ; ; and denotes the value of function at point and state.

For the integral term, the trapezoidal formula is used to discrete it at grid point (see [9]) as follows:where

In addition, the left-hand Riemann–Liouville fractional derivative can be approximated by the first-order Grünwald–Letnikov formula as follows (see [12, 24]):where are the fractional difference coefficients given as follows:

Using the discretization and applying the fully implicit method, it leads to a nonlinear algebraic equation:with the boundary and terminal conditions as follows:

For our proposed difference scheme, we have a discrete analogue form of the important property inherited by the stock loan model. Prior to presenting the proof, we need the following lemma of [24].

Lemma 1. The coefficients satisfyfor .

Lemma 2. Both the coefficients in equation (24) and in equation (25) are bounded, and

Proof. From the fact that is the density function of hyperexponential random variables , it follows thatIn addition, and or , so it is straightforward to obtain

Theorem 1. If the time step andthen the approximate stock loan values generated by the scheme (28) satisfyfor , , and .

Proof. First, we prove . Letfor all , and . It is straightforward to obtain . Hence, by substituting into (28) and after simplification, it yieldswhereFrom and , we have thatFrom the fact and according to [14], when . Therefore, we obtain thatTo sum up, combining Lemma 2, it is straightforward to obtainDefineand let be a pair of indices, such that . So, from (28), it follows thatTo reduce the form of inequality (43), we obtain thatOn the other hand,Hence, it is straightforward to obtain thatLetDefine a functionSo, if , we can obtainFromit follows thatBased on the condition , we can obtain . Naturally, . Consequently,for all , and .
Next, we prove that for all , and . Following the idea above, defineand let be a pair of indices, such that . It follows from (28) thatBy simple calculation, the following equality holds:From above, the inequalities for all are verified. So,Notice that , thenIn addition, for all . Namely, . Therefore, it yields thatfor all , and , which completes the proof.

3.3. Implementation of Numerical Method

In practice, to implement the numerical scheme (28), the semi-infinite domain must be truncated into a finite domain:

Here, we take . Then, we can check that . In such a case, we redefine the space step , and

Letthen the matrix form of numerical scheme (28) can be written as follows:where , , and

denotes the identity matrix.

, , and

In fact, system (62) is a coupling nonlinear equation so that it is very hard to obtain the numerical solution. Hence, we letwhere

Then, (62) can be written as the following nonlinear system:

Note that system (68) cannot be solved directly since the penalty function is nonlinear with respect to . In the following, Newton iteration approach is employed to solvewhere and is a damping parameter. The initial value for each time level is the given initial guess and . is a Jacobian matrix with column vector . We choose . The condition is the stopping criterion, where the positive control is a sufficiently small tolerance number. In this paper, we take , and .

4. Numerical Examples and Discussion

4.1. Discussion on the Numerical Method

In this part, we provide some numerical examples to support the theoretical results that we have obtained in this paper.

The model parameters are set to be , , , , , , , , and . All contracts of stock loan have maturity years and principal . First of all, to verify the conclusion that inequality holds in Theorem 1, we show the mesh surface of in Figure 1 for different time and stock price in two states. As shown in Figure 1, it is straightforward to conclude that the present difference scheme conserves the inequality for all , and . In Figure 2, the values of as a function of both stock price and time over the rectangular domain under difference economic state are displayed. It can be observed from Figure 2 that our numerical method produces smooth and stable approximation solutions.

(a)

(b)

(a)

(b)

In fact, the function in models (18)–(21) can be recognized as an American call option with strike and optimal exercise boundary, . Then, the values of are not less than the payoff function , and the function is increasing with respect to the duration in every state. Figures 3 and 4 show these two facts, respectively.

In Ref. [9], for solving the coupled equations resulted by regime switch, the authors proposed the iteration algorithm. Hence, we use this idea to solve system (62), and the consumed CPU time to run our method and the iteration algorithm for different and fixed is displayed in Table 1. All numerical computations were carried out by using MATLAB on the Lenovo S5 laptop with configuration as follows: Intel(R) Core(TM) i5-6300HQ, 2.30 GHz. We solve the linear systems by Matlab operator . In Table 1, and denote consumed CPU time of the iteration algorithm and our method, respectively. According to the results in Table 1, we can obtain that our method is more effective, and the advantages of our method become more obvious with the increase in .

4.2. Impact of Parameters

In this subsection, we discuss some interesting implications by considering the two-state regime-switching model in pricing stock loans. Moreover, we compare the result with the optimal redemption prices obtained by using the framework without regime switching.

In Figure 5, there are four different curves of optimal redemption prices under different parameters, and the horizontal axis denotes the time to expiration. The black and blue curves are obtained by the models without regime switching as . In the case of or , the optimal redemption prices calculated by our model are bigger than those obtained by the model without regime switching, which means the investors should choose a higher price to redeem stock when there is uncertainty in the financial market. However, as the maturity closes to zero, the uncertainty should be disappeared such that the holder should exercise the stock loan contract at the price level of as shown in Figure 5.

As shown in Figure 6, the introduction of regime switching produces higher optimal redemption prices than those calculated by the framework without regime switching . This phenomenon is the same as an American option, of which the optimal exercise values are monotonically increasing functions of volatility. The difference between blue line and green line (red line and black line) reflects the extra values of having a certain positive probability that the stock should stay at the economy state with higher risk. However, as the maturity of stock loan contract approaches, the expected amount of time spent on other economy state also decreases, which can result in the extra value due to decreasing regime switching. In addition, as time approaches to the expiration, the investor should redeem stock at price level of as shown in Figure 6.

5. Conclusions

In this paper, we investigate the pricing model of stock loan under the FMLS model with jump diffusion process in a -state regime-switching economy, which was formulated as a -dimension coupling fractional partial integro-differential equations. We first introduce a penalty term to transform the original model into one with fixed domain in every state. Secondly, a fully nonlinear implicit difference scheme is proposed, and we prove the numerical solution following the constraint inherited by the stock loan. Moreover, our numerical examples show an important phenomenon that the uncertainty coming from the regime switching should raise the optimal redemption prices of stock loans. And the results in this paper can be used to do research studies on other types of American Contingent Claims.

In fact, during the process of solving the coupling linear equation, the values of objective function in different economical state are calculated simultaneously rather than obtained through iteration, respectively; therefore, our method is more efficient. Under the same accuracy, our method is more effective than the iteration algorithm as shown the results of simulation in Table 1.

In future research, first of all, one could expand the number of underlying assets, which is also meaningful. Secondly, the investor can terminate the contract by paying liquidated damages. Finally, the fractional derivatives and the coupled term usually result in a full or dense coefficient matrix, which has significant computational and storage requirements; therefore, one could improve the computing speed of our method.

Appendix

Ifthen

Let , then

Now, substituting equations (A.1)–(A.5) into equation (6) and boundary conditions (8), we can obtain models (11)–(15).

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work on this paper was partially supported by the Basic Research Projects under the Department of Science and Technology in Guizhou with grantno. 1175 [2019], Natural Science Research Project of Guizhou Education Department with grant no. KY[2021]276, and Humanities and Social Sciences Projects of Guizhou University of Commerce with grant no. 2020YJZD005.

References

J. Xia and X. Y. Zhou, “Stock loans,” Mathematical Finance, vol. 17, no. 2, pp. 307–317, 2007.
View at: Publisher Site | Google Scholar
W. Chen, L. Xu, and S.-P. Zhu, “Stock loan valuation under a stochastic interest rate model,” Computers & Mathematics with Applications, vol. 70, no. 8, pp. 1757–1771, 2015.
View at: Publisher Site | Google Scholar
N. Cai and L. Sun, “Valuation of stock loans with jump risk,” Journal of Economic Dynamics and Control, vol. 40, pp. 213–241, 2014.
View at: Publisher Site | Google Scholar
J. R. Liang, J. Wang, W. J. Zhang, W. Y. Qiu, and F. Y. Ren, “Option pricing of a bi-fractional Black-Merton-Scholes model with the Hurst exponent H in [12, 1],” Applied Mathematics Letters, vol. 23, pp. 859–863, 2010.
View at: Publisher Site | Google Scholar
D. Prager and Q. Zhang, “Stock loan valuation under a regime-switching model with mean-reverting and finite maturity,” Journal of Systems Science and Complexity, vol. 23, no. 3, pp. 572–583, 2010.
View at: Publisher Site | Google Scholar
P. Carr and L. Wu, “The finite moment log stable process and option pricing,” The Journal of Finance, vol. 58, no. 2, pp. 753–777, 2003.
View at: Publisher Site | Google Scholar
E. F. Fama, “The behavior of stock-market prices,” The Journal of Business, vol. 38, no. 1, pp. 34–105, 1965.
View at: Publisher Site | Google Scholar
B. B. Mandelbrot, Fractals and Scaling in Finance: Discontinuity, Concentration, Risk, Springer-Verlag, Berlin, Germany, 1997.
Z. Zhou, J. T. Ma, and X. M. Guo, “Convergence of iterative laplace transform methods for a system of fractional PDEs and PIDEs arising in option pricing,” East Asian Journal on Applied Mathematics, vol. 8, no. 4, pp. 782–808, 2018.
View at: Publisher Site | Google Scholar
G. Jumarie, “Derivation and solutions of some fractional Black Scholes equations in coarse-grained space and time. Application to Mertons optimal portfolio,” Computers & Mathematics with Applications, vol. 3, no. 59, pp. 1142–1164, 2010.
View at: Publisher Site | Google Scholar
S. Kumar, A. Yildirim, Y. Khan, H. Jafari, K. Sayevand, and L. Wei, “Analytical solution of fractional Black-Scholes European option pricing equation by using laplace transform,” Journal of Fractional Calculus and Applications, vol. 2, no. 8, pp. 1–9, 2012.
View at: Google Scholar
S. Chen, F. Liu, X. Jiang, I. Turner, and V. Anh, “A fast semi-implicit difference method for a nonlinear two-sided space-fractional diffusion equation with variable diffusivity coefficients,” Applied Mathematics and Computation, vol. 257, pp. 591–601, 2015.
View at: Publisher Site | Google Scholar
A. Cartea and D. del Castillo-Negrete, “Fractional diffusion models of option prices in markets with jumps,” Physica A: Statistical Mechanics and Its Applications, vol. 2, no. 374, pp. 749–763, 2007.
View at: Publisher Site | Google Scholar
O. Marom and E. Momoniat, “A comparison of numerical solutions of fractional diffusion models in finance,” Nonlinear Analysis: Real World Applications, vol. 10, no. 6, pp. 3435–3442, 2009.
View at: Publisher Site | Google Scholar
W. Li, The Numerical Solution of Fractional Order Equation in Financial Models and its Application, Hangzhou University of Electronic Science and Technology, Hangzhou, China, 2009, Masters thesis.
G. Hariharan, “An efficient wavelet based approximation method to time fractional black-scholes european option pricing problem arising in financial market,” Applied Mathematical Sciences, vol. 69, no. 7, pp. 3445–3456, 2013.
View at: Publisher Site | Google Scholar
W. Chen, Numerical Methods for Fractional Black-Scholes Equations and Variational Inequalities Governing Option Pricing, The University of Western Australia, Crawley, Australia, 2013, Ph.D. thesis.
H. Zhang, F. Liu, I. Turner, and S. Chen, “The numerical simulation of the tempered fractional Black-Scholes equation for European double barrier option,” Applied Mathematical Modelling, vol. 40, no. 11-12, pp. 5819–5834, 2016.
View at: Publisher Site | Google Scholar
C. Y. Fan, K. L. Xiang, and S. Z. Chen, “Stock loan valuation based on the finite moment log stable process,” Computers & Mathematics with Applications, vol. 75, pp. 374–387, 2018.
View at: Publisher Site | Google Scholar
W. Wang and P. Chen, “Valuation of stock loan under uncertain stock model with floating interest rate,” Soft Computing, vol. 24, no. 3, pp. 1803–1814, 2020.
View at: Publisher Site | Google Scholar
Z. Q. Xu and F. H. Yi, “Optimal redeeming strategy of stock loans under drift uncertainty,” Mathematics of Operations Research, vol. 45, no. 1, 2020.
View at: Publisher Site | Google Scholar
A. Q. M. Khaliq and R. H. Liu, “New numerical scheme for pricing American option with regime-switching,” International Journal of Theoretical and Applied Finance, vol. 12, no. 03, pp. 319–340, 2009.
View at: Publisher Site | Google Scholar
J. D. Kandilarov and R. L. Valkov, “A numerical approach for the American call option pricing model,” Numerical Methods and Applications, vol. 6046, pp. 453–460, 2011.
View at: Publisher Site | Google Scholar
H. Hejazi, T. Moroney, and F. Liu, “Stability and convergence of a finite volume method for the space fractional advection-dispersion equation,” Journal of Computational and Applied Mathematics, vol. 255, pp. 684–697, 2014.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Congyin Fan and Peimin Chen. 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.

PDF Download Citation

Download other formats

Order printed copies