วิวิธโหราศาสตร์

3 พฤษภาคม 2559

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

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

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

day = วันที่

month = เดือน

year = ปี ค.ศ.

hr = ชั่วโมง

mn = นาที

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

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

Function neptune(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.000103311 * days + 5.233615854 - 0.125768426, 6.283185307)

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

L = modulo(Ta + 0.125768426, 6.283185307)

RV = 30.26664 * (1 - 0.0102981 ^ 2) / (1 + 0.0102981 * 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 - 2.300213054), Sin(L - 2.300213054) * Cos(0.030859167)) + 2.300213054

Hlat = Asin(Sin(L - 2.300213054) * Sin(0.030859167))

Dis = RV * Cos(Hlat)

'ผลลัพธ์

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

End Function

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

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

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

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

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

ตำแหน่งดาวเนปจูน(นิรายนะวิธี) 341.4010258 - 24.0852731799126 = 317.3157526

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

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

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

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

สรุป ตำแหน่งดาวเนปจูน(นิรายนะ) = ราศีกุมภ์ 17 องศา 18 ลิปดา 56 ฟิลิปดา

เปรียบเทียบผลกับโมดูล swiss คือ ราศีกุมภ์ 17 องศา 27 ลิปดา 6 ฟิลิปดา ต่างกันอยู่ 8 ลิปดา 10 ฟิลิปดา

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

3 พฤษภาคม 2559

(ต้องการขอรับไฟล์งาน Excel ติดต่อได้ที่ tepar2009@gmail.com)

2 พฤษภาคม 2559

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

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

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

day = วันที่

month = เดือน

year = ปี ค.ศ.

hr = ชั่วโมง

mn = นาที

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

Function uranus(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.000202859 * days + 5.291658 - 3.066206647, 6.283185307)

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

L = modulo(Ta + 3.066206647, 6.283185307)

RV = 19.30181 * (1 - 0.0428959 ^ 2) / (1 + 0.0428959 * 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.293208691), Sin(L - 1.293208691) * Cos(0.0134989)) + 1.293208691

Hlat = Asin(Sin(L - 1.293208691) * Sin(0.0134989))

Dis = RV * Cos(Hlat)

'ผลลัพธ์

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

uranus = lambda * (180 / 3.141592654)

End Function

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

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

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

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

ตำแหน่งดาวมฤตยู(นิรายนะวิธี) 21.42093018 - 24.0852731799126 = -2.664343

กรณีค่าที่ได้ติดลบ ให้นำ 360 มาบวกเข้า = 360 + -2.664343 = 357.335657

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

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

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

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

สรุป ตำแหน่งดาวมฤตยู(นิรายนะ) = ราศีมีน 27 องศา 20 ลิปดา 8 ฟิลิปดา

เปรียบเทียบผลกับโมดูล swiss คือ ราศีมีน 27 องศา 35 ลิปดา 0 ฟิลิปดา ต่างกันอยู่ 14 ลิปดา 52 ฟิลิปดา

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

2 พฤษภาคม 2559

(ต้องการขอรับไฟล์งาน Excel ติดต่อได้ที่ tepar2009@gmail.com)

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

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

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

day = วันที่

month = เดือน

year = ปี ค.ศ.

hr = ชั่วโมง

mn = นาที

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

Function saturn(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.00058096 * days + 0.365778949 - 1.550951933, 6.283185307)

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

L = modulo(Ta + 1.550951933, 6.283185307)

RV = 9.5719 * (1 - 0.0531651 ^ 2) / (1 + 0.0531651 * 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) 'Mean anomaly

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.983318858), Sin(L - 1.983318858) * Cos(0.043375621)) + 1.983318858

Hlat = Asin(Sin(L - 1.983318858) * Sin(0.043375621))

Dis = RV * Cos(Hlat)

'ผลลัพธ์

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

saturn = lambda * (180 / 3.141592654)

End Function

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

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

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

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

ตำแหน่งดาวเสาร์(นิรายนะวิธี) 253.4010088 - 24.0852731799126 = 229.3157356

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

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

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

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

สรุป ตำแหน่งดาวเสาร์(นิรายนะ) = ราศีพิจิก 19 องศา 18 ลิปดา 56 ฟิลิปดา

เปรียบเทียบผลกับโมดูล swiss คือ ราศีพิจิก 21 องศา 14 ลิปดา 37 ฟิลิปดา ต่างกันอยู่ 1 องศา 55 ลิปดา 41 ฟิลิปดา

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

2 พฤษภาคม 2559

(ต้องการขอรับไฟล์งาน Excel ติดต่อได้ที่ tepar2009@gmail.com)

1 พฤษภาคม 2559

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

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

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

day = วันที่

month = เดือน

year = ปี ค.ศ.

hr = ชั่วโมง

mn = นาที

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

Function jupiter(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.001450302 * days + 5.629731068 - 0.273978295, 6.283185307)

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

L = modulo(Ta + 0.273978295, 6.283185307)

RV = 5.202597 * (1 - 0.0484646 ^ 2) / (1 + 0.0484646 * 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.753554989), Sin(L - 1.753554989) * Cos(0.022770089)) + 1.753554989

Hlat = Asin(Sin(L - 1.753554989) * Sin(0.022770089))

Dis = RV * Cos(Hlat)

'ผลลัพธ์

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

jupiter = lambda * (180 / 3.141592654)

End Function

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

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

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

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

ตำแหน่งดาวพฤหัสบดี(นิรายนะวิธี) 163.303458 - 24.0852731799126 = 139.2181848

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

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

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

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

สรุป ตำแหน่งดาวพฤหัสบดี(นิรายนะ) = ราศีสิงห์ 19 องศา 13 ลิปดา 5 ฟิลิปดา

เปรียบเทียบผลกับโมดูล swiss คือ ราศีสิงห์ 19 องศา 16 ลิปดา 8 ฟิลิปดา ต่างกันอยู่ 56 ลิปดา 57 ฟิลิปดา

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

2 พฤษภาคม 2559

(ต้องการขอรับไฟล์งาน Excel ติดต่อได้ที่ tepar2009@gmail.com)

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

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

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

day = วันที่

month = เดือน

year = ปี ค.ศ.

hr = ชั่วโมง

mn = นาที

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

Function mar(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.009146595 * days + 4.580229857 - 5.865845667, 6.283185307)

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

L = modulo(Ta + 5.865845667, 6.283185307)

RV = 1.5236365 * (1 - 0.0934231 ^ 2) / (1 + 0.0934231 * 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 - 0.865096878), Sin(L - 0.865096878) * Cos(0.032287195)) + 0.865096878

Hlat = Asin(Sin(L - 0.865096878) * Sin(0.032287195))

Dis = RV * Cos(Hlat)

'ผลลัพธ์

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

mar = lambda * (180 / 3.141592654)

End Function

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

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

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

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

ตำแหน่งดาวอังคาร(นิรายนะวิธี) 248.279949 - 24.0852731799126 = 224.1946758

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

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

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

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

สรุป ตำแหน่งดาวอังคาร(นิรายนะ) = ราศีพิจิก 14 องศา 11 ลิปดา 40 ฟิลิปดา

เปรียบเทียบผลกับโมดูล swiss คือ ราศีพิจิก 13 องศา 38 ลิปดา 33 ฟิลิปดา ต่างกันอยู่ 33 ลิปดา 10 ฟิลิปดา

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

1 พฤษภาคม 2559

(ต้องการขอรับไฟล์งาน Excel ติดต่อได้ที่ tepar2009@gmail.com)