Computation of Limit Cycles and Their Isochrons: Fast Algorithms and Their Convergence
Computation of Limit Cycles and Their Isochrons: Fast Algorithms and Their Convergence
复制标题
极限环及其等时线的计算:快速算法及其收敛性
DOI:
10.1137/120901210
复制
发表时间:
2013
期刊:
影响因子:
--
通讯作者:
R. Llave
中科院分区:
文献类型:
--
作者:
G. Huguet;R. Llave
We present efficient algorithms to compute limit cycles and their isochrons (i.e., the sets of points with the same asymptotic phase) for planar vector fields. We formulate a functional equation for the parameterization of the invariant cycle and its isochrons, and we show that it can be solved by means of a Newton method. Using the right transformations, we can solve the equation of the Newton step efficiently. The algorithms are efficient in the sense that if we discretize the functions using $N$ points, a Newton step requires $O(N)$ storage and $O(N\log N)$ operations in Fourier discretization or $O(N)$ operations in other discretizations. We prove convergence of the algorithms and present a validation theorem in an a posteriori format. That is, we show that if there is an approximate solution of the invariance equation that satisfies some some mild nondegeneracy conditions, then there is a true solution nearby. Thus, our main theorem can be used to validate numerically computed solutions. The theorem ...