Contracting Nonlinear Observers: Convex Optimization and Learning from Data

A new approach to design of nonlinear observers (state estimators) is proposed. The main idea is to (i) construct a convex set of dynamical systems which are contracting observers for a particular system, and (ii) optimize over this set for one which minimizes a bound on state-estimation error on a simulated noisy data set. We construct convex sets of continuous-time and discrete-time observers, as well as contracting sampled-data observers for continuous-time systems. Convex bounds for learning are constructed using Lagrangian relaxation. The utility of the proposed methods are verified using numerical simulation.Comment: conference submissio

A Contraction Theory Approach to Stochastic Incremental Stability

We investigate the incremental stability properties of It\^o stochastic dynamical systems. Specifically, we derive a stochastic version of nonlinear contraction theory that provides a bound on the mean square distance between any two trajectories of a stochastically contracting system. This bound can be expressed as a function of the noise intensity and the contraction rate of the noise-free system. We illustrate these results in the contexts of stochastic nonlinear observers design and stochastic synchronization.Comment: 23 pages, 2 figure

Evaluation of stochastic effects on biomolecular networks using the generalised Nyquist stability criterion

Abstract—Stochastic differential equations are now commonly used to model biomolecular networks in systems biology, and much recent research has been devoted to the development of methods to analyse their stability properties. Stability analysis of such systems may be performed using the Laplace transform, which requires the calculation of the exponential matrix involving time symbolically. However, the calculation of the symbolic exponential matrix is not feasible for problems of even moderate size, as the required computation time increases exponentially with the matrix order. To address this issue, we present a novel method for approximating the Laplace transform which does not require the exponential matrix to be calculated explicitly. The calculation time associated with the proposed method does not increase exponentially with the size of the system, and the approximation error is shown to be of the same order as existing methods. Using this approximation method, we show how a straightforward application of the generalized Nyquist stability criterion provides necessary and sufficient conditions for the stability of stochastic biomolecular networks. The usefulness and computational efficiency of the proposed method is illustrated through its application to the problem of analysing a model for limit-cycle oscillations in cAMP during aggregation of Dictyostelium cells

Optimal control of ankle joint moment: Toward unsupported standing in paraplegia

This paper considers part of the problem of how to provide unsupported standing for paraplegics by feedback control. In this work our overall objective is to stabilize the subject by stimulation only of his ankle joints while the other joints are braced, Here, we investigate the problem of ankle joint moment control. The ankle plantarflexion muscles are first identified with pseudorandom binary sequence (PRBS) signals, periodic sinusoidal signals, and twitches. The muscle is modeled in Hammerstein form as a static recruitment nonlinearity followed by a linear transfer function. A linear-quadratic-Gaussian (LQG)-optimal controller design procedure for ankle joint moment was proposed based on the polynomial equation formulation, The approach was verified by experiments in the special Wobbler apparatus with a neurologically intact subject, and these experimental results are reported. The controller structure is formulated in such a way that there are only two scalar design parameters, each of which has a clear physical interpretation. This facilitates fast controller synthesis and tuning in the laboratory environment. Experimental results show the effects of the controller tuning parameters: the control weighting and the observer response time, which determine closed-loop properties. Using these two parameters the tradeoff between disturbance rejection and measurement noise sensitivity can be straightforwardly balanced while maintaining a desired speed of tracking. The experimentally measured reference tracking, disturbance rejection, and noise sensitivity are good and agree with theoretical expectations

Feedback control of unsupported standing in paraplegia. Part I: optimal control approach

This is the first of a pair of papers which describe an investigation into the feasibility of providing artificial balance to paraplegics using electrical stimulation of the paralyzed muscles. By bracing the body above the shanks, only stimulation of the plantarflexors is necessary. This arrangement prevents any influence from the intact neuromuscular system above the spinal cord lesion. Here, the authors extend the design of the controllers to a nested-loop LQG (linear quadratic Gaussian) stimulation controller which has ankle moment feedback (inner loops) and inverted pendulum angle feedback (outer loop). Each control loop is tuned by two parameters, the control weighting and an observer rise-time, which together determine the behavior. The nested structure was chosen because it is robust, despite changes in the muscle properties (fatigue) and interference from spasticity