开普勒轨道求解器

数值求解开普勒方程 E − e·sin(E) = M,得到偏近点角 E、真近点角 ν、轨道半径与相对速度,支持牛顿/二分/不动点三种方法、收敛历史与 SVG 轨道图

给定平近点角 M(度)、偏心率 e,以及可选的半长轴 a(AU),本工具求解超越方程

M = E − e · sin(E)

得到偏近点角 E。再由 E 推导:

  • 真近点角 ν:tan(ν/2) = √((1+e)/(1−e)) · tan(E/2)
  • 轨道半径:r = a · (1 − e·cos E)
  • 相对速度(活力方程):v/v_circ = √(2 − r/a)

提供三种求解方法:

  • 牛顿-拉夫森(推荐):使用解析导数 f′(E) = 1 − e·cos E,二次收敛。高偏心率 (e>0.8) 时以 E = π 作为初值。
  • 二分法:稳健,括号区间 [M−e, M+e],线性收敛。
  • 不动点迭代 E ← M + e·sin E:最简单,但 e → 1 时发散。

报告以 SVG 展示轨道(焦点居中的椭圆,标注近日点、远日点、焦点与当前位置及半径向量),并给出求解量、轨道类型(圆/椭圆/抛物/双曲)以及完整的迭代历史与残差,便于比较三种方法的收敛。

示例结果

1 个示例

求解中等偏心率的椭圆轨道

M=45°, e=0.3, a=1 AU,对比 E、ν、r 与牛顿迭代历史。

Kepler solution with orbit SVG and iteration history.
查看输入参数
{ "M": 45, "e": 0.3, "a": 1, "method": "newton", "maxIter": 100, "tol": "1e-12" }

关键信息

分类
地理与科学
输入类型
number, select
输出类型
html
样本覆盖
4
支持 API
Yes

概览

开普勒轨道求解器是一款专为天体力学和航天动力学设计的在线计算工具,能够数值求解超越方程 E − e·sin(E) = M。输入平近点角、偏心率及半长轴,即可快速获取偏近点角、真近点角、轨道半径与相对速度,并直观展示 SVG 轨道几何图示与算法收敛历史。

适用场景

  • 需要求解给定平近点角和偏心率下的偏近点角与真近点角,以确定天体在轨道上的具体位置时。
  • 对比牛顿法、二分法和不动点迭代法在不同偏心率(尤其是高偏心率)下的收敛速度与残差表现时。
  • 需要快速生成包含近日点、远日点及当前位置标注的 SVG 轨道几何图示用于教学或报告时。

工作原理

  • 输入平近点角 M、偏心率 e 以及半长轴 a,并选择求解算法(如牛顿法)与收敛容差。
  • 求解器根据选定算法在设定的最大迭代次数内数值求解开普勒方程 M = E − e·sin(E),计算出偏近点角 E。
  • 利用解析公式推导真近点角 ν、轨道半径 r 和相对速度,并判断轨道类型。
  • 渲染交互式 SVG 轨道图,并输出包含每一步迭代残差的收敛历史表格。

使用场景

天文学与天体力学教学中,演示开普勒方程的数值求解过程及不同数值算法的收敛特性。
航天器轨道设计与任务规划中,根据轨道根数快速计算特定时间点的航天器位置与速度。
科普与学术报告中,快速生成标准且美观的 SVG 轨道几何示意图。

用户案例

1. 求解中等偏心率椭圆轨道

天文学专业学生
背景原因
在完成天体力学作业时,需要计算一颗偏心率为 0.3 的小行星在平近点角为 45° 时的轨道位置。
解决问题
手动求解开普勒超越方程非常困难,且需要精确的偏近点角和真近点角数据。
如何使用
在求解器中输入平近点角 M = 45,偏心率 e = 0.3,半长轴 a = 1.0,选择“牛顿法”并设置收敛容差为 1e-12,点击计算。
示例配置
{
  "M": 45,
  "e": 0.3,
  "a": 1,
  "method": "newton",
  "maxIter": 100,
  "tol": "1e-12"
}
效果
成功解得偏近点角 E、真近点角 ν 和轨道半径 r,并获得清晰的 SVG 轨道图和仅需数次迭代即收敛的记录。

2. 比较高偏心率下的算法收敛性

轨道动力学研究人员
背景原因
研究高偏心率(如 e = 0.85)轨道下,不同数值求解方法的稳定性和收敛速度。
解决问题
需要对比牛顿法与不动点迭代法在极端偏心率下的表现。
如何使用
输入 M = 60,e = 0.85,分别选择“牛顿法”和“不动点迭代”运行求解,观察迭代历史。
示例配置
{
  "M": 60,
  "e": 0.85,
  "a": 1.5,
  "method": "fixedpoint",
  "maxIter": 100,
  "tol": "1e-9"
}
效果
发现不动点迭代法在高偏心率下难以收敛或迭代次数极多,而牛顿法通过合理的初值选择在极少步骤内即达到 1e-9 的精度。

用 Samples 测试

image, svg, barcode

相关专题

常见问题

为什么在高偏心率下不动点迭代法会失效?

不动点迭代的收敛条件要求导数绝对值小于 1。当偏心率 e 接近 1 时,该方法在某些区间会发散,建议使用牛顿法或二分法。

半长轴 a 设为 0 会有什么影响?

如果半长轴 a 设为 0,计算出的轨道半径 r 将以半长轴 a 为单位进行相对表示。

牛顿法在高偏心率下如何保证收敛?

当偏心率 e > 0.8 时,求解器会自动采用 E = π 作为初值,以避免牛顿法因初值不当而发散。

该工具支持哪些轨道类型?

工具支持对圆轨道、椭圆轨道、抛物线轨道和双曲线轨道进行分类与计算。

如何查看算法的收敛过程?

求解结果中会提供完整的迭代历史表格,展示每一次迭代的估算值和残差,方便对比不同算法的收敛效率。

API 文档

请求端点

POST /zh/api/tools/kepler-orbit-solver

请求参数

参数名 类型 必填 描述
M number -
e number -
a number -
method select -
maxIter number -
tol select -

响应格式

{
  "result": "
Processed HTML content
", "error": "Error message (optional)", "message": "Notification message (optional)", "metadata": { "key": "value" } }
HTML: HTML

AI MCP 文档

将此工具添加到您的 MCP 服务器配置中:

{
  "mcpServers": {
    "elysiatools-kepler-orbit-solver": {
      "name": "kepler-orbit-solver",
      "description": "数值求解开普勒方程 E − e·sin(E) = M,得到偏近点角 E、真近点角 ν、轨道半径与相对速度,支持牛顿/二分/不动点三种方法、收敛历史与 SVG 轨道图",
      "baseUrl": "https://elysiatools.com/mcp/sse?toolId=kepler-orbit-solver",
      "command": "",
      "args": [],
      "env": {},
      "isActive": true,
      "type": "sse"
    }
  }
}

你可以串联多个工具,比如:`https://elysiatools.com/mcp/sse?toolId=png-to-webp,jpg-to-webp,gif-to-webp`,最多20个。

如果遇见问题,请联系我们:[email protected]