儒略日的計算(2)

儒略日的計算(1)

http://newscienceview.blogspot.com/2009/10/1_18.html

計算原則

        以下所討論的計算方法與範例中，皆已考慮了格里曆的修正，也就是公元1582年 十月四日 的後一日為1582年 十月十五日 。而我們通常所稱的西元前一年，在數學計算式中需記為“ ０” ，西元前二年則記為“－ １” ，以此類推。

       至於數學計算式中出現的INT(X)，是指取Ｘ的整數部分，如INT(5.32)=5，INT(-4.71)=-4。

此外，文中所應用的代號敘述如下：YYYY為公曆的年份，MM為公曆月份，DD.dddd表公曆之日，其中dddd即時、分、秒換算成附加在日之後的小數部分。y為實際代入計算式中所使用的年數，m為月數。其他各符號為不同計算時所應用的參數。YYYY與y以及MM與m間的關係為：

若MM＞2，m=MM

若MM＝1或2，y=YYYY－1，m＝MM+12

公曆換算為儒略日

當計算的日期在1582年 十月十五日 之前，參數A、B以0計，但在此日期之後（含當日）者，必須另依下式計算Ａ與Ｂ參數值：

A＝INT(y/100)

B＝2－A+INT(A/4 )

根據上述原則與參數，便可計算儒略日如下：

公曆轉換為儒略日

JD=INT(365.25*y)＋INT(30.6001*(m+1))+DD.dddd+1720994.5+B

儒略日轉換為修正儒略日

MJD=JD－2400000.5

【例一】格林威治時間2001年 十二月十五日10：36：30的儒略日與修正儒略日為何？

【解】YYYY＝2001，MM=12，y=2001，m=12

DD.dddd＝15+10/24+36/1440+30/86400≒15.4420

A＝INT(2001/100)=20

B＝2－20＋INT(20/4)=－13

JD=INT(365.25*2001)+INT(30.6001*(12+1))+15.4420+1720994.5+(-13) =2452258.9420

MJD=2452258.9420－2400000.5＝52258.4420

【例二】中原標準時間2002年 一月十日15：42：15的儒略日與修正儒略日為何？

【解】YYYY=2002，MM=1，y=2001，m=13

先將中原標準時間減八小時換算成格林威治

時間2002年 一月十日7：42：15。

DD.dddd=10+7/24+42/1440+15/86400≒10.3210

A=INT(2001/100)=20

B=2－20+INT(20/4)＝－13

JD=INT(365.25*2001)+INT(30.6001*(13+1))+10.3210+172094.5+(－13)=2452284.8210

MJD=2452284.8210－2400000.5＝52284.3210

Link :
天文計算實驗室 http://blog.sina.com.tw/astro_calculator/