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

关于水星的轨道方程以及运行周期等问题的探讨 摘要: 运用万有引力定律和依〃2第二定律，在复平面坐标系中建立了水星绕太阳运动 轨迹的简单数学模型二阶微分方程组模型，然后通过解析法，求出了水星到太阳的最近距离 为匚-4.6016x101°伽），再运用K印0 第二定律求出水星绕太阳运动的周期 T - 7.6025X106（5）- 87.9919（d）。又运用基于压缩映象的求根方法求出了以及从远日点 开始的第50天（地球天）结束时水星的位置r = 4.7668xl0,0（m）o接看又运用经典的四阶 Runge-Kutte法，通过MATLAB编程，计算出取不同的步长h时所对应的水星绕太阳运动的周 期T以及从远tl点开始的第50天（地球天）结束吋水星的位置厂，比较发现结果与前面方 法所得的结果非常接近。最后利用MATLAB I田i出了水星绕太阳运动的轨迹曲线。 关键字：水星轨道周期基于压缩映象的求根方法Runge-Kutte法 一、模型建立： 设太阳屮心所在位置为复平面之原点O,在时刻/,行星位于 (4. 1)Z(r) = re (4. 1) 所表示的点P.这里r = r（r）,0 = ⑵均为/的函数，分别表示Z⑴的模和辐角。 于是行星的速度为其加速度为dZdtdtdedt 于是行星的速度为 其加速度为 dZ dt dt de dt ei0 斗理) dt dt d2Z d2Z ~d^ W[(兽“(少)+ 2兽+ 2字字)] dv dt dv dt dt (4.2) 而太阳对行星的引力依万有引力定律，大小为型笋，方向由行星位置 而太阳对行星的引力依万有引力定律，大小为 型笋，方向由行星位置P指向A阳的中 r 心O,故为一吵字 严,其屮M太阳的质量，加为行星的质量，G为万有引力常数 M = 1.989xl030(kg) G=6.672xl0~,,(N-m2/kg2) 依Newton笫二定律，我们得到 mMGd2Zdt2 mMG d2Z dt2 (4.3) 将（4.2）带入（4.3）,然后比较实部与虚部，就有d2e.drd0 nr— +2 = 0dt~ dt dtd 将（4.2）带入（4.3）,然后比较实部与虚部，就有 d2e.drd0 n r— +2 = 0 dt~ dt dt d2r W、2 MG (4.4) (4.5) 这是两个未知函数的二阶微分方程组。在确定某一行星轨道时，需要加上定解条件。我们设 当r=0时，行星正处于远日点，而远日点位于正实轴上，距原点。为/行星的线速度为心. 那么就有如下的初值条件: e ；=（）= ° dr ~dt de ~dt /=()_ — (4.6) (4.7) (4.8) (4.9) 问题（4.4） ~ （4.9）就是行星绕太阳运行的轨迹的数学模型。 将式（4.4）乘以门即得 dt\ dt (4. 10) (4. 10) 其中G为常数 这样，有向线段OP在时间△/内扫过的而积等于 C-r2-dt = C-r2-dt = 2 dt C,Ar 2 (4. 11) 显然，这正是Kepler的第二定律：从太阳指向行星的线段在单位时间内扫过的面积相等。 将式（4 10）改写后代入式（4.5） d2r C.2 _ MG 苛_ 卞一 一~ 于是我们得到了行星运行的较为简单形式的数学模型: 二、模型求解: 1?行星的轨迹 d~r C,2 _ MG — — 1dt2 r3 r2 dt r2 |m)=0 (4」2) (4.13) (4.14) . (4.16) 前面求得了行星的轨迹方程，要消去变量/,令r = -,那么式（4.13）可以改写为 U 从而 dt (4. 17) dr _ 1 du _ 1 du dO _ du dt ir dt ir d dt d d2r _ d ( du乔__5広(胁 将上式代入式（4.12）,化简后为 d2u (4. 18) 其中 MG - 1 式（4.18）是一个二阶常系数非齐次常微分方程，引进u = u-一，立即可以求出 P u = u = A cos（ — %） P 这里A和（）是待定的常数。记e = Ap,上式可写为 P 1 _￡cos(e_e(j 这就是行星的轨道方程，是一条平而二次曲线。由于行星绕太阳运行，故必有0 vevl, 这样我们得到了 K印/纣第一定律：行星的轨道是以太阳为一个焦点的椭圆上。 由于，?在f = 0时取得最大值心（远日点），这意味着此时函数cos（-（））取得最大值lo 于是就有 从而行星轨迹的方程为对水星而言，我们己知由此可得P1 从而行星轨迹的方程为 对水星而言，我们己知 由此可得 P 1 一 WCOS (4. 19) 斤)=0.6982 xl0(〃2) v() =3.886xl04(m/5) C[ = ^v0 - 2.7132xl015(m2