2020/02/18 by Jian Ren, Ren, Jian, Jinqiao Duan +1
Computer Science · Mathematics · Physics and Astronomy · #37M10 #37N40 #49K30 #62J02 #68M07 #Computational Physics (physics.comp-ph) #FOS: Computer and information sciences #FOS: Mathematics #FOS: Physical sciences #Methodology (stat.ME) #Numerical Analysis (math.NA) #cs.NA #math.NA #msc:37M10 #msc:37N40 #msc:49K30 #msc:62J02 #msc:68M07 #physics.comp-ph #stat.ME
paper · pdf · doi:10.48550/arxiv.2002.10251
23 pages, 4 figures, 7 tables
arxiv created 2020/08/20 · arxiv updated 2020/08/21
Extracting governing stochastic differential equation models from elusive data is crucial to understand and forecast dynamics for complex systems. We devise a method to extract the drift term and estimate the diffusion coefficient of a governing stochastic dynamical system, from its time-series data of the most probable transition trajectory. By the Onsager-Machlup theory, the most probable transition trajectory satisfies the corresponding Euler-Lagrange equation, which is a second order deterministic ordinary differential equation involving the drift term and diffusion coefficient. We first estimate the coefficients of the Euler-Lagrange equation based on the data of the most probable trajectory, and then we calculate the drift and diffusion coefficients of the governing stochastic dynamical system. These two steps involve sparse regression and optimization. Finally, we illustrate our method with an example and some discussions.