开普勒三大定律

在堆砌公式前，我们先建立直观的物理图像：

第一定律（椭圆定律）：想象在一张蹦床上放一个沉重的保龄球（中心天体），蹦床凹陷形成了一个"引力势阱"。当你滚出一个弹珠（航天器）时，只要速度不够逃逸，弹珠就会在势阱周围绕行，走出一条椭圆形的轨迹。保龄球不在椭圆的正中心，而是位于势阱最深的地方------即椭圆的一个焦点上。
第二定律（面积定律）：回想花样滑冰运动员在旋转时收回双臂会转得更快。航天器靠近中心天体时，引力做正功，势能转化为动能，速度加快，扫过的"扇形"变得细长；远离时动能转化为势能，速度减慢，扫过的"扇形"变得宽扁。一快一慢，刚好让这两块扇形的面积相等。
第三定律（调和定律）：想象用绳子甩一个重物。绳子越长（半长轴越大），离心力对抗引力所需的时间就越长，转一圈自然越慢。轨道越高，不仅路径变长，速度还变慢，因此周期呈非线性（3/2次方）急剧增加。

在**地心惯性坐标系（ECI）**下，我们基于二体问题假设（航天器质量 m2≪m_2 \llm2≪ 中心天体质量 m1m_1m1，且无其他摄动力），推导三大定律的数学表达。

航天器相对于中心天体的运动方程为：

r¨=−μr3r \ddot{\mathbf{r}} = -\frac{\mu}{r^3}\mathbf{r} r¨=−r3μr

通过引入比角动量 h=r×r˙\mathbf{h} = \mathbf{r} \times \dot{\mathbf{r}}h=r×r˙ 并与位置矢量 r\mathbf{r}r 叉乘积分，可得轨道的极坐标方程：

r=p1+ecos⁡ν=a(1−e2)1+ecos⁡ν r = \frac{p}{1 + e \cos\nu} = \frac{a(1-e^2)}{1 + e \cos\nu} r=1+ecosνp=1+ecosνa(1−e2)

rrr：航天器到中心天体的距离，单位：km
ppp：半通径，单位：km
aaa：轨道半长轴，描述轨道大小，单位：km
eee：轨道偏心率，描述轨道形状，无量纲（000为圆，0<e<10<e<10<e<1为椭圆）
ν\nuν：真近点角，表示航天器与近地点的角距离，单位：rad 或度
μ\muμ：引力参数（μ=G(m1+m2)≈Gm1\mu = G(m_1+m_2) \approx G m_1μ=G(m1+m2)≈Gm1），单位：km3/s2\text{km}^3/\text{s}^2km3/s2（如地球 μ≈398600.4418 km3/s2\mu \approx 398600.4418 \text{ km}^3/\text{s}^2μ≈398600.4418 km3/s2）
结论：当 0<e<10 < e < 10<e<1 时，上述方程为椭圆方程，原点位于椭圆的一个焦点上。

比角动量 h\mathbf{h}h 是守恒量（因为引力始终指向中心，无力矩）。其大小 h=r2ν˙h = r^2 \dot{\nu}h=r2ν˙。

航天器在 Δt\Delta tΔt 时间内扫过的面积微元为：

dA=12r⋅rdν=12r2ν˙dt=h2dt dA = \frac{1}{2} r \cdot r d\nu = \frac{1}{2} r^2 \dot{\nu} dt = \frac{h}{2} dt dA=21r⋅rdν=21r2ν˙dt=2hdt

因此，面积速率为：

dAdt=h2=const \frac{dA}{dt} = \frac{h}{2} = \text{const} dtdA=2h=const

对第二定律积分，轨道运行一周（周期为 TTT），扫过整个椭圆面积 A=πabA = \pi a bA=πab（其中 b=a1−e2b = a\sqrt{1-e^2}b=a1−e2 为半短轴）：

T=AdA/dt=πa21−e2h/2 T = \frac{A}{dA/dt} = \frac{\pi a^2 \sqrt{1-e^2}}{h/2} T=dA/dtA=h/2πa21−e2

由半通径定义 p=a(1−e2)=h2/μp = a(1-e^2) = h^2/\mup=a(1−e2)=h2/μ，可得 h=μa(1−e2)h = \sqrt{\mu a(1-e^2)}h=μa(1−e2) ，代入上式化简得：

T=2πa3μ T = 2\pi \sqrt{\frac{a^3}{\mu}} T=2πμa3

或者：

T2=4π2μa3 T^2 = \frac{4\pi^2}{\mu} a^3 T2=μ4π2a3

在真实航天任务中，开普勒定律是所有轨道设计的基石，但必须清楚其模型边界与工程约束：

轨道设计（如星座部署） ：利用第三定律，低地球轨道（LEO，a≈7000a \approx 7000a≈7000 km）周期约 90 分钟，地球静止轨道（GEO，a≈42164a \approx 42164a≈42164 km）周期恰好 24 小时。这是所有通信卫星选址的核心依据。
速度脉冲（ΔV\Delta VΔV）规划：利用第二定律的推论（近地点速度最大，远地点速度最小），我们在近地点进行变轨抬升远地点（如霍曼转移），能最大化奥伯特效应，最省燃料。
边界意识（极重要的误区） ：开普勒定律严格成立的前提是无摄二体问题 。在真实任务中：
- J2摄动 ：地球并非正球体（赤道隆起），这会导致轨道面进动（升交点赤经Ω\OmegaΩ漂移）和近地点旋转（近地点幅角ω\omegaω漂移）。此时轨道不再是封闭的固定椭圆，第三定律计算的周期需修正为"交点周期"或"近点周期"。
- 大气阻力 ：在低于 1000 km 的轨道，大气阻力使航天器持续失去能量，半长轴 aaa 不断减小，按第三定律周期 TTT 会变短，航天器反而越飞越快（动能增加）。
- 三体引力：深空探测或高轨卫星受日月引力影响，轨道不再是纯粹的开普勒椭圆。

进阶问题：如果考虑地球扁率（J2摄动），开普勒第三定律中的"周期"还成立吗？我们工程上常说的"交点周期"和"近点周期"与开普勒的无摄周期有何区别？

提示：J2摄动使得轨道不再是一个封闭椭圆，航天器每圈经过赤道（交点）和经过近地点的时间都会发生微小偏移，思考这二者与开普勒理论周期的定量差异，是走向高精度轨道确定的第一步。