星期、天干地支、二十八星宿計算公式

bestlong

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年有七个闰月。
但二十四个节气的十二节气和十二中气是怎么分的呢？我没有资料，估记应该是一节气一中气这样交叉。

1. unit CNYear;
   
2. 
   
3. interface
   
4. uses sysutils;
   
5.   type TCNDate = Cardinal;
   
6.   function DecodeGregToCNDate(dtGreg:TDateTime):TCNDate;
   
7.   function GetGregDateFromCN(cnYear,cnMonth,cnDay:word;bLeap:Boolean=False):TDateTime;
   
8.   function GregDateToCNStr(dtGreg:TDateTime):String;
   
9.   function isCNLeap(cnDate:TCNDate):boolean;
   
10. implementation
    
11. 
    
12. const cstDateOrg:Integer=32900; //公历1990-01-27的TDateTime表示 对应农历1990-01-01
    
13. const cstCNYearOrg=1990;
    
14. const cstCNTable:array[cstCNYearOrg..cstCNYearOrg + 60] of WORD = ( // unsigned 16-bit
    
15. 24402, 3730, 3366, 13614, 2647, 35542, 858, 1749, //1997
    
16. 23401, 1865, 1683, 19099, 1323, 2651, 10926, 1386, //2005
    
17. 32213, 2980, 2889, 23891, 2709, 1325, 17757, 2741, //2013
    
18. 39850, 1490, 3493, 61098, 3402, 3221, 19102, 1366, //2021
    
19. 2773, 10970, 1746, 26469, 1829, 1611, 22103, 3243, //2029
    
20. 1370, 13678, 2902, 48978, 2898, 2853, 60715, 2635, //2037
    
21. 1195, 21179, 1453, 2922, 11690, 3474, 32421, 3365, //2045
    
22. 2645, 55901, 1206, 1461, 14038); //2050
    
23. //建表方法：
    
24. // 0101 111101010010 高四位是闰月位置,后12位表示大小月,大月30天,小月29天,
    
25. //闰月一般算小月,但是有三个特例2017/06,2036/06,2047/05
    
26. //对于特例则高四位的闰月位置表示法中的最高为设置为1 特殊处理用wLeapNormal变量
    
27. //
    
28. //2017/06 28330->61098 2036/06 27947->60715 2047/05 23133->55901
    
29. 
    
30. //如果希望用汇编,这里有一条信息:农历不会滞后公历2个月.
    
31. //将公历转换为农历
    
32. //返回:12位年份+4位月份+5位日期
    
33. function DecodeGregToCNDate(dtGreg:TDateTime):TCNDate;
    
34. var
    
35.   iDayLeave:Integer;
    
36.   wYear,wMonth,wDay:WORD;
    
37.   i,j:integer;
    
38.   wBigSmallDist,wLeap,wCount,wLeapShift:WORD;
    
39.   label OK;
    
40. begin
    
41.   result := 0;
    
42.   iDayLeave := Trunc(dtGreg) - cstDateOrg;
    
43.   DecodeDate(IncMonth(dtGreg,-1),wYear,wMonth,wDay);
    
44.   if (iDayLeave < 0) or (iDayLeave > 22295 )then Exit;
    
45. //Raise Exception.Create('目前只能算1990-01-27以后的');
    
46. //Raise Exception.Create('目前只能算2051-02-11以前的');
    
47.   for i:=Low(cstCNTable) to High(cstCNTable) do
    
48.   begin
    
49.     wBigSmallDist := cstCNTable[i];
    
50.     wLeap := wBigSmallDist shr 12;
    
51.     if wLeap > 12 then
    
52.     begin
    
53.       wLeap := wLeap and 7;
    
54.       wLeapShift := 1;
    
55.     end else
    
56.       wLeapShift := 0;
    
57.     for j:=1 to 12 do
    
58.     begin
    
59.       wCount:=(wBigSmallDist and 1) + 29;
    
60.       if j=wLeap then
    
61.         wCount := wCount - wLeapShift;
    
62.       if iDayLeave < wCount then
    
63.       begin
    
64.         Result := (i shl 9) + (j shl 5) + iDayLeave + 1;
    
65.         Exit;
    
66.       end;
    
67.       iDayLeave := iDayLeave - wCount;
    
68.       if j=wLeap then
    
69.       begin
    
70.         wCount:=29 + wLeapShift;
    
71.         if iDayLeave < wCount then
    
72.         begin
    
73.           Result := (i shl 9) + (j shl 5) + iDayLeave + 1 + (1 shl 21);
    
74.           Exit;
    
75.         end;
    
76.         iDayLeave := iDayLeave - wCount;
    
77.       end;
    
78.       wBigSmallDist := wBigSmallDist shr 1;
    
79.     end;
    
80.   end;
    
81.   //返回值:
    
82.   // 1位闰月标志 + 12位年份+4位月份+5位日期 (共22位)
    
83. end;
    
84. 
    
85. function isCNLeap(cnDate:TCNDate):boolean;
    
86. begin
    
87.   result := (cnDate and $200000) <> 0;
    
88. end;
    
89. function GetGregDateFromCN(cnYear,cnMonth,cnDay:word;bLeap:Boolean=False):TDateTime;
    
90. var
    
91.   i,j:integer;
    
92.   DayCount:integer;
    
93.   wBigSmallDist,wLeap,wLeapShift:WORD;
    
94. begin
    
95.   // 0101 010010101111 高四位是闰月位置,后12位表示大小月,大月30天,小月29天,
    
96.   DayCount := 0;
    
97.   if (cnYear < 1990) or (cnYear >2050) then
    
98.   begin
    
99.     Result := 0;
    
100.     Exit;
     
101.   end;
     
102.   for i:= cstCNYearOrg to cnYear-1 do
     
103.   begin
     
104.     wBigSmallDist := cstCNTable[i];
     
105.     if (wBIgSmallDist and $F000) <> 0 then
     
106.       DayCount := DayCount + 29;
     
107.     DayCount := DayCount + 12 * 29;
     
108.     for j:= 1 to 12 do
     
109.     begin
     
110.       DayCount := DayCount + wBigSmallDist and 1;
     
111.       wBigSmallDist := wBigSmallDist shr 1;
     
112.     end;
     
113.   end;
     
114.   wBigSmallDist := cstCNTable[cnYear];
     
115.   wLeap := wBigSmallDist shr 12;
     
116.   if wLeap > 12 then
     
117.   begin
     
118.     wLeap := wLeap and 7;
     
119.     wLeapShift := 1; //大月在闰月.
     
120.   end else
     
121.     wLeapShift := 0;
     
122.   for j:= 1 to cnMonth-1 do
     
123.   begin
     
124.     DayCount:=DayCount + (wBigSmallDist and 1) + 29;
     
125.     if j = wLeap then
     
126.       DayCount := DayCount + 29;
     
127.     wBigSmallDist := wBigSmallDist shr 1;
     
128.   end;
     
129.   if bLeap and (cnMonth = wLeap) then //是要闰月的吗?
     
130.     DayCount := DayCount + 30 - wLeapShift;
     
131.   result := cstDateOrg + DayCount + cnDay - 1;
     
132. end;
     
133. 
     
134. //将日期显示成农历字符串.
     
135. function GregDateToCNStr(dtGreg:TDateTime):String;
     
136. const hzNumber:array[0..10] of string=('零','一','二','三','四','五','六','七','八','九','十');
     
137. 
     
138.     function ConvertYMD(Number:Word;YMD:Word):string;
     
139.     var
     
140.       wTmp:word;
     
141.     begin
     
142.       result := '';
     
143.       if YMD = 1 then
     
144.       begin //年份
     
145.         while Number > 0 do
     
146.         begin
     
147.           result := hzNumber[Number Mod 10] + result;
     
148.           Number := Number DIV 10;
     
149.         end;
     
150.         Exit;
     
151.     end;
     
152.     if Number<=10 then
     
153.     begin //可只用1位
     
154.       if YMD = 2 then //月份
     
155.         result := hzNumber[Number]
     
156.       else //天
     
157.         result := '初' + hzNumber[Number];
     
158.       Exit;
     
159.     end;
     
160.     wTmp := Number Mod 10; //个位
     
161.     if wTmp <> 0 then result := hzNumber[wTmp];
     
162.     wTmp := Number Div 10; //十位
     
163.     result := '十'+result;
     
164.     if wTmp > 1 then result := hzNumber[wTmp] + result;
     
165.   end;
     
166.   
     
167. var
     
168.   cnYear,cnMonth,cnDay:word;
     
169.   cnDate:TCNDate;
     
170.   strLeap:string;
     
171. begin
     
172.   cnDate:= DecodeGregToCNDate(dtGreg);
     
173.   if cnDate = 0 then
     
174.   begin
     
175.     result := '输入越界';
     
176.     Exit;
     
177.   end;
     
178.   cnDay := cnDate and $1F;
     
179.   cnMonth := (cnDate shr 5) and $F;
     
180.   cnYear := (cnDate shr 9) and $FFF;
     
181.   //测试第22位,为1表示闰月
     
182.   if isCNLeap(cnDate) then strLeap:='(闰)' else strLeap := '';
     
183.   result := '农历' + ConvertYMD(cnYear,1) + '年' +
     
184.     ConvertYMD(cnMonth,2) + '月' + strLeap + ConvertYMD(cnDay,3) ;
     
185. end;
     
186. end.

複製代碼