› 开普勒方程模拟器 返回
天体力学模拟器

开普勒方程模拟器 — 椭圆轨道的数值解

以牛顿迭代法实时求解超越方程 M = E − e·sin(E)，给出偏近点角、真近点角、焦点距离与公转周期。轨道形状与二次收敛过程一目了然。

参数设置
半长轴 a
AU
离心率 e
平均近点角 M
deg
中心质量 M_central
M_sun

容差 1e-10、最多 30 次迭代。牛顿法为二次收敛，通常 3〜5 次即可达到 10 位精度。

计算结果
—
偏近点角 E
—
真近点角 ν
—
焦点距离 r
—
公转周期 T
椭圆轨道与卫星位置

蓝色椭圆＝轨道／黄色圆点＝中心天体（焦点）／白色圆点＝卫星位置／绿色直线＝焦点-卫星距离 r／灰色虚线＝偏近点角 E 对应的辅助圆点

牛顿迭代收敛历史

横轴＝迭代次数 n／纵轴＝残差 log10|f(E_n)|。二次收敛使位数大致每步翻倍（容差线 1e-10）。

理论与主要公式

开普勒方程描述椭圆轨道上天体位置随时间的变化。平均近点角 M 与时间成正比，但与几何意义上的偏近点角 E 之间满足超越方程，需要数值求解。

开普勒方程与牛顿迭代：

$$M = E - e\,\sin E,\qquad E_{n+1} = E_n - \frac{E_n - e\,\sin E_n - M}{1 - e\,\cos E_n}$$

真近点角 ν 与焦点距离 r：

$$\tan\!\frac{\nu}{2} = \sqrt{\frac{1+e}{1-e}}\,\tan\!\frac{E}{2},\qquad r = a\,(1 - e\,\cos E)$$

由开普勒第三定律得到的公转周期（中心质量以太阳质量为单位）：

$$T = 2\pi\sqrt{\frac{a^{3}}{G\,M_{c}}}\;\;\Longrightarrow\;\; T_{\mathrm{yr}} = \sqrt{\frac{a_{\mathrm{AU}}^{3}}{M_{c,\odot}}}$$

$a$ 为半长轴 [AU]，$e$ 为离心率（0〜1），$M$ 为平均近点角 [rad]，$E$ 为偏近点角 [rad]，$\nu$ 为真近点角 [rad]，$M_c$ 为中心质量 [太阳质量]，$T$ 为公转周期 [年]。

什么是开普勒方程模拟器

🙋
行星沿椭圆轨道运动，那有没有一个公式，输入时刻就能直接得到位置？
🎓
差一点点——意外的是没有解析解。把与时间成正比的量记作平均近点角 M，则它与从椭圆中心测量的偏近点角 E 满足开普勒方程 $M = E - e\sin E$。这是一个超越方程，无法以初等函数对 E 求解，因此只能用牛顿迭代数值求解。本模拟器正是做这件事。
🙋
咦，连行星位置都没有解析公式？
🎓
是的，这点常让初学者吃惊。开普勒方程是「超越方程」的典型例子，至今仍出现在数值分析教材的开篇。默认 e=0.20、M=60° 时，从 E_0=60° 出发的牛顿迭代二次收敛到 E_1≈71.03°、E_2≈70.823°，3 次迭代即可达 10 位精度。请观察右下方收敛图中残差快速下降的过程。
🙋
如果把离心率提到 0.9 会怎样？
🎓
好问题。e≈0.9 时轨道极为狭长，特别是 M 接近 π 时初值 E_0=M 很差，需要更多迭代。哈雷彗星级（e=0.967）的轨道往往需要 Battin 等更巧妙的初值方法。本工具为安全起见设置容差 1e-10、最大 30 次迭代，e=0.95 也能收敛。请用滑块改变 e 并观察收敛图，能直观看到二次收敛何时生效。
🙋
右上角的真近点角 ν 与 E 有什么不同？
🎓
这正是关键。E 是从椭圆中心测量的「几何角」，ν 是从焦点（太阳）实际观测的角。默认 e=0.20、M=60° 时 E=70.82°、ν=82.10° 略有偏差。圆轨道（e=0）时 M=E=ν 全部一致；椭圆中三者的差异可视为轨道本身的「身份」。把位置转换为可观测量时使用 ν。
🙋
公转周期是开普勒第三定律吧？
🎓
没错，$T^2 \propto a^3$。以太阳质量与 AU 为单位即可简化为 $T_{\mathrm{yr}} = \sqrt{a^3/M_c}$。地球（a=1, M_c=1）为 T=1 年；哈雷彗星（a=17.8）则约 75 年。本模拟器允许调节中心质量，可立即查看白矮星（0.6 M_sun）周围或巨星（10 M_sun）周围的轨道周期差异。

常见问题

当离心率 e 较小时，开普勒方程 M = E − e·sin(E) 近似于 M ≈ E，因此 E_0 = M 是良好初值。牛顿法二次收敛意味着误差 ε 经一次迭代变为 O(ε²)，两次后变为 O(ε⁴)。例如 e=0.2 时初始误差约 0.18 rad，两次迭代后降至 1e-4 rad，三次后降至 1e-10 rad。当 e > 0.6 或 M 接近 π 时初值变差，可在收敛图中清楚看到迭代次数增加。
这正是开普勒第一定律的核心：行星沿椭圆运动，太阳位于其中一个焦点。从牛顿力学求解中心力（万有引力）二体问题可得，轨道为圆锥曲线（椭圆、抛物线或双曲线），引力源恰好位于焦点。几何上看，焦点距离 r 随轨道位置变化，因此天体在近日点附近运动较快、远日点附近运动较慢（开普勒第二定律：等面积定律）。让 M 在 0〜360° 间扫描即可看到这种速度差异。
不能直接处理。当 e ≥ 1 时轨道不闭合、周期发散，需要使用专门的「双曲线开普勒方程」或「抛物线方程」（如 Barker 方程）。穿越太阳系的彗星通常采用 $M_h = e\sinh F - F$（双曲线）或 $D + D^3/3 = 2 M_p$（Barker，抛物线）。本模拟器为安全起见将 e 上限设为 0.95，仅处理椭圆情形。
是的，GPS 接收机内部实时求解开普勒方程。GPS 卫星的轨道半径约 26,600 km、离心率约 0.01，每颗卫星通过导航电文广播其轨道根数（半长轴、离心率、近地点幅角、升交点经度、轨道倾角、过近点时刻等）。接收机据此根据当前时刻计算 M，再用牛顿法求解 E 得到卫星位置，并通过多颗卫星的距离三角定位求出自身位置。由于 e 较小，通常 2〜3 次迭代即可。

实际应用

卫星轨道规划与运行：所有人造卫星、空间站和太空望远镜的轨道预测都以开普勒方程为核心。NASA 的 SGP4/SDP4 与 ESA 的 Orekit 等轨道传播库通过迭代求解开普勒方程，由初始轨道根数预测未来位置。火星探测器运行、月球与小行星着陆序列、轨道碎片碰撞预警等都基于与本模拟器相同的计算原理。

系外行星探测与表征：凌日法与视向速度法的数据分析将观测到的亮度或速度变化拟合到「开普勒轨道模型」，反演出行星质量、离心率、公转周期等参数。开普勒空间望远镜（恰好同名）发现的两千余颗系外行星，都通过 MCMC 或最小二乘反复求解本方程才得以表征。

GPS 与 GNSS 导航：GPS、Galileo、北斗、GLONASS 等所有 GNSS 系统都依赖接收机实时求解开普勒方程：从广播星历中得到轨道根数后据此计算卫星位置。手机的 GPS 功能每秒都在执行这个计算，虽然计算量不大，但其精度是实现米级定位的基础。

小行星与彗星轨道确定：隼鸟 2 号、OSIRIS-REx 等小行星采样返回任务的目标轨道，都通过对地面观测数据进行最小二乘拟合得到轨道根数，再由开普勒方程预测未来位置。近地天体（NEO）监测网络也使用同一框架，对潜在威胁的小行星轨道进行数十年尺度的持续推算。

常见误区与注意事项

最常见的误区是「平均近点角 M 与偏近点角 E 是同一个角度」。两者只在圆轨道（e = 0）时才重合，椭圆下一般不同。M 是与时间成正比的虚拟角度，仅圆轨道时才等同于实际位置角；椭圆中 M 仅是时间的参数化。真正的几何位置由 E 给出，而从太阳观测的角度由 ν 给出。当输入 M=60° 时，模拟器返回的 E=70.82° 表示卫星在通过近日点后经过 60°/360° = 1/6 周期的时间，「从中心看处于 70.82°、从太阳看处于 82.10°」。

第二个常见误区是「牛顿法总能快速收敛」。椭圆开普勒方程的 f'(E) = 1 − e·cos(E) > 0 保证牛顿单调收敛，但 e 接近 1 且 M 接近 π 时初值 E_0 = M 较差，需要多次迭代才能进入二次收敛。本工具的最大迭代次数为 30；正式的轨道确定代码常用 Halley 法（三次收敛）或 Laguerre 法。试着调到 e = 0.95、M = 170°，可观察到迭代次数明显增加。

最后，请勿过度解读开普勒方程为「椭圆轨道的本质」。它仅对纯二体问题（太阳 + 一颗行星）严格成立。真实太阳系中存在其他行星的微弱引力摄动以及相对论效应（水星近日点进动是经典案例）。GPS 与深空探测中，开普勒轨道作为基线，再通过 Cowell 法或 Encke 法逐项加入摄动项。请把本模拟器视为「纯开普勒轨道」的入门起点，而非完整轨道模型。

相关工具

开普勒定律模拟器 — 行星轨道与周期 可视化开普勒三大定律：椭圆轨道、等面积定律和周期-轨道关系。调节轨道参数，观察行星运动和扫过面积的变化。 开普勒轨道模拟器 交互式开普勒轨道模拟器。调整半长轴与离心率，实时动画演示椭圆轨道。可视化等面积定律（第二定律），验证 T²=a³（第三定律）。含地球、火星、水星、哈雷彗星预设。 谢乐方程模拟器 — 由 X 射线衍射估算晶粒尺寸 谢乐方程 D = Kλ/(β cosθ) 实时模拟器：由 X 射线波长 λ、布拉格角 θ、实测 FWHM β_total、仪器 FWHM β_inst 计算晶粒尺寸 D、校正后 FW... 轨道力学模拟器 — 开普勒轨道与霍曼转移 模拟开普勒轨道，调节半长轴和离心率，可视化霍曼转移轨道、轨道周期和比轨道能量，适用于航天工程和天体力学教学。 轨道力学与开普勒轨道模拟器 免费在线轨道力学模拟器，可视化圆、椭圆、抛物线及双曲线轨道。提供LEO、ISS、GEO、月球、火星等预设，自动计算轨道周期、速度与能量。探索开普勒定律，助力航天学习与模拟。 抛体与轨道模拟器 通过逐步提高初速度，直观演示抛体如何从弹道飞行连续过渡到圆轨道、椭圆轨道，直至双曲线逃逸。本模拟器清晰揭示了开普勒定律与轨道力学的核心原理，并展示了圆轨道速度、逃逸速度及轨道周期的关...