农历计算程序

zopey

1. Function getjq_12(yy, jda)
   
2. k0 = Round((moon_L(jda) - earth_L(jda)) / (pai * 2), 0)
   
3. 
   
4. jd1 = jda
   
5. W1 = (k0 - 1) * pai * 2
   
6. Do
   
7.    jd0 = jd1
   
8.    stDegree = moon_L(jd0) - earth_L(jd0) - W1
   
9.    stDegreep = (moon_L(jd0 + 0.000005) - earth_L(jd0 + 0.000005) - moon_L(jd0 - 0.000005) + earth_L(jd0 - 0.000005)) / 0.00001
   
10.    jd1 = jd0 - stDegree / stDegreep
    
11. Loop Until Abs(jd1 - jd0) < 0.0000001
    
12. getjq_12a = jd1 + 8 / 24 - deltatT(yy) / 86400
    
13. 
    
14. jd2 = jda
    
15. W2 = k0 * pai * 2
    
16. Do
    
17.    jd0 = jd2
    
18.    stDegree = moon_L(jd0) - earth_L(jd0) - W2
    
19.    stDegreep = (moon_L(jd0 + 0.000005) - earth_L(jd0 + 0.000005) - moon_L(jd0 - 0.000005) + earth_L(jd0 - 0.000005)) / 0.00001
    
20.    jd2 = jd0 - stDegree / stDegreep
    
21. Loop Until Abs(jd2 - jd0) < 0.0000001
    
22. getjq_12b = jd2 + 8 / 24 - deltatT(yy) / 86400
    
23. 
    
24. If Round(getjq_12b, 0) = Round(jda, 0) Then
    
25.    getjq_12 = getjq_12b
    
26. ElseIf getjq_12b <= jda Then
    
27.    getjq_12 = getjq_12b
    
28. ElseIf getjq_12a <= jda Then
    
29.    getjq_12 = getjq_12a
    
30. ElseIf Round(getjq_12a, 0) = Round(jda, 0) Then
    
31.    getjq_12 = getjq_12a
    
32. End If
    
33. 
    
34. End Function

复制代码

zopey

cyfx2288 发表于 2019-2-17 23:50
各种网站，我都查过，农历网、台湾中央研究院的网页。各版本之间总的算来部分年度相关一两天，部分年度润 ...

甲子纪日 是不会错的，初一、十五的日期也 不会差太多（观月可得），不同朝代有差别的就是 定年首 和定 闰月。1楼 V2.0版本 已经修复了常见的bug，按第一个“无中气”定闰月 比较准确了。

zopey

YZC51 发表于 2019-2-16 22:47
古代万年历节选
明庄烈帝崇祯十五年·清太宗崇德七年壬午(公元1642年)

V2.0  生成农历（含闰11月）与之 一致。

YZC51

zopey 发表于 2019-2-18 17:08
V2.0  生成农历（含闰11月）与之 一致。

非常佩服楼主，农历函数做得精益求精！
2891 R10
3988 R12

zopey

YZC51 发表于 2019-2-18 19:17
非常佩服楼主，农历函数做得精益求精！
2891 R10
3988 R12

代入2891年， 某月初一日期 计算结果 2777302.495，四舍五入后为   2777302。
与网页计算的 2777303 有一天误差，使置闰结果不同。 这个原因未知，无法修改。

部分代码 可优化：

If jda >= getjq_12b Or Round(getjq_12b, 0) = Round(jda, 0) Then
   getjq_12 = getjq_12b
ElseIf jda >= getjq_12a Or Round(getjq_12a, 0) = Round(jda, 0) Then
   getjq_12 = getjq_12a
End If

YZC51

不会，瞎搞的
If yy = 2891 Then yzc = 21 / 4320
getjq_12 = getjq_12 + yzc

If yy = 3988 Or yy = 3989 Then yzc = 43 / 4320
getjq_24 = getjq_24 - yzc

zopey

农历的编算和颁行》国家标准：

（1）以北京时间为标准时间；
     北京时间属于协调世界时，协调世界时以原子时为基准，但通过闰秒的方式与世界时时刻的偏差不超过0.9秒，因此计算精度不计及编算时尚未正式发布的闰秒。

（2）朔日为农历月的第一个农历日；
     为了保证农历日期编排结果的唯一性，朔和节气时间的计算必须精确到在0时附近也能准确判断它们所在的农历日，在用于日期判定时对它们的计算精度要求应达到1秒。 （达不到这个精度）

（3）包含节气冬至在内的农历月为农历十一月；

（4）若从某个农历十一月开始到下一个农历十一月（不含）之间有13个农历月，则需要置闰。置闰规则为：取其中最先出现的一个不包含中气的农历月为农历闰月。

（5）农历十一月之后第2个（不计闰月）农历月为农历年的起始月。

YZC51

zopey 发表于 2019-2-19 09:45
农历的编算和颁行》国家标准：

（1）以北京时间为标准时间；

几千年到1万多年的时间跨度，精度很难达到1秒！
楼主能做到这样，就很不错啦！

zopey

中国古代 也没有“北京时间“的 说法吧，各个地方时不同，那农历初一是怎么 精确定位的呢？

WYS67

zopey 发表于 2019-2-16 12:44
http://club.excelhome.net/thread-1457643-1-1.html  3楼 计算初一（合塑时刻）的代码，应该是正确的。
...

四柱预测更重视农历节气的安排！它是严格按照立春时间为新年开始的时间--进入寅月。然后是春夏秋冬四季轮回。