Learning Nonlinear Dynamics Using Kalman Smoothing

被引:1
|
作者
Stevens-Haas, Jacob M. [1 ]
Bhangale, Yash [2 ]
Kutz, J. Nathan [1 ,3 ]
Aravkin, Aleksandr [1 ]
机构
[1] Univ Washington, Dept Appl Math, Seattle, WA 98195 USA
[2] Univ Washington, Dept Mech Engn, Seattle, WA 98195 USA
[3] Univ Washington, Dept Elect & Comp Engn, Seattle, WA 98195 USA
来源
IEEE ACCESS | 2024年 / 12卷
基金
美国国家科学基金会;
关键词
Kalman filters; Smoothing methods; Noise measurement; Mathematical models; Filtering; Nonlinear dynamical systems; Dynamical systems; Differential equations; machine learning; sparse regression; optimization; Kalman smoothing; SINDy; differential equations; FRAMEWORK; ROBUST;
D O I
10.1109/ACCESS.2024.3465390
中图分类号
TP [自动化技术、计算机技术];
学科分类号
0812 ;
摘要
Identifying Ordinary Differential Equations (ODEs) from measurement data requires both fitting the dynamics and assimilating, either implicitly or explicitly, the measurement data. The Sparse Identification of Nonlinear Dynamics (SINDy) method does so in two steps: a derivative estimation and smoothing step and a sparse regression step on a library of ODE terms. Previously, the derivative step in SINDy and its python package, pysindy, used finite difference, L1 total variation minimization, or local filters like Savitzky-Golay. We have incorporated Kalman smoothing, along with hyperparameter optimization, into the existing pysindy architecture, allowing for rapid adoption of the method. Kalman smoothing is a classical framework for assimilating the measurement data with known noise statistics. As a first SINDy step, it denoises the data by applying a prior belief that the system is an instance of integrated Brownian motion. We conduct numerical experiments on eight dynamical systems show Kalman smoothing to be the best SINDy differentiation/smoothing option in the presence of noise on four of those systems, and tied for three of them. It has particular advantage at preserving problem structure in simulation. The addition of hyperparameter optimization further makes it the most amenable method for generic data. In doing so, it is the first SINDy method for noisy data that requires only a single hyperparameter, and it gives viable results in half of the systems we test.
引用
收藏
页码:138564 / 138574
页数:11
相关论文
共 50 条
  • [41] ON USE OF THE KALMAN FILTER FOR SPATIAL SMOOTHING
    KASHIWAGI, N
    ANNALS OF THE INSTITUTE OF STATISTICAL MATHEMATICS, 1993, 45 (01) : 21 - 34
  • [42] Ensemble-based Kalman Filters in Strongly Nonlinear Dynamics
    Joshua HACKER
    Advances in Atmospheric Sciences, 2009, 26 (03) : 373 - 380
  • [43] Ensemble-based Kalman filters in strongly nonlinear dynamics
    Zhaoxia Pu
    Joshua Hacker
    Advances in Atmospheric Sciences, 2009, 26 : 373 - 380
  • [44] Ensemble-based Kalman filters in strongly nonlinear dynamics
    Pu, Zhaoxia
    Hacker, Joshua
    ADVANCES IN ATMOSPHERIC SCIENCES, 2009, 26 (03) : 373 - 380
  • [45] Reinforcement Learning using Kalman Filters
    Takahata, Kei
    Miura, Takao
    PROCEEDINGS OF THE 2019 IEEE 18TH INTERNATIONAL CONFERENCE ON COGNITIVE INFORMATICS & COGNITIVE COMPUTING (ICCI*CC 2019), 2019, : 136 - 143
  • [46] Generative learning for nonlinear dynamics
    William Gilpin
    Nature Reviews Physics, 2024, 6 : 194 - 206
  • [47] Nonlinear Dynamics of Motor Learning
    Mayer-Kress, Gottfried
    Newell, Karl M.
    Liu, Yeou-Teh
    NONLINEAR DYNAMICS PSYCHOLOGY AND LIFE SCIENCES, 2009, 13 (01) : 3 - 26
  • [48] Generative learning for nonlinear dynamics
    Gilpin, William
    NATURE REVIEWS PHYSICS, 2024, 6 (03) : 194 - 206
  • [49] Estimating model-error covariances in nonlinear state-space models using Kalman smoothing and the expectation-maximization algorithm
    Dreano, D.
    Tandeo, P.
    Pulido, M.
    Ait-El-Fquih, B.
    Chonavel, T.
    Hoteit, I.
    QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY, 2017, 143 (705) : 1877 - 1885
  • [50] Algorithm 1051: UltimateKalman, Flexible Kalman Filtering and Smoothing Using Orthogonal Transformations
    Toledo, Sivan
    ACM Transactions on Mathematical Software, 2024, 50 (04):