/Informed-HySINDy

Informed-HySINDy ( Informed Hybrid Sparse Identification of Nonlinear Dynamics )

Primary LanguageMATLAB

Informed-HySINDy

Informed-HySINDy ( Informed Hybrid Sparse Identification of Nonlinear Dynamics )

Related article: Boosting the Model Discovery of Hybrid Dynamical Systems in An Informed Sparse Regression Approach

Abstract

We present an efficient data-driven sparse identification of dynamical systems. The work aims at reconstructing the different sets of governing equations and identifying discontinuity surfaces in hybrid systems when the number of discontinuities is known a priori. In a two-stages approach, we first locate the switches between separate vector fields. Then, the dynamics among the manifolds are regressed, in this case by making use of the existing algorithm developed by Brunton et al. (SINDy algorithm). The reconstruction of the discontinuity surfaces comes as the outcome of a statistical analysis implemented via symbolic regression with small clusters (micro-clusters) and a rigid library of models. These allow to classify all the feasible discontinuities that are clustered and to reduce them into the actual discontinuity surfaces. The performances of the sparse regression hybrid model discovery are tested on two numerical examples, namely, a canonical spring–mass hopper and a free/impact electromagnetic energy harvester, engineering archetypes characterized by the presence of a single and double discontinuity, respectively. Results show that a supervised approach, i.e. where the number of discontinuities is preassigned, is computationally efficient and it determines accurately both discontinuities and set of governing equations. A large improvement in the time of computation is found with the maximum achievable reliability. Informed regression-based identification offers the prospect to outperform existing data-driven identification approaches for hybrid systems at the expense of instructing the algorithm for expected discontinuities.

Code

The repository is composed of two folders containing the implementation of the algorithm for the two numerical examples presented in the abstract ("Hopper" and "FIEH"). Each folder contain its specific readme file.