24节气公式

zopey · 发表于 2018-1-8 11:04

Public Function getjq(yy As Integer, mm As Integer)
v0 = 628.3319653318
'第1步迭代
t = 0
L0 = (48650621.66 + 6283319653.318 * t) / 10 ^ 7
'W指的是太阳黄经。1999年春分对应W=0，以后每W每增加15度对应下一个节气。
W = (mm - 5 + (yy - 1999) * 24) * 15 * 3.1415926 / 180
'第2步迭代
t = t + (W - L0) / v0
t2 = t * t
l1 = (48950621.66 + 6283319653.318 * t + 53 * t2 _
+ 334116 * Cos(4.67 + 628.307585 * t) + 2061 * Cos(2.678 + 628.3076 * t) * t) / 10 ^ 7
v1 = 628.332 + 21 * Sin(1.527 + 628.307585 * t)
'第3步迭代
t = t + (W - l1) / v1
t2 = t * t
t3 = t2 * t
t4 = t3 * t
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)) / 10 ^ 7
'第4步迭代
t = t + (W - L2) / v1
J2000 = 2451545
'地球自转修正项需完善
JD = J2000 + t * 36525 - (64.7 + (yy - 2005) * 0.4) / 86400 + 8 / 24
'转换日期
Z = Int(JD + 0.5)
F = JD + 0.5 - Z
a0 = Int((Z - 1867216.25) / 36524.25)
If Z < 2299161 Then
A = Z
Else
A = Z + 1 + a0 - Int(a0 / 4)
End If
B = A + 1524
C = Int((B - 122.1) / 365.25)
d = Int(365.25 * C)
E = Int((B - d) / 30.6001)
''''
d1 = B - d - Int(30.6001 * E) + F
If E < 14 Then
m1 = E - 1
Else
m1 = E - 13
End If
If m1 > 2 Then
y1 = C - 4716
Else
y1 = C - 4715
End If
d2 = Int((d1 - Int(d1)) * 86400)
hh1 = Int(d2 / 3600)
mm1 = Round(((d2 - hh1 * 3600) / 60), 0)
'ss1 = d2 - hh1 * 3600 - mm1 * 60
If mm1 < 10 Then mm1 = "0" & mm1
getjq = y1 & "-" & m1 & "-" & Int(d1) & " " & hh1 & ":" & mm1
End Function

复制代码

计算节气时间是制作一份合格日历的必要条件，找到如上代码，简单调试后可以正常运行了。
自定义函数 =getjq (年份，节气序号) ，精度为分钟，详见附件。

24节气公式.zip (18.94 KB, 下载次数: 906)

YZC51 · 发表于 2018-1-8 12:55

本帖最后由 YZC51 于 2018-1-8 14:58 编辑

优秀作品，简单实用！1908到2090年期间有效！

ivccav · 发表于 2018-1-8 13:14

计算得精确没有意义，权威才是根本。据说寿星天文历比日梭万年历准确，但是节气是按日梭的，人家是接近紫金天文台的。24节气计算差1秒，可能就差一天！还是用天文台的数据查表吧，因为所有现行日历都以之为标准。天文台的数据到2100年，够用了。

zopey · 发表于 2018-1-8 13:48

YZC51 发表于 2018-1-8 12:55
优秀作品，简单实用！1098到2090年期间有效！

你收藏的excel函数版的节气公式，下载看了，看不懂（公式中很多的定义名称），估计是二分逼近法。

数学上通常采用的迭代求解方法有二分逼近法和牛顿迭代法。

zopey · 发表于 2018-1-8 14:27

1908、2090 报错是什么原因呢？它们的平均值是1999 年，而1999出现在如下代码

'W指的是太阳黄经。1999年春分对应W=0，以后每W每增加15度对应下一个节气。
W = (mm - 5 + (yy - 1999) * 24) * 15 * 3.1415926 / 180

zopey · 发表于 2018-1-8 15:01

试图在计算中使用一个数字，并且该数字被强制为一个整数，但是结果超过了整数的范围，例如：
Dim x As Long
x = 2000 * 365 ' 错误: 溢出
要解决这个问题，请按如下方式键入该数字：
Dim x As Long
x = CLng(2000) * 365

根据以上帮助文件，修改变量 yy为 clng(yy), 溢出问题解决。

YZC51 · 发表于 2018-1-8 17:21

zopey 发表于 2018-1-8 13:48
你收藏的excel函数版的节气公式，下载看了，看不懂（公式中很多的定义名称），估计是二分逼近法。

数 ...

收藏的，供参考

可以汇总的24节气时刻计算30节气-.rar (75.8 KB, 下载次数: 474)

YZC51 · 发表于 2018-1-8 19:40

zopey 发表于 2018-1-8 15:01
试图在计算中使用一个数字，并且该数字被强制为一个整数，但是结果超过了整数的范围，例如：
Dim x A ...

谢谢老师分享！
这样也可以
Public Function getjq(yy, mm As Integer)

YZC51 · 发表于 2018-1-9 07:52

zopey 发表于 2018-1-8 15:01
试图在计算中使用一个数字，并且该数字被强制为一个整数，但是结果超过了整数的范围，例如：
Dim x A ...

向老师学习，在原来的基础上的文件

24节气自定义函数.rar (35.24 KB, 下载次数: 437)

02761752696 · 发表于 2018-1-9 08:27

学习一下谢谢

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

[分享] 24节气公式

评分

评分