李任飞
- 作品数：1 被引量：1H指数：1
- 供职机构：中山大学工学院更多>>
- 发文基金：中央高校基本科研业务费专项资金国家自然科学基金更多>>
- 相关领域：理学更多>>
- 一种基于加权残值法的高阶辛算法被引量：1
- 2015年
- 提出了利用加权残值法构造高阶辛算法的一种新途径。首先根据加权残值法的思想,在时间子域内给出了哈密顿正则方程伽辽金法所对应的积分方程,然后在该时间子域内采用相同的拉氏插值作为位移和动量的试函数,并将这些试函数代入到积分方程中,通过数值积分,将原动力学初值问题转为以插值点位移和动量为未知量的代数方程组。对于非线性问题,给出了一种能显著提高牛顿迭代法计算效率的初值选取方案。最后,对算法的保辛性和性能进行了详细的讨论。通过与同阶辛RK法相比较,两种方法精度几乎完全相同,但文中方法更简便,计算量更小。数值算例结果表明该法在计算精度和效率上均具有良好的性能。
- 陆克浪富明慧李纬华李任飞
- 关键词：哈密顿系统加权残值法非线性动力学伽辽金法辛算法
