第 3 章 · 开普勒问题的朗道式解法
开普勒问题在朗道笔下不是例题而是方法展示:用守恒量改写、用巧积分求解、用极限行为分类。这一节把它当作"解力学问题的一套手法"来学。
轨道方程的推导流水线
从能量方程 $E = \frac{\mu\dot r^2}{2} + U_\text{eff}(r)$ 出发:
代入 $U = -\kappa/r$,积分结果即圆锥曲线 $r = \frac{p}{1+e\cos\phi}$($p = M^2/\mu\kappa$)。学习重点不是背结果,而是流水线本身:守恒量消时间 → 径向积分 → 反解轨道。这个流水线对任何中心势通用(换势换积分难度而已)。
能量的符号分类
Runge-Lenz 矢量 $\mathbf{A} = \mathbf{v}\times\mathbf{M} - \kappa\hat r$ 的守恒性解释了"轨道不进动 + 周期只依赖 E"的反常丰富对称性——它是简并的几何化身,量子化后正对应氢原子的"意外简并"。
- $E < 0$:椭圆($e < 1$),周期 $T = \pi\kappa\sqrt{\mu/2|E|^3}$——周期只由能量决定(与角动量无关),这是反平方力的隐藏对称性;
- $E = 0$:抛物线(逃逸临界);
- $E > 0$:双曲线(散射轨道),散射角与瞄准距离的公式直接可读。
频率、进动与微扰视角
一般中心势的径向振荡频率与角向频率之比 $\omega_r/\omega_\phi$ 决定轨道是否闭合——不闭合的轨道呈玫瑰花瓣状缓慢进动。这个视角把"反平方力多特殊"量化:只有它让两频率严格可通约(Bertrand 定理的另一半)。
老师的提醒
- 三条关系要背熟:$p = M^2/\mu\kappa$、$e^2 = 1 + 2EM^2/\mu\kappa^2$、$a = -\kappa/2E$——考题与天体力学的换算全靠它们。
- 对比学习:谐振子中心力场($U \propto r^2$)同样闭轨道但周期只依赖"劲度"不依赖能量——两个 Bertrand 特例对照着记。
自检
- 推导轨道方程的积分流水线。
- 周期为什么只依赖能量?(Runge-Lenz 矢量与此的关系)
- 椭圆、抛物、双曲轨道分别对应什么能量符号?
更深一层:开普勒三参数的"速查卡"
开普勒轨道的参数换算卡(工程与天体力学高频使用,值得随身):
换算卡的天文学应用:给一组轨道观测(位置+速度),三步换算出轨道全貌($E$ → $a$ → 周期;$M$ → $p$;合出 $e$)——人造卫星入轨参数的快速校核流程。
- 半长轴 $a = -\frac{\alpha}{2E}$(能量唯一决定半长轴——与角动量无关,反平方力的特权);
- 半通径 $p = \frac{M^2}{\mu\alpha}$(角动量决定"胖瘦");
- 偏心率 $e = \sqrt{1 + \frac{2EM^2}{\mu\alpha^2}}$($e = 0$ 圆、$< 1$ 椭圆、$= 1$ 抛物线、$> 1$ 双曲);
- 周期 $T = \pi\alpha\sqrt{\frac{\mu}{2|E|^3}}$(椭圆,只依赖能量);
- 近日点速度 $v_p = \sqrt{\frac{\mu}{a}\frac{1+e}{1-e}}$。
Runge-Lenz 矢量 $\mathbf{A} = \mathbf{v}\times\mathbf{M} - \alpha\hat{\mathbf{r}}$ 是隐藏的第六参数(守恒、指向近日点),它的存在解释了"周期不依赖角动量"与"轨道不进动"这两件反平方力的特权——额外的守恒律 = 额外的对称性 = 意外的简并(量子版即氢原子的 $n$ 简并)。
在真实世界里:近日点进动——"反平方之外"的信号
真实行星轨道会缓慢进动(水星每百年 43 角秒)——因为太阳非完美球体 + 广义相对论修正。进动是"偏离反平方"的精密探针:任何微小非平方分量都让椭圆不再闭合、近日点缓慢漂移。广义相对论的一阶修正恰好补齐水星的 43 角秒,成为 GR 的三大经典验证之一。"轨道闭合性"是检验力律的最灵敏仪器——这个思想从天体力学延伸到核物理(壳模型中的轨道简并检验)。
专题深入:开普勒三参数的"速查卡"完整版
开普勒轨道的参数换算卡(工程与天体力学高频使用,值得随身):
- 半长轴 $a = -\frac{\alpha}{2E}$(能量唯一决定半长轴——与角动量无关,反平方力的特权);
- 半通径 $p = \frac{M^2}{\mu\alpha}$(角动量决定"胖瘦");
- 偏心率 $e = \sqrt{1 + \frac{2EM^2}{\mu\alpha^2}}$($e = 0$ 圆、$< 1$ 椭圆、$= 1$ 抛物线、$> 1$ 双曲);
- 周期 $T = \pi\alpha\sqrt{\frac{\mu}{2|E|^3}}$(椭圆,只依赖能量);
- 近日点速度 $v_p = \sqrt{\frac{\mu}{a}\frac{1+e}{1-e}}$。
换算卡的天文学应用:给一组轨道观测(位置+速度),三步换算出轨道全貌($E$ → $a$ → 周期;$M$ → $p$;合出 $e$)——人造卫星入轨参数的快速校核流程。
Runge-Lenz 矢量 $\mathbf{A} = \mathbf{v}\times\mathbf{M} - \alpha\hat{\mathbf{r}}$ 是隐藏的第六参数(守恒、指向近日点),它的存在解释了"周期不依赖角动量"与"轨道不进动"这两件反平方力的特权——额外的守恒律 = 额外的对称性 = 意外的简并(量子版即氢原子的 $n$ 简并)。
数学细节展开:谐振子中心力场的"另一个 Bertrand 特例"
谐振子中心力场($U \propto r^2$)同样闭轨道但周期只依赖"劲度"不依赖能量——两个 Bertrand 特例对照着记:反平方力(开普勒)与谐振子是仅有的两种"全部轨道闭合"的中心势,它们的周期-能量关系恰好相反(开普勒周期只依赖能量、谐振子周期只依赖劲度)。这个对照揭示了"反平方力多特殊"的定量版:只有它让径向与角向频率严格可通约(Bertrand 定理的另一半)。
工程案例细节:开普勒公式的"工程换算"实战
开普勒公式的工程换算实战:轨道设计(霍曼转移的两条椭圆 $a_1 < a_2$、两次变轨脉冲 $\Delta v$ 由能量差给出);引力弹弓(行星飞掠的能量交换——日心系里 $E$ 变化但行星轨道近似不变);卫星入轨(入轨点位置+速度 → $E$ → $a$ → 周期 → $M$ → $p$ → $e$——三步换算出轨道全貌)。三条关系要背熟:$p = M^2/\mu\alpha$、$e^2 = 1 + 2EM^2/\mu\alpha^2$、$a = -\alpha/2E$——考题与天体力学的换算全靠它们。
例题
开普勒轨道:某卫星 $E = -2\times10^7\ \mathrm{J/kg}$、$M = 6\times10^{10}\ \mathrm{m^2/s}$($\mu = GM_\odot = 1.33\times10^{20}$)。求半长轴 $a$ 与偏心率 $e$。
解答
半长轴:$E = -\frac{\mu}{2a}$
偏心率:$e^2 = 1 + \frac{2EM^2}{\mu^2}$:
近抛物线的狭长椭圆($a$ 大、$e\to1$)——像一颗长周期彗星轨道。
- 三条公式($a = -\mu/2E$、$p = M^2/\mu$、$e^2 = 1 + 2EM^2/\mu^2$)把能量角动量一键翻译成轨道形状。