几类辛方法的数值稳定性研究

• 作者: 伍歆 刘福窑 陆本魁
• 来源期刊: 天文学报
• 学科: 地球科学
• 关键词: 天体力学:辛算法 线性稳定性
• 年，卷（期）: 2006,（4）
• 页码范围: 418-431
• 页数: 14页
• 分类号: P135
• 字数: 9597字
• 语种: 中文
• DOI: 10.3321/j.issn:0001-5245.2006.04.009

摘要

主要对一阶隐式Euler辛方法M1、二阶隐式Euler中点辛方法M2、一阶显辛Euler方法M3和二阶leapfrog显辛积分器M4共4种辛方法及一些组合算法进行了通常意义下的线性稳定性分析.针对线性哈密顿系统,理论上找到每个数值方法的稳定区,然后用数值方法检验其正确性.对于哈密顿函数为实对称二次型的情况,为了理论推导便利,特推荐采用相似变换将二次型的矩阵对角化来研究辛方法的线性稳定性.当哈密顿分解为一个主要部分和一个小摄动次要部分且二者皆可积时,无论是线性系统还是非线性系统,这种主次分解与哈密顿具有动势能分解相比,明显扩大了辛方法的稳定步长范围.

作者信息

序号	姓名	单位	发文数	被引次数	H指数	G指数
1	伍歆	南昌大学理学院	11	42	4.0	6.0