Mathematics – Numerical Analysis
Scientific paper
2006-12-21
Mathematics
Numerical Analysis
34 pages, 7 figures
Scientific paper
Two-step hybrid methods specially adapted to the numerical integration of perturbed oscillators are obtained. The formulation of the methods is based on a refinement of classical Taylor expansions due to Scheifele [{\em Z. Angew. Math. Phys.}, {\bf 22}, 186--210 (1971)]. The key property is that those algorithms are able to integrate exactly harmonic oscillators with frequency $\omega$ and that, for perturbed oscillators, the local error contains the (small) perturbation parameter as a factor. The methods depend on a parameter $\nu=\omega h$, where $h$ is the stepsize. Based on the B2-series theory of Coleman [{\em IMA J. Numer. Anal.}, {\bf 23}, 197--220 (2003)] we derive the order conditions of this new type of methods. The linear stability and phase properties are examined. The theory is illustrated with some fourth- and fifth-order explicit schemes. Numerical results carried out on an assortment of test problems (such as the integration of the orbital motion of earth satellites) show the relevance of the theory.
No associations
LandOfFree
Two-step hybrid methods adapted to the numerical integration of perturbed oscillators does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Two-step hybrid methods adapted to the numerical integration of perturbed oscillators, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Two-step hybrid methods adapted to the numerical integration of perturbed oscillators will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-22546