This paper develops a method for obtaining guaranteed outer approximations for global attractors of continuous and discrete time nonlinear dynamical systems. The method is based on a hierarchy of semidefinite programming problems of increasing size with guaranteed convergence to the global attractor. The approach taken follows an established line of reasoning, where we first characterize the global attractor via an infinite dimensional linear programming problem (LP) in the space of Borel measures. The dual to this LP is in the space of continuous functions and its feasible solutions provide guaranteed outer approximations to the global attractor. For systems with polynomial dynamics, a hierarchy of finite-dimensional sum-of-squares tightenings of the dual LP provides a sequence of outer approximations to the global attractor with guaranteed convergence in the sense of volume discrepancy tending to zero. The method is very simple to use and based purely on convex optimization. Numerical examples with the code available online demonstrate the method.

This paper develops sparse moment-sum-of-squares approximations for three problems from nonlinear dynamical systems: region of attraction, maximum positively invariant set and global attractor. We prove general results allowing for a decomposition of these sets provided that the vector field and constraint set posses certain structure. We combine these decompositions with the methods from [10], [12] and [20] based on infinite-dimensional linear programming. For polynomial dynamics, we show that these problems admit a sparse sum-of-squares (SOS) approximation with guaranteed convergence such that the number of variables in the largest SOS multiplier is given by the dimension of the largest subsystem appearing in the decomposition. The dimension of such subsystems depends on the sparse structure of the vector field and the constraint set and can allow for a significant reduction of the size of the semidefinite program (SDP) relaxations, thereby allowing to address far larger problems without compromising convergence guarantees. The method is simple to use and based on convex optimization. Numerical examples demonstrate the approach.

Koopman and Perron–Frobenius operators for dynamical systems are becoming popular in a number of fields in science recently. Properties of the Koopman operator essentially depend on the choice of function spaces where it acts. Particularly, the case of reproducing kernel Hilbert spaces (RKHSs) is drawing increasing attention in data science. In this paper, we give a general framework for Koopman and Perron–Frobenius operators on reproducing kernel Banach spaces (RKBSs). More precisely, we extend basic known properties of these operators from RKHSs to RKBSs and state new results, including symmetry and sparsity concepts, on these operators on RKBS for discrete and continuous time systems.

No abstract