B. MENGHITUNG DATA EPHEMERIS MATAHARI
DENGAN
ALGHORITMA JEAN MEEUS (HIGH ACCURACY)
Setelah pada
postingan sebelumnya kita telah mencoba menghitung data ephemeris Matahari
dengan low accuracy, sekarang penulis akan melanjutkan dengan
perhitungan data ephemeris Matahari dengan menggunakan Algoritma Jean Meeus high
accuracy pada buku Astronomical Algorithms Chapter 24, Solar Coordinates.
Perhitungan
data ephemeris Matahari high accuracy dengan akurasi lebih
tinggi dari 0,01 detik busur bisa kita dapatkan dengan menggunakan VSOP87
Theory, dengan total jumlah koreksi sebanyak 2425 buah. 1080 koreksi untuk
bujur ekliptika, 348 koreksi untuk lintang ekliptika dan 997 koreksi untuk
jarak Matahari ke Bumi. (ada yang mau nyoba? ^_^). Sedangkan koreksi pada
perhitungan high accuracy dengan Algoritma Jean Meeus sebenarnya
merupakan reduksi dari VSOP87 theory dengan mengambil koreksi-koreksi yang
penting. Total koreksi pada Algoritma Jean Meeus sebanyak 159 koreksi, dengan
kesalahan tidak lebih dari 1 detik untuk periode tahun -2000 sampai 6000
(keren..^_^).
Seperti
biasanya, buku ini cukup membingungkan. Karena tidak semua cara perhitungannya
tersedia di chapter tersebut. Kita harus buka bolak-balik beberapa chapter
sebelum dan sesudahnya, seperti Chapter 21, Chapter 53, Chapter 11 dan sampai harus buka-buka juga bagian Appendix
II di akhir buku.dll.. But, jangan
berkecil hati, di sini rumus-rumus yang tersebar di Chapter-Chapter
tersebut sudah penulis kumpulkan di sini.^_^
Pada postingan
ini, kita akan mencoba dulu menghitung data ephemeris Matahari high accuracy.
Menghitung data
ephemeris Matahari High Accuracy
Perhitungan
dengan metode ini sebenarnya langkah-langkahnya cukup simple. Pertama kita buka
tabel koreksi pada Appendix II di bagian akhir buku. Kemudian cara menggunakan
koreksi tabel tersebut menggunakan rumus yang ada pada Chapter 31. Tapi
sebelumnya harus menghitung Julian Day UT dan JDE pada Chapter 7 dan 9. Setelah
itu menghitung koreksi nutasi, dan true obliquity pada Chapter 21, kemudian
menghitung koreksi Aberasi dan mengubah bujur dan lintang geosentrik Matahari
menjadi Apparent pada Chapter 24 dll..(simple tapi yach lumayan mumet.>,<)
Agar lebih
mudah dalam memahami perhitungannya, mari kita coba dengan contoh.
# Menghitung
Data Ephemeris Matahari
Tanggal 26 Desember 2019
Jam = 7 : 00 : 00 WIB
A. Menghitung
Julian Day (Chapter 7)
Sebelum masuk
ke rumus Julian Day, Tanggal dan Jamnya harus dirubah ke UT atau GMT
7:00:00
WIB = 0:00:00 GMT
- Jika bulan > 2, maka M = bulan Y = Tahun
- Jika bulan = 1 atau 2, maka M = bulan + 12 Y = Tahun - 1
Maka M = 12, dan Y = 2019
A = INT(Y/100)
= 20
B = 2 – A + INT
(A/4)
= -13
Control B (tgl)
|
control1 ==> Jumlah hari sampai 05-10-1582 =
578130
|
control2 ==> Jumlah hari sampai 26-12 - 2019 =
737825
|
C = 13
|
JD = INT(365,25 * (Y + 4716)) +
INT(30,6001
* (M + 1)) + C + (Jam + Menit/60
+Detik/3600)
/ 24 – 1524,5
= INT(365.25 *
(2019 + 4716)) + INT(30.6001 * (12 + 1)) + 13 + (0 + 0/ 60 + 0/ 3600) / 24 -
1524.5
|
= 2458843.5
|
B. Menghitung
Julian Day Ephemeris (Chapter 9)
∆T = -15 + ((JD –
2382148)2 / 41048480
= 118,4196099 detik
= 0,001370597 hari
JDE
= JD + ∆T
= 128.298843714798 detik
|
= 0.00148494032077312 hari
|
C. Menghitung
Periodic Terms pada Appendix III
Pada Appendix
III, terdapat tabel-tabel koreksi orbit planet-planet. Tiap data orbit
planet-planet tersebut ada 3 macam data, yaitu L untuk bujur heliosentris
planet, B untuk lintang heliosentris planet dan R untuk jarak planet ke
Matahari. Kita akan menggunakan tabel Earth, karena kita akan menghitung posisi
Matahari dari Bumi. Untuk data Bumi, L terbagi menjadi 6 kelompok (L0,
L1,L2,L3,L4, dan L5), B terbagi menjadi 2 kelompok (B0 dan B2) dan R terbagi
menjadi 5 kelompok (R0, R1,R2,R3,dan R4).
Tapi sebelumnya
kita harus menghitung T (TD) dan Ʈ
T (TD) = (JDE –
2452545,0) / 36525
= 0.0199822080354286
Ʈ
= T / 10
= 0.0199822080354286
Penggunaan
tabel ini dengan rumus:
A * Cos(B + C *
Ʈ)
Konstanta B dan
C sudah tersaji dalam satuan Radian. Karena banyaknya suku koreksi tersebut,
penulis merasa tangel mencantumkannya di sini (gomena
sai..^_^). Tapi tak
apalah, namanya juga belajar. OK, kita hitung.
L0
|
A * COS(B + C * Ʈ)
|
1
|
=175347046 * COS(0 + 0 * Ʈ)
|
2
|
=3341656 * COS(4.6692568 + 6283.07585 * Ʈ)
|
3
|
=34894 * COS(4.6261 + 12566.1517 * Ʈ)
|
4
|
=3497 * COS(2.7441 + 5753.3849 * Ʈ)
|
5
|
=3418 * COS(2.8289 + 3.5231 * Ʈ)
|
6
|
=3136 * COS(3.6277 + 77713.7715 * Ʈ)
|
7
|
=2676 * COS(4.4181 + 7860.4194 * Ʈ)
|
8
|
=2343 * COS(6.1352 + 3930.2097 * Ʈ)
|
9
|
=1324 * COS(0.7425 + 11506.7698 * Ʈ)
|
10
|
=1273 * COS(2.0371 + 529.691 * Ʈ)
|
11
|
=1199 * COS(1.1096 + 1577.3435 * Ʈ)
|
12
|
=990 * COS(5.233 + 5884.927 * Ʈ)
|
13
|
=902 * COS(2.045 + 26.298 * Ʈ)
|
14
|
=857 * COS(3.508 + 398.149 * Ʈ)
|
15
|
=780 * COS(1.179 + 5223.694 * Ʈ)
|
16
|
=753 * COS(2.533 + 5507.553 * Ʈ)
|
17
|
=505 * COS(4.583 + 18849.228 * Ʈ)
|
18
|
=492 * COS(4.205 + 775.523 * Ʈ)
|
19
|
=357 * COS(2.92 + 0.067 * Ʈ)
|
20
|
=317 * COS(5.849 + 11790.629 * Ʈ)
|
21
|
=284 * COS(1.899 + 796.298 * Ʈ)
|
22
|
=271 * COS(0.315 + 10977.079 * Ʈ)
|
23
|
=243 * COS(0.345 + 5486.778 * Ʈ)
|
24
|
=206 * COS(4.806 + 2544.314 * Ʈ)
|
25
|
=205 * COS(1.869 + 5573.143 * Ʈ)
|
26
|
=202 * COS(2.458 + 6069.777 * Ʈ)
|
27
|
=156 * COS(0.833 + 213.299 * Ʈ)
|
28
|
=132 * COS(3.411 + 2942.463 * Ʈ)
|
29
|
=126 * COS(1.083 + 20.775 * Ʈ)
|
30
|
=115 * COS(0.645 + 0.98 * Ʈ)
|
31
|
=103 * COS(0.636 + 4694.003 * Ʈ)
|
32
|
=102 * COS(0.976 + 15720.839 * Ʈ)
|
33
|
=102 * COS(4.267 + 7.114 * Ʈ)
|
34
|
=99 * COS(6.21 + 2146.17 * Ʈ)
|
35
|
=98 * COS(0.68 + 155.42 * Ʈ)
|
36
|
=86 * COS(5.98 + 161000.69 * Ʈ)
|
37
|
=85 * COS(1.3 + 6275.96 * Ʈ)
|
38
|
=85 * COS(3.67 + 71430.7 * Ʈ)
|
39
|
=80 * COS(1.81 + 17260.15 * Ʈ)
|
40
|
=79 * COS(3.04 + 12036.46 * Ʈ)
|
41
|
=75 * COS(1.76 + 5088.63 * Ʈ)
|
42
|
=74 * COS(3.5 + 3154.69 * Ʈ)
|
43
|
=74 * COS(4.68 + 801.82 * Ʈ)
|
44
|
=70 * COS(0.83 + 9437.76 * Ʈ)
|
45
|
=62 * COS(3.98 + 8827.39 * Ʈ)
|
46
|
=61 * COS(1.82 + 7084.9 * Ʈ)
|
47
|
=57 * COS(2.78 + 6286.6 * Ʈ)
|
48
|
=56 * COS(4.39 + 14143.5 * Ʈ)
|
49
|
=56 * COS(3.47 + 6279.55 * Ʈ)
|
50
|
=52 * COS(0.19 + 12139.55 * Ʈ)
|
51
|
=52 * COS(1.33 + 1748.02 * Ʈ)
|
52
|
=51 * COS(0.28 + 5856.48 * Ʈ)
|
53
|
=49 * COS(0.49 + 1194.45 * Ʈ)
|
54
|
=41 * COS(5.37 + 8429.24 * Ʈ)
|
55
|
=41 * COS(2.4 + 19651.05 * Ʈ)
|
56
|
=39 * COS(6.17 + 10447.39 * Ʈ)
|
57
|
=37 * COS(6.04 + 10213.29 * Ʈ)
|
58
|
=37 * COS(2.57 + 1059.38 * Ʈ)
|
59
|
=36 * COS(1.71 + 2352.87 * Ʈ)
|
60
|
=36 * COS(1.78 + 6812.77 * Ʈ)
|
61
|
=33 * COS(0.59 + 17789.85 * Ʈ)
|
62
|
=30 * COS(0.44 + 83996.85 * Ʈ)
|
63
|
=30 * COS(2.74 + 1349.87 * Ʈ)
|
64
|
=25 * COS(3.16 + 4690.48 * Ʈ)
|
|
|
L1
|
A * COS(B + C * Ʈ)
|
1
|
= 628331966747 * Cos(0 + 0 * Ʈ)
|
2
|
= 206059 * Cos(2.678235 + 6283.07585 * Ʈ)
|
3
|
= 4303 * Cos(2.6351 + 12566.1517 * Ʈ)
|
4
|
= 425 * Cos(1.59 + 3.523 * Ʈ)
|
5
|
= 119 * Cos(5.796 + 26.298 * Ʈ)
|
6
|
= 109 * Cos(2.966 + 1577.344 * Ʈ)
|
7
|
= 93 * Cos(2.59 + 18849.23 * Ʈ)
|
8
|
= 72 * Cos(1.14 + 529.69 * Ʈ)
|
9
|
= 68 * Cos(1.87 + 398.15 * Ʈ)
|
10
|
= 67 * Cos(4.41 + 5507.55 * Ʈ)
|
11
|
= 59 * Cos(2.89 + 5223.69 * Ʈ)
|
12
|
= 56 * Cos(2.17 + 155.42 * Ʈ)
|
13
|
= 45 * Cos(0.4 + 796.3 * Ʈ)
|
14
|
= 36 * Cos(0.47 + 775.52 * Ʈ)
|
15
|
= 29 * Cos(2.65 + 7.11 * Ʈ)
|
16
|
= 21 * Cos(5.34 + 0.98 * Ʈ)
|
17
|
= 19 * Cos(1.85 + 5486.78 * Ʈ)
|
18
|
= 19 * Cos(4.97 + 213.3 * Ʈ)
|
19
|
= 17 * Cos(2.99 + 6275.96 * Ʈ)
|
20
|
= 16 * Cos(0.03 + 2544.31 * Ʈ)
|
21
|
= 16 * Cos(1.43 + 2146.17 * Ʈ)
|
22
|
= 15 * Cos(1.21 + 10977.08 * Ʈ)
|
23
|
= 12 * Cos(2.83 + 1748.02 * Ʈ)
|
24
|
= 12 * Cos(3.26 + 5088.63 * Ʈ)
|
25
|
= 12 * Cos(5.27 + 1194.45 * Ʈ)
|
26
|
= 12 * Cos(2.08 + 4694 * Ʈ)
|
27
|
= 11 * Cos(0.77 + 553.57 * Ʈ)
|
28
|
= 10 * Cos(1.3 + 6286.6 * Ʈ)
|
29
|
= 10 * Cos(4.24 + 1349.87 * Ʈ)
|
30
|
= 9 * Cos(2.7 + 242.73 * Ʈ)
|
31
|
= 9 * Cos(5.64 + 951.72 * Ʈ)
|
32
|
= 8 * Cos(5.3 + 2352.87 * Ʈ)
|
33
|
= 6 * Cos(2.65 + 9437.76 * Ʈ)
|
34
|
= 6 * Cos(4.67 + 4690.48 * Ʈ)
|
L2
|
A * COS(B + C * Ʈ)
|
1
|
= 52919 * Cos(0 + 0 * Ʈ)
|
2
|
= 8720 * Cos(1.0721 + 6283.0758 * Ʈ)
|
3
|
= 309 * Cos(0.867 + 12566.152 * Ʈ)
|
4
|
= 27 * Cos(0.05 + 3.52 * Ʈ)
|
5
|
= 16 * Cos(5.19 + 26.3 * Ʈ)
|
6
|
= 16 * Cos(3.68 + 155.42 * Ʈ)
|
7
|
= 10 * Cos(0.76 + 18849.23 * Ʈ)
|
8
|
= 9 * Cos(2.06 + 77713.77 * Ʈ)
|
9
|
= 7 * Cos(0.83 + 775.52 * Ʈ)
|
10
|
= 5 * Cos(4.66 + 1577.34 * Ʈ)
|
11
|
= 4 * Cos(1.03 + 7.11 * Ʈ)
|
12
|
= 4 * Cos(3.44 + 5573.14 * Ʈ)
|
13
|
= 3 * Cos(5.14 + 796.3 * Ʈ)
|
14
|
= 3 * Cos(6.05 + 5507.55 * Ʈ)
|
15
|
= 3 * Cos(1.19 + 242.73 * Ʈ)
|
16
|
= 3 * Cos(6.12 + 529.69 * Ʈ)
|
17
|
= 3 * Cos(0.31 + 398.15 * Ʈ)
|
18
|
= 3 * Cos(2.28 + 553.57 * Ʈ)
|
19
|
= 2 * Cos(4.38 + 5223.69 * Ʈ)
|
20
|
= 2 * Cos(3.75 + 0.98 * Ʈ)
|
L3
|
A * COS(B + C * Ʈ)
|
1
|
= 289 * Cos(5.844 + 6283.076 * Ʈ)
|
2
|
= 35 * Cos(0 + 0 * Ʈ)
|
3
|
= 17 * Cos(5.49 + 12566.15 * Ʈ)
|
4
|
= 3 * Cos(5.2 + 155.42 * Ʈ)
|
5
|
= 1 * Cos(4.72 + 3.52 * Ʈ)
|
6
|
= 1 * Cos(5.3 + 18849.23 * Ʈ)
|
7
|
= 1 * Cos(5.97 + 242.73 * Ʈ)
|
L4
|
A * COS(B + C * Ʈ)
|
1
|
= 114 * Cos(3.142 + 0 * Ʈ)
|
2
|
= 8 * Cos(4.13 + 6283.08 * Ʈ)
|
3
|
= 1 * Cos(3.84 + 12566.15 * Ʈ)
|
L5
|
A * COS(B + C * Ʈ)
|
1
|
= 1 * Cos(3.14 + 0 * Ʈ)
|
Setelah itu hitung hasil tiap
kelompok koreksi. Maka:
|
|
L0 = 174807309.018875
|
|
L1 = 628331790826.155
|
|
L2 = 58226.7548097241
|
|
L3 = 288.618643763524
|
|
L4 = -120.022334618986
|
|
L5 = -0.99999873172754
|
Untuk
mendapatkan bujur heliosentris Bumi menggunakan rumus:
L = (L0 + L1 *
Ʈ+ L2 * Ʈ2+ L3 * Ʈ3+ L4 * Ʈ4+ L5 * Ʈ5)
/ 108
L = 127.302638918321 radian
|
= 93.9039308976317 derajat
|
Hasil ini merupakan bujur ekliptika bumi diukur dari
matahari. Untuk mendapatkan bujur ekliptika matahari, hasil tersebut
ditambahkan 180o, karena posisi Bumi diukur dari matahari
merupakan kebalikan dari posisi Matahari diukur dari Bumi.
|
|
= L + 180
|
|
= 273.903930897632 derajat
|
|
Setelah itu menghitung koreksi untuk dengan rumus
|
|
t = 273.903930897632° -
0.09033”
|
|
= 273.903905805965 derajat
|
|
Hasil ini merupakan true geometric longitude matahari
|
|
Kita beralih ke tabel lintang heliosentris Bumi.
|
|
B0
|
A * COS (B + C *
Ʈ)
|
1
|
= 280 * Cos(3.199 + 84334.662 * Ʈ)
|
2
|
= 102 * Cos(5.422 + 5507.553 * Ʈ)
|
3
|
= 80 * Cos(3.88 + 5223.69 * Ʈ)
|
4
|
= 44 * Cos(3.7 + 2352.87 * Ʈ)
|
5
|
= 32 * Cos(4 + 1577.34 * Ʈ)
|
B1
|
A * COS (B + C *
Ʈ)
|
1
|
= 9 *
Cos(3.9 + 5507.55 * Ʈ)
|
2
|
= 6 * Cos(1.73 + 5223.69 * Ʈ)
|
B0 = -102.005477710248
|
|
B1 = 10.4776974612464
|
|
Untuk mendapatkan lintang geometrik Matahari menggunakan
rumus:
|
|
B = (B0 + B1 * Ʈ) / 108
|
|
= -1.01796110179845E-06 radian
|
|
= 0.209969549429538 detik busur
(hasilnya dinegatifkan)
|
|
=
-0.209969549429538
|
|
Setelah itu menghitung koreksi B dengan menghitung
|
|
koreksi bujur ekliptika matahari terlebih dahulu:
|
|
λ’
= - 1,397 * T - 0,00031 *T2
|
|
=
273.624767073429
|
|
baru kemudian hitung koreksi B dengan rumus:
|
|
∆B
= 0,03916 * (COS λ’ - SIN λ’)
|
|
= 0.0415574308293185
|
|
ß = B + ∆B
|
|
=
-0.168412118600219
|
|
|
|
Kita beralih ke tabel R jarak Bumi ke Matahari
|
|
R0
|
A * COS (B + C *
Ʈ)
|
1
|
= 100013989 * Cos(0
+ 0 * Ʈ)
|
2
|
= 1670700 *
Cos(3.0984635 + 6283.07585 * Ʈ)
|
3
|
= 13956 *
Cos(3.05525 + 12566.1517 * Ʈ)
|
4
|
= 3084 * Cos(5.1985
+ 77713.7715 * Ʈ)
|
5
|
= 1628 * Cos(1.1739
+ 5753.3849 * Ʈ)
|
6
|
= 1576 * Cos(2.8469
+ 7860.4194 * Ʈ)
|
7
|
= 925 * Cos(5.453 +
11506.77 * Ʈ)
|
8
|
= 542 * Cos(4.564 +
3930.21 * Ʈ)
|
9
|
= 472 * Cos(3.661 +
5884.927 * Ʈ)
|
10
|
= 346 * Cos(0.964 +
5507.553 * Ʈ)
|
11
|
= 329 * Cos(5.9 +
5223.694 * Ʈ)
|
12
|
= 307 * Cos(0.299 +
5573.143 * Ʈ)
|
13
|
= 243 * Cos(4.273 +
11790.629 * Ʈ)
|
14
|
= 212 * Cos(5.847 +
1577.344 * Ʈ)
|
15
|
= 186 * Cos(5.022 +
10977.079 * Ʈ)
|
16
|
= 175 * Cos(3.012 +
18849.228 * Ʈ)
|
17
|
= 110 * Cos(5.055 +
5486.778 * Ʈ)
|
18
|
= 98 * Cos(0.89 +
6069.78 * Ʈ)
|
19
|
= 86 * Cos(5.69 +
15720.84 * Ʈ)
|
20
|
= 86 * Cos(1.27 +
161000.69 * Ʈ)
|
21
|
= 65 * Cos(0.27 +
17260.15 * Ʈ)
|
22
|
= 63 * Cos(0.92 +
529.69 * Ʈ)
|
23
|
= 57 * Cos(2.01 +
83996.85 * Ʈ)
|
24
|
= 56 * Cos(5.24 +
71430.7 * Ʈ)
|
25
|
= 49 * Cos(3.25 +
2544.31 * Ʈ)
|
26
|
= 47 * Cos(2.58 +
775.52 * Ʈ)
|
27
|
= 45 * Cos(5.54 +
9437.76 * Ʈ)
|
28
|
= 43 * Cos(6.01 +
6275.96 * Ʈ)
|
29
|
= 39 * Cos(5.36 +
4694 * Ʈ)
|
30
|
= 38 * Cos(2.39 +
8827.39 * Ʈ)
|
31
|
= 37 * Cos(0.83 +
19651.05 * Ʈ)
|
32
|
= 37 * Cos(4.9 +
12139.55 * Ʈ)
|
33
|
= 36 * Cos(1.67 +
12036.46 * Ʈ)
|
34
|
= 35 * Cos(1.84 +
2942.46 * Ʈ)
|
35
|
= 33 * Cos(0.24 +
7084.9 * Ʈ)
|
36
|
= 32 * Cos(0.18 +
5088.63 * Ʈ)
|
37
|
= 32 * Cos(1.78 +
398.15 * Ʈ)
|
38
|
= 28 * Cos(1.21 +
6286.6 * Ʈ)
|
39
|
= 28 * Cos(1.9 +
6279.55 * Ʈ)
|
40
|
= 26 * Cos(4.59 +
10447.39 * Ʈ)
|
R1
|
A * COS (B + C *
Ʈ)
|
1
|
= 103019 *
Cos(1.10749 + 6283.07585 * Ʈ)
|
2
|
= 1721 * Cos(1.0644
+ 12566.1517 * Ʈ)
|
3
|
= 702 * Cos(3.142 +
0 * Ʈ)
|
4
|
= 32 * Cos(1.02 +
18849.23 * Ʈ)
|
5
|
= 31 * Cos(2.84 +
5507.55 * Ʈ)
|
6
|
= 25 * Cos(1.32 +
5223.69 * Ʈ)
|
7
|
= 18 * Cos(1.42 +
1577.34 * Ʈ)
|
8
|
= 10 * Cos(5.91 +
10977.08 * Ʈ)
|
9
|
= 9 * Cos(1.42 +
6275.96 * Ʈ)
|
10
|
= 9 * Cos(0.27 +
5486.78 * Ʈ)
|
R2
|
A * COS (B + C *
Ʈ)
|
= 4359 * Cos(5.7846
+ 6283.0758 * Ʈ)
|
|
= 124 * Cos(5.579 +
12566.152 * Ʈ)
|
|
= 12 * Cos(3.14 + 0
* Ʈ)
|
|
= 9 * Cos(3.63 +
77713.77 * Ʈ)
|
|
= 6 * Cos(1.87 +
5573.14 * Ʈ)
|
|
= 3 * Cos(5.47 +
18849.23 * Ʈ)
|
|
R3
|
A * COS (B + C *
Ʈ)
|
= 145 * Cos(4.273 +
6283.076 * Ʈ)
|
|
= 7 * Cos(3.92 +
12566.15 * Ʈ)
|
|
R4
|
A * COS (B + C *
Ʈ)
|
= 4 * Cos(2.56 +
6283.08 * Ʈ)
|
|
R0 = 98349222.1127999
|
|
R1 = 56741.5222021953
|
|
R2 = 3634.28293602439
|
|
R3 = -82.171127369278
|
|
R4 = -3.07096662876891
|
|
Maka true geocentric distance bisa diketahui dengan rumus:
|
|
R = (R0 + R1 * Ʈ+ R2 * Ʈ2+ R3 * Ʈ3+
R4 * Ʈ4) / 108
|
|
= 0.983503573841725
|
|
D. Menghitung Nutasi dan True Obliquity
|
|
Pertama menghitung mean obliquity dengan rumus:
|
|
U = T/100
|
|
= 0.00199822080354286
|
|
Εo = 23o 26’ 21,448” – ( 4680.93 * U –
1,55 * U2 +1999,25 * U3 -51,38 * U4
-249,67 * U5
|
|
-39,05 * U6
+7,12 * U7 +27,87 * U8 +5,79 * U9
+2,45 * U10) / 3600
|
|
= 23.4366928652695 derajat
|
Sebelum
menghitung True obliquity, kita harus menghitung koreksi ∆ε dengan menggunakan tabel di bawah ini yang merupakan
tabel terms of the 1980 IAU Theory of Nutations:
D
|
M
|
M'
|
F
|
Ω
|
∆ψ
Coefficients of The sine |
∆ε
Coefficients
of The Cosine
|
|||
0
|
0
|
0
|
0
|
1
|
-171996
|
-174,2
|
92025
|
8,9
|
|
-2
|
0
|
0
|
2
|
2
|
-13187
|
-1,6
|
5736
|
-3,1
|
|
0
|
0
|
0
|
2
|
2
|
-2274
|
-0,2
|
977
|
-0,5
|
|
0
|
0
|
0
|
0
|
2
|
2062
|
0,2
|
-895
|
0,5
|
|
0
|
1
|
0
|
0
|
0
|
1426
|
-3,4
|
54
|
-0,1
|
|
0
|
0
|
1
|
0
|
0
|
712
|
0,1
|
-7
|
0
|
|
-2
|
1
|
0
|
2
|
2
|
-517
|
1,2
|
224
|
-0,6
|
|
0
|
0
|
0
|
2
|
1
|
-386
|
-0,4
|
200
|
0
|
|
-2
|
-1
|
0
|
2
|
2
|
217
|
-0,5
|
129
|
-0,1
|
|
-2
|
0
|
0
|
2
|
1
|
129
|
0,1
|
-95
|
0,3
|
|
0
|
0
|
1
|
0
|
1
|
63
|
0,1
|
|||
0
|
0
|
-1
|
0
|
1
|
-58
|
-0,1
|
-70
|
0
|
|
0
|
2
|
0
|
0
|
0
|
17
|
-0,1
|
-53
|
0
|
|
-2
|
2
|
0
|
2
|
2
|
-16
|
0,1
|
|||
0
|
0
|
1
|
2
|
2
|
-301
|
0
|
-33
|
0
|
|
-2
|
0
|
1
|
0
|
0
|
-158
|
0
|
26
|
0
|
|
0
|
0
|
-1
|
2
|
2
|
123
|
0
|
32
|
0
|
|
2
|
0
|
0
|
0
|
0
|
63
|
0
|
27
|
0
|
|
2
|
0
|
-1
|
2
|
2
|
-59
|
0
|
|||
0
|
0
|
1
|
2
|
1
|
-51
|
0
|
-24
|
0
|
|
-2
|
0
|
2
|
0
|
0
|
48
|
0
|
16
|
0
|
|
0
|
0
|
-2
|
2
|
1
|
46
|
0
|
13
|
0
|
|
2
|
0
|
0
|
2
|
2
|
-38
|
0
|
|||
0
|
0
|
2
|
2
|
2
|
-31
|
0
|
-12
|
0
|
|
0
|
0
|
2
|
0
|
0
|
29
|
0
|
|||
-2
|
0
|
1
|
2
|
2
|
29
|
0
|
|||
0
|
0
|
0
|
2
|
0
|
26
|
0
|
-10
|
0
|
|
-2
|
0
|
0
|
2
|
0
|
-22
|
0
|
|||
0
|
0
|
-1
|
2
|
1
|
21
|
0
|
-8
|
0
|
|
2
|
0
|
-1
|
0
|
1
|
16
|
0
|
7
|
0
|
|
0
|
1
|
0
|
0
|
1
|
-15
|
0
|
9
|
0
|
|
-2
|
0
|
1
|
0
|
1
|
-13
|
0
|
7
|
0
|
|
0
|
-1
|
0
|
0
|
1
|
-12
|
0
|
6
|
0
|
|
0
|
0
|
2
|
-2
|
0
|
11
|
0
|
|||
2
|
0
|
-1
|
2
|
1
|
-10
|
0
|
5
|
0
|
|
2
|
0
|
1
|
2
|
2
|
-8
|
0
|
3
|
0
|
|
0
|
1
|
0
|
2
|
2
|
7
|
0
|
-3
|
0
|
|
-2
|
1
|
1
|
0
|
0
|
-7
|
0
|
|||
0
|
-1
|
0
|
2
|
2
|
-7
|
0
|
3
|
0
|
|
2
|
0
|
0
|
2
|
1
|
-7
|
0
|
3
|
0
|
|
2
|
0
|
1
|
0
|
0
|
6
|
0
|
|||
-2
|
0
|
2
|
2
|
2
|
6
|
0
|
-3
|
0
|
|
-2
|
0
|
1
|
2
|
1
|
6
|
0
|
-3
|
0
|
|
2
|
0
|
-2
|
0
|
1
|
-6
|
0
|
3
|
0
|
|
2
|
0
|
0
|
0
|
1
|
-6
|
0
|
3
|
0
|
|
0
|
-1
|
1
|
0
|
0
|
5
|
0
|
|||
-2
|
-1
|
0
|
2
|
1
|
-5
|
0
|
3
|
0
|
|
-2
|
0
|
0
|
0
|
1
|
-5
|
0
|
3
|
0
|
|
0
|
0
|
2
|
2
|
1
|
-5
|
0
|
3
|
0
|
|
-2
|
0
|
2
|
0
|
1
|
4
|
0
|
|||
-2
|
1
|
0
|
2
|
1
|
4
|
0
|
|||
0
|
0
|
1
|
-2
|
0
|
4
|
0
|
|||
-1
|
0
|
1
|
0
|
0
|
-4
|
0
|
|||
-2
|
1
|
0
|
0
|
0
|
-4
|
0
|
|||
1
|
0
|
0
|
0
|
0
|
-4
|
0
|
|||
0
|
0
|
1
|
2
|
0
|
3
|
0
|
|||
0
|
0
|
-2
|
2
|
2
|
-3
|
0
|
|||
-1
|
-1
|
1
|
0
|
0
|
-3
|
0
|
|||
0
|
1
|
1
|
0
|
0
|
-3
|
0
|
|||
0
|
-1
|
1
|
2
|
2
|
-3
|
0
|
|||
2
|
-1
|
-1
|
2
|
2
|
-3
|
0
|
|||
0
|
0
|
3
|
2
|
2
|
-3
|
0
|
|||
2
|
-1
|
0
|
2
|
2
|
-3
|
0
|
Sebelum
menghitung ∆ψ dan ∆ε dengan tabel
ini, kita menghitung dulu Multiple argumentsnya:
a) D = 297,85036 + 445267,11148*T - 0,0019142*T2
+ T3/189474
|
= 297,85036 +
445267,11148 * 0.199822080354286 - 0,0019142 *
|
0.199822080354286² +
0.199822080354286³ / 189474
|
=
352.050812887654 derajat
|
b) Mo = 357,52772 + 35999,05034*T - 0,0001603*T2
- T3/300000
|
=
357,52772 + 35999,05034 * 0.199822080354286 - 0,0001603 * 0.199822080354286²
|
- 0.199822080354286³ / 300000
|
= 350.932843290273 derajat
|
c) Mc = 134,96298 + 477198,867398*T +
0,0086972*T2 + T3/56250
|
=
134,96298 + 477198,867398 *0,126147653 + 0,0086972 *0,1261476532
|
+ 0,1261476533/56250
|
= 89.8337535885512 derajat
|
d) F = 93,27191 + 483202,017538*T - 0,0036825*T2
+ T3/327270
|
= 93,27191 +
483202,017538 *0,126147653 - 0,0036825 *0,1261476532
|
+
0,1261476533/327270
|
= 167.704138817659 derajat
|
e) Ωc
= 125,04452 - 1934,136261*T +
0,0020708*T2 + T3/450000
|
=
125,04452 - 1934,136261*0,126147653 + 0,0020708*0,1261476532
|
+ 0,1261476533/450000
|
= 98.5614713407414 derajat
|
Cara
menggunakan tabel di atas untuk menghitung ∆ψ adalah dengan format:
Coefficient *
Sin (Multiple Arguments)
Baris pertama = (-171996 + -174,2 *T) *
Sin (Ωc)
|
= -170113.820297936
|
Baris kedua = …..
|
Silahkan lanjutkan sendiri sampai Baris ke-63. Setelah itu
jumlahkan semuanya. Harus diingat,
|
bahwa hasil koreksi dengan
tabel ini menggunakan format 0,0001 detik,
|
maka setelah semua koreksi dijumlahkan dibagi 10000, maka
akan ketemu:
|
∆ψ = koreksi/10000/ 3600
|
=
-0.00469131429944207 derajat
|
Baik ! Kita lanjutkan:
|
Baris ke-2 = (-13187 + -1.6 * T ) * SIN (-2 * D + 2 * F +
2 * Ωc)
|
Baris ke-3 = (-2274 + -0.2 * T ) * SIN ( 2 * F + 2 * Ωc)
|
Baris ke-4 = (2062 + 0.2 * T ) * SIN ( 2 * Ωc)
|
Baris ke-5 = (1426 + -3.4 * T ) * SIN ( Mo)
|
Baris ke-6 = (712 + 0.1 * T ) * SIN ( Mc)
|
Baris ke-7 = (-517 + 1.2 * T ) * SIN (-2 * D + Mo + 2 * F
+ 2 * Ωc)
|
Baris ke-8 = (-386 + -0.4 * T ) * SIN ( 2 * F + Ωc)
|
Baris ke-9 = (217 + -0.5 * T ) * SIN (-2 * D + -Mo + 2 * F
+ 2 * Ωc)
|
Baris ke-10 = (129 + 0.1 * T ) * SIN (-2 * D + 2 * F + Ωc)
|
Baris ke-11 = 63 + 0.1 * T ) * SIN ( Mc + Ωc)
|
Baris ke-12 = (-58 + -0.1 * T ) * SIN ( -Mc + Ωc)
|
Baris ke-13 = (17 + -0.1 * T ) * SIN ( 2 * Mo)
|
Baris ke-14 = (-16 + 0.1 * T ) * SIN (-2 * D + 2 * Mo + 2
* F + 2 * Ωc)
|
Baris ke-15 = -301 * SIN ( Mc + 2 * F + 2 * Ωc)
|
Baris ke-16 = -158 * SIN (-2 * D + Mc)
|
Baris ke-17 = 123 * SIN ( -Mc + 2 * F + 2 * Ωc)
|
Baris ke-18 = 63 * SIN (2 * D)
|
Baris ke-19 = -59 * SIN (2 * D + -Mc + 2 * F + 2 * Ωc)
|
Baris ke-20 = -51 * SIN ( Mc + 2 * F + Ωc)
|
Baris ke-21 = 48 * SIN (-2 * D + 2 * Mc)
|
Baris ke-22 = 46 * SIN ( -2 * Mc + 2 * F + Ωc)
|
Baris ke-23 = -38 * SIN (2 * D + 2 * F + 2 * Ωc)
|
Baris ke-24 = -31 * SIN ( 2 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-25 = 29 * SIN ( 2 * Mc)
|
Baris ke-26 = 29 * SIN (-2 * D + Mc + 2 * F + 2 * Ωc)
|
Baris ke-27 = 26 * SIN ( 2 * F)
|
Baris ke-28 = -22 * SIN (-2 * D + 2 * F)
|
Baris ke-29 = 21 * SIN ( -Mc + 2 * F + Ωc)
|
Baris ke-30 = 16 * SIN (2 * D + -Mc + Ωc)
|
Baris ke-31 = -15 * SIN ( Mo + Ωc)
|
Baris ke-32 = -13 * SIN (-2 * D + Mc + Ωc)
|
Baris ke-33 = -12 * SIN ( -Mo + Ωc)
|
Baris ke-34 = 11 * SIN ( 2 * Mc + -2 * F)
|
Baris ke-35 = -10 * SIN (2 * D + -Mc + 2 * F + Ωc)
|
Baris ke-36 = -8 * SIN (2 * D + Mc + 2 * F + 2 * Ωc)
|
Baris ke-37 = 7 * SIN ( Mo + 2 * F + 2 * Ωc)
|
Baris ke-38 = -7 * SIN (-2 * D + Mo + Mc)
|
Baris ke-39 = -7 * SIN ( -Mo + 2 * F + 2 * Ωc)
|
Baris ke-40 = -7 * SIN (2 * D + 2 * F + Ωc)
|
Baris ke-41 = 6 * SIN (2 * D + Mc)
|
Baris ke-42 = 6 * SIN (-2 * D + 2 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-43 = 6 * SIN (-2 * D + Mc + 2 * F + Ωc)
|
Baris ke-44 = -6 * SIN (2 * D + -2 * Mc + Ωc)
|
Baris ke-45 = -6 * SIN (2 * D + Ωc)
|
Baris ke-46 = 5 * SIN ( -Mo + Mc)
|
Baris ke-47 = -5 * SIN (-2 * D + -Mo + 2 * F + Ωc)
|
Baris ke-48 = -5 * SIN (-2 * D + Ωc)
|
Baris ke-49 = -5 * SIN ( 2 * Mc + 2 * F + Ωc)
|
Baris ke-50 = 4 * SIN (-2 * D + 2 * Mc + Ωc)
|
Baris ke-51 = 4 * SIN (-2 * D + Mo + 2 * F + Ωc)
|
Baris ke-52 = 4 * SIN ( Mc + -2 * F)
|
Baris ke-53 = -4 * SIN (-D + Mc)
|
Baris ke-54 = -4 * SIN (-2 * D + Mo)
|
Baris ke-55 = -4 * SIN (D)
|
Baris ke-56 = 3 * SIN ( Mc + 2 * F)
|
Baris ke-57 = -3 * SIN ( -2 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-58 = -3 * SIN (-D + -Mo + Mc)
|
Baris ke-59 = -3 * SIN ( Mo + Mc)
|
Baris ke-60 = -3 * SIN ( -Mo + Mc + 2 * F + 2 * Ωc)
|
Baris ke-61 = -3 * SIN (2 * D + -Mo + -Mc + 2 * F + 2 *
Ωc)
|
Baris ke-62 = -3 * SIN ( 3 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-63 = -3 * SIN (2 * D + -Mo + 2 * F + 2 * Ωc)
|
∆ψ = koreksi/10000/ 3600
|
=
-0.00469131429944207 derajat
|
Sekarang menghitung koreksi ∆ε
|
Cara menggunakan tabel di atas untuk menghitung ∆ε adalah dengan format:
|
Coefficient * Cos (Multiple Arguments)
|
Baris pertama = (92025 + 8,9 *T) * Cos (Ωc)
|
= -13700.0650176995
|
Baris kedua =
|
Silahkan lanjutkan sendiri sampai Baris ke-49. Setelah itu
jumlahkan semuanya. Harus diingat,
|
bahwa hasil koreksi dengan tabel ini menggunakan format
0,0001 detik,
|
maka setelah semua koreksi dijumlahkan dibagi 10000, maka
akan ketemu:
|
∆ε = koreksi/10000/ 3600
|
=
-0.000551148761150953 derajat
|
Baik ! Kita lanjutkan:
|
Baris ke-2 = (5736 - 3.1 * T) * COS (-2 * D + 2 * F + 2 *
Ωc)
|
Baris ke-3 = (977 - 0.5 * T) * COS (2 * F + 2 * Ωc)
|
Baris ke-4 = (-895 + 0.5 * T) * COS (2 * Ωc)
|
Baris ke-5 = (54 - 0.1 * T) * COS ( Mo)
|
Baris ke-6 = (-7 ) * COS ( Mc)
|
Baris ke-7 = (224 - 0.6 * T) * COS (-2 * D + Mo + 2 * F + 2 * Ωc)
|
Baris ke-8 = (200 ) * COS (2 * F + 1 * Ωc)
|
Baris ke-9 = (129 - 0.1 * T) * COS (-2 * D - Mo + 2 * F +
2 * Ωc)
|
Baris ke-10 = (-95 + 0.3 * T) * COS (-2 * D + 2 * F + Ωc)
|
Baris ke-11 =
|
Baris ke-12 = -70 * COS ( - Mc + Ωc)
|
Baris ke-13 = -53 * COS (2 * Mo)
|
Baris ke-14 =
|
Baris ke-15 = -33 * COS ( + 1 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-16 = 26 * COS (-2 * D + Mc)
|
Baris ke-17 = 32 * COS ( - Mc + 2 * F + 2 * Ωc)
|
Baris ke-18 = 27 * COS (2 * D)
|
Baris ke-19 =
|
Baris ke-20 = -24 * COS ( Mc + 2 * F + Ωc)
|
Baris ke-21 = 16 * COS (-2 * D + 2 * Mc)
|
Baris ke-22 = 13 * COS ( - 2 * Mc + 2 * F + Ωc)
|
Baris ke-23 =
|
Baris ke-24 = -12 * COS (2 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-25 =
|
Baris ke-26 =
|
Baris ke-27 = -10 * COS (2 * F)
|
Baris ke-28 =
|
Baris ke-29 = -8 * COS ( - Mc + 2 * F + Ωc)
|
Baris ke-30 = 7 * COS (2 * D - Mc + Ωc)
|
Baris ke-31 = 9 * COS (Mo + Ωc)
|
Baris ke-32 = 7 * COS (-2 * D + Mc + Ωc)
|
Baris ke-33 = 6 * COS ( - Mo + Ωc)
|
Baris ke-34 =
|
Baris ke-35 = 5 * COS (2 * D - Mc + 2 * F + Ωc)
|
Baris ke-36 = 3 * COS (2 * D + Mc + 2 * F + 2 * Ωc)
|
Baris ke-37 = -3 * COS ( Mo + 2 * F + 2 * Ωc)
|
Baris ke-38 =
|
Baris ke-39 = 3 * COS ( - Mo + 2 * F + 2 * Ωc)
|
Baris ke-40 = 3 * COS (2 * D + 2 * F + Ωc)
|
Baris ke-41 =
|
Baris ke-42 = -3 * COS (-2 * D + 2 * Mc + 2 * F + 2 * Ωc)
|
Baris ke-43 = -3 * COS (-2 * D + Mc + 2 * F + Ωc)
|
Baris ke-44 = 3 * COS (2 * D - 2 * Mc + Ωc)
|
Baris ke-45 = 3 * COS (2 * D + Ωc)
|
Baris ke-46 =
|
Baris ke-47 = 3 * COS (-2 * D - Mo + 2 * F + Ωc)
|
Baris ke-48 = 3 * COS (-2 * D + Ωc)
|
Baris ke-49 = 3 * COS (2 * Mc + 2 * F + Ωc)
|
Sekali lagi …
|
∆ε = koreksi/10000/ 3600
|
=
-0.000551148761150953 derajat
|
Setelah itu menghitung true obliquity dengan rumus:
|
ε = Eo + ∆ε
|
= 23.4366928652695 -0.000551148761150953
|
= 23.4361417165084 derajat
|
= 23° 26’ 10.11”
|
E. Menghitung koreksi Aberasi, dengan rumus:
|
= -20,4898” / R
|
=
-0.00578707720286033
|
F. Menghitung Apparent Longitude Matahari, dengan rumus:
|
λ = t + ∆ψ + Aberasi
|
= 273.903905805965
+ -0.00469131429944207 -0.00578707720286033
|
=
273.893427414463
|
= 273° 53’
36.34”
|
G. Menghitung Apparent Right Ascension (α)
|
Sebelum menghitung α, koreksi dulu Longitudenya:
|
Ω
= 125,04 – 1934,136 *T
|
= 125,04 – 1934,136 * 0.199822080354286
|
= -261.443079208117 derajat
|
= 98.5569207918828 derajat
|
= 1.72014276844561 radian
|
Sebelum menghitung α, ε harus dikoreksi dengan rumus:
|
+0,00256 * Cos (Ω)
|
Maka:
|
E = ε +0,00256 * Cos (Ω)
|
=
23.4361417165084 + 0,00256 * Cos (98.5569207918828)
|
= 23.435760809293 derajat
|
kontrol (λ), untuk menyamakan kuadran ARA dengan ALM
|
kontrol (E * λ), untuk menyamakan kuadran ARA dengan ALM
|
Menghitung α dengan rumus:
|
Tan α = Cos (E) * Tan (λ)
|
=
Cos (23.435760809293) * Tan (273.893427414463)
|
=
-85.757738318015
|
nilai α harus sama kuadrannya dengan , maka ditambahkan
180 atau 360
|
α = 360 - 85.757738318015
|
= 274.242261681985 derajat
|
= 274° 14’
32.14”
|
H. Menghitung Apparent Declination
|
Sin = Sin E * Sin λ
|
=
SIN (23.435760809293) * SIN (273.893427414463)
|
= -23.3784522783986 derajat
|
= -23° 22’
42.43”
|
I. Menghitung
Equation of Time (Chapter 27)
Ketika Matahari
Semu melintasi meridian pengamat, maka pada saat itu terjadi waktu istiwa
pertengahan. Waktu istiwa hakiki yaitu ketika Matahari hakiki melintasi
meridian. Equation of time adalah perbedaan antara waktu hakiki dan waktu
pertengahan. Atau dengan kata lain, equation of time adalah perbedaan antara
Sudut waktu matahari hakiki dan semu.
Sebelum
menghitung equation of time, kita harus menghitung dulu rata-rata bujur
matahari dengan rumus:
Lo =
280,4664567 + 360007,6982779 *T - 0,03032028 *T 2 +T3 /49931
-T4 /
15299 - T5 /
1988000
|
= 7474.21516599369 derajat
|
= 274.215165993686 derajat
|
Eq = Lo - 0,0057183 – α + ∆ψ *COS ε
|
=
-0.0371182875839366 * 4 (dimenitkan)
|
= -0.14847 menit
|
= -0 menit 8.908 detik
|
kalau mau dijamkan, angka derajat tadi harus dibagi 15
|
Eq = -0.0371182875839366 / 15
|
= -0.0025 jam
|
= -0 menit 8.908 detik
|
J. Menghitung Semi Diameter Matahari (Chapter 53)
|
Sd = 15’ 59,63” / R
|
= 15’ 59,63” / 0.983503573841725
|
= 0° 16’ 15.726”
|
Kesimpulan:
Data Ephemeris
Matahari pada tanggal 26 Desember 2019 jam 7: 00: 00 WIB
Apparent
Longitude = 273° 53’ 36.34”
|
Apparent
Lattitude = -0.168”
|
True Obliquity = 23° 26’ 10.11”
|
Apparent Right Ascension = 274° 14’ 32.14”
|
Apparent Declination = -23° 22’ 42.43”
|
Equation of
Time
= -0 menit 8.908 detik
|
Semi
diameter = 0° 16’ 15.726”
|
True Geocentric Distance =
0.983503573841725 AU
|
Jika hasil
tersebut dibandingkan dengan Data Ephemeris Matahari pada Eph_Syauqin 1 yang Low
Accuracy:
Apparent
Longitude = 273° 53’
48.805”
|
Apparent
Lattitude = 0
|
True
Obliquity
= 23° 26’ 14.046”
|
Apparent Right Ascension =
274° 14’ 45.844”
|
Apparent Declination
= -23° 22’ 45.986”
|
Equation of
Time
= - 0 menit 10.409 detik
|
Semi
diameter
= 0° 16’ 15.713”
|
True Geocentric Distance = 0.983516506211852 (au)
|
File: Eph_Syauqin 1
Diposting:
Tidak ada komentar:
Posting Komentar