Saved in:
Bibliographic Details
Main Authors: Stevens-Haas, Jacob, Bhangale, Yash, Aravkin, Aleksandr, Kutz, Nathan
Format: Preprint
Published: 2024
Subjects:
Online Access:https://arxiv.org/abs/2405.03154
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866909190861094912
author Stevens-Haas, Jacob
Bhangale, Yash
Aravkin, Aleksandr
Kutz, Nathan
author_facet Stevens-Haas, Jacob
Bhangale, Yash
Aravkin, Aleksandr
Kutz, Nathan
contents 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 involves a derivative estimation step (and optionally, smoothing) and a sparse regression step on a library of candidate ODE terms. Kalman smoothing is a classical framework for assimilating the measurement data with known noise statistics. Previously, derivatives in SINDy and its python package, pysindy, had been estimated by finite difference, L1 total variation minimization, or local filters like Savitzky-Golay. In contrast, Kalman allows discovering ODEs that best recreate the essential dynamics in simulation, even in cases when it does not perform as well at recovering coefficients, as measured by their F1 score and mean absolute error. We have incorporated Kalman smoothing, along with hyperparameter optimization, into the existing pysindy architecture, allowing for rapid adoption of the method. Numerical experiments on a number of dynamical systems show Kalman smoothing to be the most amenable to parameter selection and best at preserving problem structure in the presence of noise.
format Preprint
id arxiv_https___arxiv_org_abs_2405_03154
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Learning Nonlinear Dynamics Using Kalman Smoothing
Stevens-Haas, Jacob
Bhangale, Yash
Aravkin, Aleksandr
Kutz, Nathan
Dynamical Systems
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 involves a derivative estimation step (and optionally, smoothing) and a sparse regression step on a library of candidate ODE terms. Kalman smoothing is a classical framework for assimilating the measurement data with known noise statistics. Previously, derivatives in SINDy and its python package, pysindy, had been estimated by finite difference, L1 total variation minimization, or local filters like Savitzky-Golay. In contrast, Kalman allows discovering ODEs that best recreate the essential dynamics in simulation, even in cases when it does not perform as well at recovering coefficients, as measured by their F1 score and mean absolute error. We have incorporated Kalman smoothing, along with hyperparameter optimization, into the existing pysindy architecture, allowing for rapid adoption of the method. Numerical experiments on a number of dynamical systems show Kalman smoothing to be the most amenable to parameter selection and best at preserving problem structure in the presence of noise.
title Learning Nonlinear Dynamics Using Kalman Smoothing
topic Dynamical Systems
url https://arxiv.org/abs/2405.03154