Solving Fractional Stochastic Differential Equations via a Bilinear Time-Series Framework

This paper introduces a novel numerical approach for solving fractional stochastic differential equations (FSDEs) using bilinear time-series models, driven by the Caputo–Katugampola (C-K) fractional derivative. The C-K operator generalizes classical fractional derivatives by incorporating an additional parameter, enabling the enhanced modeling of memory effects and hereditary properties in stochastic systems. The primary contribution of this work is the development of an efficient numerical framework that combines bilinear time-series discretization with the C-K derivative to approximate solutions for FSDEs, which are otherwise analytically intractable due to their nonlinear and memory-dependent nature. We rigorously analyze the impact of fractional-order dynamics on system behavior. The bilinear time-series framework provides a computationally efficient alternative to traditional methods, leveraging multiplicative interactions between past observations and stochastic innovations to model complex dependencies. A key advantage of our approach is its flexibility in handling both stochasticity and fractional-order effects, making it suitable for applications in a famous nuclear physics model. To validate the method, we conduct a comparative analysis between exact solutions and numerical approximations, evaluating convergence properties under varying fractional orders and discretization steps. Our results demonstrate robust convergence, with simulations highlighting the superior accuracy of the C-K operator over classical fractional derivatives in preserving system dynamics. Additionally, we provide theoretical insights into the stability and error bounds of the discretization scheme. Using the changes in the number of simulations and the operator parameters of Caputo–Katugampola, we can extract some properties of the stochastic fractional differential model, and also note the influence of Brownian motion and its formulation on the model, the main idea posed in our contribution based on constructing the fractional solution of a proposed fractional model using known bilinear time series illustrated by application in nuclear physics models.