Pleiades问题的Runge-Kutta法求解与其平面轨迹.pdfVIP

维普资讯 第29卷 第6期 宜春学院学报 (自然科学) Vo1．29．No．6 2007年 12月 JournalofYichunUniversity(naturalflcience) Dee．20o7 Pleiades问题的Runge—Kutta法求解与其平面轨迹 陈全发 (湘南学院数学系，湖南 郴州 423000) 摘 要：该文针对天体力学中的Pleiades问题，提 出其数学模型和解决问题 的算法，选择经典四阶显式 Runge—Kutta法来获得在一定时间内的7个天体各 自的坐标近似值，利用Matlab语言编写相关程序，在平面直角 坐标系x—Y中绘制其运行轨道． 关键词：Pleiades；Runge—Kutta法 ；坐标 ；轨迹 中图分类号：015 文献标识码：A 文章编号：1671—380X (2007)06-0059—04 Runge—KuttaM ethodAndNumericalSimulationForPleiadesProblem CHENQuart—fa (MathematicsDepartmentofXiangnanUniversity，Chenzhouttunan423000China) Abstract：Accordingtotheproblem ofcelestialmechanicsofPleiades，loadingthemathematicalmodelandthealgorithm ofsolution， choosingtheRunge—Kuttamethodtogettheapproximateofcoordinatesfor7bodies，simulatingthebody’Strajectorieswithprogram— melP8ofMatlabintheplanespace． Keywords：Pleiades；RK—method；coordinates；trajectory 0 引言 X0= (3，3，一1，一3，2，一2，2) ． Pleiades问题是一个天体力学问题，它描述的是7个天 这里 3 0 体 (坐标分别为 ( ，Y)， =1，2，…，7，质量为 ) 在平面直角坐标系x—Y中的运动轨迹．其实，Pleiades问 L)，= (0，0，0，一1．25，1，0，0) 题可以看作是由14个特殊的二阶非刚性微分方程构成．若 根据力学知识可以获得该问题 的表达式 ，如：考虑物 将其写成一阶形式，就形成一个含有28个一阶微分方程的 体 i，由牛顿第二定律 ，物体 i的运动方程为： 微分方程组．该问题 的数学描述是： (见参考文献 [1] = Mipi”(t)，P= ( ，Y)． (4) [2])函数形式 另一方面，根据万有引力定律，物体 i受到物体j的引 。“(t)= (。)， 。(0) Z0， 。(0) =。 (1) 力 的作用， 其中，。ER“，0≤￡≤3．