也來談談節氣 近道作於2002年1月19日
請問各位同好有人徹底去研究過節氣的公式,或規律的麼?
一個回歸年365.242199074日,冬至歲差50.2564秒(或50.2735秒),怎樣實際運算
呢?夏至呢?春分呢?小寒,大寒呢?節氣是全球性的,只是地球繞日某種定義下的軌
道交點的時間,如台北上午10:00交節氣,同時間日本是同日上午11:00交同一節氣
,只要同一經度(南北子午線),交節氣時間一定相同.根據我的整理,自西元1985年
起往前每100年一次,至有資料的西元前1500為止,將節氣時間與儒略周日結合,
就可計算每年平均歲差,例如此約3500年間,自西元前1500年起,冬至平均歲差約
在66至50秒之間,立春約45至47秒,春分約7至15秒,立夏約-26至-21秒,夏至約
-31至-47秒,立秋約-8至-41秒,秋分約在29至-14秒之間,立冬約在60至25秒,歲
差與年差成近似線性關係,可用EXCEL多項式解之,其觀念大概如下:
如曆元取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年立春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)