BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Rodolphe Sepulchre appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We formulate and discuss the affine-invariant matrix midrange problem on the cone of $n\times n$ positive definite Hermitian matrices $\mathbb{P}(n)$, which is based on the Thompson metric. A particular computationally efficient midpoint of this metric is investigated as a highly scalable candidate for an average of two positive definite matrices within this context, before studying the $N$-point problem in the vector and matrix settings.
The Scaled Relative Graph (SRG) is a generalization of the Nyquist diagram that may be plotted for nonlinear operators, and allows nonlinear robustness margins to be defined graphically. This abstract explores techniques for shaping the SRG of an operator in order to maximize these robustness margins.
Following the seminal work of Zames, the input-output theory of the 70s acknowledged that incremental properties (e.g. incremental gain) are the relevant quantities to study in nonlinear feedback system analysis. Yet, non-incremental analysis has dominated the use of dissipativity theory in nonlinear control from the 80s. Results connecting dissipativity theory and incremental analysis are scattered and progress has been limited. This abstract investigates whether this theoretical gap is of fundamental nature and considers new avenues to circumvent it.
Given a non-contracting trajectory of a nonlinear system, we consider the question of designing an input perturbation that makes the perturbed trajectory contracting. This paper stresses the analogy of this question with the classical question of feedback stabilization. In particular, it is shown that the existence of an output variable that ensures contraction of the inverse system facilitates the design of a contracting input perturbation. We illustrate the relevance of this question in parameter estimation.
Spikes and rhythms organize control and communication in the animal world, in contrast to the bits and clocks of digital technology. As continuous-time signals that can be counted, spikes have a mixed nature. This paper reviews ongoing efforts to develop a control theory of spiking systems. The central thesis is that the mixed nature of spiking results from a mixed feedback principle, and that a control theory of mixed feedback can be grounded in the operator theoretic concept of maximal monotonicity. As a nonlinear generalization of passivity, maximal monotonicity acknowledges at once the physics of electrical circuits, the algorithmic tractability of convex optimization, and the feedback control theory of incremental passivity. We discuss the relevance of a theory of spiking control systems in the emerging age of event-based technology.
The splitting algorithms of monotone operator theory find zeros of sums of relations. This corresponds to solving series or parallel one-port electrical circuits, or the negative feedback interconnection of two subsystems. One-port circuits with series and parallel interconnections, or block diagrams with multiple forward and return paths, give rise to current-voltage relations consisting of nested sums and inverses. In this extended abstract, we present new splitting algorithms specially suited to these structures, for interconnections of monotone and anti-monotone relations.
This paper considers balanced truncation of discrete-time Hankel $k$-positive systems, characterized by Hankel matrices whose minors up to order $k$ are nonnegative. Our main result shows that if the truncated system has order $k$ or less, then it is Hankel totally positive ($\infty$-positive), meaning that it is a sum of first order lags. This result can be understood as a bridge between two known results: the property that the first-order truncation of a positive system is positive ($k=1$), and the property that balanced truncation preserves state-space symmetry. It provides a broad class of systems where balanced truncation is guaranteed to result in a minimal internally positive system.
Neuromorphic engineering is a rapidly developing field that aims to take inspiration from the biological organization of neural systems to develop novel technology for computing, sensing, and actuating. The unique properties of such systems call for new signal processing and control paradigms. The article introduces the mixed feedback organization of excitable neuronal systems, consisting of interlocked positive and negative feedback loops acting in distinct timescales. The principles of biological neuromodulation suggest a methodology for designing and controlling mixed-feedback systems neuromorphically. The proposed design consists of a parallel interconnection of elementary circuit elements that mirrors the organization of biological neurons and utilizes the hardware components of neuromorphic electronic circuits. The interconnection structure endows the neuromorphic systems with a simple control methodology that reframes the neuronal control as an input-output shaping problem. The potential of neuronal control is illustrated on elementary network examples that suggest the scalability of the mixed-feedback principles.
This paper studies the variation diminishing property of $k$-positive linear time-invariant (LTI) systems, which map inputs with $k-1$ sign changes to outputs with at most the same variation. We characterize this property for the Toeplitz and Hankel operators of finite-dimensional systems. Our main result is that these operators have a dominant approximation in the form of series or parallel interconnections of $k$ first order positive systems. This is shown by expressing the $k$-positivity of a LTI system as the external positivity (that is, $1$-positivity) of $k$ compound LTI systems. Our characterization generalizes well known properties of externally positive systems ($k=1$) and totally positive systems ($k=\infty$; also known as relaxation systems).
Motivated by neuronal models from neuroscience, we consider the system identification of simple feedback structures whose behaviors include nonlinear phenomena such as excitability, limit-cycles and chaos. We show that output feedback is sufficient to solve the identification problem in a two-step procedure. First, the nonlinear static characteristic of the system is extracted, and second, using a feedback linearizing law, a mildly nonlinear system with an approximately-finite memory is identified. In an ideal setting, the second step boils down to the identification of a LTI system. To illustrate the method in a realistic setting, we present numerical simulations of the identification of two classical systems that fit the assumed model structure.
We define geometric matrix midranges for positive definite Hermitian matrices and study the midrange problem from a number of perspectives. Special attention is given to the midrange of two positive definite matrices before considering the extension of the problem to $N > 2$ matrices. We compare matrix midrange statistics with the scalar and vector midrange problem and note the special significance of the matrix problem from a computational standpoint. We also study various aspects of geometric matrix midrange statistics from the viewpoint of linear algebra, differential geometry and convex optimization.
The paper presents a model reduction framework geared towards the analysis and design of systems that switch and oscillate. While such phenomena are ubiquitous in nature and engineering, model reduction methods are not well developed for non-equilibrium behaviors. The proposed framework addresses this need by exploiting recent advances on dominance theory. Classical balanced truncation for linear time-invariant systems is used to develop a dominance-preserving model reduction method for Lure systems, i.e. systems that can be decomposed as the feedback interconnection of a linear system and a static nonlinearity. The method is illustrated by approximating the oscillatory behavior of a discretized heat flow control system.
Morris-Lecar model is arguably the simplest dynamical model that retains both the slow-fast geometry of excitable phase portraits and the physiological interpretation of a conductance-based model. We augment this model with one slow inward current to capture the additional property of bistability between a resting state and a spiking limit cycle for a range of input current. The resulting dynamical system is a core structure for many dynamical phenomena such as slow spiking and bursting. We show how the proposed model combines physiological interpretation and mathematical tractability and we discuss the benefits of the proposed approach with respect to alternative models in the literature.
We introduce new partial orders on the set $S^+_n$ of positive-definite matrices of dimension $n$ derived from the homogeneous geometry of $S^+_n$ induced by the natural transitive action of the general linear group $GL(n)$. The orders are induced by affine-invariant cone fields, which arise naturally from a local analysis of the orders that are compatible with the homogeneous structure of $S^+_n$. We then take a geometric approach to the study of monotone functions on $S^+_n$ and establish a number of relevant results, including an extension of the well-known Löwner-Heinz theorem derived using differential positivity with respect to affine-invariant cone fields.
We study spatio-temporal bursting in a three-scale reaction diffusion equation organized by the winged cusp singularity. For large time-scale separation the model exhibits traveling bursts, whereas for large space-scale separation the model exhibits standing bursts. Both behaviors exhibit a common singular skeleton, whose geometry is fully determined by persistent bifurcation diagrams of the winged cusp. The modulation of spatio-temporal bursting in such a model naturally translates into paths in the universal unfolding of the winged-cusp.
The notion of path-complete positivity is introduced as a way to generalize the property of positivity from one LTI system to a family of switched LTI systems whose switching rule is constrained by a finite automaton. The generalization builds upon the analogy between stability and positivity, the former referring to the contraction of a norm, the latter referring to the contraction of a cone (or, equivalently, a projective norm). We motivate and investigate the potential of path-positivity and we propose an algorithm for the automatic verification of positivity.
This paper exploits a basic connection between sequential quadratic programming and Riemannian gradient optimization to address the general question of selecting a metric in Riemannian optimization, in particular when the Riemannian structure is sought on a quotient manifold. The proposed method is shown to be particularly insightful and efficient in quadratic optimization with orthogonality and/or rank constraints, which covers most current applications of Riemannian optimization in matrix manifolds.
The paper looks at a scaled variant of the stochastic gradient descent algorithm for the matrix completion problem. Specifically, we propose a novel matrix-scaling of the partial derivatives that acts as an efficient preconditioning for the standard stochastic gradient descent algorithm. This proposed matrix-scaling provides a trade-off between local and global second order information. It also resolves the issue of scale invariance that exists in matrix factorization models. The overall computational complexity is linear with the number of known entries, thereby extending to a large-scale setup. Numerical comparisons show that the proposed algorithm competes favorably with state-of-the-art algorithms on various different benchmarks.
The paper shows that normally hyperbolic one-dimensional compact attractors of smooth dynamical systems are characterized by differential positivity, that is, the pointwise infinitesimal contraction of a smooth cone field. The result is analog to the characterization of zero-dimensional hyperbolic attractors by differential stability, which is the pointwise infinitesimal contraction of a Riemannian metric.
The paper considers an extension of factor analysis to moving average processes. The problem is formulated as a rank minimization of a suitable spectral density. It is shown that it can be adequately approximated via a trace norm convex relaxation.
This paper considers the problem of identifying multivariate autoregressive (AR) sparse plus low-rank graphical models. Based on the corresponding problem formulation recently presented, we use the alternating direction method of multipliers (ADMM) to efficiently solve it and scale it to sizes encountered in neuroimaging applications. We apply this decomposition on synthetic and real neuroimaging datasets with a specific focus on the information encoded in the low-rank structure of our model. In particular, we illustrate that this information captures the spatio-temporal structure of the original data, generalizing classical component analysis approaches.
The paper highlights the role of a positive feedback gating mechanism at the cellular level in the robust- ness and modulation properties of rhythmic activities at the circuit level. The results are presented in the context of half-center oscillators, which are simple rhythmic circuits composed of two reciprocally connected inhibitory neuronal populations. Specifically, we focus on rhythms that rely on a particu- lar excitability property, the post-inhibitory rebound, an intrinsic cellular property that elicits transient membrane depolarization when released from hyperpolarization. Two distinct ionic currents can evoke this transient depolarization: a hyperpolarization-activated cation current and a low-threshold T-type calcium current. The presence of a slow activation is specific to the T-type calcium current and provides a slow-positive feedback at the cellular level that is absent in the cation current. We show that this slow- positive feedback is necessary and sufficient to endow the network rhythm with physiological modulation and robustness properties. This study thereby identifies an essential cellular property to be retained at the network level in modeling network robustness and modulation.
The paper proposes an identification procedure for autoregressive gaussian stationary stochastic processes wherein the manifest (or observed) variables are mostly related through a limited number of latent (or hidden) variables. The method exploits the sparse plus low-rank decomposition of the inverse of the manifest spectral density and the efficient convex relaxations recently proposed for such decomposition.
This paper studies the link between resting-state functional connectivity (FC), measured by the correlations of the fMRI BOLD time courses, and structural connectivity (SC), estimated through fiber tractography. Instead of a static analysis based on the correlation between SC and the FC averaged over the entire fMRI time series, we propose a dynamic analysis, based on the time evolution of the correlation between SC and a suitably windowed FC. Assessing the statistical significance of the time series against random phase permutations, our data show a pronounced peak of significance for time window widths around 20-30 TR (40-60 sec). Using the appropriate window width, we show that FC patterns oscillate between phases of high modularity, primarily shaped by anatomy, and phases of low modularity, primarily shaped by inter-network connectivity. Building upon recent results in dynamic FC, this emphasizes the potential role of SC as a transitory architecture between different highly connected resting state FC patterns. Finally, we show that networks implied in consciousness-related processes, such as the default mode network (DMN), contribute more to these brain-level fluctuations compared to other networks, such as the motor or somatosensory networks. This suggests that the fluctuations between FC and SC are capturing mind-wandering effects.
Differential analysis aims at inferring global properties of nonlinear behaviors from the local analysis of the linearized dynamics. The paper motivates and illustrates the use of differential analysis on the nonlinear pendulum model, an archetype example of nonlinear behavior. Special emphasis is put on recent work by the authors in this area, which includes a differential Lyapunov framework for contraction analysis and the concept of differential positivity.
The paper introduces and studies differentially positive systems, that is, systems whose linearization along an arbitrary trajectory is positive. A generalization of Perron Frobenius theory is developed in this differential framework to show that the property induces a (conal) order that strongly constrains the asymptotic behavior of solutions. The results illustrate that behaviors constrained by local order properties extend beyond the well-studied class of linear positive systems and monotone systems, which both require a constant cone field and a linear state space.
Borrowing the concept of organizing center from singularity theory, the paper proposes a methodology to realize nonlinear behaviors such as switches, relaxation oscillators, or bursters from core circuits that reveal the fundamental role of monotonicity and feedback in their robustness and modulation.
We use the qualitative insight of a planar neuronal phase portrait to detect an excitability switch in arbitrary conductance-based models from a simple mathematical condition. The condition expresses a balance between ion channels that provide a negative feedback at resting potential (restorative channels) and those that provide a positive feedback at resting potential (regenerative channels). Geometrically, the condition imposes a transcritical bifurcation that rules the switch of excitability through the variation of a single physiological parameter. Our analysis of six different published conductance based models always finds the transcritical bifurcation and the associated switch in excitability, which suggests that the mathematical predictions have a physiological relevance and that a same regulatory mechanism is potentially involved in the excitability and signaling of many neurons.
Lyapunov's second theorem is an essential tool for stability analysis of differential equations. The paper provides an analog theorem for incremental stability analysis by lifting the Lyapunov function to the tangent bundle. The Lyapunov function endows the state-space with a Finsler structure. Incremental stability is inferred from infinitesimal contraction of the Finsler metrics through integration along solutions curves.
We investigate performance bounds for feedback control of distributed plants where the controller can be centralized (i.e. it has access to measurements from the whole plant), but sensors only measure differences between neighboring subsystem outputs. Such "distributed sensing" can be a technological necessity in applications where system size exceeds accuracy requirements by many orders of magnitude. We formulate how distributed sensing generally limits feedback performance robust to measurement noise and to model uncertainty, without assuming any controller restrictions (among others, no "distributed control" restriction). A major practical consequence is the necessity to cut down integral action on some modes. We particularize the results to spatially invariant systems and finally illustrate implications of our developments for stabilizing the segmented primary mirror of the European Extremely Large Telescope.
Optimization on manifolds is a rapidly developing branch of nonlinear optimization. Its focus is on problems where the smooth geometry of the search space can be leveraged to design efficient numerical algorithms. In particular, optimization on manifolds is well-suited to deal with rank and orthogonality constraints. Such structured constraints appear pervasively in machine learning applications, including low-rank matrix completion, sensor network localization, camera network registration, independent component analysis, metric learning, dimensionality reduction and so on. The Manopt toolbox, available at www.manopt.org, is a user-friendly, documented piece of software dedicated to simplify experimenting with state of the art Riemannian optimization algorithms. We aim particularly at reaching practitioners outside our field.
Exploiting the specific structure of neuron conductance-based models, the paper investigates the mathematical modeling of neuronal bursting modulation. The proposed approach combines singularity theory and geometric singular perturbations to capture the geometry of multiple time-scales attractors in the neighborhood of high-codimension singularities. We detect a three-time scale bursting attractor in the universal unfolding of the winged cusp singularity and discuss the physiological relevance of the bifurcation and unfolding parameters in determining a physiological modulation of bursting. The results suggest generality and simplicity in the organizing role of the winged cusp singularity for the global dynamics of conductance based models.
Neuronal spiking exhibits an exquisite combination of modulation and robustness properties, rarely matched in artificial systems. We exploit the particular interconnection structure of conductance based models to investigate this remarkable property. We find that much of neuronal modulation and robustness can be explained by separating the total transmembrane current into three different components corresponding to the three time scales of neuronal bursting. Each equivalent current aggregates many ionic contributions into an equivalent voltage-dependent conductance, which defines a key modulation parameter. Plugging those equivalent feedback gains in a minimal abstract model recovers many experimental modulation scenarii as modulatory paths in elementary two-parameter charts. Likewise, robustness owes to the many possible physiological realizations of a same equivalent conductance, highlighting the role of equivalent conductances as prominent targets for neuromodulation and intrinsic homeostasis.
We use a computational model to propose a physiological mechanism by which transient control of beta oscillations in the indirect pathway of the basal ganglia is orchestrated at the cellular level. Our model includes a simple and robust mechanism by which a cellular switch (from bursting to tonic) almost instantaneously translates into a functional gating switch (from blocking to conducive) in an excitatory-inhibitory network. Applied to the control of beta oscillations in the basal ganglia, the model shows the modulation of beta activity under the action of a transient depolarization, for instance a dopamine signal. The model predicts, by analogy to the thalamocortical circuit, a novel gating function by which the transfer of cortical spikes through the indirect pathway is blocked under the inhibitory drive preceding movement but briefly released at the onset of movement execution.
Differential passivity is a property that allows to check with a pointwise criterion that a system is incrementally passive, a property that is relevant to study interconnected systems in the context of regulation, synchronization, and estimation. The paper investigates how restrictive is the property, focusing on a class of open gradient systems encountered in the coenergy modeling framework of physical systems, in particular the Brayton-Moser formalism for nonlinear electrical circuits.
Dissipativity is an essential concept of systems theory. The paper provides an extension of dissipativity, named differential dissipativity, by lifting storage functions and supply rates to the tangent bundle. Differential dissipativity is connected to incremental stability in the same way as dissipativity is connected to stability. It leads to a natural formulation of differential passivity when restricting to quadratic supply rates. The paper also shows that the interconnection of differentially passive systems is differentially passive, and provides preliminary examples of differentially passive electrical systems.
Motivated by population studies of Diffusion Tensor Imaging, the paper investigates the use of mean-based and dispersion-based permutation tests to define and compute the significance of a statistical test for data taking values on nonlinear manifolds. The paper proposes statistical tests that are computationally tractable and geometrically sound for Diffusion Tensor Imaging.
An important property of the Kalman filter is that the underlying Riccati flow is a contraction for the natural metric of the cone of symmetric positive definite matrices. The present paper studies the geometry of a low-rank version of the Kalman filter. The underlying Riccati flow evolves on the manifold of fixed rank symmetric positive semidefinite matrices. Contraction properties of the low-rank flow are studied by means of a suitable metric recently introduced by the authors.
Statistical analysis of Diffusion Tensor Imaging (DTI) data requires a computational framework that is both numerically tractable (to account for the high dimensional nature of the data) and geometric (to account for the nonlinear nature of diffusion tensors). Building upon earlier studies that have shown that a Riemannian framework is appropriate to address these challenges, the present paper proposes a novel metric and an accompanying computational framework for DTI data processing. The proposed metric retains the geometry and the computational tractability of earlier methods grounded in the affine invariant metric. In addition, and in contrast to earlier methods, it provides an interpolation method which preserves anisotropy, a central information carried by diffusion tensor data.
The paper studies the excitability properties of a generalized FitzHugh-Nagumo model. The model differs from the purely competitive FitzHugh-Nagumo model in that it accounts for the effect of cooperative gating variables such as activation of calcium currents. Excitability is explored by unfolding a pitchfork bifurcation that is shown to organize five different types of excitability. In addition to the three classical types of neuronal excitability, two novel types are described and distinctly associated to the presence of cooperative variables.
This paper establishes a global contraction property for networks of phase-coupled oscillators characterized by a monotone coupling function. The contraction measure is a total variation distance. The contraction property determines the asymptotic behavior of the network, which is either finite-time synchronization or asymptotic convergence to a splay state.
We consider a continuum of phase oscillators on the circle interacting through an impulsive instantaneous coupling. In contrast with previous studies on related pulse-coupled models, the stability results obtained in the continuum limit are global. For the nonlinear transport equation governing the evolution of the oscillators, we propose (under technical assumptions) a global Lyapunov function which is induced by a total variation distance between quantile densities. The monotone time evolution of the Lyapunov function completely characterizes the dichotomic behavior of the oscillators: either the oscillators converge in finite time to a synchronous state or they asymptotically converge to an asynchronous state uniformly spread on the circle. The results of the present paper apply to popular phase oscillators models (e.g. the well-known leaky integrate-and-fire model) and draw a strong parallel between the analysis of finite and infinite populations. In addition, they provide a novel approach for the (global) analysis of pulse-coupled oscillators.
The generalization of the geometric mean of positive scalars to positive definite matrices has attracted considerable attention since the seminal work of Ando. The paper generalizes this framework of matrix means by proposing the definition of a rank-preserving mean for two or an arbitrary number of positive semi-definite matrices of fixed rank. The proposed mean is shown to be geometric in that it satisfies all the expected properties of a rank-preserving geometric mean. The work is motivated by operations on low-rank approximations of positive definite matrices in high-dimensional spaces.
Fifty years ago, Fitzugh introduced a phase portrait that became famous for a twofold reason: it captured in a physiological way the qualitative behavior of Hodgkin-Huxley model and it revealed the power of simple dynamical models to unfold complex firing patterns. To date, in spite of the enormous progresses in qualitative and quantitative neural modeling, this phase portrait has remained the core picture of neuronal excitability. Yet, a major difference between the neurophysiology of 1961 and of 2011 is the recognition of the prominent role of calcium channels in firing mechanisms. We show that including this extra current in Hodgkin-Huxley dynamics leads to a revision of Fitzugh-Nagumo phase portrait that affects in a fundamental way the reduced modeling of neural excitability. The revisited model considerably enlarges the modeling power of the original one. In particular, it captures essential electrophysiological signatures that otherwise require non-physiological alteration or considerable complexication of the classical model. As a basic illustration, the new model is shown to highlight a core dynamical mechanism by which the calcium conductance controls the two distinct firing modes of thalamocortical neurons.
The paper addresses the problem of learning a regression model parameterized by a fixed-rank positive semidefinite matrix. The focus is on the nonlinear nature of the search space and on scalability to high-dimensional problems. The mathematical developments rely on the theory of gradient descent algorithms adapted to the Riemannian geometry that underlies the set of fixed-rank positive semidefinite matrices. In contrast with previous contributions in the literature, no restrictions are imposed on the range space of the learned matrix. The resulting algorithms maintain a linear complexity in the problem size and enjoy important invariance properties. We apply the proposed algorithms to the problem of learning a distance function parameterized by a positive semidefinite matrix. Good performance is observed on classical benchmarks.
Convergence analysis of consensus algorithms is revisited in the light of the Hilbert distance. Tsitsiklis Lyapunov function is shown to be the Hilbert distance to consensus in log coordinates. Birkhoff theorem, which proves contraction of the Hilbert metric for any positive homogeneous monotone map, provides an early yet general convergence result for consensus algorithms. Because Birkhoff theorem holds in arbitrary cones, we extend consensus algorithms to the cone of positive definite matrices. The proposed generalization finds applications in the convergence analysis of quantum stochastic maps, which are a generalization of stochastic maps to non-commutative probability spaces.