python - 请帮助解释我试图用 Python 重写的 Matlab 函数

Question

所以我最近一直在学习 Python，我想我会通过尝试重新创建一些我在硕士学位期间使用的 Matlab 代码来挑战自己。除了看看我是否可以之外，我没有其他任何理由这样做，而且由于我不再拥有 Matlab，因此能够使用代码也很不错。该程序用于绘制围绕月球的轨道，因此，它使用轨道力学和特定模块，如 SpiceyPy。

到目前为止，我已经能够自己转换大部分代码，但特别是我坚持使用一个功能。问题是我不完全理解代码在做什么，所以我无法重新创建它。

这是函数的链接：https ://uk.mathworks.com/matlabcentral/mlc-downloads/downloads/submissions/64408/versions/1/previews/LHPOP/prop/accelpntmasses.m/index.html

任何人都可以提供的任何帮助将不胜感激。

标签： pythonmatlaborbital-mechanics

Answer 1

看起来该函数的目的是计算来自“其他行星”（指除地球以外的所有物体）的第三个物体的加速度或扰动。其他物体被视为点质量。循环遍历每个点质量，并通过调用 JPL SPICE 文件函数之一来提取它们的位置和“mu”项（GM 变量）。它累积所有质量的总加速度并将其返回。

有关 JPL SPICE 文件功能的更多信息，请参阅JPL SPICE。