Many systems in science and engineering can be modelled as coupled or forced nonlinear oscillators, which may possess quasi-periodic or phase-locked invariant tori. Since there exist routes to chaos involving the break-down of invariant tori, these phenomena attract considerable attention. This paper presents a new algorithm for the computation and continuation of quasi-periodic invariant tori of ordinary differential equations that is based on a natural parametrisation of such tori. Since this parametrisation is uniquely defined, the proposed method requires neither the computation of a base of a transversal bundle, nor re-meshing during continuation. It is independent of the stability type of the torus and examples of attracting and saddle-type tori are given. The algorithm is robust in the sense that it can compute approximations to weakly resonant tori. The performance of the method is demonstrated with examples.

pdf copy of the paper (2.8M)