关于水星的轨道方程以及运行周期等问题的探讨.docVIP

关于水星的轨道方程以及运行周期等问题的探讨 摘要：运用万有引力定律和Newton第二定律，在复平面坐标系中建立了水星绕太阳运动轨迹的简单数学模型二阶微分方程组模型，然后通过解析法，求出了水星到太阳的最近距离为，再运用Kepler第二定律求出水星绕太阳运动的周期。又运用基于压缩映象的求根方法求出了以及从远日点开始的第50天（地球天）结束时水星的位置。接着又运用经典的四阶 Runge-Kutte法，通过MATLAB编程，计算出取不同的步长h时所对应的水星绕太阳运动的周期T以及从远日点开始的第50天（地球天）结束时水星的位置r，比较发现结果与前面方法所得的结果非常接近。最后利用MATLAB画出了水星绕太阳运动的轨迹曲线。 关键字：水星 轨道 周期 基于压缩映象的求根方法 Runge-Kutte法 一、模型建立： 设太阳中心所在位置为复平面之原点O,在时刻t,行星位于 （4.1） 所表示的点P.这里均为的函数,分别表示的模和辐角。 于是行星的速度为 其加速度为 (4.2） 而太阳对行星的引力依万有引力定律，大小为,方向由行星位置P指向太阳的中心O，故为,其中M太阳的质量，m为行星的质量，G为万有引力常数 依Newton第二定律，我们得到 （4.3） 将（4.2）带入（4.3），然后比较实部与虚部，就有 这是两个未知函数的二阶微分方程组。在确定某一行星轨道时，需要加上定解条件。我们设当t=0时，行星正处于远日点，而远日点位于正实轴上，距原点O为，行星的线速度为.那么就有如下的初值条件： 问题（4.4）~（4.9）就是行星绕太阳运行的轨迹的数学模型。 将式（4.4）乘以r，即得 从而 （4.10） 其中为常数 这样，有向线段在时间内扫过的面积等于 （4.11） 显然，这正是Kepler的第二定律：从太阳指向行星的线段在单位时间内扫过的面积相等。 将式（4.10）改写后代入式（4.5） 于是我们得到了行星运行的较为简单形式的数学模型： 二、模型求解： 1.行星的轨迹 前面求得了行星的轨迹方程，要消去变量t，令，那么式（4.13）可以改写为 （4.17） 从而 将上式代入式（4.12），化简后为 （4.18） 其中 式（4.18）是一个二阶常系数非齐次常微分方程，引进，立即可以求出 这里A和是待定的常数。记上式可写为 这就是行星的轨道方程，是一条平面二次曲线。由于行星绕太阳运行，故必有，这样我们得到了Kepler第一定律：行星的轨道是以太阳为一个焦点的椭圆上。 由于r在时取得最大值（远日点），这意味着此时函数取得最大值1。于是就有 从而行星轨迹的方程为 （4.19） 对水星而言，我们已知 由此可得 显然，在式（4.19）中时，水星到达近日点，故其到太阳的距离为 带入上面的数据，从而计算得到水星到太阳的最近距离为 2.运行周期 设行星的周期为T，那么利用Kepler第二定律，即式（4.11），我们有 （4.20） 上式左端为行星轨迹椭圆所围的面积，记为S 我们已经求得椭圆的半长轴，半短轴，从而有 将上式代入式(4.20),解得 (4.21) 对于水星，将有关数据代入，易得 3.行星的位置 Kepler第二定律，即式（4.11），而该式可改写为 （4.22） 从而可得 （4.23） 如果我们要求时相应的和，则意味着首先要解方程 （4.24） 其中 在求出了时的后，立即可以由式（4.19）得到相应的r. 但是在实际问题中，即使可以积分出初等表达式，要由方程（4.24）来求出的解析解的可能性也是很小的。当然，我们可以利用计算编程利用迭代的方法，逐步增大的值来求得最符合条件的解，但这显然不是一个较为合理的方法。我们可以利用数值方法来对该问题进行求解，基于压缩映像的求根方法是其中较为有效的一种方法。其算法如下： 我们已经求得椭圆的半长轴，半短轴，从