สูตรการคำนวณตำแหน่งดาวเคราะห์แบบดาราศาสตร์ : ศุกร์

สูตรการคำนวณตำแหน่งดาวเคราะห์แบบดาราศาสตร์ : ศุกร์

ฟังก์ชั่น venus คือฟังก์ชั่น VBA บน Excel ใช้ในการหาค่าตำแหน่งดาวศุกร์แบบดาราศาสตร์สากล หรือระบบสายนะวิธี มีตัวแปรที่รับเข้าดังต่อไปนี้

day = วันที่

month = เดือน

year = ปี ค.ศ.

hr = ชั่วโมง

mn = นาที

และฟังก์ชั่นเสริมที่ใช้ประกอบการคำนวณคือ Modulo , Atan2 , Asin , Acos (ดูรายละเอียดฟังก์ชั่นเสริมได้ที่การคำนวณดาวพุธ)

สูตรการคำนวณมีดังนี้

Function venus(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.027962932 * days + 4.135390984 - 2.299663276, 6.283185307)

Ta = Ma + (2 * 0.0067933 - 0.0067933 ^ 3 / 4 + 5 / 96 * 0.0067933 ^ 5) * Sin(Ma) + (5 * 0.0067933 ^ 2 / 4 - 11 / 24 * 0.0067933 ^ 4) * Sin(2 * Ma) + (13 * 0.0067933 ^ 3 / 12 - 43 / 64 * 0.0067933 ^ 5) * Sin(3 * Ma) + 103 / 96 * 0.0067933 ^ 4 * Sin(4 * Ma) + 1097 / 960 * 0.0067933 ^ 5 * Sin(5 * Ma)

L = modulo(Ta + 2.299663276, 6.283185307)

RV = 0.7233238 * (1 - 0.0067933 ^ 2) / (1 + 0.0067933 * 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.338473805), Sin(L - 1.338473805) * Cos(0.059249041)) + 1.338473805

Hlat = Asin(Sin(L - 1.338473805) * Sin(0.059249041))

Dis = RV * Cos(Hlat)

'ผลลัพธ์

lambda = modulo(Atn(Dis * Sin(Le - Hlong) / (RVe - Dis * Cos(Le - Hlong))) + 3.141592654 + Le, 6.283185307)

venus = lambda * 180 / 3.141592654

End Function

ตัวอย่างการคำนวณ 1 พฤษภาคม 2016 เวลา 12.00 น.(กรีนิช)

รูปแบบสูตร venus(1,5,2016,12,0)

ค่าที่ได้คือ 31.55983249 องศา

(เปรียบเทียบกับค่าที่ได้จากโมดูล swiss คือ 31.8173858814644 ต่างกันอยู่ 0.257553389 องศา)

ในกรณีที่จะทำเป็นค่าตำแหน่งดาวในระบบนิรายนะวิธี ให้นำค่าอายนางศมาลบออกจากค่าที่คำนวณได้ กรณีนี้ใช้ค่าอายนางศแบบลาหิรีที่คำนวณได้ในวันเวลาตามตัวอย่าง คือ 24.0852731799126

ตำแหน่งดาวศุกร์(นิรายนะวิธี)  31.55983249 - 24.0852731799126 = 7.47455931

1.ราศี 7.47455931 หาร 30 ลัพธ์เป็นราศี = 0

2.องศา จำนวนเต็มของเศษจากการหาร = 7.47455931 = 7

3.ลิปดา ทศนิยมที่เหลือคูณด้วย 60 = 0.47455931 X 60 = 28.4735586 = 28

4.ฟิลิปดา ทศนิยมที่เหลือในข้อ 3. คูณด้วย 60 = 0.4735586 X 60 = 28.413516= 28

สรุป ตำแหน่งดาวศุกร์(นิรายนะ) = ราศีเมษ 7 องศา 28 ลิปดา 28 ฟิลิปดา

เปรียบเทียบผลกับโมดูล swiss คือ ราศีเมษ 7 องศา 22 ลิปดา 22 ฟิลิปดา ต่างกันอยู่ 6 ลิปดา 6 ฟิลิปดา

ธีรพร  เพชรกำแพง

1 พฤษภาคม 2559
(ต้องการขอรับไฟล์งาน Excel ติดต่อได้ที่ tepar2009@gmail.com)