Algebraic-numerical method for the slightly perturbed harmonic oscillator

Nadeau, A.; Guyard, J.; Feix, M. R.

doi:10.1090/S0025-5718-1974-0349020-9

Algebraic-numerical method for the slightly perturbed harmonic oscillator
HTML articles powered by AMS MathViewer

by A. Nadeau, J. Guyard and M. R. Feix PDF

Math. Comp. 28 (1974), 1057-1066 Request permission

Abstract:

The solution of slightly perturbed harmonic oscillators can easily be obtained in the form of a series given by Poisson’s method. However, this perturbation method leads to secular terms unbounded for large time (the time unit being the fundamental period of the harmonic oscillator), which prevent the use of finite series. The analytical elimination of such terms was first solved by Poincaré and, more recently, generalized by Krylov and Bogoliubov. Unfortunately, these methods are very difficult to handle and are not easily carried out for high orders. A numerical reinitialization method is combined here with the Poisson perturbation treatment to avoid the growth of secular terms and therefore to get the solution at any time. The advantages of such a method is that the analytical work can be carried to high orders keeping the step of numerical integration to a relatively large value (compared to a purely numerical method). This algorithm has been tested on the Mathieu equation. A method for the computation of the eigenvalues of this equation is given. By properly selecting the order of the perturbation and the time step of reinitialization, we can recover, at any order, all the effects of the slight perturbation (including all the unstable zones). Consequently, such a method is a useful intermediate between purely analytical and purely numerical algorithms.

References

N. N. Bogoliubov and Y. A. Mitropolsky, Asymptotic methods in the theory of non-linear oscillations, Translated from the second revised Russian edition, International Monographs on Advanced Mathematics and Physics, Hindustan Publishing Corp., Delhi, Gordon and Breach Science Publishers, Inc., New York, 1961. MR 0141845
The plasma in a magnetic field: A symposium on magnetohydrodynamics. , Stanford University Press, Stanford, Calif., 1958. MR 0115576

J. Mathematical Phys.

Jacques Guyard, André Nadeau, Germain Baumann, and Marc R. Feix, Eigenvalues of the Hill equation to any order in the adiabatic limit, J. Mathematical Phys. 12 (1971), 488–492. MR 280104, DOI 10.1063/1.1665611