5D2 Our Moon: Modelling Orbital Motion with Differential Equations | 5D2 我们的月球：微分方程轨道建模

📚 5D2 Our Moon: Modelling Orbital Motion with Differential Equations | 5D2 我们的月球：微分方程轨道建模

In Further Mathematics, the celestial dance of our Moon around the Earth offers a classic application of second‑order differential equations and central force mechanics. The problem ‘5D2 Our Moon’ invites students to move beyond simple circular approximations and derive the actual elliptical orbit, connecting Newton’s law of gravitation, conservation of angular momentum and the polar form of an orbit equation. By modelling the Moon as a particle in a central force field, we can predict its shape, speed, period and energy with satisfying rigour.

在进阶数学中，月球绕地球运行的“天体之舞”为二阶微分方程和中心力力学提供了经典的应用场景。“5D2 我们的月球”这道题目要求学生突破简单的圆形近似，推导出真实的椭圆轨道，并将牛顿万有引力定律、角动量守恒与极坐标轨道方程联系在一起。通过把月球视为中心力场中的质点，我们能够以严谨的方式预测其轨道形状、速度、周期和能量，获得令人满足的物理图景。

1. The Central Force Problem | 中心力问题

Any motion under a central force, where the force is always directed towards a fixed point (the centre), can be analysed using polar coordinates (r, θ). The gravitational pull of the Earth on the Moon is a central force, directed towards the Earth’s centre. Because the force has zero transverse component, the angular momentum of the Moon about the Earth is conserved—a fact that dramatically simplifies the equations of motion.

所有中心力运动——即力始终指向一固定点（中心）的运动——都可以用极坐标 (r, θ) 进行分析。地球对月球的引力是中心力，总是指向地心。由于该力的横向分量为零，月球相对于地球的角动量守恒，这一事实极大地简化了运动方程。

2. Radial and Transverse Components | 径向与横向分量

In polar coordinates, the acceleration of a particle can be resolved into a radial component r̈ − rθ̇² and a transverse component (1/r) d/dt (r²θ̇). For a central force of magnitude F(r) directed towards the origin, Newton’s second law gives us two differential equations. The transverse equation immediately yields r²θ̇ = constant, which we denote by h, the specific angular momentum.

在极坐标中，质点的加速度可分解为径向分量 r̈ − rθ̇² 和横向分量 (1/r) d/dt (r²θ̇)。对于指向原点、大小为 F(r) 的中心力，牛顿第二定律给出两个微分方程。横向方程直接给出 r²θ̇ = 常数，记作 h，即单位质量的角动量。

3. Newton’s Law of Gravitation | 牛顿万有引力定律

The gravitational force between the Earth (mass M) and the Moon (mass m) is attractive and obeys an inverse‑square law. In vector form its magnitude is GMm/r², directed towards the Earth. For the Moon’s relative motion, the effective central acceleration is −GM/r², where the minus sign indicates attraction towards the centre. Using μ = GM simplifies the radial equation.

地球（质量 M）与月球（质量 m）之间的万有引力呈平方反比规律，大小为 GMm/r²，方向指向地心。在月球的相对运动中，有效的中心加速度为 −GM/r²，负号表示指向中心的吸引。令 μ = GM 可以简化径向方程。

4. Setting Up the Differential Equations | 建立微分方程

The radial equation for a body of unit mass moving under gravity is r̈ − rθ̇² = −μ/r². Together with the constant angular momentum condition r²θ̇ = h, we can eliminate θ̇ to obtain a single differential equation for r(t). However, it is far more revealing to transform this into an equation that relates the orbit shape r(θ).

单位质量天体在引力作用下的径向方程为 r̈ − rθ̇² = −μ/r²。结合角动量守恒条件 r²θ̇ = h，消去 θ̇ 后可以得到关于 r(t) 的单一微分方程。然而，将其转换成描述轨道形状 r(θ) 的方程会更加清晰地揭示轨道的性质。

r̈ − h²/r³ = −μ/r²

5. Conservation of Angular Momentum | 角动量守恒

The constancy of h = r²θ̇ has both geometric and physical meaning: it represents twice the areal velocity swept out by the radius vector. This is Kepler’s second law. Because h is constant, we can use it to change the independent variable from time t to angle θ via the chain rule, arriving at an elegant differential equation for u = 1/r.

h = r²θ̇ 的恒定兼具几何和物理意义：它代表矢径扫过面积速率的两倍。这正是开普勒第二定律。由于 h 为常数，可利用链式法则将自变量从时间 t 换成角度 θ，从而得到一个关于 u = 1/r 的简洁微分方程。

d/dt = (h/r²) d/dθ = h u² d/dθ

6. Transforming to an Orbit Equation | 转换为轨道方程

By substituting r = 1/u and differentiating with respect to θ, the radial equation simplifies dramatically. After applying the chain rule twice, the acceleration term r̈ transforms into −h²u² d²u/dθ². The resulting equation is linear in u, making it solvable by standard methods for second‑order linear ODEs.

代入 r = 1/u 并对 θ 求导，径向方程得到极大的简化。两次应用链式法则后，加速度项 r̈ 转化为 −h²u² d²u/dθ²。得到的方程关于 u 是线性的，因此可以用标准的二阶线性常微分方程方法求解。

d²u/dθ² + u = μ/h²

7. Solving for the Orbital Shape | 求解轨道形状

The homogeneous equation d²u/dθ² + u = 0 has solutions of the form A cos θ + B sin θ. A particular solution is simply the constant μ/h². Thus the general solution is u(θ) = μ/h² + C cos(θ − θ₀), which can be rewritten as u = (μ/h²)(1 + e cos(θ − θ₀)). The constant e is the eccentricity of the orbit, and θ₀ defines the orientation of the major axis.

齐次方程 d²u/dθ² + u = 0 的解具有 A cos θ + B sin θ 的形式。一个特解就是常数 μ/h²。因此通解为 u(θ) = μ/h² + C cos(θ − θ₀)，可改写为 u = (μ/h²)(1 + e cos(θ − θ₀))。其中常数 e 是轨道的偏心率，θ₀ 决定了长轴的取向。

r = (h²/μ) / (1 + e cos(θ − θ₀))

8. Eccentricity and Energy | 偏心率与能量

The eccentricity e is determined by the initial conditions—essentially the total mechanical energy of the Moon. For 0 ≤ e < 1, the orbit is an ellipse; e = 0 gives a circle, and e = 1 a parabola. The Moon's orbital eccentricity is about 0.0549, indicating an almost circular ellipse. The semi‑major axis a and semi‑latus rectum L are related to h and μ.

偏心率 e 由初始条件——本质上是月球的总机械能——决定。当 0 ≤ e < 1 时，轨道为椭圆；e = 0 对应圆，e = 1 为抛物线。月球轨道的偏心率约 0.0549，表明它是一个近圆的椭圆。半长轴 a 和半通径 L 与 h 和 μ 存在直接关系。

e = √(1 + 2Eh²/μ²) , L = h²/μ = a(1 − e²)

9. Applying to Our Moon | 应用到我们的月球

Substituting measured values for the Earth–Moon distance (perigee ≈ 363 300 km, apogee ≈ 405 500 km) yields a semi‑major axis of about 384 400 km and the small eccentricity above. The angular momentum h for the Moon can be inferred from the orbital period and geometry, confirming that the central‑force model matches observation remarkably well, even before considering perturbations from the Sun.

代入地月距离的实测值（近地点 ≈ 363 300 km，远地点 ≈ 405 500 km），可算出半长轴约 384 400 km 以及上述的小偏心率。利用轨道周期和几何关系可以反推出月球的角动量 h，证实中心力模型与观测高度吻合，甚至在未考虑太阳摄动之前就已十分精确。

10. Period and Kepler’s Third Law | 周期与开普勒第三定律

Integrating the areal velocity over one full orbit gives the orbital period T. The result is Kepler’s third law: T² = (4π²/μ) a³. For the Earth–Moon system, using μ = GM_earth ≈ 3.986×10⁵ km³/s² and a ≈ 384 400 km, the predicted sidereal period is about 27.3 days, matching the actual sidereal month.

对一整轨道的面积速率进行积分可以得到轨道周期 T。结果就是开普勒第三定律：T² = (4π²/μ) a³。对地月系统，使用 μ = GM_earth ≈ 3.986×10⁵ km³/s² 和 a ≈ 384 400 km，预测出的恒星周期约为 27.3 天，与实际恒星月吻合。

T = 2π√(a³/μ)

11. Perturbations and Realism | 摄动与真实复杂性

While the simple two‑body central‑force model explains the core elliptical orbit, a real Moon is subject to solar perturbations, Earth’s oblateness and tidal interactions. These cause the slow precession of the lunar perigee (apsidal precession) and oscillations in eccentricity. In Further Mathematics, such effects are often modelled by adding small perturbing terms to the differential equation, offering a glimpse into more advanced celestial mechanics.

尽管简单的二体中心力模型成功解释了基本的椭圆轨道，但真实的月球还受到太阳摄动、地球扁率和潮汐相互作用的影响。这些因素导致月球近地点缓慢进动（拱线进动）以及偏心率的振荡。在进阶数学中，这类效应通常通过在微分方程中增加小的摄动项来建模，从而一窥更高级的天体力学。

12. Summary and Exam Tips | 总结与应试技巧

The ‘Our Moon’ problem highlights the power of differential equations to describe seemingly complex celestial motion in a clean, deductive manner. When tackling such questions, always start by writing down the radial and transverse equations, state the conserved angular momentum, use the substitution u = 1/r, and carefully differentiate. Sketch the polar orbit, label peri‑ and apo‑centres, and check that your derived eccentricity and period are physically reasonable. Memorising the final solution u = μ/h² (1 + e cos θ) saves time but should always be supported by methodical derivation.

“我们的月球”问题突显了微分方程以清晰、演绎的方式描述复杂天体运动的威力。在应对此类题目时，务必先写出径向和横向方程，明确角动量守恒，使用 u = 1/r 代换，并仔细求导。绘制极坐标轨道草图，标出近地点和远地点，并检查所推导的偏心率和周期在物理上是否合理。熟记最终解 u = μ/h² (1 + e cos θ) 可以节省时间，但始终应辅以条理清晰的推导过程。

Published by TutorHao | Further Mathematics Revision Series | aleveler.com

更多咨询请联系16621398022（同微信）