Local stability analysis using simulations and sum-of-squares programming

Ufuk Topcu, Andrew Packard, Peter Seiler

Research output: Contribution to journalArticlepeer-review

190 Scopus citations


The problem of computing bounds on the region-of-attraction for systems with polynomial vector fields is considered. Invariant subsets of the region-of-attraction are characterized as sublevel sets of Lyapunov functions. Finite-dimensional polynomial parametrizations for Lyapunov functions are used. A methodology utilizing information from simulations to generate Lyapunov function candidates satisfying necessary conditions for bilinear constraints is proposed. The suitability of Lyapunov function candidates is assessed solving linear sum-of-squares optimization problems. Qualified candidates are used to compute invariant subsets of the region-of-attraction and to initialize various bilinear search strategies for further optimization. We illustrate the method on small examples from the literature and several control oriented systems.

Original languageEnglish (US)
Pages (from-to)2669-2675
Number of pages7
Issue number10
StatePublished - Oct 2008

Bibliographical note

Funding Information:
This work was sponsored by the Air Force Office of Scientific Research, USAF, under grant/contract number FA9550-05-1-0266. The views and conclusions contained herein are those of the authors and should not be interpreted as necessarily representing the official policies or endorsements, either expressed or implied, of the AFOSR or the US Government.


  • Local stability
  • Nonlinear dynamics
  • Region-of-attraction
  • Simulations
  • Sum-of-squares programming


Dive into the research topics of 'Local stability analysis using simulations and sum-of-squares programming'. Together they form a unique fingerprint.

Cite this