中华农历论坛知识讨论区历法知识 → 请教de406行星历表的结构


  共有24019人关注过本帖树形打印

主题:请教de406行星历表的结构

美女呀,离线,留言给我吧!
蔡越
  1楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
请教de406行星历表的结构  发帖心情 Post By:2009/4/27 13:34:00

1  728
......

475  728

1是什么意思? 728又是什么意思?

哪位有 JPL星历表文件结构示意图 可以让我学习学习?

谢谢!

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
xjw01
  2楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3804 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2009/4/27 18:24:00

结合上次发的c++代码看
图片点击可在新窗口打开查看此主题相关图片如下:结构图.png
图片点击可在新窗口打开查看

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
美女呀,离线,留言给我吧!
蔡越
  3楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
  发帖心情 Post By:2009/4/27 19:06:00

谢谢xjw01, 虽然我还是看不懂这个1 和475究竟是什么, 但毕竟终于有了一张图可以研究了, 谢谢!

1 是金星 那么 0 在哪里? 12是天平动, 13以后是什么? 13以后是其他行星吗?


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
美女呀,离线,留言给我吧!
蔡越
  4楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
  发帖心情 Post By:2009/4/27 20:03:00

以下 header.405

KSIZE=  2036    NCOEFF=  1018



GROUP  1010



JPL Planetary Ephemeris DE405/DE405
Start Epoch: JED=  2305424.5 1599 DEC 09 00:00:00
Final Epoch: JED=  2525008.5 2201 FEB 20 00:00:00



GROUP  1030



  2305424.50  2525008.50        32.



GROUP  1040



  156
  DENUM  LENUM  TDATEF  TDATEB  CENTER  CLIGHT  AU      EMRAT  GM1    GM2 
  GMB    GM4    GM5    GM6    GM7    GM8    GM9    GMS    RAD1    RAD2 
  RAD4    JDEPOC  X1      Y1      Z1      XD1    YD1    ZD1    X2      Y2   
  Z2      XD2    YD2    ZD2    XB      YB      ZB      XDB    YDB    ZDB 
  X4      Y4      Z4      XD4    YD4    ZD4    X5      Y5      Z5      XD5 
  YD5    ZD5    X6      Y6      Z6      XD6    YD6    ZD6    X7      Y7   
  Z7      XD7    YD7    ZD7    X8      Y8      Z8      XD8    YD8    ZD8 
  X9      Y9      Z9      XD9    YD9    ZD9    XM      YM      ZM      XDM 
  YDM    ZDM    XS      YS      ZS      XDS    YDS    ZDS    BETA    GAMMA
  J2SUN  GDOT    MA0001  MA0002  MA0004  MAD1    MAD2    MAD3    RE      ASUN 
  PHI    THT    PSI    OMEGAX  OMEGAY  OMEGAZ  AM      J2M    J3M    J4M 
  C22M    C31M    C32M    C33M    S31M    S32M    S33M    C41M    C42M    C43M 
  C44M    S41M    S42M    S43M    S44M    LBET    LGAM    K2M    TAUM    AE   
  J2E    J3E    J4E    K2E0    K2E1    K2E2    TAUE0  TAUE1  TAUE2  DROTEX
  DROTEY  GMAST1  GMAST2  GMAST3  KVC    IFAC    PHIC    THTC    PSIC    OMGCX
  OMGCY  OMGCZ  PSIDOT  MGMIS  ROTEX  ROTEY                               



GROUP  1041



  156
  0.405000000000000000D+03  0.405000000000000000D+03  0.000000000000000000D+00
  0.119970525194723000D+17  0.000000000000000000D+00  0.299792457999999984D+06
.........
  0.646682543384255465D-13  0.127748118910414607D-13  0.333405877296029502D-14
  0.000000000000000000D+00  0.299999999999999974D-03 -0.425951830000000000D-02
  0.408844299999999994D+00 -0.171450900000000006D+01  0.000000000000000000D+00
-0.158167070000000005D-05  0.229888000000000009D+00  0.000000000000000000D+00
  0.100000000000000000D+01  0.000000000000000000D+00  0.000000000000000000D+00



GROUP  1050



    3  171  231  309  342  366  387  405  423  441  753  819  899
    14    10    13    11    8    7    6    6    6    13    11    10    10
    4    2    2    1    1    1    1    1    1    8    2    4    4



GROUP  1070



==============================================================
以下ascp2000.405


    1  1018                                                                   
  0.245153650000000000D+07  0.245156850000000000D+07 -0.338008787742210925D+08 
  0.113078768223014772D+08  0.397860220581820933D+06 -0.255631650982068168D+05 
...............
  0.329526795987821133D-06 -0.361299724751331767D-07 -0.101434726607598293D-07 
  0.148495325313020708D-08  0.000000000000000000D+00  0.000000000000000000D+00 
    2  1018                                                                   
  0.245156850000000000D+07  0.245160050000000000D+07  0.512033995950585008D+08 
  0.200245113866126607D+07 -0.115597941240896890D+07 -0.422458541757940184D+05
................
  229  1018                                                                   
  0.245883250000000000D+07  0.245886450000000000D+07 -0.468225142464447618D+08 
  0.855287673857185431D+07  0.612484375662173959D+06 -0.271197032404459242D+05 
............
  0.171576773599179644D-06 -0.401686155703562924D-07 -0.519119817867234511D-08 
  0.105036480471432147D-08  0.000000000000000000D+00  0.000000000000000000D+00 

还是搞不清楚......

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
xjw01
  5楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3804 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2009/4/27 22:16:00

KSIZE= 2036 NCOEFF= 1018

nCOEFF为本表系数的个数

KSIZE=2倍nCOEFF

====================

JPL Planetary Ephemeris DE405/DE405
Start Epoch: JED= 2305424.5 1599 DEC 09 00:00:00
Final Epoch: JED= 2525008.5 2201 FEB 20 00:00:00

指DE405/LE405适用历元

======================

2305424.50 2525008.50 32.

指适用历元,32指每小段系数适用的天数,一个小段可能有14个或11个或6个数据等

=======================

156
DENUM LENUM TDATEF TDATEB CENTER CLIGHT AU EMRAT GM1 GM2
GMB GM4 GM5 GM6 GM7 GM8 GM9 GMS RAD1 RAD2
RAD4 JDEPOC X1 Y1 Z1 XD1 YD1 ZD1 X2 Y2
Z2 XD2 YD2 ZD2 XB YB ZB XDB YDB ZDB
X4 Y4 Z4 XD4 YD4 ZD4 X5 Y5 Z5 XD5
……

其中156指数积分常数有156个,包括每个星体的初始坐标、初始速度、天文学数等。这些常数在进行切比雪夫计算时基本不用,只用了两三个,如光速、AU长度等

DENUM LENUM……指常数名称,具体含义在说明书中

==========================

156
0.405000000000000000D+03 0.405000000000000000D+03 0.000000000000000000D+00
0.119970525194723000D+17 0.000000000000000000D+00 0.299792457999999984D+06

……

指常数值

==========================

GROUP 1050

3 171 231 309 342 366 387 405 423 441 753 819 899
14 10 13 11 8 7 6 6 6 13 11 10 10
4 2 2 1 1 1 1 1 1 8 2 4 4

非常重要。指每个星体所用系数表的位置索引

3 水星系数位置

171 金星系数位置

231 地球系数位置

……

最后一个星体(其实不是,是天平动)位置在899

14指水星系数个数(每32/4天用了14个系数),32天分4段

10指金星系数个数(每32/2天用了10个系数),32天分2段

13指地球星系数个数(每32/2天用了13个系数)

11指火星系数个数(每32/1天用了11个系数)

如果系数个数为0,说明该版本的DE星历表不含该星体的系数表

……

其实表中数字关系为,171=3+14*4*3=171,因为有三个坐标所以乘3,不过应注意,章动的坐标数是2,其它是3

所有系数加起来之后再加上2就是1018

 

====================
1 1018
0.245153650000000000D+07 0.245156850000000000D+07 -0.338008787742210925D+08
0.113078768223014772D+08 0.397860220581820933D+06 -0.255631650982068168D+05
……

1表示第1块,1018表示第1块有1018个系数

0.245153650000000000D+07 0.245156850000000000D+07表示本块的时间范围

即2451536.5儒略日到2451568.5儒略日,二者相减正好是32天

====================

-0.338008787742210925D+08
……

金星数据从-0.338008787742210925D+08 开始


 具体细节详节源代码,我没有找到JPL对数据表的详细解释,也是他们的C++源代码阅读得到的。他们提供的源代码很长很长,看了简值要吐血。可以这么说,写这段代码的人程序设计功力不太好!以至于我们阅读起来很麻烦。必要的话也可以参考我的代码,不过我的注释太少。

 有了DE405,并不是说就可以算出与别人相同的数据,你最好用DE405自带的校验数据进行校验,确保无误。之后,你还要理解很多坐标变换的东西,才可能得到与《天文年历》或JPL网站或swiss星历表一模一样的星历。

DE405至少可以压缩20——80倍,不然太大。主要方法是重新拟合,这需要了解切比雪夫多项式的性质。

DE406是DE405的压缩版,适用范围也大一些,但还是可以压缩20倍以上,原来有600多M,压缩后可以达至30M左右或更少

 

 


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
美女呀,离线,留言给我吧!
蔡越
  6楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
  发帖心情 Post By:2009/4/27 22:30:00

敬礼! 感谢xjw01! 这样明白一些了.



支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
xjw01
  7楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:蝙蝠侠 帖子:866 积分:3804 威望:3 精华:20 注册:2008/3/20 22:14:00
  发帖心情 Post By:2009/4/27 22:39:00

还是使用VSOP87吧,精度高,数据量小。压缩后只用几十KB,少量修正后精度确高达0.1至0.01角秒(具体看压缩量),不比《天文年历》差多少。

VSOP87本身拟合了DE200


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
美女呀,离线,留言给我吧!
蔡越
  8楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
  发帖心情 Post By:2009/4/27 22:45:00

以下是引用xjw01在2009-4-27 22:39:00的发言:

还是使用VSOP87吧,精度高,数据量小。压缩后只用几十KB,少量修正后精度确高达0.1至0.01角秒,不比《天文年历》差多少。


VSOP87本身拟合了DE200



VSOP87我找到的资料不全, 没有月亮和冥王星的轨道计算, 所以耽搁了好几年.


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
帅哥哟,离线,有人找我吗?
ymy111
  9楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游侠 帖子:190 积分:1585 威望:0 精华:1 注册:2008/7/6 15:34:00
  发帖心情 Post By:2009/4/29 10:51:00


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
美女呀,离线,留言给我吧!
蔡越
  10楼 个性首页 | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
  发帖心情 Post By:2009/5/4 4:41:00

以下是引用ymy111在2009-4-29 10:51:00的发言:

http://bbs.nongli.com/dispbbs_2_18614.html


这里应该有---你是否能改成js程序


http://www.namipan.com/downfile/VSOP87Functions.zip/bc473299b4eb015ea16ad4ad1a925381a2d5b33bba671200


这个事VSOP87源代码



第1个是北京天文同好的源码

第2个是Jay Tanner 杰伊唐纳的代码

杰伊唐纳放出来的代码就是缺了月亮的轨道计算.


[此贴子已经被作者于2009-5-4 5:53:50编辑过]

支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部

返回版面帖子列表

请教de406行星历表的结构








签名