de406行星历表的结构

以下 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

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

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左右或更少