TY - JOUR
T1 - Numerical strategies for filtering partially observed stiff stochastic differential equations
AU - Harlim, John
N1 - Funding Information:
This work is motivated from discussions at the SAMSI working group for the interaction between deterministic and stochastic dynamics. The author thanks Peter Kramer for sharing his knowledge on the HMM, Andrew J. Majda for his comments on the manuscript, and Emily L. Kang for her thoughts after testing the MMF algorithm on different systems. The research of the author is partially supported by the NC State University startup grant and the SAMSI teaching buyout for the Fall 2009 semester.
PY - 2011/2/1
Y1 - 2011/2/1
N2 - In this paper, we present a fast numerical strategy for filtering stochastic differential equations with multiscale features. This method is designed such that it does not violate the practical linear observability condition and, more importantly, it does not require the computationally expensive cross correlation statistics between multiscale variables that are typically needed in standard filtering approach. The proposed filtering algorithm comprises of a " macro-filter" that borrows ideas from the Heterogeneous Multiscale Methods and a " micro-filter" that reinitializes the fast microscopic variables to statistically reflect the unbiased slow macroscopic estimate obtained from the macro-filter and macroscopic observations at asynchronous times. We will show that the proposed micro-filter is equivalent to solving an inverse problem for parameterizing differential equations. Numerically, we will show that this microscopic reinitialization is an important novel feature for accurate filtered solutions, especially when the microscopic dynamics is not mixing at all.
AB - In this paper, we present a fast numerical strategy for filtering stochastic differential equations with multiscale features. This method is designed such that it does not violate the practical linear observability condition and, more importantly, it does not require the computationally expensive cross correlation statistics between multiscale variables that are typically needed in standard filtering approach. The proposed filtering algorithm comprises of a " macro-filter" that borrows ideas from the Heterogeneous Multiscale Methods and a " micro-filter" that reinitializes the fast microscopic variables to statistically reflect the unbiased slow macroscopic estimate obtained from the macro-filter and macroscopic observations at asynchronous times. We will show that the proposed micro-filter is equivalent to solving an inverse problem for parameterizing differential equations. Numerically, we will show that this microscopic reinitialization is an important novel feature for accurate filtered solutions, especially when the microscopic dynamics is not mixing at all.
UR - http://www.scopus.com/inward/record.url?scp=78649325777&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=78649325777&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2010.10.016
DO - 10.1016/j.jcp.2010.10.016
M3 - Article
AN - SCOPUS:78649325777
SN - 0021-9991
VL - 230
SP - 744
EP - 762
JO - Journal of Computational Physics
JF - Journal of Computational Physics
IS - 3
ER -