你好许老师！getCoordOne这个函数有两个问题请教您。
void getCoordOne(int a, int bCenter, double* r) { //在pv中提取一个星体的坐标
  int i,j;
  double *pa=pv[a], *pb=pv[bCenter];
  for (i=0; i<6; i++) {
    if(a==13||a==14) r[i] = pa[i]; // 提取章动(13)和天平动(14)
    else r[i] = pa[i]-pb[i];
  }
  if(a==14) r[4]=r[5]=0;
}
一是15被定义成月球的地心坐标，这里做坐标平移时并没有把15排除掉。
二是14是月球天平动，包括三个欧拉角及变化率，不知为何把14的后两个数置0;而13是章动，包括黄经与交角章动及变化率，应为4个数。
还有想请教：计算中首先算的是SSB下的星体位置，SSB是不是就是和J2000平行的，要得到J2000下的值是否还要通过所谓ICRF进行转换。

可能是我弄错了。也许应改为if(a==13) r[4]=r[5]=0;很久没有使用DE星表了，当时编写程序的过程忘了许多。

不过DE406中没有章动和天平动。所以我上次写的程序在DE406中是没有问题的，实际上我利用testpro.406验证过。

 

你说到：“SSB是不是就是和J2000平行的”，我不太明白您的意思。

 

确定一个坐标系，须知道三个坐标系的方向及坐标原点。SSB指的是坐标原点建在太阳系质心。J2000坐标系，指的是坐标轴的方向，它由历元J2000的黄道、赤道确定。所以SSB和J2000坐标之前不能比较。

感谢您的答复，看来不同的版本程序还是要稍作改动的。不知道DE200,或者更新的DE410是否也通用。
第二个我的意思是问DE405/6所求的坐标是在ICRF下的描述吗。SSB是指太阳系质心坐标系，是在卫星轨道控制领域里常用的一个标记。ICRF是ICRS(International Earth Rotation and Reference Systems Service)定义的坐标系。J2000.0自然很常用了。这些坐标系的转动关系又是如何呢。

ICRF不是基于黄道赤道定义的坐标系，它是利用遥远的星系定义的坐标系。可以认为它是一个无旋的坐标系，即随时间推移，坐标轴的方向可以认为没有变化。
ICRF与J2000赤道坐标系还是很接近的，实际上它于FK5更接近。

在直角坐标系中，ICRF旋转到J2000坐标系，可以通过旋转矩阵来解决问题，不过我没有找到这个矩阵，也许使用FK5的那个就可以了。我有一个球面旋转使用的参数，利用它也可以得到直角坐标旋转用的矩阵。通过这个矩阵旋转后，与NASA星历表比对，也没有明显误差，误差大约在毫角秒级以下。

旋转方法参考如下：

function X406to2000(a){
 var J=new Array(
   1.0000000000000,-0.0000000742925,0.0000000734750,
   0.0000000389354, 0.9174821361227,0.3977769851257,
  -0.0000000969638,-0.3977769851257,0.9174821361227);
 var b=new Array();
 b[0] = a[0]*J[0] + a[1]*J[1] + a[2]*J[2];
 b[1] = a[0]*J[3] + a[1]*J[4] + a[2]*J[5];
 b[2] = a[0]*J[6] + a[1]*J[7] + a[2]*J[8];
 return b;
}
//如果不使用X406to2000，可在球面坐标中参考以下3行
//b[0]+=0.05028/rad;
//b=llrConv(b,-(23*3600+26*60+21.4096)/rad);
//b[0]-=0.0381/rad;

谢谢楼主的帖子！有参考价值。

严重同意.......谢谢。

许老师，我在VC6.0中实现了您给的星历算法，可程序在链接过程中总出错：
Linking...
sxwnl_eph.obj : error LNK2005: "void __cdecl d2e(char *,int)" (?d2e@@YAXPADH@Z) already defined in EphemUtil.obj
sxwnl_eph.obj : error LNK2005: "void __cdecl gotoNextGroup(char *,struct _iobuf *)" (?gotoNextGroup@@YAXPADPAU_iobuf@@@Z) already defined in EphemUtil.obj.......

不知道是什么原因，您要是有时间的话帮我看一下吧！
谢谢了~

希望许老师能写一个基于瑞士星历表的实用软件

精品好贴，向上顶。同时再次谢谢许兄的奉献和劳动。