Anda setuju dengan Blog ini ? Kalau ya, silahkan Klik !

Sabtu, 12 Oktober 2019

MENGHITUNG DATA EPHEMERIS MATAHARI DENGAN ALGHORITMA JEAN MEEUS (HIGH ACCURACY)


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