June  2016, 3(2): 191-210. doi: 10.3934/jcd.2016010

Determination of the basin of attraction of a periodic orbit in two dimensions using meshless collocation

1. 

Department of Mathematics, University of Sussex, Falmer BN1 9QH, United Kingdom

Received  December 2016 Revised  February 2017 Published  April 2017

A contraction metric for an autonomous ordinary differential equation is a Riemannian metric such that the distance between adjacent solutions contracts over time. A contraction metric can be used to determine the basin of attraction of a periodic orbit without requiring information about its position or stability. Moreover, it is robust to small perturbations of the system.
    In two-dimensional systems, a contraction metric can be characterised by a scalar-valued function. In [9], the function was constructed as solution of a first-order linear Partial Differential Equation (PDE), and numerically constructed using meshless collocation. However, information about the periodic orbit was required, which needed to be approximated.
    In this paper, we overcome this requirement by studying a second-order PDE, which does not require any information about the periodic orbit. We show that the second-order PDE has a solution, which defines a contraction metric. We use meshless collocation to approximate the solution and prove error estimates. In particular, we show that the approximation itself is a contraction metric, if the collocation points are dense enough. The method is applied to two examples.
Citation: Peter Giesl, James McMichen. Determination of the basin of attraction of a periodic orbit in two dimensions using meshless collocation. Journal of Computational Dynamics, 2016, 3 (2) : 191-210. doi: 10.3934/jcd.2016010
References:
[1]

D. Angeli, A Lyapunov approach to incremental stability properties,, IEEE Trans. Automat. Contr., 47 (2002), 410. doi: 10.1109/9.989067. Google Scholar

[2]

E. Aylward, P. Parrilo and J.-J. Slotine, Stability and robustness analysis of nonlinear systems via contraction metrics and SOS programming,, Automatica, 44 (2008), 2163. doi: 10.1016/j.automatica.2007.12.012. Google Scholar

[3]

V. Boichenko, G. Leonov and V. Reitmann, Dimension Theory for Ordinary Differential Equations,, volume 141 of Teubner-Texte zur Mathematik [Teubner Texts in Mathematics]. B. G. Teubner Verlagsgesellschaft mbH, (2005). doi: 10.1007/978-3-322-80055-8. Google Scholar

[4]

G. Borg, A Condition for the Existence of Orbitally Stable Solutions of Dynamical Systems,, Kungliga Tekniska Högskolan Handlingar Stockholm 153, (1960). Google Scholar

[5]

M. Buhmann, Radial Basis Functions: Theory and Implementations,, volume 12 of Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, (2003). doi: 10.1017/CBO9780511543241. Google Scholar

[6]

F. Forni and R. Sepulchre, A differential Lyapunov framework for Contraction Analysis,, IEEE Trans. Automat. Control, 59 (2014), 614. doi: 10.1109/TAC.2013.2285771. Google Scholar

[7]

P. Giesl, Necessary conditions for a limit cycle and its basin of attraction,, Nonlinear Anal., 56 (2004), 643. doi: 10.1016/j.na.2003.07.020. Google Scholar

[8]

P. Giesl, Construction of Global Lyapunov Functions Using Radial Basis Functions,, volume 1904 of Lecture Notes in Mathematics. Springer, (1904). Google Scholar

[9]

P. Giesl, On the determination of the basin of attraction of a periodic orbit in two-dimensional systems,, Journal of Mathematical Analysis and Applications, 335 (2007), 461. doi: 10.1016/j.jmaa.2007.01.069. Google Scholar

[10]

P. Giesl, On the determination of the basin of attraction of periodic orbits in three- and higher-dimensional systems,, J. Math. Anal. Appl., 354 (2009), 606. doi: 10.1016/j.jmaa.2009.01.027. Google Scholar

[11]

P. Giesl and S. Hafstein, Construction of a CPA contraction metric for periodic orbits using semidefinite optimization,, Nonlinear Anal., 86 (2013), 114. doi: 10.1016/j.na.2013.03.012. Google Scholar

[12]

P. Giesl and S. Hafstein, Review on computational methods for Lyapunov functions,, Discrete Contin. Dyn. Syst. Ser. B, 20 (2015), 2291. doi: 10.3934/dcdsb.2015.20.2291. Google Scholar

[13]

P. Giesl and H. Wendland, Meshless collocation: Error estimates with application to dynamical systems,, SIAM J. Numer. Anal., 45 (2007), 1723. doi: 10.1137/060658813. Google Scholar

[14]

P. Hartman, Ordinary Differential Equations,, Wiley, (1964). Google Scholar

[15]

P. Hartman and C. Olech, On global asymptotic stability of solutions of differential equations,, Trans. Amer. Math. Soc., 104 (1962), 154. doi: 10.2307/1993939. Google Scholar

[16]

A. Iske, Perfect Centre Placement for Radial Basis Function Methods,, Technical report, (1999). Google Scholar

[17]

G. Leonov, I. Burkin and A. Shepelyavyi, Frequency Methods in Oscillation Theory,, Ser. Math. and its Appl., (1996). doi: 10.1007/978-94-009-0193-3. Google Scholar

[18]

D. Lewis, Metric properties of differential equations,, Amer. J. Math., 71 (1949), 294. doi: 10.2307/2372245. Google Scholar

[19]

W. Lohmiller and J.-J. Slotine, On contraction analysis for non-linear systems,, Automatica, 34 (1998), 683. doi: 10.1016/S0005-1098(98)00019-3. Google Scholar

[20]

I. Manchester and J.-J. Slotine, Transverse contraction criteria for existence, stability, and robustness of a limit cycle,, Systems Control Lett., 63 (2014), 32. doi: 10.1016/j.sysconle.2013.10.005. Google Scholar

[21]

J. McMichen, Determination of Areas and Basins of Attraction in Planar Dynamical Systems using Meshless Collocation,, PhD thesis, (2016). Google Scholar

[22]

B. Stenström, Dynamical systems with a certain local contraction property,, Math. Scand., 11 (1962), 151. doi: 10.7146/math.scand.a-10661. Google Scholar

[23]

H. Wendland, Error estimates for interpolation by compactly supported radial basis functions of minimal degree,, Journal of Approximation Theory, 93 (1998), 258. doi: 10.1006/jath.1997.3137. Google Scholar

[24]

H. Wendland, Scattered Data Approximation,, volume 17 of Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, (2005). Google Scholar

show all references

References:
[1]

D. Angeli, A Lyapunov approach to incremental stability properties,, IEEE Trans. Automat. Contr., 47 (2002), 410. doi: 10.1109/9.989067. Google Scholar

[2]

E. Aylward, P. Parrilo and J.-J. Slotine, Stability and robustness analysis of nonlinear systems via contraction metrics and SOS programming,, Automatica, 44 (2008), 2163. doi: 10.1016/j.automatica.2007.12.012. Google Scholar

[3]

V. Boichenko, G. Leonov and V. Reitmann, Dimension Theory for Ordinary Differential Equations,, volume 141 of Teubner-Texte zur Mathematik [Teubner Texts in Mathematics]. B. G. Teubner Verlagsgesellschaft mbH, (2005). doi: 10.1007/978-3-322-80055-8. Google Scholar

[4]

G. Borg, A Condition for the Existence of Orbitally Stable Solutions of Dynamical Systems,, Kungliga Tekniska Högskolan Handlingar Stockholm 153, (1960). Google Scholar

[5]

M. Buhmann, Radial Basis Functions: Theory and Implementations,, volume 12 of Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, (2003). doi: 10.1017/CBO9780511543241. Google Scholar

[6]

F. Forni and R. Sepulchre, A differential Lyapunov framework for Contraction Analysis,, IEEE Trans. Automat. Control, 59 (2014), 614. doi: 10.1109/TAC.2013.2285771. Google Scholar

[7]

P. Giesl, Necessary conditions for a limit cycle and its basin of attraction,, Nonlinear Anal., 56 (2004), 643. doi: 10.1016/j.na.2003.07.020. Google Scholar

[8]

P. Giesl, Construction of Global Lyapunov Functions Using Radial Basis Functions,, volume 1904 of Lecture Notes in Mathematics. Springer, (1904). Google Scholar

[9]

P. Giesl, On the determination of the basin of attraction of a periodic orbit in two-dimensional systems,, Journal of Mathematical Analysis and Applications, 335 (2007), 461. doi: 10.1016/j.jmaa.2007.01.069. Google Scholar

[10]

P. Giesl, On the determination of the basin of attraction of periodic orbits in three- and higher-dimensional systems,, J. Math. Anal. Appl., 354 (2009), 606. doi: 10.1016/j.jmaa.2009.01.027. Google Scholar

[11]

P. Giesl and S. Hafstein, Construction of a CPA contraction metric for periodic orbits using semidefinite optimization,, Nonlinear Anal., 86 (2013), 114. doi: 10.1016/j.na.2013.03.012. Google Scholar

[12]

P. Giesl and S. Hafstein, Review on computational methods for Lyapunov functions,, Discrete Contin. Dyn. Syst. Ser. B, 20 (2015), 2291. doi: 10.3934/dcdsb.2015.20.2291. Google Scholar

[13]

P. Giesl and H. Wendland, Meshless collocation: Error estimates with application to dynamical systems,, SIAM J. Numer. Anal., 45 (2007), 1723. doi: 10.1137/060658813. Google Scholar

[14]

P. Hartman, Ordinary Differential Equations,, Wiley, (1964). Google Scholar

[15]

P. Hartman and C. Olech, On global asymptotic stability of solutions of differential equations,, Trans. Amer. Math. Soc., 104 (1962), 154. doi: 10.2307/1993939. Google Scholar

[16]

A. Iske, Perfect Centre Placement for Radial Basis Function Methods,, Technical report, (1999). Google Scholar

[17]

G. Leonov, I. Burkin and A. Shepelyavyi, Frequency Methods in Oscillation Theory,, Ser. Math. and its Appl., (1996). doi: 10.1007/978-94-009-0193-3. Google Scholar

[18]

D. Lewis, Metric properties of differential equations,, Amer. J. Math., 71 (1949), 294. doi: 10.2307/2372245. Google Scholar

[19]

W. Lohmiller and J.-J. Slotine, On contraction analysis for non-linear systems,, Automatica, 34 (1998), 683. doi: 10.1016/S0005-1098(98)00019-3. Google Scholar

[20]

I. Manchester and J.-J. Slotine, Transverse contraction criteria for existence, stability, and robustness of a limit cycle,, Systems Control Lett., 63 (2014), 32. doi: 10.1016/j.sysconle.2013.10.005. Google Scholar

[21]

J. McMichen, Determination of Areas and Basins of Attraction in Planar Dynamical Systems using Meshless Collocation,, PhD thesis, (2016). Google Scholar

[22]

B. Stenström, Dynamical systems with a certain local contraction property,, Math. Scand., 11 (1962), 151. doi: 10.7146/math.scand.a-10661. Google Scholar

[23]

H. Wendland, Error estimates for interpolation by compactly supported radial basis functions of minimal degree,, Journal of Approximation Theory, 93 (1998), 258. doi: 10.1006/jath.1997.3137. Google Scholar

[24]

H. Wendland, Scattered Data Approximation,, volume 17 of Cambridge Monographs on Applied and Computational Mathematics. Cambridge University Press, (2005). Google Scholar

[1]

Peter Giesl, Holger Wendland. Construction of a contraction metric by meshless collocation. Discrete & Continuous Dynamical Systems - B, 2019, 24 (8) : 3843-3863. doi: 10.3934/dcdsb.2018333

[2]

Peter Giesl, Holger Wendland. Approximating the basin of attraction of time-periodic ODEs by meshless collocation. Discrete & Continuous Dynamical Systems - A, 2009, 25 (4) : 1249-1274. doi: 10.3934/dcds.2009.25.1249

[3]

Peter Giesl, Holger Wendland. Approximating the basin of attraction of time-periodic ODEs by meshless collocation of a Cauchy problem. Conference Publications, 2009, 2009 (Special) : 259-268. doi: 10.3934/proc.2009.2009.259

[4]

Peter Giesl. Converse theorem on a global contraction metric for a periodic orbit. Discrete & Continuous Dynamical Systems - A, 2019, 39 (9) : 5339-5363. doi: 10.3934/dcds.2019218

[5]

Peter Giesl. Necessary condition for the basin of attraction of a periodic orbit in non-smooth periodic systems. Discrete & Continuous Dynamical Systems - A, 2007, 18 (2&3) : 355-373. doi: 10.3934/dcds.2007.18.355

[6]

Peter Giesl, James McMichen. Determination of the area of exponential attraction in one-dimensional finite-time systems using meshless collocation. Discrete & Continuous Dynamical Systems - B, 2018, 23 (4) : 1835-1850. doi: 10.3934/dcdsb.2018094

[7]

Gerard Gómez, Josep–Maria Mondelo, Carles Simó. A collocation method for the numerical Fourier analysis of quasi-periodic functions. II: Analytical error estimates. Discrete & Continuous Dynamical Systems - B, 2010, 14 (1) : 75-109. doi: 10.3934/dcdsb.2010.14.75

[8]

Peter Giesl. Construction of a finite-time Lyapunov function by meshless collocation. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2387-2412. doi: 10.3934/dcdsb.2012.17.2387

[9]

Hjörtur Björnsson, Sigurdur Hafstein, Peter Giesl, Enrico Scalas, Skuli Gudmundsson. Computation of the stochastic basin of attraction by rigorous construction of a Lyapunov function. Discrete & Continuous Dynamical Systems - B, 2019, 24 (8) : 4247-4269. doi: 10.3934/dcdsb.2019080

[10]

Helge Dietert, Josephine Evans, Thomas Holding. Contraction in the Wasserstein metric for the kinetic Fokker-Planck equation on the torus. Kinetic & Related Models, 2018, 11 (6) : 1427-1441. doi: 10.3934/krm.2018056

[11]

Li-Ming Yeh. Pointwise estimate for elliptic equations in periodic perforated domains. Communications on Pure & Applied Analysis, 2015, 14 (5) : 1961-1986. doi: 10.3934/cpaa.2015.14.1961

[12]

Giuseppe Savaré. Self-improvement of the Bakry-Émery condition and Wasserstein contraction of the heat flow in $RCD (K, \infty)$ metric measure spaces. Discrete & Continuous Dynamical Systems - A, 2014, 34 (4) : 1641-1661. doi: 10.3934/dcds.2014.34.1641

[13]

Alexander Pankov. Nonlinear Schrödinger Equations on Periodic Metric Graphs. Discrete & Continuous Dynamical Systems - A, 2018, 38 (2) : 697-714. doi: 10.3934/dcds.2018030

[14]

Huan-Zhen Chen, Zhao-Jie Zhou, Hong Wang, Hong-Ying Man. An optimal-order error estimate for a family of characteristic-mixed methods to transient convection-diffusion problems. Discrete & Continuous Dynamical Systems - B, 2011, 15 (2) : 325-341. doi: 10.3934/dcdsb.2011.15.325

[15]

Xingwen Hao, Yachun Li, Qin Wang. A kinetic approach to error estimate for nonautonomous anisotropic degenerate parabolic-hyperbolic equations. Kinetic & Related Models, 2014, 7 (3) : 477-492. doi: 10.3934/krm.2014.7.477

[16]

Patrick Henning, Mario Ohlberger. A-posteriori error estimate for a heterogeneous multiscale approximation of advection-diffusion problems with large expected drift. Discrete & Continuous Dynamical Systems - S, 2016, 9 (5) : 1393-1420. doi: 10.3934/dcdss.2016056

[17]

Anete S. Cavalcanti. An existence proof of a symmetric periodic orbit in the octahedral six-body problem. Discrete & Continuous Dynamical Systems - A, 2017, 37 (4) : 1903-1922. doi: 10.3934/dcds.2017080

[18]

Tatiane C. Batista, Juliano S. Gonschorowski, Fábio A. Tal. Density of the set of endomorphisms with a maximizing measure supported on a periodic orbit. Discrete & Continuous Dynamical Systems - A, 2015, 35 (8) : 3315-3326. doi: 10.3934/dcds.2015.35.3315

[19]

Xueting Tian, Shirou Wang, Xiaodong Wang. Intermediate Lyapunov exponents for systems with periodic orbit gluing property. Discrete & Continuous Dynamical Systems - A, 2019, 39 (2) : 1019-1032. doi: 10.3934/dcds.2019042

[20]

Gerard Gómez, Josep–Maria Mondelo, Carles Simó. A collocation method for the numerical Fourier analysis of quasi-periodic functions. I: Numerical tests and examples. Discrete & Continuous Dynamical Systems - B, 2010, 14 (1) : 41-74. doi: 10.3934/dcdsb.2010.14.41

 Impact Factor: 

Metrics

  • PDF downloads (13)
  • HTML views (0)
  • Cited by (0)

Other articles
by authors

[Back to Top]