星期、干支、二十八宿计算公式
(按阳历月日计算)
1. 求星期公式
星期=[5+A(实际天数)] mod 7
2. 干支计算公式
六十甲子干支序号,从1->59->0。
六十甲子干支序号=[23+A(实际天数)] mod 60
3. 二十八宿计算公式
二十八宿序号=[23+A(实际天数)] mod 28
4. 实际天数A的计算
A=B(基本天数)+C(闰日天数)
B=(计算年-1)*365+(要计算到年的月日天数)
例:1984年2月1日的基本天数B=(1984-1)*365+(31+1)=723827(天),其中,31是1月为31天,1为2月1日为1天。
公元308年8月28日的基本天数B=(308-1)*365+(31+28+31+30+31+30+31+27)=112055+239=112294(天)
这里的(要计算到年的月日天数),用的是公历,月日天数的规则我好象小学就学过了。哈哈……
C=(计算年-1) div 4 -误差修正值 + fixValue2
fixValue2为0或者1。常值为0,当年数为闰年(公历闰年法)之中的3月1日之后的为1。
误差修正值推算:
公元元年1月1日至1582年10月14日为0。
1582年10月15日至1699年12月31日为10。
从1701年1月1日起每增加一个世纪累加1,但能被400除尽的世纪不累加1。此方法推算即可。
——有一个问题,1700年这一年的修正值应为多少呢?算法中正好没有讲到,但看来应该是10。
例1701年1月1日起误差值为11,而1801年1月1日起误差修正值为12,而1901年1月1日起误差修正值为13,
但2001年误差修正值仍为13,因为2000年能被400整除,故不累加。而2101年1月1日起误差修正值为14。
5. 实例:1998.3.15的星期、干支与二十八宿
B=(1998-1)*365+(31+28+15)=728979
C=(1998-1) div 4 - 13 + 0 = 486
A=B+C=728979+486=729465
星期序号=(5+729465) mod 7=0,即为星期日
干支序号=(13+729465) mod 60=58,即为辛酉
二十八宿序号=(23+729465) mod 28=4,即为房
===================================================
好可怕!还有一些其它公式……但好象有些参数不知道怎么得到:
二十四节交节日算法:
用已知年的交接时辰加上22个小时35分,超过24要减去24,分数足60进1个小时,即得到8年后的各节交节时辰。
如2000年雨水交节时辰为16时22分,则2008年雨水交节时辰为14时52分。
因为16时22分+22时35分=38时57分。38-24=14时。
谁知道公元元年到公元八年的交节日,这个算法就可以实现了。——好象逆算法可以解决这个问题。谁试试?
农历闰月算法:
农历中,二十四节气(十二节气和十二中气)的中气落在月末的话,下个月就没有中气。农历将这种有节(节气)无气(中气)的月份规定为闰月。平均计算,19年有七个闰月。
但二十四个节气的十二节气和十二中气是怎么分的呢?我没有资料,估记应该是一节气一中气这样交叉。
====================================
啊,差不多了,好象弄明白就可以做万年历的算法程序了。
不过我没有时间,大家将代码弄出来了的话,别忘了给我一份儿。哈哈……