星期、干支、二十八宿計算公式 (2000年12月6日) 

本站更新  分類:   作者:來自未名碼頭  推薦:   閱讀次數:765  
(http://www.codesky.net)  

--------------------------------------------------------------------------------

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 
1195, 21179, 1453, 2922, 11690, 3474, 32421, 3365, //2045 
2645, 55901, 1206, 1461, 14038); //2050 
//建表方法： 
// 0101 111101010010 高四位是閏月位置,后12位表示大小月,大月30
天,小月29天, 
//閏月一般算小月,但是有三個特例2017/06,2036/06,2047/05 
//對于特例則高四位的閏月位置表示法中的最高為設置為1 特殊處理
用wLeapNormal變量 
// //2017/06 28330->61098 2036/06 27947->60715 2047/05 
23133->55901 

//如果希望用匯編,這里有一條信息:農歷不會滯后公歷2個月. 
//將公歷轉換為農歷 
//返回:12位年份+4位月份+5位日期 
function DecodeGregToCNDate(dtGreg:TDateTime):TCNDate; 
var 
iDayLeave:Integer; 
wYear,wMonth,wDay:WORD; 
i,j:integer; 
wBigSmallDist,wLeap,wCount,wLeapShift:WORD; 
label OK; 
begin 
result := 0; 
iDayLeave := Trunc(dtGreg) - cstDateOrg; 
DecodeDate(IncMonth(dtGreg,-1),wYear,wMonth,wDay); 
if (iDayLeave < 0) or (iDayLeave > 22295 )then Exit; 
//Raise Exception.Create('目前只能算1990-01-27以后的'); 
//Raise Exception.Create('目前只能算2051-02-11以前的'); 
for i:=Low(cstCNTable) to High(cstCNTable) do begin 
wBigSmallDist := cstCNTable[i]; 
wLeap := wBigSmallDist shr 12; 
if wLeap > 12 then begin 
wLeap := wLeap and 7; 
wLeapShift := 1; 
end else 
wLeapShift := 0; 
for j:=1 to 12 do begin 
wCount:=(wBigSmallDist and 1) + 29; 
if j=wLeap then wCount := wCount - wLeapShift; 
if iDayLeave < wCount then begin 
Result := (i shl 9) + (j shl 5) + iDayLeave + 1; 
Exit; 
end; 
iDayLeave := iDayLeave - wCount; 
if j=wLeap then begin 
wCount:=29 + wLeapShift; 
if iDayLeave < wCount then begin 
Result := (i shl 9) + (j shl 5) + iDayLeave + 1 + (1 shl 21); 
Exit; 
end; 
iDayLeave := iDayLeave - wCount; 
end; 
wBigSmallDist := wBigSmallDist shr 1; 
end; 
end; 
//返回值: 
// 1位閏月標志 + 12位年份+4位月份+5位日期 (共22位) 
end; 
function isCNLeap(cnDate:TCNDate):boolean; 
begin 
result := (cnDate and $200000) <> 0; 
end; 
function 
GetGregDateFromCN(cnYear,cnMonth,cnDay:word;bLeap:Boolean=Fal
se):TDateTime; 
var 
i,j:integer; 
DayCount:integer; 
wBigSmallDist,wLeap,wLeapShift:WORD; 
begin 
// 0101 010010101111 高四位是閏月位置,后12位表示大小月,大月30
天,小月29天, 
DayCount := 0; 
if (cnYear < 1990) or (cnYear >2050) then begin 
Result := 0; 
Exit; 
end; 
for i:= cstCNYearOrg to cnYear-1 do begin 
wBigSmallDist := cstCNTable[i]; 
if (wBIgSmallDist and $F000) <> 0 then DayCount := DayCount + 
29; 
DayCount := DayCount + 12 * 29; 
for j:= 1 to 12 do begin 
DayCount := DayCount + wBigSmallDist and 1; 
wBigSmallDist := wBigSmallDist shr 1; 
end; 
end; 
wBigSmallDist := cstCNTable[cnYear]; 
wLeap := wBigSmallDist shr 12; 
if wLeap > 12 then begin 
wLeap := wLeap and 7; 
wLeapShift := 1; //大月在閏月. 
end else 
wLeapShift := 0; 
for j:= 1 to cnMonth-1 do begin 
DayCount:=DayCount + (wBigSmallDist and 1) + 29; 
if j=wLeap then DayCount := DayCount + 29; 
wBigSmallDist := wBigSmallDist shr 1; 
end; 
if bLeap and (cnMonth = wLeap) then //是要閏月的嗎? 
DayCount := DayCount + 30 - wLeapShift; 
result := cstDateOrg + DayCount + cnDay - 1; 
end; 

//將日期顯示成農歷字符串. 
function GregDateToCNStr(dtGreg:TDateTime):String; 
const hzNumber:array[0..10] of string=('零','一','二','三','
四','五','六','七','八','九','十'); 
function ConvertYMD(Number:Word;YMD:Word):string; 
var 
wTmp:word; 
begin 
result := ''; 
if YMD = 1 then begin //年份 
while Number > 0 do begin 
result := hzNumber[Number Mod 10] + result; 
Number := Number DIV 10; 
end; 
Exit; 
end; 
if Number<=10 then begin //可只用1位 
if YMD = 2 then //月份 
result := hzNumber[Number] 
else //天 
result := '初' + hzNumber[Number]; 
Exit; 
end; 
wTmp := Number Mod 10; //個位 
if wTmp <> 0 then result := hzNumber[wTmp]; 
wTmp := Number Div 10; //十位 
result:='十'+result; 
if wTmp > 1 then result := hzNumber[wTmp] + result; 
end; 
var 
cnYear,cnMonth,cnDay:word; 
cnDate:TCNDate; 
strLeap:string; 
begin 
cnDate:= DecodeGregToCNDate(dtGreg); 
if cnDate = 0 then begin 
result := '輸入越界'; 
Exit; 
end; 
cnDay := cnDate and $1F; 
cnMonth := (cnDate shr 5) and $F; 
cnYear := (cnDate shr 9) and $FFF; 
//測試第22位,為1表示閏月 
if isCNLeap(cnDate) then strLeap:='(閏)' else strLeap := ''; 
result := '農歷' + ConvertYMD(cnYear,1) + '年' + 
ConvertYMD(cnMonth,2) + '月' 
+ strLeap + ConvertYMD(cnDay,3) ; 
end; 
end. 


 
 


--------------------------------------------------------------------------------