Reconstruction of diffusions using spectral data from timeseries

Abstract

A numerical technique for the reconstruction of diffusion processes (diffusions, in short) from data is presented. The drift and diffusion coeffcients of the generator of the diffusion are found by minimizing an object function which measures the difference between the eigenspectrum of the operator and a reference eigenspectrum. The reference spectrum can be obtained, in discretized form, from time-series through the construction of a discrete-time Markov chain. Discretization of the Fokker-Planck operator turns minimization of the object function into a quadratic programming problem on a convex domain, for which well-established solution methods exist. The technique is a generalization of a reconstruction procedure for continuous-time Markov chain generators, recently developed by the authors. The technique also allows us to derive the coeffcients in the homogenized diffusion for the slow variables in systems with multiple timescales.