- 中华农历论坛 (http://bbs.nongli.net/index.asp)
-- 历法知识 (http://bbs.nongli.net/list.asp?boardid=2)
---- 请教de406行星历表的结构 (http://bbs.nongli.net/dispbbs.asp?boardid=2&id=18710)
-- 发布时间：2009/4/27 13:34:00
-- 请教de406行星历表的结构
1 728
......
475 728
1是什么意思? 728又是什么意思?
哪位有 JPL星历表文件结构示意图 可以让我学习学习?
谢谢!
-- 发布时间：2009/4/27 18:24:00
--
结合上次发的c++代码看
此主题相关图片如下：结构图.png
-- 发布时间：2009/4/27 19:06:00
--
谢谢xjw01, 虽然我还是看不懂这个1 和475究竟是什么, 但毕竟终于有了一张图可以研究了, 谢谢!
1 是金星 那么 0 在哪里? 12是天平动, 13以后是什么? 13以后是其他行星吗?
-- 发布时间：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
还是搞不清楚......
-- 发布时间：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左右或更少
-- 发布时间：2009/4/27 22:30:00
--
敬礼! 感谢xjw01! 这样明白一些了.
-- 发布时间：2009/4/27 22:39:00
--
还是使用VSOP87吧，精度高，数据量小。压缩后只用几十KB，少量修正后精度确高达0.1至0.01角秒（具体看压缩量），不比《天文年历》差多少。
VSOP87本身拟合了DE200
-- 发布时间：2009/4/27 22:45:00
--
还是使用VSOP87吧，精度高，数据量小。压缩后只用几十KB，少量修正后精度确高达0.1至0.01角秒，不比《天文年历》差多少。
VSOP87本身拟合了DE200
VSOP87我找到的资料不全, 没有月亮和冥王星的轨道计算, 所以耽搁了好几年.
-- 发布时间：2009/4/29 10:51:00
--
http://bbs.nongli.com/dispbbs_2_18614.html
这里应该有---你是否能改成js程序
http://www.namipan.com/downfile/VSOP87Functions.zip/bc473299b4eb015ea16ad4ad1a925381a2d5b33bba671200
这个事VSOP87源代码
-- 发布时间：2009/5/4 4:41:00
--
http://bbs.nongli.com/dispbbs_2_18614.html
这里应该有---你是否能改成js程序
http://www.namipan.com/downfile/VSOP87Functions.zip/bc473299b4eb015ea16ad4ad1a925381a2d5b33bba671200
这个事VSOP87源代码
第1个是北京天文同好的源码
第2个是Jay Tanner 杰伊唐纳的代码
杰伊唐纳放出来的代码就是缺了月亮的轨道计算.