中华农历论坛知识讨论区历法知识 → 大家帮忙看这样求de405星历的结果是否正确? 怎么转换为日心黄经?


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

主题:大家帮忙看这样求de405星历的结果是否正确? 怎么转换为日心黄经?

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


加好友 发短信
等级:论坛游民 帖子:46 积分:535 威望:0 精华:0 注册:2009/4/27 13:24:00
大家帮忙看这样求de405星历的结果是否正确? 怎么转换为日心黄经?  发帖心情 Post By:2009/5/5 0:21:00

首先:
按照天文算法上的例子

1992年12月20日0时(力学时)
jd=2448976.5

金星的日心黄经:
    L = -68.6592582弧度 = -3933°.88572 = +26°.11428

金星的日心黄纬:
    B = -0.0457399(弧度) = -2°.62070,
距离:
R = 0.724603 (AU)

其次:
//--------------------------------
以下是通过喷气动力实验室的de405 中的2个文件 ascp1980.405 header.405计算出来的:

位置:x、y、z,单位:天文单位 速度:vx、vy、vz,单位:天文单位/天
  金星位置x:  0.65169409182828708  金星速度vx: -0.00899973297509975
  金星位置y: 0.31050718862610466  金星速度yx: 0.01626860234545932
  金星位置z: 0.09843463973061885  金星速度zx: 0.00788808009503498
  金星d      0.2888953              0.72856659394408141
  d = sqrt(x*x + y*y + z*z)
//================================================================

请教:

1. 上述位置:x、y、z 速度:vx、vy、vz 是否是准确的?
2. 利用太阳系质心坐标x、y、z vx、vy、vz 怎么变换为日心黄经?
  就是说 0.65169409182828708这些怎么变成 L = -68.6592582弧度 = -3933°.88572 = +26°.11428
3. 太阳系质心坐标 和 日心坐标 是否相同? 百度搜索 太阳系质心坐标 资料实在很少.

谢谢!



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


加好友 发短信
等级:论坛游侠 帖子:159 积分:583 威望:0 精华:1 注册:2008/9/3 16:55:00
  发帖心情 Post By:2009/5/5 8:31:00


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


加好友 发短信
等级:版主 帖子:2068 积分:4263 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2009/5/5 11:39:00

楼主计算的数据与我查JPL网站的数据应当是一致的.

我按太阳系质心为参考点查得的数据如下:

2448976.500000000 = A.D. 1992-Dec-20 00:00:00.0000 (CT)
 X = 6.516940918282871E-01 Y = 3.240398267452110E-01 Z =-3.320065014902096E-02
 VX=-8.999732975099746E-03 VY= 1.806384889286169E-02 VZ= 7.658936193954949E-04
 LT= 4.207849718743506E-03 RG= 7.285665939440813E-01 RR=-5.091436901910274E-05


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


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



以下是引用浪-淘-沙在2009-5-5 11:39:00的发言:

楼主计算的数据与我查JPL网站的数据应当是一致的.


我按太阳系质心为参考点查得的数据如下:


2448976.500000000 = A.D. 1992-Dec-20 00:00:00.0000 (CT)
 X = 6.516940918282871E-01 Y = 3.240398267452110E-01 Z =-3.320065014902096E-02
 VX=-8.999732975099746E-03 VY= 1.806384889286169E-02 VZ= 7.658936193954949E-04
 LT= 4.207849718743506E-03 RG= 7.285665939440813E-01 RR=-5.091436901910274E-05



谢谢热心的两位老师!

x、y、z vx、vy、vz LT RG RR 究竟代表什么意思呢? 读了提供的资料我还是很迷糊. 正在消化中.

是否有一C语言函数能够演示将 x、y、z vx、vy、vz 进行运算, 得到的结果是 L = -68.6592582弧度 = -3933°.88572 = +26°.11428 ?

论坛有没有既懂轨道运算 又是熟悉delphi编程语言的老师? 可否指点一下, 我的QQ是: 379873

感谢!


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


加好友 发短信
等级:版主 帖子:2068 积分:4263 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2009/5/5 12:57:00

Coordinate system description:

  Ecliptic and Mean Equinox of Reference Epoch

    Reference epoch: J2000.0
    xy-plane: plane of the Earth's orbit at the reference epoch
    x-axis  : out along ascending node of instantaneous plane of the Earth's
              orbit and the Earth's mean equator at the reference epoch
    z-axis  : perpendicular to the xy-plane in the directional (+ or -) sense
              of Earth's north pole at the reference epoch.

Symbol meaning [1 AU=149597870.691 km, 1 day=86400.0 s]:

    JDCT     Epoch Julian Date, Coordinate Time
      X      x-component of position vector (AU)                              
      Y      y-component of position vector (AU)                              
      Z      z-component of position vector (AU)                              
      VX     x-component of velocity vector (AU/day)                          
      VY     y-component of velocity vector (AU/day)                          
      VZ     z-component of velocity vector (AU/day)                          
      LT     One-way down-leg Newtonian light-time (day)  (就是光行时了)                    
      RG     Range; distance from coordinate center (AU)   (这个就是星体之间的距离)                    
      RR     Range-rate; radial velocity wrt coord. center (AU/day)           

States have been corrected for Newtonian light-time.
***************

 

编算星历表,请找"许剑伟"先生(xjw01),他是这个论坛的高人.

寿星万年历,就是一个例子,请楼主参考一下,是开放源程序的.

 

delphi编程语言,似乎用的人较少了.但楼主可以到专门的程序论坛里找同好一块研究啊.


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


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

以下是引用浪-淘-沙在2009-5-5 12:57:00的发言:

Coordinate system description:

  Ecliptic and Mean Equinox of Reference Epoch

    Reference epoch: J2000.0
    xy-plane: plane of the Earth's orbit at the reference epoch
    x-axis  : out along ascending node of instantaneous plane of the Earth's
              orbit and the Earth's mean equator at the reference epoch
    z-axis  : perpendicular to the xy-plane in the directional (+ or -) sense
              of Earth's north pole at the reference epoch.

Symbol meaning [1 AU=149597870.691 km, 1 day=86400.0 s]:

    JDCT     Epoch Julian Date, Coordinate Time
      X      x-component of position vector (AU)                              
      Y      y-component of position vector (AU)                              
      Z      z-component of position vector (AU)                              
      VX     x-component of velocity vector (AU/day)                          
      VY     y-component of velocity vector (AU/day)                          
      VZ     z-component of velocity vector (AU/day)                          
      LT     One-way down-leg Newtonian light-time (day)  (就是光行时了)                    
      RG     Range; distance from coordinate center (AU)   (这个就是星体之间的距离)                    
      RR     Range-rate; radial velocity wrt coord. center (AU/day)           

States have been corrected for Newtonian light-time.
***************


 


编算星历表,请找"许剑伟"先生(xjw01),他是这个论坛的高人.


寿星万年历,就是一个例子,请楼主参考一下,是开放源程序的.


 


delphi编程语言,似乎用的人较少了.但楼主可以到专门的程序论坛里找同好一块研究啊.



谢谢老师细心回复!


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


加好友 发短信
等级:版主 帖子:2068 积分:4263 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2009/5/5 15:19:00

楼主太谦虚了.

以后回复请尽量少用"老师"二个字.至少我不敢当.

大家在这里相聚,都算是朋友.

 

引用一句我从小说里看到的:"出家有早迟,得道无先后".(金庸是说:师弟比师兄得道更早.师弟就用这句话回复他,说师兄还是师兄,这个名份不必改的.)

 

也许我比你虚长几岁,有些知识我比你先学.

但这并不表示我比你有知识.后学先进的例子多了去.

 

大家一块努力吧.

有问题大家一块来讨论解决.


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


加好友 发短信
等级:版主 帖子:2068 积分:4263 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2009/5/5 18:02:00

以下是引用蔡越在2009-5-5 0:21:00的发言:

1. 上述位置:x、y、z 速度:vx、vy、vz 是否是准确的?
2. 利用太阳系质心坐标x、y、z vx、vy、vz 怎么变换为日心黄经?
就是说 0.65169409182828708这些怎么变成 L = -68.6592582弧度 = -3933°.88572 = +26°.11428

3. 太阳系质心坐标 和 日心坐标 是否相同? 百度搜索 太阳系质心坐标 资料实在很少.

谢谢!


关于坐标变成经度,很容易的:

有了XY坐标,用反正切函数就可以计算出角度.

按我查的JPL坐标数据,计算得到的角度为:

角度 X Y
26.21236794 0.6494062274379880 0.3197216206752680
26°12'44.52"

而实际上,这一时刻的黄经为26°.11428(26°6'51.40")

差了5'53"(353角秒)左右.

原因在哪?因为我查到的坐标是以JD2000.0的春分点为参考点的.

而26°6'51.40"这个数据是以当日春分点为参考点的.

按平均每年50.29"的春分西移来推算,1992年12月20日与2000年1月1日有7年多些,春分点的位置就差了352多角秒.

 

至于春分点具体的每年西移数值,这个计算公式,请楼主向XJW01朋友请教.


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


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

以下是引用浪-淘-沙在2009-5-5 18:02:00的发言:

关于坐标变成经度,很容易的:


有了XY坐标,用反正切函数就可以计算出角度.


按我查的JPL坐标数据,计算得到的角度为:




















角度XY
26.212367940.6494062274379880 0.3197216206752680
26°12'44.52"


而实际上,这一时刻的黄经为26°.11428(26°6'51.40")


差了5'53"(353角秒)左右.


原因在哪?因为我查到的坐标是以JD2000.0的春分点为参考点的.


而26°6'51.40"这个数据是以当日春分点为参考点的.


按平均每年50.29"的春分西移来推算,1992年12月20日与2000年1月1日有7年多些,春分点的位置就差了352多角秒.


 


至于春分点具体的每年西移数值,这个计算公式,请楼主向XJW01朋友请教.



谢谢!

感谢前辈的详细解释, 这样明白多一些了.

努力学习中, 我写的程序得到的结果已经是这样了, 看来不久就要打破几年来的裹步不前了.

太阳: a = -3967662.8282221 b = -134983002.5208531 c = -58525295.0310535
太阳弧转角 a = -227330334.58806377600 b = -7733956350.44863415000

水星: a = 3.5317554 b = 0.0537852 c = 0.4199397
水星弧转角 a = 202.35467768563 b = 3.08166313593

金星: a = 0.4574925 b = -0.0457297 c = 0.7246017
金星弧转角 a = 26.21238891682 b = -2.62011670415

火星: a = 1.72 b = 0.02 c = 1.59
火星弧转角 a = 98.80405888417 b = 1.40151419228

木星: a = 3.18 b = 0.02 c = 5.45
木星弧转角 a = 182.28348827101 b = 1.29135608234

土星: a = 5.57 b = -0.02 c = 9.87
土星弧转角 a = 319.41662394165 b = -1.08066958717

天王: a = 5.03 b = -0.01 c = 19.57
天王弧转角 a = 287.98912875096 b = -0.43151623040

海王: a = 5.04 b = 0.01 c = 30.19
海王弧转角 a = 288.63671853272 b = 0.69589289246

冥王: a = 4.07 b = 0.25 c = 29.72
冥王弧转角 a = 233.21289830268 b = 14.51415862193

月 : a = 3.83 b = -0.06 c = 0.00
月 弧转角 a = 219.46431177864 b = -3.54544680713

再次感谢 浪-淘-沙 老师!


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


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

你接下来要学习实践一下:

1、坐标旋转、平移

2、岁差计算

3、章动、光行差计算

 

计算出的数据如果与《寿星万年历》相同,那么在星历计算就成功一半多了。更精密的找JPL

 

之后就是星历算法的优化问题,这也是一个大难题。

再之后是天象的计算,这又是一个大难题。


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部
总数 13 1 2 下一页

返回版面帖子列表

大家帮忙看这样求de405星历的结果是否正确? 怎么转换为日心黄经?








签名