开普勒方程求解器
M = E − e·sinE（牛顿迭代求 E）
Kepler Equation
Kepler Equation
📐 计算公式
M = E − e·sinE；牛顿迭代收敛快。 e=0.1、M=90° 时 E≈95.7°。
由平近点角 M 与偏心率 e 数值求解偏近点角 E。
📌 计算说明
- M = E − e·sinE；牛顿迭代收敛快。
- e=0.1、M=90° 时 E≈95.7°。
📚 深度解析：开普勒方程求解器
💡 常见使用场景
- 给定平近点角 M 与偏心率 e，数值迭代求解偏近点角 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。
❓ 常见问题（FAQ）
- 为什么需要迭代而不是直接解？
- E 同时出现在线性项与三角函数内，无初等闭式解，通常用牛顿—拉夫森迭代快速收敛。
- 偏心率 e 很大时还准吗？
- 迭代法仍适用，但初值要好（如高 e 用 M+0.85e）；近抛物轨道需特殊初值避免收敛慢。