Symplectic integrators for the matrix Hill equation

详细信息    查看全文

• 作者：Philipp Bader^a ; ^{p.bader@latrobe.edu.au} ; Sergio Blanes^b ; ^{serblaza@imm.upv.es} ; Enrique Ponsoda^b ; ^{eponsoda@imm.upv.es} ; Muaz Seydaoğlu^c ; ^{muasey@upvnet.upv.es}
• 关键词：65L07 ; 65L05 ; 65Z05
• 刊名：Journal of Computational and Applied Mathematics
• 出版年：2017
• 出版时间：15 May 2017
• 年：2017
• 卷：316
• 期：Complete
• 页码：47-59
• 全文大小：616 K
• 卷排序：316

文摘

We consider the numerical integration of the matrix Hill equation. Parametric resonances can appear and this property is of great interest in many different physical applications. Usually, Hill’s equations originate from a Hamiltonian function and the fundamental matrix solution is a symplectic matrix. This is a very important property to be preserved by the numerical integrators. In this work we present new sixth-and eighth-order symplectic exponential integrators that are tailored to Hill’s equation. The methods are based on an efficient symplectic approximation to the exponential of high dimensional coupled autonomous harmonic oscillators and yield accurate results for oscillatory problems at a low computational cost. The proposed methods can also be used for solving general second order linear differential equations where their performance will depend on how the methods are finally adapted to each particular problem or the qualitative properties one is interested to preserve. Several numerical examples illustrate the performance of the new methods.