LAUSR

In this paper, we study the numerical schemes for the two-dimensional Fokker-Planck equation governing the probability density function of the tempered fractional Brownian motion. The main challenges of the numerical schemes come from the singularity in the time direction. When 0 < H < 0.5, a change of variables ∂t2H=2Ht2H−1∂t\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\partial \left (t^{2H}\right )=2Ht^{2H-1}\partial t$\end{document} avoids the singularity of numerical computation at t = 0, which naturally results in nonuniform time discretization and greatly improves the computational efficiency. For H > 0.5, the time span dependent numerical scheme and nonuniform time discretization are introduced to ensure the effectiveness of the calculation and the computational efficiency. The stability and convergence of the numerical schemes are demonstrated by using Fourier method. By numerically solving the corresponding Fokker-Planck equation, we obtain the mean squared displacement of stochastic processes, which conforms to the characteristics of the tempered fractional Brownian motion.