Constructing Linear Operators Using Classical Perturbation Theory

This work introduces a methodology for generating linear operators that approximately represent nonlinear systems of perturbed ordinary differential equations. This is done through the application of classical perturbation theory via the Lindstedt-Poincar\'e expansion, followed by an extension of the space of configuration that guarantees the linear representation of the expanded system of differential equations. To ensure that such a linear representation exists, this paper uses polynomial basis functions. Pseudo-code describing the implementation of the proposed method is listed. The method is applied to the Duffing oscillator as well as to the J2 problem, with and without atmospheric drag, both analyzed using an osculating formulation. Additionally, conditions on the osculating Keplerian elements that produce low-eccentricity frozen orbits are presented, and a modification of the Lindstedt-Poincar\'e method is proposed to enable the generation of linear operators that dynamically adapt to changes in the frequency of the motion. Finally, the proposed method is compared with alternatives in the literature.