中华农历论坛知识讨论区历法知识 → 农历24节气计算的问题,请许剑伟先生不吝赐教!


  共有9311人关注过本帖平板打印

主题:农历24节气计算的问题,请许剑伟先生不吝赐教!

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


加好友 发短信
等级:新手上路 帖子:3 积分:244 威望:0 精华:0 注册:2021/7/21 17:17:00
农历24节气计算的问题,请许剑伟先生不吝赐教!  发帖心情 Post By:2021/8/9 17:55:00

 尊敬的许剑伟先生,在网上看了你的一篇《农历24节气算法》,我按照里面的方法试着算2007年的24节气,但怎么也算不出你的结果,现把春分和清明的计算结果展示如下,不知哪步错了,请指教。牛顿迭代算法步骤:第1步迭代:t = 0第2步迭代:t = t + ( W – L0(t) ) / v0第3步迭代:t = t + ( W – L1(t) ) / v1(t)第4步迭代:t = t + ( W – L2(t) ) / v1(t) (一)春分第1步迭代:W =0 t = 0第2步迭代:v0 = 628.3319653318L0 = 48950621.66 + 6283319653.318 * t= 48950621.66t = t + (W - L0) / v0 =0+(0-48950621.66)/628.3319653318=-77905.6676731926 第3步迭代:      L1 = (48950621.66 + 6283319653.318 * t + 53 * t * t + 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10000000 = -4891845.50859252       '弧度   V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =647.9743   t = t + (W - L1) /V1 =-77905.6676731926+(0+4891845.50859252)/647.9743=-70356.2240040676第4步迭代:    t2 = t * t: t3 = t2 * t: t4 = t3 * t t2=4949998256.11053  :t3=-348263186126657 :t4=2.45024827354973E+19 L2 = (48950621.66 + 6283319653.318 * t + 52.9674 * t2 + 0.00432 * t3 - 0.001124 * t4 + 334166 * Cos(4.669257 + 628.307585 * t) + 3489 * Cos(4.6261 + 1256.61517 * t) + 350 * Cos(2.744 + 575.3385 * t) + 342 * Cos(2.829 + 0.3523 * t) + 314 * Cos(3.628 + 7771.3771 * t) + 268 * Cos(4.418 + 786.0419 * t) + 234 * Cos(6.135 + 393.021 * t) + 132 * Cos(0.742 + 1150.677 * t) + 127 * Cos(2.037 + 52.9691 * t) + 120 * Cos(1.11 + 157.7344 * t) + 99 * Cos(5.23 + 588.493 * t) + 90 * Cos(2.05 + 2.63 * t) + 86 * Cos(3.51 + 39.815 * t) + 78 * Cos(1.18 + 522.369 * t) + 75 * Cos(2.53 + 550.755 * t) + 51 * Cos(4.58 + 1884.923 * t) + 49 * Cos(4.21 + 77.552 * t) + 36 * Cos(2.92 + 0.07 * t) + 32 * Cos(5.85 + 1179.063 * t) + 28 * Cos(1.9 + 79.63 * t) + 27 * Cos(0.31 + 1097.71 * t) + 2060.6 * Cos(2.67823 + 628.307585 * t) * t + 43 * Cos(2.635 + 1256.6152 * t) * t + 8.72 * Cos(1.072 + 628.3076 * t) * t2 - 994 - 834 * Sin(2.1824 - 33.75705 * t) - 64 * Sin(3.5069 + 1256.66393 * t)) / 10000000  =-2.798408E+08   V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t)=628.332 + 21 * Sin(1.527 + 628.307585 * (-70356.2240040676))=639.8343   t = t + (W - L2) / V1    =-70356.2240040676+(0-(-2.798408E+08))/639.8343= 367008.307245932JD = 2451545 + t * 36525 - TD_V / 86400 + 8 / 24     (2007年TD_V=65.5)  =2451545 +367008.307245932 * 36525 -65.5/86400 + 8 / 24=13407429967.4903春分日期 -0-12.99026    (实际春分 2007-03-21 08:06:59) (二)清明第1步迭代:W =15 t = 0第2步迭代:v0 = 628.3319653318L0 = 48950621.66 + 6283319653.318 * t= 48950621.66t = t + (W - L0) / v0 =0+(15-48950621.66)/628.3319653318=-77905.6438004618第3步迭代:      L1 = (48950621.66 + 6283319653.318 * t + 53 * t * t + 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10000000 = -4891841.99029775       '弧度   V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t) =618.2513   t = t + (W - L1) /V1 =-77905.6438004618+(15+4891841.99029775)/618.2513=-69993.2348256777第4步迭代:    t2 = t * t: t3 = t2 * t: t4 = t3 * t t2=4899052921.36246:t3=-342900561548345 :t4=2.40007195263101E+19   L2 = (48950621.66 + 6283319653.318 * t + 52.9674 * t2 + 0.00432 * t3 - 0.001124 * t4 + 334166 * Cos(4.669257 + 628.307585 * t) + 3489 * Cos(4.6261 + 1256.61517 * t) + 350 * Cos(2.744 + 575.3385 * t) + 342 * Cos(2.829 + 0.3523 * t) + 314 * Cos(3.628 + 7771.3771 * t) + 268 * Cos(4.418 + 786.0419 * t) + 234 * Cos(6.135 + 393.021 * t) + 132 * Cos(0.742 + 1150.677 * t) + 127 * Cos(2.037 + 52.9691 * t) + 120 * Cos(1.11 + 157.7344 * t) + 99 * Cos(5.23 + 588.493 * t) + 90 * Cos(2.05 + 2.63 * t) + 86 * Cos(3.51 + 39.815 * t) + 78 * Cos(1.18 + 522.369 * t) + 75 * Cos(2.53 + 550.755 * t) + 51 * Cos(4.58 + 1884.923 * t) + 49 * Cos(4.21 + 77.552 * t) + 36 * Cos(2.92 + 0.07 * t) + 32 * Cos(5.85 + 1179.063 * t) + 28 * Cos(1.9 + 79.63 * t) + 27 * Cos(0.31 + 1097.71 * t) + 2060.6 * Cos(2.67823 + 628.307585 * t) * t + 43 * Cos(2.635 + 1256.6152 * t) * t + 8.72 * Cos(1.072 + 628.3076 * t) * t2 - 994 - 834 * Sin(2.1824 - 33.75705 * t) - 64 * Sin(3.5069 + 1256.66393 * t)) / 10000000 =-2.741786E+08   V1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t)=628.332 + 21 * Sin(1.527 + 628.307585 * (-69993.2348256777))=608.6536   t = t + (W - L2) / V1    =-69993.2348256777+(15-(-2.741786E+08))/608.6536= 380474.202674322   JD = 2451545 + t * 36525 - TD_V / 86400 + 8 / 24    (2007年TD_V=65.5)  =2451545 +380474.202674322 * 36525 -65.5/86400 + 8 / 24=13899271798.0122 清明日期 --9-30.5122 (实际清明 2007-03-21 08:06:59) 算出的结果相差十万八千里,不知哪里错了,请许剑伟先生不吝赐教!在此谢谢了!

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

返回版面帖子列表

农历24节气计算的问题,请许剑伟先生不吝赐教!








签名