This work generalizes the subdiffusive Black-Scholes model by introducing the variable exponent in order to provide adequate descriptions for the option pricing, where the variable exponent may account for the variation of the memory property. In addition to standard nonlinear-to-linear transformation, we apply a further spatial-temporal transformation to convert the model to a more tractable form in order to circumvent the difficulties caused by the "non-positive, non-monotonic" variable-exponent memory kernel. An interesting phenomenon is that the spatial transformation not only eliminates the advection term but naturally turns the original noncoercive spatial operator into a coercive one due to the specific structure of the Black-Scholes model, which thus avoids imposing constraints on coefficients. Then we perform numerical analysis for both the semi-discrete and fully discrete schemes to support numerical simulation. Numerical experiments are carried out to substantiate the theoretical results.
Meihui ZHANG
,
Yaxue LIU
,
Mengmeng LIU
,
Wenlin QIU
,
Xiangcheng ZHENG
. GENERALIZING SUBDIFFUSIVE BLACK-SCHOLES MODEL BY VARIABLE EXPONENT: MODEL TRANSFORMATION AND NUMERICAL APPROXIMATION[J]. Acta mathematica scientia, Series B, 2026
, 46(2)
: 993
-1010
.
DOI: 10.1007/s10473-026-0224-0
[1] Adams R, Fournier J. Sobolev Spaces. San Diego: Elsevier, 2003
[2] Alvarez-Ramirez J, Alvarez J, Rodriguez E, Fernandez-Anaya G. Time-varying Hurst exponent for US stock markets. Physica A, 2008, 387: 6159-6169
[3] Araneda A. A multifractional option pricing formula. Fluct Noise Lett, 2024, 23: Art 2450060
[4] Askey R, Roy R.Gamma function, NIST Handbook of Mathematical Functions. Cambridge: Cambridge University Press, 2010
[5] Black F, Scholes M. The pricing of options and corporate liabilities. J Political Economy, 1973, 81: 637-654
[6] Cao G, Han Y, Cui W. Time-varying long memories of the Chinese currency and stock markets based on the Hurst exponent. Fluct Noise Lett, 2014, 13: Art 1450007
[7] Capiński M Kopp E. The Black-Scholes model. Cambridge: Cambridge University Press, UK, 2012
[8] Carr P, Wu L. The finite moment log stable process and option pricing. J Finance, 2003, 2: 597-626
[9] Chen W, Xu X, Zhu S. A predictor-corrector approach for pricing American options under the finite moment log-stable model. Appl Numer Math, 2015, 97: 15-29
[10] De Staelen R, Hendy A. Numerically pricing double barrier options in a time-fractional Black-Scholes model. Comput Math Appl, 2017, 74: 1166-1175
[11] Duan J, Lu L, Chen L, An Y. Fractional model and solution for the Black-Scholes equation. Math Meth Appl Sci, 2018, 41: 697-704
[12] Farrell P, Hegarty A, Miller J, O'Riordan E, Shishkin G. Robust Computational Techniques for Boundary Layers. Boca Raton, FL: Chapman & Hall/CRC, 2000
[13] Golbabai A, Nikan O, Nikazad T. Numerical analysis of time fractional Black-Scholes European option pricing model arising in financial market. Comp Appl Math, 2019, 38: Article 173
[14] Guo C, Fang S, He Y. Derivation and application of some fractional Black-Scholes equations driven by fractional G-Brownian motion. Comput Econ, 2023, 61: 1681-1705
[15] Jin B.Fractional Differential Equations: An Approach via Fractional Derivatives. Cham: Springer, 2021
[16] Jumarie G. Stock exchange fractional dynamics defined as fractional exponential growth driven by (usual) Gaussian white noise. Application to fractional Black-Scholes equations. Insur Math Econ, 2008, 42, 271-287
[17] Jumarie G. Derivation and solutions of some fractional Black-Scholes equations in coarse-grained space and time. Application to merton's optimal portfolio. Comput Math Appl, 2010, 59: 1142-1164
[18] Korbel J, Luchko Y. Modeling of financial processes with a space-time fractional diffusion equation of varying order. Fract Calc Appl Anal, 2016, 19: 1414-1433
[19] Krzyżanowski G, Magdziarz M. A tempered subdiffusive Black-Scholes model. Fract Calc Appl Anal,2024, 27: 1800-1834
[20] Lorenzo C, Hartley T. Variable order and distributed order fractional operators. Nonlinear Dynamics, 2002, 29: 57-98
[21] McLean W, Thomée V. Numerical solution of an evolution equation with a positive-type memory term. J Austral Math Soc Ser B,1993, 35: 23-70
[22] Nikan O, Avazzadeh Z, Machado J. Numerical study of the nonlinear anomalous reaction-subdiffusion process arising in the electroanalytical chemistry. J Comput Sci, 2021, 53: Art 101394
[23] Qiu W, Zheng X. Numerical analysis for high-order methods for variable-exponent fractional diffusion-wave equation. J Sci Comput, 2025, 105: Article 53
[24] Shirzadi M, Dehghan M, Bastani A. On the pricing of multi-asset options under jump-diffusion processes using meshfree moving least-squares approximation. Commun Nonlinear Sci Numer Simul, 2020, 84: Art 105160
[25] kovränek T, Podlubny I, Peträš I. Modeling of the national economies in state-space: A fractional calculus approach. Econ Model,2012, 29: 1322-1327
[26] Soleymani F, Zhu S. Error and stability estimates of a time-fractional option pricing model under fully spatial-temporal graded meshes. J Comput Appl Math, 2023, 425: Art 115075
[27] Thomee V.Galerkin Finite Element Methods for Parabolic Problems. Berlin: Springer, 2007
[28] Wang Y. A compact finite difference method for solving a class of time fractional convection-subdiffusion equations. BIT, 2015, 55: 1187-1217
[29] Wyss W. The fractional Black-Scholes equation. Frac Calc Appl Anal, 2017, 3: 51-62
[30] Zeng Z, Chen Y. A local discontinuous Galerkin method for time-fractional diffusion equations. Acta Math Sci, 2023, 43: 839-854
[31] Zhang H, Liu F, Turner I, Yang Q. Numerical solution of the time fractional Black-Scholes model governing European options. Comput Math Appl, 2016, 71: 1772-1783
[32] Zhang M, Zheng X. Numerical approximation to a variable-order time-fractional Black-Scholes model with applications in option pricing. Comput Econ, 2023, 62: 1155-1175
[33] Zheng X. Two methods addressing variable-exponent fractional initial and boundary value problems and Abel integral equation. CSIAM Trans Appl Math, 2025, 6: 666-710