Suriyayart : Thai Classical Astronomy

สุริยยาตร์ : ตำราคำนวณตำแหน่งดาวนพเคราะห์ ที่ตกทอดมาสู่สยาม

เป็นหลักในการเทียบทำปฏิทินจันทรคติไทย และ ปฏิทินโหร ตลอดจนการผูกดวง

วันเสาร์, พฤษภาคม 21, 2548

หรคุณ และ ปฏิทินสากล Julian

เทคนิคการแปลง หรคุณ กับ ปฏิทินสากลนี้
จะใช้สำหรับปฏิทินสากลแบบ Julian
ซึ่งเป็นปฏิทินสากล ที่ใช้ก่อนวันที่ 15 ตุลาคม ค.ศ. 1582

Julian --> HD
------------------


ในวันที่ d เดือน m ปีค.ศ. y เวลา hr:mn:sc
ตามปฏิทินสากลแบบ Julian
สามารถ หา hd ได้ตามขั้นตอนดังนี้


1. ถ้า m > 2 แล้ว jy = y และ jm = m + 1
ถ้า m <= 2 แล้ว jy = y - 1 และ jm = m + 13

2. intgr = Floor[ Floor[365.25* jy] + Floor[30.6001* jm] + d + 1720995 ]

3. dayfrac = hr/24.0-0.5
ถ้า dayfrac < 0 แล้ว dayfrac = dayfrac + 1 และ intgr = intgr - 1

4. frac = dayfrac + (mn + sc/60)/60/24

5. hd = Round[ (intgr + frac - 1954167.5)*100000] / 100000



HD --> Julian
-------------------


การหาวันที่ d เดือน m ค.ศ. y เวลา hr:mn:sc
ตามปฏิทินสากลแบบ Julian
จาก hd ที่ต้องการ


1. intgr = Floor[hd + 1954167.5]

2. frac = (hd + 1954167.5) - intgr

3. j1 = intgr

4. dayfrac = frac + 0.5

5. ถ้า ( dayfrac > 1 ) แล้ว dayfrac = dayfrac - 1 และ j1 = j1+1

6. j2 = j1 + 1524

7. j3 = Floor[ 6680.0 + ( (j2 - 2439870) - 122.1 )/365.25 ]

8. j4 = Floor[ j3*365.25 ]

9. j5 = Floor[ (j2 - j4)/30.6001 ]

10. d = Floor[ j2 - j4 - Floor[j5*30.6001] ]

11. m = Floor[ j5 - 1]

12. ถ้า ( m > 12 ) แล้ว m = m - 12

13. y = Floor[ j3 - 4715 ]

14. ถ้า ( m > 2 ) แล้ว y = y - 1

15. ถ้า ( y <= 0 ) แล้ว y = y - 1

16. hr = Floor[dayfrac*24]

17. mn = Floor[(dayfrac*24 - hr)*60];

18. sc = Round[((dayfrac*24 - hr)*60.0 - mn)*60.0];

19. ถ้า sc >= 60 แล้ว sc = sc - 60 และ mn = mn+1

แปลงจาก hd จะได้วันที่ตามปฏิทินสากลแบบ Julian
คือ วันที่ d เดือน m ค.ศ. y เวลา hr:mn:sc

0 ความคิดเห็น:

แสดงความคิดเห็น

สมัครสมาชิก ส่งความคิดเห็น [Atom]

<< หน้าแรก