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] Ying Lv, Yan-Fang Xue, Chun-Lei Tang. Ground state homoclinic orbits for a class of asymptotically periodic second-order Hamiltonian systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (3) : 1627-1652. doi: 10.3934/dcdsb.2020176 [2] Jiangtao Yang. Permanence, extinction and periodic solution of a stochastic single-species model with Lévy noises. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020371 [3] Peizhao Yu, Guoshan Zhang, Yi Zhang. Decoupling of cubic polynomial matrix systems. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 13-26. doi: 10.3934/naco.2020012 [4] Xin Guo, Lexin Li, Qiang Wu. Modeling interactive components by coordinate kernel polynomial models. Mathematical Foundations of Computing, 2020, 3 (4) : 263-277. doi: 10.3934/mfc.2020010 [5] Guojie Zheng, Dihong Xu, Taige Wang. A unique continuation property for a class of parabolic differential inequalities in a bounded domain. Communications on Pure & Applied Analysis, , () : -. doi: 10.3934/cpaa.2020280 [6] Jing Zhou, Cheng Lu, Ye Tian, Xiaoying Tang. A SOCP relaxation based branch-and-bound method for generalized trust-region subproblem. Journal of Industrial & Management Optimization, 2021, 17 (1) : 151-168. doi: 10.3934/jimo.2019104 [7] Mengyu Cheng, Zhenxin Liu. Periodic, almost periodic and almost automorphic solutions for SPDEs with monotone coefficients. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021026 [8] Julian Tugaut. Captivity of the solution to the granular media equation. Kinetic & Related Models, , () : -. doi: 10.3934/krm.2021002 [9] Shuxing Chen, Jianzhong Min, Yongqian Zhang. Weak shock solution in supersonic flow past a wedge. Discrete & Continuous Dynamical Systems - A, 2009, 23 (1&2) : 115-132. doi: 10.3934/dcds.2009.23.115 [10] Yukihiko Nakata. Existence of a period two solution of a delay differential equation. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1103-1110. doi: 10.3934/dcdss.2020392 [11] Yuanfen Xiao. Mean Li-Yorke chaotic set along polynomial sequence with full Hausdorff dimension for $\beta$-transformation. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 525-536. doi: 10.3934/dcds.2020267 [12] Feimin Zhong, Jinxing Xie, Yuwei Shen. Bargaining in a multi-echelon supply chain with power structure: KS solution vs. Nash solution. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2020172 [13] Yi An, Bo Li, Lei Wang, Chao Zhang, Xiaoli Zhou. Calibration of a 3D laser rangefinder and a camera based on optimization solution. Journal of Industrial & Management Optimization, 2021, 17 (1) : 427-445. doi: 10.3934/jimo.2019119 [14] Rong Chen, Shihang Pan, Baoshuai Zhang. Global conservative solutions for a modified periodic coupled Camassa-Holm system. Electronic Research Archive, 2021, 29 (1) : 1691-1708. doi: 10.3934/era.2020087 [15] Dong-Ho Tsai, Chia-Hsing Nien. On space-time periodic solutions of the one-dimensional heat equation. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3997-4017. doi: 10.3934/dcds.2020037 [16] Vaibhav Mehandiratta, Mani Mehra, Günter Leugering. Fractional optimal control problems on a star graph: Optimality system and numerical solution. Mathematical Control & Related Fields, 2021, 11 (1) : 189-209. doi: 10.3934/mcrf.2020033 [17] Tinghua Hu, Yang Yang, Zhengchun Zhou. Golay complementary sets with large zero odd-periodic correlation zones. Advances in Mathematics of Communications, 2021, 15 (1) : 23-33. doi: 10.3934/amc.2020040 [18] Yicheng Liu, Yipeng Chen, Jun Wu, Xiao Wang. Periodic consensus in network systems with general distributed processing delays. Networks & Heterogeneous Media, 2020  doi: 10.3934/nhm.2021002 [19] Zhihua Liu, Yayun Wu, Xiangming Zhang. Existence of periodic wave trains for an age-structured model with diffusion. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021009 [20] Kai Zhang, Xiaoqi Yang, Song Wang. Solution method for discrete double obstacle problems based on a power penalty approach. Journal of Industrial & Management Optimization, 2020  doi: 10.3934/jimo.2021018

Impact Factor: