Region of attraction analysis of nonlinear stochastic systems using Polynomial Chaos Expansion

A method is presented to estimate the region of attraction (ROA) of stochastic systems with finite second moment and uncertainty-dependent equilibria. The approach employs Polynomial Chaos (PC) expansions to represent the stochastic system by a higher-dimensional set of deterministic equations. We first show how the equilibrium point of the deterministic formulation provides the stochastic moments of an uncertainty-dependent equilibrium point of the stochastic system. A connection between the boundedness of the moments of the stochastic system and the Lyapunov stability of its PC expansion is then derived. Defining corresponding notions of a ROA for both system representations, we show how this connection can be leveraged to recover an estimate of the ROA of the stochastic system from the ROA of the PC expanded system. Two optimization programs, obtained from sum-of-squares programming techniques, are provided to compute inner estimates of the ROA. The first optimization program uses the Lyapunov stability arguments to return an estimate of the ROA of the PC expansion. Based on this result and user specifications on the moments for the initial conditions, the second one employs the shown connection to provide the corresponding ROA of the stochastic system. The method is demonstrated by two examples.

[1]  James Anderson,et al.  Region of Attraction Estimation Using Invariant Sets and Rational Lyapunov Functions , 2016, Autom..

[2]  Kwang-Ki K. Kim,et al.  Wiener's Polynomial Chaos for the Analysis and Control of Nonlinear Dynamical Systems with Probabilistic Uncertainties [Historical Perspectives] , 2013, IEEE Control Systems.

[3]  Mircea Grigoriu,et al.  On the accuracy of the polynomial chaos approximation for random variables and stationary stochastic processes. , 2003 .

[4]  Franz S. Hover,et al.  Application of polynomial chaos in stability and control , 2006, Autom..

[5]  Graziano Chesi Estimating the domain of attraction for uncertain polynomial systems , 2004, Autom..

[6]  Shu-jinWu,et al.  Boundedness of Nonlinear Differential Systems with Impulsive Effect on Random Moments , 2004 .

[7]  T. Sullivan Introduction to Uncertainty Quantification , 2015 .

[8]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[9]  Andrés Marcos,et al.  An equilibrium-independent region of attraction formulation for systems with uncertainty-dependent equilibria , 2018, 2018 IEEE Conference on Decision and Control (CDC).

[10]  A. Garulli,et al.  LMI‐based computation of optimal quadratic Lyapunov functions for odd polynomial systems , 2005 .

[11]  Johan Löfberg,et al.  Pre- and Post-Processing Sum-of-Squares Programs in Practice , 2009, IEEE Transactions on Automatic Control.

[12]  Lorenzo Fagiano,et al.  On the guaranteed accuracy of Polynomial Chaos Expansions , 2011, IEEE Conference on Decision and Control and European Control Conference.

[13]  Omar M. Knio,et al.  Spectral Methods for Uncertainty Quantification , 2010 .

[14]  R. Khasminskii Stochastic Stability of Differential Equations , 1980 .

[15]  A. Packard,et al.  Stability region analysis using sum of squares programming , 2006, 2006 American Control Conference.

[16]  Raktim Bhattacharya,et al.  Linear quadratic regulation of systems with stochastic parameter uncertainties , 2009, Autom..

[17]  G. Stengle A nullstellensatz and a positivstellensatz in semialgebraic geometry , 1974 .

[18]  Ufuk Topcu,et al.  Robust Region-of-Attraction Estimation , 2010, IEEE Transactions on Automatic Control.

[19]  J. Fisher,et al.  Stability analysis of stochastic systems using polynomial chaos , 2008, 2008 American Control Conference.

[20]  P. Olver Nonlinear Systems , 2013 .

[21]  Veit Hagenmeyer,et al.  Comments on Truncation Errors for Polynomial Chaos Expansions , 2017, IEEE Control Systems Letters.

[22]  P. Parrilo Structured semidefinite programs and semialgebraic geometry methods in robustness and optimization , 2000 .

[23]  Andrew Packard,et al.  Control Applications of Sum of Squares Programming , 2005 .

[24]  Amir Ali Ahmadi,et al.  DSOS and SDSOS Optimization: More Tractable Alternatives to Sum of Squares and Semidefinite Optimization , 2017, SIAM J. Appl. Algebra Geom..

[25]  Tony A. Wood,et al.  Stability Verification for Periodic Trajectories of Autonomous Kite Power Systems , 2018, 2018 European Control Conference (ECC).

[26]  Richard D. Braatz,et al.  On stability of stochastic linear systems via polynomial chaos expansions , 2017, 2017 American Control Conference (ACC).

[27]  Tony A. Wood,et al.  Transverse Contraction-Based Stability Analysis for Periodic Trajectories of Controlled Power Kites with Model Uncertainty , 2018, 2018 IEEE Conference on Decision and Control (CDC).

[28]  Ryan Feeley,et al.  Some controls applications of sum of squares programming , 2003, 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475).

[29]  Frank Kozin,et al.  A survey of stability of stochastic systems , 1969, Autom..

[30]  Ufuk Topcu,et al.  Local stability analysis using simulations and sum-of-squares programming , 2008, Autom..

[31]  Skuli Gudmundsson,et al.  Probabilistic Basin of Attraction and Its Estimation Using Two Lyapunov Functions , 2018, Complex..

[32]  Jean-Jacques Slotine,et al.  Wasserstein Contraction of Stochastic Nonlinear Systems , 2019, 1902.08567.

[33]  D. Xiu,et al.  Modeling uncertainty in flow simulations via generalized polynomial chaos , 2003 .

[34]  Russ Tedrake,et al.  Finite-time regional verification of stochastic non-linear systems , 2011, Int. J. Robotics Res..

[35]  Peter Seiler,et al.  Region of attraction analysis with Integral Quadratic Constraints , 2019, Autom..

[36]  Nicolas Tabareau,et al.  A Contraction Theory Approach to Stochastic Incremental Stability , 2007, IEEE Transactions on Automatic Control.