Constructing Linear Operators Using Classical Perturbation Theory

TL;DR

This paper presents a novel method using classical perturbation theory to construct linear operators that approximate nonlinear perturbed differential systems, demonstrated on orbital models and with adaptive frequency features.

Contribution

It introduces a new approach combining Lindstedt-Poincaré expansion and polynomial basis functions to generate linear representations of nonlinear systems, including adaptive frequency adjustments.

Findings

Successfully applied to Duffing oscillator and J2 problem

Achieved low-eccentricity frozen orbit conditions

Compared favorably with existing methods in literature

Abstract

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.