Source author record

Igor Mezic

Igor Mezic appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2026arXiv

Analytic Extended Dynamic Mode Decomposition

We develop a novel EDMD-type algorithm that captures the spectrum of the Koopman operator defined on a reproducing kernel Hilbert space of analytic functions. This method, which we call analytic EDMD, relies on an orthogonal projection on polynomial subspaces, which is equivalent to a data-driven Taylor approximation. In the case of dynamics with a hyperbolic equilibrium, analytic EDMD demonstrates excellent performance to capture the lattice-structured Koopman spectrum based on the eigenvalues of the linearized system at the equilibrium. Moreover, it yields the Taylor approximation of associated principal eigenfunctions. Since the method preserves the triangular structure of the operator, it does not suffer from spectral pollution and, moreover, arbitrary accuracy on the spectrum can be reached with a fixed finite dimension of the approximation and with a (possibly non-uniform) sampling over an arbitrary set of nonzero measure. The performance of analytic EDMD is illustrated with numerical examples and is assessed through a comparative study with related methods. Finally, the method is complemented with theoretical results, proving strong convergence of the eigenfunctions and providing error bounds on the spectrum estimation.

preprint2022arXiv

An Operator Theoretic View on Pruning Deep Neural Networks

The discovery of sparse subnetworks that are able to perform as well as full models has found broad applied and theoretical interest. While many pruning methods have been developed to this end, the naïve approach of removing parameters based on their magnitude has been found to be as robust as more complex, state-of-the-art algorithms. The lack of theory behind magnitude pruning's success, especially pre-convergence, and its relation to other pruning methods, such as gradient based pruning, are outstanding open questions in the field that are in need of being addressed. We make use of recent advances in dynamical systems theory, namely Koopman operator theory, to define a new class of theoretically motivated pruning algorithms. We show that these algorithms can be equivalent to magnitude and gradient based pruning, unifying these seemingly disparate methods, and find that they can be used to shed light on magnitude pruning's performance during the early part of training.

preprint2020arXiv

Convex computation of extremal invariant measures of nonlinear dynamical systems and Markov processes

We propose a convex-optimization-based framework for computation of invariant measures of polynomial dynamical systems and Markov processes, in discrete and continuous time. The set of all invariant measures is characterized as the feasible set of an infinite-dimensional linear program (LP). The objective functional of this LP is then used to single-out a specific measure (or a class of measures) extremal with respect to the selected functional such as physical measures, ergodic measures, atomic measures (corresponding to, e.g., periodic orbits) or measures absolutely continuous w.r.t. to a given measure. The infinite-dimensional LP is then approximated using a standard hierarchy of finite-dimensional semidefinite programming problems (SDPs), the solutions of which are truncated moment sequences, which are then used to reconstruct the measure. In particular, we show how to approximate the support of the measure as well as how to construct a sequence of weakly converging absolutely continuous approximations. The presented framework, where a convex functional is minimized or maximized among all invariant measures, can be seen as a generalization of and a computational method to carry out the so called ergodic optimization, where linear functionals are optimized over the set of invariant measures. Finally, we also describe how the presented framework can be adapted to compute eigenmeasures of the Perron-Frobenius operator.

preprint2020arXiv

Extended Dynamic Mode Decomposition with Learned Koopman Eigenfunctions for Prediction and Control

This paper presents a novel learning framework to construct Koopman eigenfunctions for unknown, nonlinear dynamics using data gathered from experiments. The learning framework can extract spectral information from the full nonlinear dynamics by learning the eigenvalues and eigenfunctions of the associated Koopman operator. We then exploit the learned Koopman eigenfunctions to learn a lifted linear state-space model. To the best of our knowledge, our method is the first to utilize Koopman eigenfunctions as lifting functions for EDMD-based methods. We demonstrate the performance of the framework in state prediction and closed loop trajectory tracking of a simulated cart pole system. Our method is able to significantly improve the controller performance while relying on linear control methods to do nonlinear control.

preprint2020arXiv

On Numerical Approximations of the Koopman Operator

We study numerical approaches to computation of spectral properties of composition operators. We provide a characterization of Koopman Modes in Banach spaces using Generalized Laplace Analysis. We cast the Dynamic Mode-Decomposition type methods in the context of Finite Section theory of infinite dimensional operators, and provide an example of a mixing map for which the finite section method fails. Under assumptions on the underlying dynamics, we provide the first result on the convergence rate under sample size increase in the finite-section approximation. We study the error in the Krylov subspace version of the finite section method and prove convergence in pseudospectral sense for operators with pure point spectrum. This result indicates that Krylov sequence-based approximations can have low error without an exponential-in-dimension increase in the number of functions needed for approximation.

preprint2020arXiv

Spectral analysis of mixing in 2D high-Reynolds flows

We use spectral analysis of Eulerian and Lagrangian dynamics to study the advective mixing in an incompressible 2D bounded cavity flow. A significant property of such a rotational flow at high Reynolds numbers is that mixing in its core is slower than wall-adjacent areas and corner eddies. We explain this property by considering the resonance between frequencies of unsteady motion -- captured by the Koopman spectral analysis of the velocity field -- and the circulation frequency of Lagrangian tracers in the mean flow. In high-Reynolds rotational 2D flows, the vorticity in the rotational core is uniformly distributed, which leads to uniform distribution of circulation periods in the mean flow, i.e., the kinematics in the core of mean flow is like rigid-body rotation. When this ``rigid" core is exposed to velocity fluctuations arising from bifurcations at high Reynolds, it shows more resilience toward resonance in Lagrangian motion and hence mixes more slowly compared to other areas. We also show how our qualitative resonance argument extends to chaotic flows where the classical tools of dynamical systems are not applicable.

preprint2016arXiv

Finite-time Partitions for Lagrangian Structure Identification in Gulf Stream Eddy Transport

We develop a methodology to identify finite-time Lagrangian structures from data and models using an extension of the Koopman operator-theoretic methods developed for velocity fields with simple (periodic, quasi-periodic) time-dependence. To achieve this, the notion of the Finite Time Ergodic (FiTER) partition is developed and rigorously justified. In combination with a clustering-based approach, the methodology enables identification of the temporal evolution of Lagrangian structures in a classic, benchmark, oceanographic transport problem, namely the cross-stream flux induced by the interaction of a meso- scale Gulf Stream Ring eddy with the main jet. We focus on a single mixing event driven by the interaction between an energetic cold core ring (a cyclone), the strong jet, and a number of smaller scale cyclones and anticyclones. The new methodology enab les reconstruction of Lagrangian structures in three dimensions and analysis of their time-evolution.

preprint2015arXiv

Extreme phase sensitivity in systems with fractal isochrons

Sensitivity to initial conditions is usually associated with chaotic dynamics and strange attractors. However, even systems with (quasi)periodic dynamics can exhibit it. In this context we report on the fractal properties of the isochrons of some continuous-time asymptotically periodic systems. We define a global measure of phase sensitivity that we call the phase sensitivity coefficient and show that it is an invariant of the system related to the capacity dimension of the isochrons. Similar results are also obtained with discrete-time systems. As an illustration of the framework, we compute the phase sensitivity coefficient for popular models of bursting neurons, suggesting that some elliptic bursting neurons are characterized by isochrons of high fractal dimensions and exhibit a very sensitive (unreliable) phase response.

preprint2015arXiv

Global stability analysis using the eigenfunctions of the Koopman operator

We propose a novel operator-theoretic framework to study global stability of nonlinear systems. Based on the spectral properties of the so-called Koopman operator, our approach can be regarded as a natural extension of classic linear stability analysis to nonlinear systems. The main results establish the (necessary and sufficient) relationship between the existence of specific eigenfunctions of the Koopman operator and the global stability property of hyperbolic fixed points and limit cycles. These results are complemented with numerical methods which are used to estimate the region of attraction of the fixed point or to prove in a systematic way global stability of the attractor within a given region of the state space.

preprint2014arXiv

Existence of invariant tori in three dimensional maps with degeneracy

We prove a KAM-type result for the persistence of two-dimensional invariant tori in perturbations of integrable action-angle-angle maps with degeneracy, satisfying the intersection property. Such degenerate action-angle-angle maps arise upon generic perturbation of three-dimensional volume-preserving vector fields, which are invariant under volume-preserving action of $S^1$ when there is no motion in the group action direction for the unperturbed map. This situation is analogous to degeneracy in Hamiltonian systems. The degenerate nature of the map and the unequal number of action and angle variables make the persistence proof non-standard. The persistence of the invariant tori as predicted by our result has implications for the existence of barriers to transport in three-dimensional incompressible fluid flows. Simulation results indicating existence of two-dimensional tori in a perturbation of swirling Hill's spherical vortex flow are presented.

preprint2013arXiv

Isostables, isochrons, and Koopman spectrum for the action-angle representation of stable fixed point dynamics

For asymptotically periodic systems, a powerful (phase) reduction of the dynamics is obtained by computing the so-called isochrons, i.e. the sets of points that converge toward the same trajectory on the limit cycle. Motivated by the analysis of excitable systems, a similar reduction has been attempted for non-periodic systems admitting a stable fixed point. In this case, the isochrons can still be defined but they do not capture the asymptotic behavior of the trajectories. Instead, the sets of interest-that we call isostables-are defined in literature as the sets of points that converge toward the same trajectory on the stable slow manifold of the fixed point. However, it turns out that this definition of the isostables holds only for systems with slow-fast dynamics. Also, efficient methods for computing the isostables are missing. The present paper provides a general framework for the definition and the computation of the isostables of stable fixed points, which is based on the spectral properties of the so-called Koopman operator. More precisely, the isostables are defined as the level sets of a particular eigenfunction of the Koopman operator. Through this approach, the isostables are unique and well-defined objects related to the asymptotic properties of the system. Also, the framework reveals that the isostables and the isochrons are two different but complementary notions which define a set of action-angle coordinates for the dynamics. In addition, an efficient algorithm for computing the isostables is obtained, which relies on the evaluation of Laplace averages along the trajectories. The method is illustrated with the excitable FitzHugh-Nagumo model and with the Lorenz model. Finally, we discuss how these methods based on the Koopman operator framework relate to the global linearization of the system and to the derivation of special Lyapunov functions.

preprint2011arXiv

Scalable Approach to Uncertainty Quantification and Robust Design of Interconnected Dynamical Systems

Development of robust dynamical systems and networks such as autonomous aircraft systems capable of accomplishing complex missions faces challenges due to the dynamically evolving uncertainties coming from model uncertainties, necessity to operate in a hostile cluttered urban environment, and the distributed and dynamic nature of the communication and computation resources. Model-based robust design is difficult because of the complexity of the hybrid dynamic models including continuous vehicle dynamics, the discrete models of computations and communications, and the size of the problem. We will overview recent advances in methodology and tools to model, analyze, and design robust autonomous aerospace systems operating in uncertain environment, with stress on efficient uncertainty quantification and robust design using the case studies of the mission including model-based target tracking and search, and trajectory planning in uncertain urban environment. To show that the methodology is generally applicable to uncertain dynamical systems, we will also show examples of application of the new methods to efficient uncertainty quantification of energy usage in buildings, and stability assessment of interconnected power networks.