Numerical Analysis of a Class of Fractional Langevin Equations With the Block-by-Block Method

Abstract

The fractional Langevin equation is of great scientific significance and engineering application value. Based on the classical block-by-block method, the numerical solution of a class of fractional Langevin equations with Caputo derivatives was obtained. Through introduction of the quadratic Lagrange basis function interpolation, the block-by-block convergent nonlinear equations were constructed, and the numerical solution of the Langevin equation was obtained by coupling in each block. Under the condition of 0<α<1, the stochastic Taylor expansion was used to prove that the block-by-block method is (3+α)-order convergent. Numerical experiments show that, the block-by-block method is stable and convergent under different values of α and time step <i>h</i>,and overcomes the existing methods’ disadvantages of slow speed and poor accuracy for solving fractional Langevin equations.