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

Ufuk Topcu, Andrew Packard, Peter Seiler

Research output: Contribution to journalArticlepeer-review

152 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 1 2008


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

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

Cite this