สูตรการคำนวณตำแหน่งดาวเคราะห์แบบดาราศาสตร์
: พลูโต
ฟังก์ชั่น pluto คือฟังก์ชั่น VBA บน Excel ใช้ในการหาค่าตำแหน่งดาวพลูโตแบบดาราศาสตร์สากล
หรือระบบสายนะวิธี มีตัวแปรที่รับเข้าดังต่อไปนี้
day = วันที่
month = เดือน
year = ปี ค.ศ.
hr = ชั่วโมง
mn = นาที
และฟังก์ชั่นเสริมที่ใช้ประกอบการคำนวณคือ Modulo , Atan2 , Asin , Acos (ดูรายละเอียดฟังก์ชั่นเสริมได้ที่การคำนวณดาวพุธ)
สูตรการคำนวณมีดังนี้
Function pluto(day, month, year, hr, mn As Double) As Double
Dim dj2000, days, Ma, Ta, L, RV As Double
dj2000 = 367 * year - Int(7 * (year + Int((month + 9) / 12)) / 4) +
Int(275 * month / 9) + day - 730531.5 + ((hr + (mn / 60)) / 24)
days = dj2000 + 864.5
'องค์ประกอบของดาวเคราะห์
Ma = modulo(0.0000690814 * days + 4.114885983 - 3.923543792,
6.283185307)
Ta = Ma + (2 * 0.2501272 - 0.2501272 ^ 3 / 4 + 5 / 96 * 0.2501272 ^
5) * Sin(Ma) + (5 * 0.2501272 ^ 2 / 4 - 11 / 24 * 0.2501272 ^ 4) * Sin(2 * Ma)
+ (13 * 0.2501272 ^ 3 / 12 - 43 / 64 * 0.2501272 ^ 5) * Sin(3 * Ma) + 103 / 96
* 0.2501272 ^ 4 * Sin(4 * Ma) + 1097 / 960 * 0.2501272 ^ 5 * Sin(5 * Ma)
L = modulo(Ta + 3.923543792, 6.283185307)
RV = 39.5804 * (1 - 0.2501272 ^ 2) / (1 + 0.2501272 * Cos(Ta))
'องค์ประกอบของโลก
Dim Mae, Tae, Le, RVe, Hlong, Hlat, Dis, lambda, beta, alpha, delta
As Double
Mae = modulo(0.017201609 * days + 5.731722874 - 1.795100806,
6.283185307)
Tae = Mae + (2 * 0.0166967 - 0.0166967 ^ 3 / 4) * Sin(Mae) + 5 *
0.0166967 ^ 2 / 4 * Sin(2 * Mae) + 13 * 0.0166967 ^ 3 / 12 * Sin(3 * Mae)
Le = modulo(Tae + 1.795100806, 6.283185307)
RVe = 1 * (1 - 0.0166967 ^ 2) / (1 + 0.0166967 * Cos(Tae))
'มุมมอง Helio
Hlong = Atan2(Cos(L - 1.926552024), Sin(L - 1.926552024) *
Cos(0.298824279)) + 1.926552024
Hlat = Asin(Sin(L - 1.926552024) * Sin(0.298824279))
Dis = RV * Cos(Hlat)
'ผลลัพธ์
lambda = modulo(Atn(RVe * Sin(Hlong - Le) / (Dis - RVe * Cos(Hlong
- Le))) + Hlong, 6.283185307)
pluto = lambda * (180 / 3.141592654)
End Function
ตัวอย่างการคำนวณ
1 พฤษภาคม 2016 เวลา 12.00 น.(กรีนิช)
รูปแบบสูตร pluto(1,5,2016,12,0)
ค่าที่ได้คือ 287.3219655 องศา
(เปรียบเทียบกับค่าที่ได้จากโมดูล swiss คือ 287.442411271609 ต่างกันอยู่ 0.120445772 องศา)
ในกรณีที่จะทำเป็นค่าตำแหน่งดาวในระบบนิรายนะวิธี
ให้นำค่าอายนางศมาลบออกจากค่าที่คำนวณได้
กรณีนี้ใช้ค่าอายนางศแบบลาหิรีที่คำนวณได้ในวันเวลาตามตัวอย่าง คือ
24.0852731799126
ตำแหน่งดาวพลูโต(นิรายนะวิธี) 287.3219655 - 24.0852731799126 = 263.2366923
1.ราศี 263.2366923 หาร 30
ลัพธ์เป็นราศี = 8
2.องศา จำนวนเต็มของเศษจากการหาร = 23.23669232 = 23
3.ลิปดา ทศนิยมที่เหลือคูณด้วย 60 = 0.23669232 X 60 = 14.20153921 = 14
4.ฟิลิปดา
ทศนิยมที่เหลือในข้อ 3. คูณด้วย 60 = 0.20153921 X 60 = 12.09235231 = 12
สรุป ตำแหน่งดาวพลูโต(นิรายนะ) = ราศีธนู 23 องศา
14 ลิปดา 12 ฟิลิปดา
เปรียบเทียบผลกับโมดูล swiss คือ ราศีธนู 23
องศา 21 ลิปดา 25 ฟิลิปดา ต่างกันอยู่ 7 ลิปดา 13 ฟิลิปดา
ธีรพร เพชรกำแพง
3 พฤษภาคม 2559
(ต้องการขอรับไฟล์งาน
Excel
ติดต่อได้ที่ tepar2009@gmail.com)