为什么需要迭代而不是直接解？

E 同时出现在线性项与三角函数内，无初等闭式解，通常用牛顿—拉夫森迭代快速收敛。

偏心率 e 很大时还准吗？

迭代法仍适用，但初值要好（如高 e 用 M+0.85e）；近抛物轨道需特殊初值避免收敛慢。

Kepler Equation

📐 计算公式

M = E − e·sinE；牛顿迭代收敛快。 e=0.1、M=90° 时 E≈95.7°。

由平近点角 M 与偏心率 e 数值求解偏近点角 E。

平近点角 M (°)

偏心率 e

📌 计算说明

示例：M=1.0 rad、e=0.3 求 E

解 M=E−e·sinE，即 E−0.3·sinE=1.0。用牛顿迭代：取初值 E₀=1.0，E₁=1.0−(1.0−0.3·sin1.0−1.0)/(1−0.3·cos1.0)=1.0+0.2521/1.1620≈1.2170；再迭代一次得 E≈1.2189 rad。