Fictitious time integration method for seeking periodic orbits of nonlinear dynamical systems

摘要

In the paper, we compute the period and periodic orbit of an n-dimensional nonlinear dynamical system by developing two iterative algorithms based on the fictitious time integration method (FTIM). Periodicity condition, also known as the Poincare map, is a necessary condition for the existence of a periodic motion in the state space, which consists of n implicit nonlinear algebraic equations (NAEs). Instead of solving the NAEs by the Newton-Raphson method, we derive an equivalent nonlinear scalar equation to determine unknown period by using the FTIM. The resulting sequence of the iterated periods monotonically converges to a desired period, wherein the unknown initial point on the periodic orbit is determined simultaneously. We find that the second iterative algorithm using the matrix shape function to convert the periodic problem into a corresponding initial value problem is convergent faster than the first iterative algorithm. A key point is that the periodicity condition is satisfied automatically by the second iterative algorithm when the terminal values of the new variables are convergent. Numerical examples exhibit some major advantages of these two iterative algorithms.