我正在研究可能是一个独特的用例 - 我想使用Skyfield对假想的星系进行一些计算。我会通过创建我自己的星历,并使用它而不是实际的星历来做到这一点。我发现的问题是我找不到API的文档来替换我自己的星历表。
有文件吗? Skyfield是否有足够的灵活性来做我正在尝试的事情?
编辑: 为了澄清我的要求,我明白我将不得不做一些引力模型(我完全愿意配置这个房子里的每台电脑,平板电脑,有线电视盒和烤面包机,以便在几天内处理这些数字:),但在我真正深入研究之前,我想知道数据是什么样的。如果它只是一个带有许多名为numpy 2d数组的模块......这使得它相当容易,但我没有在任何地方看到这个。
答案 0 :(得分:1)
Skyfield使用的JPL发布的星历表,如DE405和DE406以及DE421,只是为每个星球提供了一个大数字表。例如,海王星的位置可以以7天的增量指定,其中对于从星历的开始到结尾的每个7天周期,该表提供了一组多项式系数,可用于估计海王星在任何时刻的位置。从那个7天开始到结束。如果我理解的话,多项式是经过设计的,这样它们的一阶和二阶导数就可以在一个结束和下一个结束的那一刻与之前和之后的7天多项式平滑地进行网格划分。
JPL通过记录行星的位置来生成这些巨大的表格,就像我们在人类历史上记录它们一样,采用我们认为理想行星将根据引力理论移动的规则,太阳风的阻力,行星的自己的旋转和动力学,它的卫星等,并试图选择一个符合理论的行星的“真实路径”,同时尽可能接近实际观察到的位置。
这是一个很大的计算问题,我认为它需要相当多的技巧。如果你无法完美地匹配所有的观察结果 - 你永远无法做到 - 那么你必须决定哪些是优先考虑的,哪些可能不是那么准确。
对于一个假想的系统,你将不得不从头做起(可能是?)一个引力动力学模拟。如果我理解正确的话,有关这个主题的各种教科书中记录的几种可能的方法。无论您选择哪一个,都可以让您为假想的行星生成 x , y , z 位置,并且您可能会在Skyfield中将这些位置实例化为ICRS如果你想使用Skyfield来计算距离,观察或绘制图表,那就是位置。
虽然我没有自己使用它,但我看到了很好的评论:
http://www.amazon.com/Solar-System-Dynamics-Carl-Murray/dp/0521575974