matlab ode45三体问题,三 体 问 题 初 探

摘 要：三体问题(three-body

problem)是天体力学中的基本力学模型，分为一般三体问题、圆型限制性三体问题和椭圆型限制性三体问题.本文就圆型限制性三体问题进行研究并绘出轨迹图，分析运动的混沌现象.

对圆型限制性三体问题的研究方法也适用于一般三体问题.

关键词：圆型限制性三体问题 天体动力学 混沌

一、引言

三体问题是研究三个可视为质点的天体在相互之间的万有引力作用下的运动规律问题，分为一般三体问题、圆型限制性三体问题和椭圆型限制性三体问题.

在一般三体问题中，每一个天体在其他两个天体的万有引力作用下的运动方程都可以表示成3个二阶的常微分方程，或6个一阶的常微分方程。因此，一般三体问题的运动方程为十八阶方程，必须得到18个积分才能得到完全解。然而，目前还只能得到三体问题的10个初积分，还远不能解决三体问题。

研究三体问题的方法大致可分为三类：第一类是分析方法，其基本原理是把天体的坐标和速度展开为时间或其他小参数的级数形式的近似分析表达式，从而讨论天体的坐标或轨道要素随时间的变化；第二类是定性方法，采用微分方程的定性理论来研究长时间内三体运动的宏观规律和全局性质；第三类是数值方法，这是直接根据微分方程的计算方法得出天体在某些时刻的具体位置和速度。这三类方法各有利弊，对新积分的探索和各类方法的改进是研究三体问题中很重要的课题。三体问题中有两个特例：早在1772年，拉格朗日已发现三体问题中存在两个定型解(如

)，即等边三角形定型解和直线形定型解.本文只讨论三体问题中最简单的一种情况：圆型限制性三体问题，即研究一个很小的质量体在两个有限质量体(围绕公共质心作圆周运动)的万有引力作用下运动的规律问题.讨论这个问题是有实际意义