Polynomial Optimization with Applications to Stability Analysis and Control - Alternatives to Sum of Squares

In this paper, we explore the merits of various algorithms for polynomial optimization problems, focusing on alternatives to sum of squares programming. While we refer to advantages and disadvantages of Quantifier Elimination, Reformulation Linear Techniques, Blossoming and Groebner basis methods, our main focus is on algorithms defined by Polya's theorem, Bernstein's theorem and Handelman's theorem. We first formulate polynomial optimization problems as verifying the feasibility of semi-algebraic sets. Then, we discuss how Polya's algorithm, Bernstein's algorithm and Handelman's algorithm reduce the intractable problem of feasibility of semi-algebraic sets to linear and/or semi-definite programming. We apply these algorithms to different problems in robust stability analysis and stability of nonlinear dynamical systems. As one contribution of this paper, we apply Polya's algorithm to the problem of H_infinity control of systems with parametric uncertainty. Numerical examples are provided to compare the accuracy of these algorithms with other polynomial optimization algorithms in the literature.Comment: AIMS Journal of Discrete and Continuous Dynamical Systems - Series

A robust LMI approach on nonlinear feedback stabilization of continuous state-delay systems with Lipschitzian nonlinearities : experimental validation

This paper suggests a novel nonlinear state-fe edback stabilization control law using linear matrix inequalities for a class oftime-delayed nonlinear dynamic systems with Lipschitz nonlinearity conditions. Based on the Lyapunov–Krasovskiistability theory, the asymptotic stabilization criterion is derived in the linear matrix inequality form and the coef¿cients ofthe nonlinear state-feedback controller are determined. Meanwhile, an appropriate criterion to ¿nd the proper feedbackgain matrix F is also provided. The robustness purpose against nonlinear functions and time delays is guaranteed in thisscheme. Moreover , the problem of robust H!performance analysis for a class of nonlinear time-delayed system s withexternal disturbance is studied in this paper. Simulations are presented to demonstrate the pro¿ciency of the offeredtechnique. For this purpos e, an unstable nonlinear numerical system and a rotary inverted pendulum system have beenstudied in the simulation section. Moreover, an experimental study of the practical rotary inverted pendul um system isprovided. These results con¿rm the expected satisfactory performance of the suggested method.Peer ReviewedPostprint (author's final draft

An improved stability criterion for discrete-time time-delayed Lur’e systemwith sector-bounded nonlinearities

The absolute stability problem of discrete-time time-delayed Lur\u27e systems with sector bounded nonlinearities is investigated in this paper. Firstly, a modified Lyapunov-Krasovskii functional (LKF) is designed with augmenting additional double summation terms, which complements more coupling information between the delay intervals and other system state variables than some previous LKFs. Secondly, some improved delay-dependent absolute stability criteria based on linear matrix inequality form (LMI) are proposed via the modified LKF and the relaxed free-matrix-based summation inequality technique application. The stability criteria are less conservative than some results previously proposed. The reduction of the conservatism mainly relies on the full use of the relaxed summation inequality technique based on the modified LKF. Finally, two common numerical examples are presented to show the effectiveness of the proposed approach

Stability analysis of linear ODE-PDE interconnected systems

Les systèmes de dimension infinie permettent de modéliser un large spectre de phénomènes physiques pour lesquels les variables d'états évoluent temporellement et spatialement. Ce manuscrit s'intéresse à l'évaluation de la stabilité de leur point d'équilibre. Deux études de cas seront en particulier traitées : l'analyse de stabilité des systèmes interconnectés à une équation de transport, et à une équation de réaction-diffusion. Des outils théoriques existent pour l'analyse de stabilité de ces systèmes linéaires de dimension infinie et s'appuient sur une algèbre d'opérateurs plutôt que matricielle. Cependant, ces résultats d'existence soulèvent un problème de constructibilité numérique. Lors de l'implémentation, une approximation est réalisée et les résultats sont conservatifs. La conception d'outils numériques menant à des garanties de stabilité pour lesquelles le degré de conservatisme est évalué et maîtrisé est alors un enjeu majeur. Comment développer des critères numériques fiables permettant de statuer sur la stabilité ou l'instabilité des systèmes linéaires de dimension infinie ? Afin de répondre à cette question, nous proposons ici une nouvelle méthode générique qui se décompose en deux temps. D'abord, sous l'angle de l'approximation sur les polynômes de Legendre, des modèles augmentés sont construits et découpent le système original en deux blocs : d'une part, un système de dimension finie approximant est isolé, d'autre part, l'erreur de troncature de dimension infinie est conservée et modélisée. Ensuite, des outils fréquentiels et temporels de dimension finie sont déployés afin de proposer des critères de stabilité plus ou moins coûteux numériquement en fonction de l'ordre d'approximation choisi. En fréquentiel, à l'aide du théorème du petit gain, des conditions suffisantes de stabilité sont obtenues. En temporel, à l'aide du théorème de Lyapunov, une sous-estimation des régions de stabilité est proposée sous forme d'inégalité matricielle linéaire et une sur-estimation sous forme de test de positivité. Nos deux études de cas ont ainsi été traitées à l'aide de cette méthodologie générale. Le principal résultat obtenu concerne le cas des systèmes EDO-transport interconnectés, pour lequel l'approximation et l'analyse de stabilité à l'aide des polynômes de Legendre mène à des estimations des régions de stabilité qui convergent exponentiellement vite. La méthode développée dans ce manuscrit peut être adaptée à d'autres types d'approximations et exportée à d'autres systèmes linéaires de dimension infinie. Ce travail ouvre ainsi la voie à l'obtention de conditions nécessaires et suffisantes de stabilité de dimension finie pour les systèmes de dimension infinie.Infinite dimensional systems allow to model a large panel of physical phenomena for which the state variables evolve both temporally and spatially. This manuscript deals with the evaluation of the stability of their equilibrium point. Two case studies are treated in particular: the stability analysis of ODE-transport, and ODE-reaction-diffusion interconnected systems. Theoretical tools exist for the stability analysis of these infinite-dimensional linear systems and are based on an operator algebra rather than a matrix algebra. However, these existence results raise a problem of numerical constructibility. During implementation, an approximation is performed and the results are conservative. The design of numerical tools leading to stability guarantees for which the degree of conservatism is evaluated and controlled is then a major issue. How can we develop reliable numerical criteria to rule on the stability or instability of infinite-dimensional linear systems? In order to answer this question, one proposes here a new generic method, which is decomposed in two steps. First, from the perspective of Legendre polynomials approximation, augmented models are built and split the original system into two blocks: on the one hand, a finite-dimensional approximated system is isolated, on the other hand, the infinite-dimensional truncation error is preserved and modeled. Then, frequency and time tools of finite dimension are deployed in order to propose stability criteria that have high or low numerical load depending on the approximated order. In frequencies, with the aid of the small gain theorem, sufficient stability conditions are obtained. In temporal, with the aid of the Lyapunov theorem, an under estimate of the stability regions is proposed as a linear matrix inequality and an over estimate as a positivity test. Our two case studies have been treated with this general methodology. The main result concerns the case of ODE-transport interconnected systems, for which the approximation and stability analysis using Legendre polynomials leads to exponentially fast converging estimates of stability regions. The method developed in this manuscript can be adapted to other types of approximations and exported to other infinite-dimensional linear systems. Thus, this work opens the way to obtain necessary and sufficient finite-dimensional conditions of stability for infinite-dimensional systems

Distributed H-infinity filtering for polynomial nonlinear stochastic systems in sensor networks

Copyright [2010] IEEE. This material is posted here with permission of the IEEE. Such permission of the IEEE does not in any way imply IEEE endorsement of any of Brunel University's products or services. Internal or personal use of this material is permitted. However, permission to reprint/republish this material for advertising or promotional purposes or for creating new collective works for resale or redistribution must be obtained from the IEEE by writing to [email protected]. By choosing to view this document, you agree to all provisions of the copyright laws protecting it.In this paper, the distributed H1 filtering problem is addressed for a class of polynomial nonlinear stochastic systems in sensor networks. For a Lyapunov function candidate whose entries are polynomials, we calculate its first- and second-order derivatives in order to facilitate the use of Itos differential role. Then, a sufficient condition for the existence of a feasible solution to the addressed distributed H1 filtering problem is derived in terms of parameter-dependent linear matrix inequalities (PDLMIs). For computational convenience, these PDLMIs are further converted into a set of sums of squares (SOSs) that can be solved effectively by using the semidefinite programming technique. Finally, a numerical simulation example is provided to demonstrate the effectiveness and applicability of the proposed design approach.This work was supported in part by the Engineering and Physical Sciences Research Council (EPSRC) of the U.K. under Grant GR/S27658/01, the Royal Society of the U.K., the National 973 Program of China under Grant 2009CB320600, the National Natural Science Foundation of China under Grant 60974030 and the Alexander von Humboldt Foundation of Germany

On the necessity of sufficient LMI conditions for time-delay systems arising from Legendre approximation

This work is dedicated to the stability analysis of time-delay systems with a single constant delay using the Lyapunov-Krasovskii theorem. This approach has been widely used in the literature and numerous sufficient conditions of stability have been proposed and expressed as linear matrix inequalities (LMI). The main criticism of the method that is often pointed out is that these LMI conditions are only sufficient, and there is a lack of information regarding the reduction of the conservatism. Recently, scalable methods have been investigated using Bessel-Legendre inequality or orthogonal polynomial-based inequalities. The interest of these methods relies on their hierarchical structure with a guarantee of reduction of the level of conservatism. However, the convergence is still an open question that will be answered for the first time in this paper. The objective is to prove that the stability of a time-delay system implies the feasibility of these scalable LMI, at a sufficiently large order of the Legendre polynomials. Moreover, the proposed contribution is even able to provide an analytic estimation of this order, giving rise to a necessary and sufficient LMI for the stability of time-delay systems

Finite-region stability of 2-D singular Roesser systems with directional delays

In this paper, the problem of finite-region stability is studied for a class of two-dimensional (2-D) singular systems described by using the Roesser model with directional delays. Based on the regularity, we first decompose the underlying singular 2-D systems into fast and slow subsystems corresponding to dynamic and algebraic parts. Then, with the Lyapunov-like 2-D functional method, we construct a weighted 2-D functional candidate and utilize zero-type free matrix equations to derive delay-dependent stability conditions in terms of linear matrix inequalities (LMIs). More specifically, the derived conditions ensure that all state trajectories of the system do not exceed a prescribed threshold over a pre-specified finite region of time for any initial state sequences when energy-norms of dynamic parts do not exceed given bounds