Journal
MULTISCALE MODELING & SIMULATION
Volume 9, Issue 4, Pages 1588-1623Publisher
SIAM PUBLICATIONS
DOI: 10.1137/100795917
Keywords
parameter estimation; diffusion process; stochastic differential equation; generator; discrete sampling; multiscale analysis; homogenization; subsampling
Funding
- Netherlands Organization for Scientific Research through the research cluster Nonlinear Dynamics of Natural Systems
- NSF [DMS02-39625, DMS07-08140]
- ONR [N00014-04-1-0565]
Ask authors/readers for more resources
In this paper we present a new procedure for the estimation of diffusion processes from discretely sampled data. It is based on the close relation between eigenpairs of the diffusion operator L and those of the conditional expectation operator P-t, a relation stemming from the semigroup structure P-t = exp(tL) for t >= 0. It allows for estimation without making time discretization errors, an aspect that is particularly advantageous in the case of data with low sampling frequency. After estimating eigenpairs of L via eigenpairs of P-t, we infer the drift and diffusion functions that determine L by fitting L to the estimated eigenpairs using a convex optimization procedure. We present numerical examples in which we apply the procedure to one- and two-dimensional diffusions, reversible as well as nonreversible. In the second part of the paper, we consider estimation of coarse-grained (homogenized) diffusion processes from multiscale data. We show that eigenpairs of the homogenized diffusion operator are asymptotically close to eigenpairs of the underlying multiscale diffusion operator. This implies that we can infer the correct homogenized process from data of the multiscale process, using the estimation procedure discussed in the first part of the paper. This is illustrated with numerical examples.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available