以文本方式查看主题

-  中华农历论坛  (http://bbs.nongli.net/index.asp)
--  历法知识  (http://bbs.nongli.net/list.asp?boardid=2)
----  大家帮忙看这样求de405星历的结果是否正确? 怎么转换为日心黄经?  (http://bbs.nongli.net/dispbbs.asp?boardid=2&id=18840)
--  作者：蔡越
--  发布时间：2009/5/5 0:21:00
--  大家帮忙看这样求de405星历的结果是否正确? 怎么转换为日心黄经?
首先：
按照天文算法上的例子

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. 太阳系质心坐标 和 日心坐标 是否相同? 百度搜索 太阳系质心坐标 资料实在很少.

谢谢！
--  作者：ldlcau
--  发布时间：2009/5/5 8:31:00
--  

http://my.nthu.edu.tw/~res9202/intro/club/wi/can/commit/no5.pdf

 

http://www.astro.uio.no/~bgranslo/aares/calculate.html

 

 

 JPL行星精密历表的使用

http://astrodynamics.googlecode.com/files/JPL%E7%B2%BE%E5%AF%86%E5%8E%86%E8%A1%A8%E7%9A%84%E4%B8%8B%E8%BD%BD%E4%BD%BF%E7%94%A8%E8%AF%B4%E6%98%8E.pdf

 

 
--  作者：浪-淘-沙
--  发布时间：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
--  作者：蔡越
--  发布时间：2009/5/5 12:26:00
--  
以下是引用ldlcau在2009-5-5 8:31:00的发言：

http://my.nthu.edu.tw/~res9202/intro/club/wi/can/commit/no5.pdf


 


http://www.astro.uio.no/~bgranslo/aares/calculate.html


 


 


 JPL行星精密历表的使用


http://astrodynamics.googlecode.com/files/JPL%E7%B2%BE%E5%AF%86%E5%8E%86%E8%A1%A8%E7%9A%84%E4%B8%8B%E8%BD%BD%E4%BD%BF%E7%94%A8%E8%AF%B4%E6%98%8E.pdf


 


 



以下是引用浪-淘-沙在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

感谢!
--  作者：浪-淘-沙
--  发布时间：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编程语言,似乎用的人较少了.但楼主可以到专门的程序论坛里找同好一块研究啊.
--  作者：蔡越
--  发布时间：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编程语言,似乎用的人较少了.但楼主可以到专门的程序论坛里找同好一块研究啊.



谢谢老师细心回复!
--  作者：浪-淘-沙
--  发布时间：2009/5/5 15:19:00
--  

楼主太谦虚了.

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

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

 

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

 

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

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

 

大家一块努力吧.

有问题大家一块来讨论解决.
--  作者：浪-淘-沙
--  发布时间：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朋友请教.
--  作者：蔡越
--  发布时间：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

再次感谢 浪-淘-沙 老师!
--  作者：xjw01
--  发布时间：2009/5/7 9:40:00
--  

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

1、坐标旋转、平移

2、岁差计算

3、章动、光行差计算

 

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

 

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

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