辛算法的纠飘研究

辛算法较RK(Runge-Kutta)方法,保持辛结构不变或保持哈密顿系统规律性不变是突出的优点,但点态数值精度并不理想.推导出了三阶、四阶辛算法的漂移量计算公式,通过补偿漂移量就能提高点态数值精度,既保辛结构又保证点态数值高精度,适合于工程应用.建立了分数步对称辛算法(即FSJS算法)的纠漂公式,制定了漂移的约束标准.相关算例的数值结果表明:三阶FSJS算法漂移量最小,点态数值精度更高.