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

关于水星的轨道方程以及运行周期等问题的探讨 理学院07级计算2班**数学实验论文 关于水星的轨道方程以及运行周期等问题的探讨 摘要:运用万有引力定律和Newton第二定律，在复平面坐标系中建立了水星绕太阳运动轨迹的简单数学模型二阶微分方程组模型，然后通过解析法，求出了水星到太阳的最近距离 10r,4.6016，10(m)为，再运用Kepler第二定律求出水星绕太阳运动的周期m 6T,7.6025，10(s),87.9919(d)。又运用基于压缩映象的求根方法求出了以及从远日点 10rm,，4.766810()开始的第50天(地球天)结束时水星的位置。接着又运用经典的四阶 Runge-Kutte法，通过MATLAB编程，计算出取不同的步长h时所对应的水星绕太阳运动的周期T以及从远日点开始的第50天(地球天)结束时水星的位置r，比较发现结果与前面方法所得的结果非常接近。最后利用MATLAB画出了水星绕太阳运动的轨迹曲线。 关键字:水星 轨道 周期 基于压缩映象的求根方法 Runge-Kutte法 一、模型建立: 设太阳中心所在位置为复平面之原点O,在时刻t,行星位于 i,,(t),re (4.1) 所表示的点P.这里均为的函数,分别表示的模和辐角。 r,r(t),,,,(t)tZ(t) 于是行星的速度为 dZdrddrd,,iii,,, ,，,，()eireeirdtdtdtdtdt 其加速度为 222dZdrdddrd,,,2i,erir,[(,())，(，2)] (4.2) 222dtdtdtdtdtdt mMG 而太阳对行星的引力依万有引力定律，大小为,方向由行星位置P指向太阳的中2r mMGi,心O，故为,其中M太阳的质量，m为行星的质量，G为万有引力常数 e,2r 30M1.98910kg,，() ,1122G6.67210,，,(N/kgm) 依Newton第二定律，我们得到 理学院07级计算2班**数学实验论文 2mMGdZ,i,, (4.3) em22rdt 将(4.2)带入(4.3)，然后比较实部与虚部，就有 2,ddrd,,(4.4)r20 ，,,2,dtdtdt ,2drdMG,2,(4.5)r(),,,22,dtdtr, 这是两个未知函数的二阶微分方程组。在确定某一行星轨道时，需要加上定解条件。我们设 rv当t=0时，行星正处于远日点，而远日点位于正实轴上，距原点O为，行星的线速度为.00 那么就有如下的初值条件: ,rr, (4.6) t,00,,,0 (4.7),t,0 ,,dr,,0 (4.8) t,0dt, ,vd,0 (4.9) ,,t,0rdt,0, 问题(4.4)~(4.9)就是行星绕太阳运行的轨迹的数学模型。 将式(4.4)乘以r，即得 dd,,,2r,0 ,,dtdt,, 从而 d,2 (4.10) r,C1dt C其中为常数 1 C,rv 100 OP这样，有向线段在时间内扫过的面积等于 ,t t，,t1C,t,d21 (4.11) rdt,,t22dt 显然，这正是Kepler的第二定律:从太阳指向行星的线段在单位时间内扫过的面积相等。 将式(4.10)改写后代入式(4.5) 22CdrMG1,,, 232dtrr于是我们得到了行星运行的较为简单形式的数学模型: 理学院07级计算2班**数学实验论文 22,drMGC1,,, (4.12),232dtrr,Cd,,1, (4.13)2,dtr, rr, (4.14),00t, ,dr,,0 (4.15).t,0dt, ,,,0 (4.16)t,0, , 二、模型求解: 1.行星的轨迹 1前面求得了行星的轨迹方程，要消去变量t，令，那么式(4.13)可以改写为 r,u d,2 (4.17) ,Cu1dt 从而 ,drdududdu11C ,,,,,,122dtdtd,dtd,uu 22drddudu,,22,,C,,Cu ,,1122dtd,dtd,,, 将上式代入式(4.12)，化简后为 2du1，,u (4.18) 2pd, 2C1,p其中 MG 1u,u,式(4.18)是一个二阶常系数非齐次常微分方程，引进，立即可以求出 p 1u,,u,Acos(,),, 0p ,这里A和是待定的常数。记上式可写为 e,Ap,0 pr, 1,ecos(,,,)0这就是行星的轨道方程，是一条平面二次曲线。由于行星绕太阳运行，故必有，0,e,1 这样我们得到了Kepler第一定律:行星的轨道是以太阳为一个焦点的椭圆上。 理学院07级计算2班**数学实验论文 rcos(,,,) 由于r在时取得最大值(远日点)，这意味着此时函数取得最大值1。t,000 于是就有 p ,,0,e,1,0r0从而行星轨迹的