摘? 要： 基于農(nóng)歷節(jié)氣的內(nèi)在規(guī)律提出了一種新型實(shí)用列表算法——按位列表法(BSMA)，利用單片機(jī)Tm8706實(shí)現(xiàn)了電子萬(wàn)年歷的研究和開(kāi)發(fā)。
關(guān)鍵詞： 節(jié)氣? 列表法? 微處理器? 單片機(jī)? 算法
?
農(nóng)歷24節(jié)氣是我國(guó)民間傳統(tǒng)節(jié)令,對(duì)于農(nóng)事活動(dòng)有廣泛的應(yīng)用價(jià)值。因而帶有24節(jié)氣的電子萬(wàn)年歷在農(nóng)村有廣泛的市場(chǎng)。
本文通過(guò)分析24節(jié)氣,找出了其中內(nèi)在的隱含規(guī)律,并且根據(jù)這種規(guī)律提出了一種最新的查表方法。由于這種查表方法是基于位查詢(xún)的,所以稱(chēng)之為按位列表法(Bit Schedule Method Arithmetic,BSMA)。
1? 陰歷節(jié)氣的隱含規(guī)律
為了便于考察50年中無(wú)規(guī)律的陰歷節(jié)氣分布時(shí)間,本文利用C語(yǔ)言在計(jì)算機(jī)中實(shí)現(xiàn)了50年(2000～2049年)的陽(yáng)歷和陰歷節(jié)氣分布時(shí)間表。由于數(shù)據(jù)龐大,所以?xún)H列出其中部分陽(yáng)歷年對(duì)應(yīng)的前6個(gè)節(jié)氣日期,見(jiàn)表1。
?
表1中,每年的節(jié)氣具體日期雖然都在一定的范圍內(nèi)變動(dòng),但是變動(dòng)的趨勢(shì)毫無(wú)規(guī)律可尋,每個(gè)月的具體節(jié)氣日期有時(shí)候推前一、二天,有時(shí)候推后一、二天。但是,從這張復(fù)雜表格中的眾多數(shù)據(jù)中可以發(fā)現(xiàn)一個(gè)規(guī)律:雖然,無(wú)法獲得這些節(jié)氣變化的具體日期,但是,這些日期的變化總是在3天之內(nèi)發(fā)生。例如:每年的第1個(gè)節(jié)氣(小寒)出現(xiàn)在1月的5日、6日2天中的某一天;每年的第2個(gè)節(jié)氣(大寒)出現(xiàn)在1月的19日、20日、21日3天中的某一天。這樣,只要固定每個(gè)節(jié)氣的基本日期,然后專(zhuān)門(mén)針對(duì)這二、三天編制一個(gè)相對(duì)日期的表格就可以了。考慮到最多變化范圍是3天,所以這個(gè)相對(duì)日期可以取值為1、0和-1。
選定的基本日期是:1月到3月的節(jié)氣基本日期:6/20、4/19、6/21;4月到6月的節(jié)氣基本日期:5/20、6/21、6/21;6月到9月的節(jié)氣基本日期:7/23、8/23、8/23;10月到12月的節(jié)氣基本日期:8/24、8/22、7/22。
根據(jù)以上的思路,可得到表2所示的相對(duì)值表。
仔細(xì)觀察表2,可以得到以下關(guān)于節(jié)氣日期相對(duì)值的規(guī)律性結(jié)論:
(1)只有4個(gè)節(jié)氣(大寒、夏至、寒露、小雪)是3個(gè)相對(duì)值(1,-1和0)同時(shí)出現(xiàn)。
(2)其他的20個(gè)節(jié)氣中僅僅出現(xiàn)了0和-1這2個(gè)相對(duì)值。
(3)在有3個(gè)相對(duì)值的4個(gè)節(jié)氣(大寒、夏至、寒露、小雪)中,2044年以前,出現(xiàn)的非零相對(duì)值都是1;而2044年以及2044年以后,出現(xiàn)的非零相對(duì)值都是-1;
?
也就是說(shuō)在50年內(nèi),大寒(2)、夏至(12)、寒露(19)和小雪(22)這4個(gè)節(jié)氣日期變化相對(duì)較多,一共有3天,它們的日期相對(duì)值可能需要使用3個(gè)狀態(tài)值來(lái)表示。
圖1是以上4條農(nóng)歷節(jié)氣日期規(guī)律的相對(duì)值分布簡(jiǎn)圖。
?
根據(jù)表2和圖1所示的農(nóng)歷節(jié)氣日期排布規(guī)律和形式,就可以構(gòu)造一個(gè)表格,并且結(jié)合相應(yīng)的算法來(lái)計(jì)算某個(gè)陽(yáng)歷日期所對(duì)應(yīng)的陰歷節(jié)氣了。
2?按位列表法(BSMA)及其編碼方法
利用上節(jié)分析結(jié)果,設(shè)想使用最少的字節(jié)來(lái)編制1個(gè)表格,使其包含了50年內(nèi)的所有相對(duì)值分布信息。由于1年包含農(nóng)歷的24個(gè)節(jié)氣,也就是說(shuō)共有24個(gè)日期相對(duì)值狀態(tài)量。
為了能夠盡量減少存儲(chǔ)空間,可以使用1個(gè)位來(lái)代表1個(gè)農(nóng)歷的節(jié)氣日期相對(duì)值。這樣一年的24個(gè)節(jié)氣就可以用24個(gè)位來(lái)表示,也就是3個(gè)字節(jié)(每個(gè)字節(jié)8個(gè)位)表示。這就是所說(shuō)的按位列表法。
由于農(nóng)歷節(jié)氣日期相對(duì)值最多有3種狀態(tài):-1、0和1,而1個(gè)位僅僅有2種狀態(tài):0和1。為了能夠把3種狀態(tài)融合進(jìn)2種狀態(tài)中,本文根據(jù)上節(jié)中節(jié)氣日期相對(duì)值的規(guī)律性結(jié)論,采取了具有條件的二值狀態(tài)來(lái)表示三值狀態(tài),制定了如表3所示的編碼規(guī)則。
?
利用以上的規(guī)則就可對(duì)50年的24個(gè)節(jié)氣進(jìn)行編碼,得到如表4所示的編碼表。
?
由此看出,使用這種按位編碼的方式對(duì)農(nóng)歷日期相對(duì)值進(jìn)行編碼時(shí),一年24個(gè)節(jié)氣的日期只需要3個(gè)字節(jié)就可以明確表示了。50年的節(jié)氣表,只需要150個(gè)字節(jié)的存儲(chǔ)空間。年份越多,這種編碼發(fā)方式能夠節(jié)約的存儲(chǔ)空間越大,優(yōu)勢(shì)越明顯。
3?BSMA法在單片機(jī)中的實(shí)現(xiàn)
本文使用的單片機(jī)是十速公司的4位機(jī),型號(hào)為T(mén)m8706。這種單片機(jī)的結(jié)構(gòu)與8位機(jī)類(lèi)似,但是功能遠(yuǎn)遜于8位計(jì)算機(jī)。選此機(jī)型的原因是它具有結(jié)構(gòu)簡(jiǎn)單、成本低的優(yōu)勢(shì),同時(shí)還包括了普通單片機(jī)的基本功能。它包含了256個(gè)字節(jié)的RAM、2KB的ROM、定時(shí)器、外部中斷、內(nèi)部中斷,此外,還有省電功能和驅(qū)動(dòng)LED的功能。
由于其他的功能程序已經(jīng)占用了較大的存儲(chǔ)空間。所以能夠用于農(nóng)歷節(jié)氣查詢(xún)的程序的空間不大,僅僅有0.5KB的容量。也就是說(shuō)表格和查表程序一共只有0.5KB的空間。表格占用的空間越大,可寫(xiě)程序的空間就越小。
如果使用表1所示的50年農(nóng)歷24節(jié)氣日期表,僅僅表格就需要使用50×24個(gè)字節(jié)=1200個(gè)字節(jié)的存儲(chǔ)空間,其要求空間遠(yuǎn)遠(yuǎn)超出了0.5KB的容量。而使用LBA編碼方法得到的表格,只需要150個(gè)字節(jié),存儲(chǔ)空間僅僅是表1所需存儲(chǔ)空間的1/8,其余部分可以供查表程序使用。
在單片機(jī)中,查表程序是使用匯編程序完成的,實(shí)現(xiàn)起來(lái)也很容易。為了便于理解,這里把查表程序用C語(yǔ)言實(shí)現(xiàn)的部分介紹如下:
輸入?yún)?shù)是陽(yáng)歷的年月日:day、month、year,要求查表尋找它附近的農(nóng)歷節(jié)氣和農(nóng)歷節(jié)氣的日期。
if(day>=15) half_mon=1; //判斷是否是上半月
else half_mon=0;
col=(month-1)*2+half_mon; //獲得節(jié)氣順序標(biāo)號(hào)(0～23)
kc_byte=col/8; //獲得節(jié)氣日期的相對(duì)值所在的
//字節(jié)位置(0～2)
kc_bit=col%8; //獲得節(jié)氣日期的相對(duì)值在字節(jié)中的
//位位置(0～7)
tmp_modi=modi_datas[year-2000][kc_byte];
//獲得節(jié)氣日期相對(duì)值所在字節(jié)
value=( (tmp_modi< //日期相對(duì)值狀態(tài) if(value!=0){ //判斷年份,以決定節(jié)氣相對(duì)值1代表1 //還是-1。 if( (col==1 || col==11 || col==18 || col==21) && year<2044) ????????????? J_qi=days[col]+1; //年份<2044的4個(gè)特定 //月份中,取+1。 else ????????????? J_qi=days[col]-1; //其他情況下,取-1。 } else J_qi=days[col]; //狀態(tài)為0,節(jié)氣日期相對(duì)值 //為節(jié)氣基本日期值。 由此程序即可獲得所輸入陽(yáng)歷日期(年、月、日)附近的農(nóng)歷節(jié)氣準(zhǔn)確日期值。 通過(guò)上面的BSMA列表編碼和查詢(xún)的方法,實(shí)現(xiàn)了基于四位單片機(jī)的萬(wàn)年歷的研制和實(shí)現(xiàn)。它用匯編程序方便地實(shí)現(xiàn)了陰歷24節(jié)氣功能,同時(shí)還可以方便地算出某年的某天距離下一個(gè)節(jié)氣的天數(shù)。 從上面可以看出,按位列表法(BSMA)思路清晰明了,實(shí)現(xiàn)起來(lái)非常方便,同時(shí)節(jié)省了大量的ROM空間。所做萬(wàn)年歷的時(shí)間跨度越大,LBA算法在節(jié)省ROM空間的優(yōu)勢(shì)就越明顯。 ? 參考文獻(xiàn) 1? 李勇,許邦信.太陽(yáng)黃經(jīng)和節(jié)氣時(shí)間的近似推算與Newcomb太陽(yáng)表.南京大學(xué)學(xué)報(bào),1995;31(3) 2? 張培瑜,黃洪峰.中歷及二十四節(jié)氣時(shí)刻計(jì)算.廣西科學(xué),1994;1(3)