# American Institute of Mathematical Sciences

January  2021, 8(1): 59-97. doi: 10.3934/jcd.2021004

## A general framework for validated continuation of periodic orbits in systems of polynomial ODEs

 1 VU Amsterdam, Department of Mathematics, De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 2 Rutgers University, Department of Mathematics, 110 Frelinghusen Rd, Piscataway, NJ 08854-8019, USA

Received  February 2019 Revised  November 2019 Published  October 2020

Fund Project: The authors are partially supported by NWO-VICI grant 639033109

In this paper a parametrized Newton-Kantorovich approach is applied to continuation of periodic orbits in arbitrary polynomial vector fields. This allows us to rigorously validate numerically computed branches of periodic solutions. We derive the estimates in full generality and present sample continuation proofs obtained using an implementation in Matlab. The presented approach is applicable to any polynomial vector field of any order and dimension. A variety of examples is presented to illustrate the efficacy of the method.

Citation: Jan Bouwe van den Berg, Elena Queirolo. A general framework for validated continuation of periodic orbits in systems of polynomial ODEs. Journal of Computational Dynamics, 2021, 8 (1) : 59-97. doi: 10.3934/jcd.2021004
##### References:

show all references

##### References:
Representation of the matrix introduced in Equation (25)
The shape of the operator $Q$, as defined in Lemma 5.1. The vertical columns represent the vectors $y^{[k]}$, where $k$ denotes the column index
The boundaries of the validation interval $(r_{\min},r_{\max})$ versus the number of modes $K$ for various values of $\nu$
On the right: the minimal number of modes $K$ necessary for the validation of a periodic orbit of (52) depending on $\mu$. The stepsize in $\mu$ is not chosen uniformly, because the problem is more sensitive to changes in $\mu$ for $\mu$ bigger, therefore requesting a steeper change in the number of modes. One the left: several of the validated periodic orbits
The computation time versus the number of modes $K$. The reference line has a slope corresponding to $K^3$
The computation time versus the dimension $N$ of the system (53). The reference lines have slopes corresponding to quadratic and cubic dependence on $N$
The computation time versus the order $D$ of the system (54). The reference lines have slopes corresponding to quadratic and cubic dependence on $D$
Validation of a periodic solution for the Lorenz system (55) with the classical parameter values $\sigma = 10$, $\beta = 8/3$ and $\rho = 28$. The solution has a period close to 25. This solution was validated with $\nu = 1+10^{-6}$ and $K = 800$
Step size versus $\mu$. Some oscillation is noticeable in the step size. This is due to the fact that the stepsize is decreased when the number of modes is changed, see Section 9. During this validation, the number of modes has been increased automatically from 4 to about 400
Validated continuation for the Lorenz system (55) from $\rho = 28$ to $\rho = 14.8$. In the left graph, the lower orbit still has two distinct swirls in the left lobe, but they are too close to each other to be seen in this graph
Validated continuation for $10$ coupled Lorenz systems (77) from $\epsilon = 0.14205$ to $\epsilon = 0.14252$. In the left graph we have depicted in blue the orbit of the first three component, in green the orbit of the last three component, with the biggest validated $\epsilon$, in the right graph we depicted the values of $\epsilon$, where we can notice the adaptation of the stepsize
Continuation through the fold of the Rychkov system (78). The validation started at the top, where a small initial stepsize was chosen. One can see how the stepsize was automatically adjusted along the validated curve. As can be seen from the values along the horizontal axis, we only depict a validated continuation for parameter values very close to saddle-node
 [1] Rongchang Liu, Jiangyuan Li, Duokui Yan. New periodic orbits in the planar equal-mass three-body problem. Discrete & Continuous Dynamical Systems - A, 2018, 38 (4) : 2187-2206. doi: 10.3934/dcds.2018090 [2] Zaihong Wang, Jin Li, Tiantian Ma. An erratum note on the paper: Positive periodic solution for Brillouin electron beam focusing system. Discrete & Continuous Dynamical Systems - B, 2013, 18 (7) : 1995-1997. doi: 10.3934/dcdsb.2013.18.1995 [3] Tobias Geiger, Daniel Wachsmuth, Gerd Wachsmuth. Optimal control of ODEs with state suprema. Mathematical Control & Related Fields, 2021  doi: 10.3934/mcrf.2021012 [4] Francisco Braun, Jaume Llibre, Ana Cristina Mereu. Isochronicity for trivial quintic and septic planar polynomial Hamiltonian systems. Discrete & Continuous Dynamical Systems - A, 2016, 36 (10) : 5245-5255. doi: 10.3934/dcds.2016029 [5] Jérôme Ducoat, Frédérique Oggier. On skew polynomial codes and lattices from quotients of cyclic division algebras. Advances in Mathematics of Communications, 2016, 10 (1) : 79-94. doi: 10.3934/amc.2016.10.79 [6] V. V. Zhikov, S. E. Pastukhova. Korn inequalities on thin periodic structures. Networks & Heterogeneous Media, 2009, 4 (1) : 153-175. doi: 10.3934/nhm.2009.4.153 [7] Giovanni Cimatti. Forced periodic solutions for piezoelectric crystals. Communications on Pure & Applied Analysis, 2005, 4 (2) : 475-485. doi: 10.3934/cpaa.2005.4.475 [8] Qigang Yuan, Jingli Ren. Periodic forcing on degenerate hopf bifurcation. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2857-2877. doi: 10.3934/dcdsb.2020208 [9] Shanjian Tang, Fu Zhang. Path-dependent optimal stochastic control and viscosity solution of associated Bellman equations. Discrete & Continuous Dynamical Systems - A, 2015, 35 (11) : 5521-5553. doi: 10.3934/dcds.2015.35.5521 [10] Changpin Li, Zhiqiang Li. Asymptotic behaviors of solution to partial differential equation with Caputo–Hadamard derivative and fractional Laplacian: Hyperbolic case. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021023 [11] Rafael Luís, Sandra Mendonça. A note on global stability in the periodic logistic map. Discrete & Continuous Dynamical Systems - B, 2020, 25 (11) : 4211-4220. doi: 10.3934/dcdsb.2020094 [12] Jaume Llibre, Luci Any Roberto. On the periodic solutions of a class of Duffing differential equations. Discrete & Continuous Dynamical Systems - A, 2013, 33 (1) : 277-282. doi: 10.3934/dcds.2013.33.277 [13] Zhiming Guo, Zhi-Chun Yang, Xingfu Zou. Existence and uniqueness of positive solution to a non-local differential equation with homogeneous Dirichlet boundary condition---A non-monotone case. Communications on Pure & Applied Analysis, 2012, 11 (5) : 1825-1838. doi: 10.3934/cpaa.2012.11.1825 [14] Jianping Gao, Shangjiang Guo, Wenxian Shen. Persistence and time periodic positive solutions of doubly nonlocal Fisher-KPP equations in time periodic and space heterogeneous media. Discrete & Continuous Dynamical Systems - B, 2021, 26 (5) : 2645-2676. doi: 10.3934/dcdsb.2020199 [15] Cécile Carrère, Grégoire Nadin. Influence of mutations in phenotypically-structured populations in time periodic environment. Discrete & Continuous Dynamical Systems - B, 2020, 25 (9) : 3609-3630. doi: 10.3934/dcdsb.2020075 [16] Marcelo Messias. Periodic perturbation of quadratic systems with two infinite heteroclinic cycles. Discrete & Continuous Dynamical Systems - A, 2012, 32 (5) : 1881-1899. doi: 10.3934/dcds.2012.32.1881 [17] Wei-Jian Bo, Guo Lin, Shigui Ruan. Traveling wave solutions for time periodic reaction-diffusion systems. Discrete & Continuous Dynamical Systems - A, 2018, 38 (9) : 4329-4351. doi: 10.3934/dcds.2018189 [18] Luke Finlay, Vladimir Gaitsgory, Ivan Lebedev. Linear programming solutions of periodic optimization problems: approximation of the optimal control. Journal of Industrial & Management Optimization, 2007, 3 (2) : 399-413. doi: 10.3934/jimo.2007.3.399 [19] José Raúl Quintero, Juan Carlos Muñoz Grajales. On the existence and computation of periodic travelling waves for a 2D water wave model. Communications on Pure & Applied Analysis, 2018, 17 (2) : 557-578. doi: 10.3934/cpaa.2018030 [20] Xiaoming Wang. Quasi-periodic solutions for a class of second order differential equations with a nonlinear damping term. Discrete & Continuous Dynamical Systems - S, 2017, 10 (3) : 543-556. doi: 10.3934/dcdss.2017027

Impact Factor: