24节气公式

baigangliao · 发表于 2018-1-9 10:44

提示: 作者被禁止或删除内容自动屏蔽

zopey · 发表于 2018-1-11 09:11

'地球自转修正项需完善
JD = J2000 + t * 36525 - (64.7 + (yy - 2005) * 0.4) / 86400 + 8 / 24修改为
JD = J2000 + t * 36525 - deltaT(yy) / 86400 + 8 / 24

Public Function deltaT(year As Integer)
If year >= 2005 And year < 2014 Then
deltaT = 64.7 + (year - 2005) * 0.4
ElseIf year >= 2014 And year < 2114 Then
deltaT = -20 + 31 * ((year - 1820) / 100) ^ 2 + (year - 2114) * ((-20 + 31 * ((2014 - 1820) / 100) ^ 2) - (64.7 + (2014 - 2005) * 0.4)) / 100
ElseIf year > 2114 Then
deltaT = -20 + 31 * ((year - 1820) / 100) ^ 2
ElseIf year >= -4000 And year < -500 Then
y1 = -4000: y2 = -500
t = (year - y1) / (y2 - y1) * 10
A = 108371.7: B = -13036.8: C = 392: d = 0
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= -500 And year < -150 Then
y1 = -500: y2 = -150
t = (year - y1) / (y2 - y1) * 10
A = 17201: B = -627.82: C = 16.17: d = -0.3413
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= -150 And year < 150 Then
y1 = -150: y2 = 150
t = (year - y1) / (y2 - y1) * 10
A = 12200.6: B = -346.41: C = 5.403: d = -0.1593
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 150 And year < 500 Then
y1 = 150: y2 = 500
t = (year - y1) / (y2 - y1) * 10
A = 9113.8: B = -328.13: C = -1.647: d = 0.0377
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 500 And year < 900 Then
y1 = 500: y2 = 900
t = (year - y1) / (y2 - y1) * 10
A = 5707.5: B = -391.41: C = 0.915: d = 0.3145
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 900 And year < 1300 Then
y1 = 900: y2 = 1300
t = (year - y1) / (y2 - y1) * 10
A = 2203.4: B = -283.45: C = 13.034: d = -0.1778
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1300 And year < 1600 Then
y1 = 1300: y2 = 1600
t = (year - y1) / (y2 - y1) * 10
A = 490.1: B = -57.35: C = 2.085: d = -0.0072
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1600 And year < 1700 Then
y1 = 1600: y2 = 1700
t = (year - y1) / (y2 - y1) * 10
A = 120: B = -9.81: C = -1.532: d = 0.1403
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1700 And year < 1800 Then
y1 = 1700: y2 = 1800
t = (year - y1) / (y2 - y1) * 10
A = 10.2: B = -0.91: C = 0.51: d = -0.037
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1800 And year < 1830 Then
y1 = 1800: y2 = 1830
t = (year - y1) / (y2 - y1) * 10
A = 13.4: B = -0.72: C = 0.202: d = -0.0193
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1830 And year < 1860 Then
y1 = 1830: y2 = 1860
t = (year - y1) / (y2 - y1) * 10
A = 7.8: B = -1.81: C = 0.416: d = -0.247
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1860 And year < 1880 Then
y1 = 1860: y2 = 1880
t = (year - y1) / (y2 - y1) * 10
A = 8.3: B = -0.13: C = -0.406: d = 0.0292
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1880 And year < 1900 Then
y1 = 1880: y2 = 1900
t = (year - y1) / (y2 - y1) * 10
A = -5.4: B = 0.32: C = -0.183: d = 0.0173
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1900 And year < 1920 Then
y1 = 1900: y2 = 1920
t = (year - y1) / (y2 - y1) * 10
A = -2.3: B = 2.06: C = 0.169: d = -0.0135
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1920 And year < 1940 Then
y1 = 1920: y2 = 1940
t = (year - y1) / (y2 - y1) * 10
A = 21.2: B = 1.69: C = -0.304: d = 0.0167
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1940 And year < 1960 Then
y1 = 1940: y2 = 1960
t = (year - y1) / (y2 - y1) * 10
A = 24.2: B = 1.22: C = -0.064: d = 0.0031
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1960 And year < 1980 Then
y1 = 1960: y2 = 1980
t = (year - y1) / (y2 - y1) * 10
A = 33.2: B = 0.51: C = 0.231: d = -0.0109
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 1980 And year < 2000 Then
y1 = 1980: y2 = 2000
t = (year - y1) / (y2 - y1) * 10
A = 51: B = 1.29: C = -0.026: d = 0.0032
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
ElseIf year >= 2000 And year < 2005 Then
y1 = 2000: y2 = 2005
t = (year - y1) / (y2 - y1) * 10
A = 63.87: B = 0.1: C = 0: d = 0
deltaT = A + B * t + C * t ^ 2 + d * t ^ 3
Else
deltaT = 108372
End If
End Function

复制代码

zopey · 发表于 2018-1-11 09:23

比对了911查询网的节气数据，基本一致。

24节气公式(年份修订).rar (18.62 KB, 下载次数: 380)

circler · 发表于 2018-1-11 10:35

学习一下，谢谢！

YZC51 · 发表于 2018-1-12 10:17

zopey 发表于 2018-1-11 09:23
比对了911查询网的节气数据，基本一致。

谢谢老师分享！已经很完善啦！
白玉微瑕的是：
1900——2100间只有此一不足！
1923-2-20
应为
1923-2-19
还请老师修正为盼！谢谢！

zopey · 发表于 2018-1-12 11:36

这个问题很难说清雨水: 1923年02月20日00时00分16秒1923-2-19 23:59:40 计算误差总是有的，究竟谁对呢？

YZC51 · 发表于 2018-1-12 11:43

zopey 发表于 2018-1-12 11:36
这个问题很难说清雨水: 1923年02月20日00时00分16秒1923-2-19 23:59:40 计算误差总是有的，究竟谁对呢？ ...

老师说的有道理！
感觉还是以天文台为准好些！
请看附件

1923.rar (46.92 KB, 下载次数: 224)

大朋 · 发表于 2019-1-16 08:33

excel学习永无止境啊。谢谢啦。

zopey · 发表于 2019-1-16 15:02

'第1步迭代
t = 0
L0 = (48650621.66 + 6283319653.318 * t) / 10 ^ 7

'第2步迭代
t = t + (W - L0) / v0
t2 = t * t
l1 = (48950621.66 + 6283319653.318 * t + 53 * t2 _

===============

48650621.66 和 48950621.66 是同一个参数，所以这里有个笔误。
似乎对最后的迭代结果也没影响。

haidaoc · 发表于 2019-8-9 08:15

谢谢楼主分享！

		自动登录	找回密码
密码			免费注册

[分享] 24节气公式

评分

最新热点 /1

baigangliao baigangliao 当前离线性别保密 UID 1254158 阅读权限 0 注册时间 2010-11-9 在线时间小时最后登录 1970-1-1 积分 2437 财富分享技术精华头像被屏蔽 TA的精华主题 TA的得分主题	发表于 2018-1-9 10:44 \| 显示全部楼层 [广告] Excel易用宝 - 提升Excel的操作效率 · Excel / WPS表格插件 ★免费下载 ★ ★ 使用帮助★ 提示: 作者被禁止或删除内容自动屏蔽

	回复举报