Question

我最近下载了pyephem。我非常喜欢用户界面。

然而，我得到了非常不稳定的结果，特别是当我试图在过去的某个时刻获得对恒星衰退的估计时。我已经使用了与http://rhodesmill.org/brandon/projects/pyephem-manual.html#FixedObjects,Precession,andEpochs中描述的Thuban / Polaris示例关联的代码，我得到了完全不同的结果。我想知道是否有其他东西需要添加到我的脚本中才能使事情正常工作，或者如果pyephem中有新的错误。

我的脚本如下所示：

import ephem
import ephem.stars

polaris = ephem.readdb("Polaris,f|M|F7,2:31:48.704,89:15:50.72,2.02,2000")
polaris.compute()    

thuban = ephem.readdb("Thuban,f|V|A0,14:4:23.3,64:22:33,3.65,2000")
thuban.compute()

print 'my result for current polaris dec',polaris.dec
print 'webpage result                    89:15:50.73\n'

print 'my result for current thuban dec',thuban.dec
print 'webpage result                   64:22:32.99\n'

polaris.compute(epoch='-2800')
print 'my result for polaris dec in -2800',polaris.dec
print 'webpage result                     63:33:17.63\n'

thuban.compute(epoch='-2800')
print 'my result for thuban dec in -2800',thuban.dec
print 'webpage result                    89:54:34.97\n'

运行它会产生以下输出。我可以理解为什么前两个比较略有不同（网页上的例子是在最近创建的），但你可以看到，对于任何一个恒星，当前日期和时期-2800之间的偏差没有变化，但是在网页上的例子，偏差在两个时代之间对两颗星都有很大影响：

my result for current polaris dec 89:19:19.5
webpage result                    89:15:50.73

my result for current thuban dec 64:18:39.5
webpage result                   64:22:32.99

my result for polaris dec in -2800 89:19:19.5
webpage result                     63:33:17.63

my result for thuban dec in -2800 64:18:39.5
webpage result                    89:54:34.97

Answer 1

.dec属性会以多种方式自动调整，这会干扰教程尝试显示进动如何工作，所以 - 如果你仔细查看教程的代码 - 我想你会发现它总是要求Polaris和Thuban的.a_dec“天体测量偏差”。如果您通过在属性名称的前面添加a_来对上面的代码进行轻微的拼写更正，那么我认为您会发现您的打印答案与教程中打印的答案一致！

Pyephem：在时代计算中没有正确应用岁差？

1 个答案: