1

我最近下载了 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
4

1 回答 1

0

该.dec属性以多种方式自动调整，这会干扰教程试图显示岁差如何工作，所以——如果你仔细查看教程的代码——我想你会发现它总是要求.a_dec“天文偏角” ”北极星和图班。如果您通过添加到属性名称的前面来对上面自己的代码进行轻微的拼写更正a_，那么我想您会发现您的打印答案与教程中打印的答案一致！

于 2013-04-20T16:32:58.147 回答