论文部分内容阅读
与Runge-Kutta（RK）方法相比，辛算法具有保持相空间辛结构不变或保哈密顿函数不变的突出优点。但是，在时域上，同阶的辛算法与Runge-Kutta法具有相同的数值精度，即辛算法在计算过程中也存在相位误差，导致解的数值精度不高。为了提高辛算法在时域上解的精度，首先根据哈密顿函数的特点将哈密顿系统归结为2种类型，然后建立了不同类型下辛算法的相位误差公式，归纳出各类相位漂移的特点，进而提出了一种适当的纠漂方法，使得辛算法在时域上获得了很高的数值精度。相关算例的数值结果验证了纠漂理论的有效性和可靠性。