說到節氣公式,只要您不要求十分精確,實在也不是什麼難事,我曾經用每100年取一
個節氣時間,約3000多年30多個數據,就可以造出一個近似公式來,雖然沒有命理網節氣計
算的神準功能,卻可以讓我逍遙於千百年干支領域之中.目前網路上最為久遠的萬年曆當
推陳達明先生的萬年曆,也只服務到西元前104年而已,雖然大陸作出了西元前1500年之後
的陰陽對照曆,也有干支,但是他們用儒略曆,換算費時,且只有表可查,並沒有程式化,節氣
與朔日分開不同的表,不方便應用.
我在兩年前曾經列出這些公式,如果您有興趣,不妨看看能不能拿來濫竽充數一番 :
如曆元取1985年,西元前1500年與1985年差為-3484年(因為沒有西元0年)
1985年冬至為2446422.256,-3484年的歲差為66.6秒,則西元前1500年冬至日
為2446422.256+(-3484)*365.242199074+(-3484)*66.6/24/3600=1173915.749
時間0.749*23.93446959=17.92691772=17時56分
(補註:我先用每個選用年節氣時間與1985年的節氣 時間求出平均歲差,再造歲差的公式.)
1985年立春2446101.217,西元500年與1985年差為-1485年
-1485年的歲差為45.9秒,則西元500年立春日
為2446101.217+(-1485)*365.242199074+(-1485)*45.9/24/3600=1903715.762
時間:0.762*23.93446959=18.23806583=18時14分
X=Y-1985 (Y若為西元前年數先減1再用負數,如西元前1500年,以-1499年去算)
JD=JDo+(aX^2+bX+c)/24/3600X+365.242199074X ...公式1
a b c
冬至 -5E-07 -0.0074 47.406
小寒 -5E-07 -0.0058 46.09
大寒 -5E-07 -0.0041 44.783
立春 -6E-07 -0.0026 43.362
雨水 -4E-07 -0.001 33.909
驚蟄 -2E-07 0.0004 24.316
春分 -1E-07 0.0017 14.619
清明 2E-07 0.0023 2.3576
榖雨 5E-07 0.0029 -10.001
立夏 9E-07 0.004 -21.772
小滿 9E-07 0.0019 -31.021
芒種 1E-06 0.0003 -39.611
夏至 1E-06 -0.0012 -48.22
小暑 1E-06 -0.0029 -46.559
大暑 1E-06 -0.0049 -45.21
立秋 1E-06 -0.0069 -43.86
處暑 7E-07 -0.0089 -34.344
白露 4E-07 -0.0109 -24.886
秋分 1E-07 -0.0127 -15.438
寒露 -1E-07 -0.00128 -2.5353
霜降 -3E-07 -0.0124 10.621
立冬 -4E-07 -0.0121 23.519
小雪 -5E-07 -0.0106 31.524
大雪 -5E-07 -0.009 39.484
西元1985年各節氣資料如下(JDo)
小寒2446071.735
大寒2446086.457
立春2446101.217
雨水2446116.047
驚蟄2446130.972
春分2446146.010
清明2446161.177
榖雨2446176.478
立夏2446194.907
小滿2446207.448
芒種2446223.083
夏至2446238.783
小暑2446254.515
大暑2446270.235
立秋2446285.922
處暑2446301.526
白露2446317.037
秋分2446332.423
寒露2446347.685
霜降2446362.809
立冬2446377.814
小雪2446392.704
大雪2446407.513
冬至2446422.256
以上公式每節氣大約只用了34個點,能有如此成果已算不錯.而因為書上只有二分二至
與四立八個節氣的資料(指西元前部分),不足的其他的16個節氣資料是我以某種假設造出
來的,有時間的人把每一年各節氣時間輸入建檔,要用時交給電腦去搜尋取出,又快又準,但
是要讓他推西元前1600年的節氣,就超出能力範圍了.懶人要找公式,由電腦去計算,雖有小
誤差,如30分,15分,又有何妨?電腦高手可以寫個程式算出來比較看看,我相信一定比史學教
授推薦的近似公式更好用.
附註:
史學大師丁驌先生冬至公式(曆元1938年2429255.843055戊子冬至)
JD=1721416.45945+365.24219879Y ...公式2(參考[商周定氣冬至儒日表]一書)
或
冬至公式(曆元1900年2415376.238冬至)
JD=1721414.961+365.2427775Y ...公式3
比較(冬至)
西元前1500年公式1=1173915.728(17:25) X=-1499-1985=-3484
公式2=1173918.403(相差2日以上) Y=-1499
公式3=1173916.038(相差6小時以上)
書上=1173915.749(17:56)(參考大陸張培瑜先生著[三千五百年曆日天象]一書,國家圖
書館有藏)
西元2002年公式1=2452631.383(9:10) X=2002-1985=17
公式2=2452631.341(8:10) Y=2002
公式3=2452631.002(0:03)
氣象局=2452631.385775557(9:14)