万年历的算法

摘自: 星期、干支、二十八宿计算公式打印本页关闭本窗口

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年有七个闰月。

但二十四个节气的十二节气和十二中气是怎么分的呢？我没有资料，

估记应该是一节气一中气这样交叉。:(

unit CNYear;

interface

uses sysutils;

type TCNDate = Cardinal;

function DecodeGregToCNDate(dtGreg:TDateTime):TCNDate;

function

GetGregDateFromCN(cnYear,cnMonth,cnDay:word;bLeap:Boolean=Fal se):TDateTime;

function GregDateToCNStr(dtGreg:TDateTime):String;

function isCNLeap(cnDate:TCNDate):boolean;

implementation

const cstDateOrg:Integer=32900; //公历1990-01-27的TDateTime表

示对应农历1990-01-01

const cstCNYearOrg=1990;

const cstCNTable:array[cstCNYearOrg..cstCNYearOrg + 60] of

WORD=( // unsigned 16-bit

24402, 3730, 3366, 13614, 2647, 35542, 858, 1749, //1997

23401, 1865, 1683, 19099, 1323, 2651, 10926, 1386, //2005

32213, 2980, 2889, 23891, 2709, 1325, 17757, 2741, //2013

39850, 1490, 3493, 61098, 3402, 3221, 19102, 1366, //2021

2773, 10970, 1746, 26469, 1829, 1611, 22103, 3243, //2029

1370, 13678, 2902, 48978, 2898, 2853, 60715, 2635, //2037