保持Runge-Lenz向量的数值方法   总被引：2，自引：2，他引：0  

刘福窑  伍歆  陆本魁 《天文学报》2005,46(3):294-306

对孤立积分和能够保持Runge-Lenz向量的梯形公式进行详尽讨论．孤立积分就是限制粒子运动区域的不变量，具有n个自由度的自治可积哈密顿系统且只有n个互相对合的独立孤立积分，并且其他孤立积分的存在对粒子的运动是有意义的，Kepler二体系统存在能量积分、角动量积分和Runge-Lenz向量．对于平面运动情况，这三类积分中只有3个独立孤立积分；而对于三维空间情形，该三类积分仅有5个是独立的．就前者而言，Kepler二体平面运动积分构成该系统中的对称群SO（3），经过Levi-Civita变换，它可以转化为二维各向同性谐振子系统中的对称群，而该对称群能够被梯形公式准确保持，另一方面，对于后者梯形公式对这三类积分的严格保持还可以在5个Kepler轨道根数n、e、i、Ω和w上得到体现。  相似文献   

有摄圆型限制性三体问题平动点稳定的充要条件及应用     

舒斯会  陆本魁  陈务深  刘福窑 《天文学报》2005,46(1):37-45

摘要给出了一个判断有摄圆型限制性三体问题平动点稳定性的充要条件．该条件只依赖于平动点变分方程的特征方程系数的一个简单关系，使用很方便．用所得到的条件，讨论了任意外力摄动对经典圆型限制性三体问题三角平动点稳定性的影响和惯性阻力摄动对Robe圆型限制性三体问题主要平动点的稳定性的影响．  相似文献   

几类辛方法的数值稳定性研究   总被引：1，自引：0，他引：1  

刘福窑  伍歆  陆本魁 《天文学报》2006,47(4):418-431

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

一个膺三阶辛积分器   总被引：1，自引：1，他引：0  

刘福窑  伍歆  陆本魁 《天文学报》2004,45(4):402-412

在太阳系动力学中,辛积分器已成为研究哈密顿系统的长期定性演化的最佳工具.对于可积分离的哈密顿系统H=H0+∑i=1N∈iHi(∈≤1),构造了一个膺三阶辛积分器.它大约相当于Wisdom-Holman二阶辛积分器的一次校正或Forest-Ruth四阶辛算法的精度.此外,含力梯度的辛算法也适合处理哈密顿系统H=Ho(q,P)+∈H1(q),其精度好于原辛积分器,但不优越于相应膺高阶辛积分器.  相似文献   

平动点线性稳定条件及阻力对限制性三体问题三角平动点线性稳定性的影响     

舒斯会  陆本魁  刘福窑  陈务深 《天文学报》2004,45(3):301-309

给出了受摄限制性三体问题平动点线性稳定性的一些判断条件，条件只与相应的平动点切映像的特征方程系数有关，使用方便，用这些判断条件，讨论了一些阻力对经典限制性三体问题三角平动点线性稳定性的影响，改进了Murray等的一些结果。  相似文献