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

开普勒方程 数值解模拟器 — Newton 法与迭代收敛

使用 Newton-Raphson 法求解 M = E − e sin E 的中高级数值模拟器。每一帧实时计算离心近点角 E 与真近点角 ν，可观察迭代次数与残差收敛过程。法则概念理解请参考《开普勒法则》，轨道整体行为请参考《开普勒轨道》。

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

容差 1e-10，最多 30 次迭代进行收敛判定。Newton 法为二次收敛，通常 3～5 次迭代即可达到 10 位精度。

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

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

Newton 迭代收敛履历

横轴=迭代次数 n／纵轴=残差 log10|f(E_n)|／二次收敛的特点是每次迭代后有效数字数量大约翻倍增加（容差线 1e-10）。

理论与主要公式

开普勒方程是从时刻求出椭圆轨道天体位置的超越方程。平均近点角 M 与时间成正比，但实际位置对应的离心近点角 E 与 M 有非线性关系，需要数值求解。

开普勒方程与 Newton 迭代：

$$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，真近点角 ν 和焦点距离 r 就立刻能算出来，这样就知道卫星在哪儿。但问题是 E 关于 E 的这个方程属于超越方程，解析解不存在，只能用 Newton 法数值求解。本模拟器就是演示这个过程的。
🙋
连行星位置都不能用公式直接算出来吗？
🎓
这确实让人惊讶，但这就是现实。所以超越方程一直是数值计算教科书开篇的经典例子。默认参数 e=0.20, M=60° 时，从 E_0 = 60° 开始，Newton 迭代后 E_1 ≈ 71.03°, E_2 ≈ 70.823°，二次收敛，3 次迭代就精确到 10 位数。右下角的收敛图能看到残差一下子掉下来，就是这个效果。
🙋
离心率增大的话会怎么样？
🎓
好眼光。离心率 e 表示轨道有多"拉长"，0 是圆形，越接近 1 越细长。e 增大，Newton 法收敛变慢，特别是当 e > 0.6 并且 M 接近 π 时，初值 E_0 = M 不再是好近似，需要更多迭代。e=0.967 的哈雷彗星那水平，本模拟器的上限 e=0.95 也接近了。调滑块把 e 改成 0.95，看收敛图，迭代次数会明显增加，能直观感受这一点。
🙋
那个真近点角 ν，和 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）T≈75 年。本模拟器还能改中心质量，白矮星（0.6 M_sun）或巨星（10 M_sun）周围的周期马上就能算。

常见问题

当离心率 e 较小时，开普勒方程 M = E − e·sin(E) 近似为 M ≈ E，所以 E_0 = M 是不错的近似值。Newton 法是二次收敛，初值偏离真值 ε，一次迭代后偏离变成 O(ε²)，两次迭代后 O(ε⁴)。e=0.2 的情况下初值偏离约 0.18 rad，经 2 次迭代变成 10⁻⁴ rad，3 次迭代变成 10⁻¹⁰ rad，快速收敛。但 e > 0.6 或 M 接近 π 时，初值质量下降，迭代次数会增加，本模拟器能观察到这一点。
这是开普勒第一定律的核心：行星运行于椭圆轨道，太阳位于椭圆的一个焦点。这是从牛顿力学导出的中心力（万有引力）二体问题解的结论——轨道必为圆锥曲线，引力源必位于焦点。直观上讲，焦点处太阳到卫星的距离 r 随位置变化（椭圆的几何性质），卫星在近日点速度快、远日点速度慢（面积速度不变 = 开普勒第二定律）。模拟器把 M 从 0～360° 扫一遍，就能看到卫星在近日点附近走得快、远日点走得慢。
不能直接用。e ≥ 1 轨道不再闭合，周期趋向无穷，需要用专门的双曲线或抛物线方程（如 Barker 方程等）。彗星横贯太阳系时走的是抛物线或双曲线，对应的是不同的超越方程。本模拟器为安全起见上限 e=0.95，保持在椭圆轨道范围内。
是的，GPS 接收机内部实时求解开普勒方程。GPS 卫星绕地球的轨道半径约 26,600 km、离心率约 0.01（近乎圆形）。每颗卫星广播自己的轨道根数（半长轴、离心率、近点角、升交点经度、轨道倾斜、近点通过时刻等）。接收机据此计算当前 M，用 Newton 法求 E，进而求出卫星位置，结合多颗卫星三角定位出自己的位置。e 很小，通常 2～3 次迭代就够了。

现实应用

卫星轨道规划与运营：所有人造卫星、空间站、航天望远镜的轨道预报都以开普勒方程为核心。NASA 的 SGP4/SDP4、ESA 的 Orekit 等轨道库，将初始轨道根数传播到将来的任意时刻时，都要反复求解开普勒方程。火星探测器运营、月球与小行星着陆序列、轨道碎片碰撞预报，无不需要这个计算。

系外行星探测与表征：凌星法和径向速度法的数据分析中，观测的亮度与速度曲线要拟合到"开普勒轨道模型"，用最小二乘或 MCMC 反复求解开普勒方程来推断行星参数（质量、离心率、公转周期）。Kepler 空间望远镜（巧合同名）发现的 2000+ 系外行星都是这样表征的。

GNSS 测位系统：GPS、Galileo、BeiDou、GLONASS 等所有全球导航系统，接收机都需实时从广播的星历根数求解开普勒方程，算出卫星位置。这是每秒进行的计算，虽然计算量小，但精度直接影响米级定位的可靠性。

小行星与彗星轨道确定：小行星探测任务（隼鸟2号、OSIRIS-REx 等）的标的天体轨道，由地面望远镜观测数据通过最小二乘拟合轨道根数，再用开普勒方程预报将来位置。近地天体（NEO）监测网络也靠这个持续跟踪、预测数十年后的潜在危险小行星位置。

常见误区与注意

最常见的误解是以为"M（平均近点角）和 E（离心近点角）是同一个角"。只有在圆轨道（e=0）时两者才相等，椭圆轨道下它们是不同的。M 是时间的代数表示（按时间线性增长），E 则是中心那边的几何角，两者通过超越方程相连。真实位置角是 ν（太阳视角）。当输入 M=60° 时输出 E=70.82° 的含义是：卫星出发后经过 60°/360°=1/6 个周期，此时中心角 70.82°、太阳视角 82.10°。

另一个常见误解是"Newton 法总是快速收敛"。虽然椭圆开普勒方程 f'(E) = 1 − e·cos(E) > 0 保证了单调收敛，但当 e 接近 1、M 接近 π 时，初值 E_0 = M 的质量变差，二次收敛前可能需多次迭代。本工具限制最多 30 次迭代，容差 1e-10。若尝试 e=0.95、M=170°，会看到迭代次数明显增多（收敛图折线段数增加）。精确轨道计算有时采用 Halley 法（三阶收敛）或 Laguerre 法等高阶方法。

最后，切勿"过度夸大开普勒方程的通用性"。它严格适用于二体问题（太阳+单个行星），实际太阳系中其他行星的重力扰动虽小但存在，相对论效应也有（水星近日点进动为典型例子）。卫星轨道计算和深空探测采用 Cowell 法、Encke 法等在开普勒轨道基础上加入摄动项的方法。本模拟器代表的是"纯开普勒轨道"的起点，是理解更复杂模型的基础。

关联工具

开普勒法则模拟器 通过轨道动画和面积速度，直观确认开普勒第 1、2、3 法则的模拟器。可交互式探索离心率、轨道周期、行星间距离的关系，天文学教学工具。 开普勒轨道模拟器 开普勒轨道模拟器实时可视化行星运动。调整离心率和长半轴，在实时动画中体验椭圆轨道和速度变化。万有引力作用下哈雷彗星、地球、火星轨道的动画演示... 舍勒公式 模拟器 — X 射线衍射结晶尺度估算 舍勒公式 D = Kλ/(β cosθ) 实时模拟器。从 X 射线波长 λ、Bragg 角 θ、实测 FWHM β_total、仪器 FWHM β_inst 推算结晶子尺... 轨道力学模拟器 基于开普勒法则的椭圆轨道实时动画。可交互式探索离心率、轨道长半轴、周期相互关系的天文学工具。 轨道力学模拟器 轨道力学模拟器，调整长轴、离心率、倾斜角，实时可视化开普勒椭圆轨道。Newton-Raphson 法数值求解开普勒方程，精确计算卫星位置和运动动画... 轨道·抛体模拟器 轨道·抛体模拟器"NovaSolver"，仅需调整初速度，即可从抛体飞行延伸到圆形轨道、椭圆轨道，乃至双曲线脱出轨道，实时可视化。开普勒法则、轨道周期、离心...

使用指南

  1. 用滑块设置轨道半长轴 a [AU] 和离心率 e（例如：a=6878 km、e=0.0015）
  2. 输入平均近点角 M [度]，启动 Newton-Raphson 法迭代计算
  3. 离心近点角 E 和真近点角 ν 实时收敛，同时输出焦点距离 r 和周期 T

计算示例

国际空间站轨道（a=6878 km、e=0.0015、M=45°）的情况下，Newton 法从初值 E₀=M 开始，3～4 次迭代收敛到 E≈45.04°。焦点距离 r=6877.6 km，公转周期 T≈90.3 分钟。地球轨道（a=1.496×10⁸ km、e=0.0167、M=100°）时 E≈100.98°、周期 T=365.25 天，可用于轨道传播初始条件。

实务建议