中华农历论坛知识讨论区历法知识 → [求助] 《天文计算》第49章 月亮的升降交点 误差


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

主题:[求助] 《天文计算》第49章 月亮的升降交点 误差

帅哥哟,离线,有人找我吗?
shtugis
  1楼 个性首页 | QQ | 信息 | 搜索 | 邮箱 | 主页 | UC


加好友 发短信
等级:新手上路 帖子:2 积分:223 威望:0 精华:0 注册:2010/6/17 10:10:00
[求助] 《天文计算》第49章 月亮的升降交点 误差  发帖心情 Post By:2012/11/29 10:26:00

大家好!我在测试《天文计算》第49章  月亮的升降交点    p260页的习题(49.a  计算1987年5月,月亮经过升交点的时刻)时,我将该书中的计算方法整理如下:
    year=1987.37
    #k=(year-2000.05)*13.4223
    k=-170   
    T=k/1342.23
    D=183.6380+331.73735691*k+0.0015057*T**2+0.00000209*T**3-0.000000010*T**4
    M=17.4006+26.82037250*k+0.0000999*T**2+0.00000006*T**3
    M1=38.3776+355.52747322*k+0.0123577*T**2+0.000014628*T**3-0.000000069*T**4   
    O=123.9767-1.44098949*k+0.0020625*T**2+0.00000214*T**3-0.000000016*T**4   
    V=299.75+132.85*T-0.009173*T**2   
    P=O+272.75-2.3*T   
    E=1-0.002516*T-0.0000074*T**2   
#以下6行把度数归到0-360
    D=D%360
    M=M%360
    M1=M1%360
    O=O%360
    V=V%360
    P=P%360 

    #角度中含M或-M时,还须乘上E  (原书中的提醒) 
    JDE=2451565.1619+27.212220817*k
    +0.0002572*T**2
    +0.000000021*T**3
    -0.000000000088*T**4
    -0.4721*sin(M1)
    -0.1649*sin(2*D)
    -0.0868*sin(2*D-M1)
    +0.0084*sin(2*D+M1)
    -E*0.0083*sin(2*D-M)
    -E*0.0039*sin(2*D-M-M1)
    +0.0034*sin(2*M1)
    -0.0031*sin(2*D-2*M1)
    +E*0.0030*sin(2*D+M)
    +E*0.0028*sin(M-M1)
    +E*0.0026*sin(M)
    +0.0025*sin(4*D)
    +0.0024*sin(D)
    +E*0.0022*sin(M+M1)
    +0.0017*sin(O)
    +0.0014*sin(4*D-M1)
    +E*0.0005*sin(2*D+M-M1)
    +E*0.0004*sin(2*D-M+M1)
    -(E**2)*0.0003*sin(2*D-2*M)
    +0.0003*E*sin(4*D-M)
    +0.0003*sin(V)
    +0.0003*sin(P)

        上述程序(python)计算结果:JDE = 2446939.08436
    《天文计算》中的结果:        JDE = 2446938.76803
    结果不一致,特此请教各位!
    我将程序与《天文计算》原书对照了N遍,仍然找不到问题所在!


    原书中的公式会不会在哪有问题?
[此贴子已经被作者于2012-11-29 13:38:34编辑过]


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


加好友 发短信
等级:版主 帖子:2058 积分:4273 威望:5 精华:4 注册:2008/11/13 21:03:00
  发帖心情 Post By:2012/11/29 12:04:00

《寿星天文历》计算的结果是2446938.7670833333与《天文计算》提供的答案基本一致。

 

与你计算的结果相差约7.5小时。

 

你要查一下,计算过程中,有没漏写了哪个参数,或者哪个参数输入时有笔误?


支持(0中立(0反对(0单帖管理 | 引用 | 回复 回到顶部

返回版面帖子列表

[求助] 《天文计算》第49章 月亮的升降交点 误差








签名