천문 계산
2009.08.12 02:31

IAU2000 장동 모델

조회 수 5567 댓글 0
?

단축키

Prev이전 문서

Next다음 문서

크게 작게 위로 아래로 댓글로 가기 인쇄 첨부
?

단축키

Prev이전 문서

Next다음 문서

크게 작게 위로 아래로 댓글로 가기 인쇄 첨부

IAU 2000 장동 모델

 

IAU2000 장동 모델은 국제천문연맹에서 정한 장동 계산식으로 MHB2000 모델에 기반을 두고 있습니다. 수 세기의 범위에서 마이크로초(microarcsecond, μas) 단위의 정밀도를 가지는 IAU2000A 모델과 이보다 조금 낮은 정밀도를 가지지만 항의 수를 대폭 줄여 계산을 편리하게 만든 IAU2000B 모델이 있습니다. IAU2000B 모델은 77개의 항을 가지며 1995년부터 2050년 사이의 기간에 한하여 1밀리초(milliarcsecond, mas) 이하의 정밀도를 가집니다. 원래의 IAU2000A 모델은 1365개의 항을 갖고 있으며 계산 정밀도와 계산항의 수를 제외하면 IAU2000A와 IAU2000B는 근본적으로 동일합니다.

 

이 글에서는 IAU2000A 모델의 계산법을 설명합니다. IAU2000A 모델에서 계산한 결과는 IAU2000 세차 모델, IAU2006 세차 모델과 함께 쓸 수 있습니다(IAU2006에서 쓰려면 약간의 보정이 필요합니다).

 

기존에 쓰던 IAU1980 장동 모형에서는 태양과 달의 영향만을 고려했지만, 새로운 IAU2000 장동 모형은 행성의 영향도 반영하고 있습니다. 그만큼 계산식이 길어지긴 했지만 정밀도가 훨씬 높아졌습니다.

 

(1) 기본 식

IAU2000 장동 모형에서 시간은 율리우스 적일로 나타난 TDB를 기본으로 하고 있지만 매우 높은 정밀도를 요구하지 않는다면 TT로 입력을 해도 충분합니다. 장동 계산에 사용하는 t는 다음과 같이 계산합니다.

 

t = (JDTT - 2451545.0) / 36525


 

계산에 필요한 14개의 값은 다음과 같이 계산합니다.

 

먼저 태양과 달의 운동에 관한 식입니다. 단위는 모두 초(arcsecond, as)입니다.

 

l = 달의 평균 근점이각 = 485868.249036 + 1717915923.2178 t + 31.8792 t2 + 0.051635 t3 - 0.00024470 t4

l' = 해의 평균 근점이각 = 1287104.79305 + 129596581.0481 t - 0.5532 t2 + 0.000136 t3 - 0.00001149 t4

F = L - Ω = 335779.526232 + 1739527262.8478 t - 12.7512 t2 - 0.001037 t3 + 0.00000417 t4

D = 해와 달의 평균 이각 = 1072260.70369 + 1602961601.2090 t - 6.3706 t2 + 0.006593 t3 - 0.00003169 t4

Ω = 달의 승교점 경도 = 450160.398036 - 6962890.5431 t + 7.4722 t2 + 0.007702 t3 - 0.00005939 t4

(L = 달의 평균 황경)


 

행성의 운동에 관한 식입니다.

 

Lm = 수성의 평균 황경 = 908103.259872 + 538101628.688982 t

Lv = 금성의 평균 황경 = 655127.283060 + 210664136.433548 t

Le = 지구의 평균 황경 = 361679.244588 + 129597742.283429 t

Lma = 화성의 평균 황경 = 1279558.798488 + 68905077.493988 t

Lj = 목성의 평균 황경 = 123665.467464 + 10925660.377991 t

Ls = 토성의 평균 황경 = 180278.799480 + 4399609.855732 t

Lu = 천왕성의 평균 황경 = 1130598.018396 + 1542481.193933 t

Ln = 해왕성의 평균 황경 = 1095655.195728 + 786550.320774 t


 

행성과 달의 위치에 관한 식 외에도 지구의 일반세차량을 계산할 필요가 있습니다.

 

Pa = 5028.8200 t+ 1.112022 t2


 

 

(2) 장동 운동 계산

기본 수치의 값을 계산했으면 따로 첨부해 놓은 표에 수록된 수치를 이용해 장동 값을 계산합니다.

황경 방향의 장동 운동량 Δφ와 황위 방향의 장동 운동량 Δε는

 

Δφ = ΣN i=0  (A + A't)sin(Argument) + (A'' + A'''t)cos(Argument)

Δε = ΣN i=0  (B + B't)sin(Argument) + (B'' + B'''t)cos(Argument)


 

로 계산할 수 있습니다.

여기에서 'Argument'로 표시해놓은 값은 첨부해 놓은 표에서 l, l', D, F, Om, Lm, Lv, Le, Lma, Lj, Ls, Lu, Ln, Pa 열 아레에 적힌 값을 (1)에서 계산한 해와 달, 행성의 위치 자료와 곱해서 구할 수 있습니다. 예를 들어 아래 표와 같이 되어 있다면 Argument는 이렇게 계산합니다.

 

행번호

l

l'

F

D

Om

Lm

Lv

Le

Lma

Lj

Ls

Lu

Ln

Pa

A

A'

B

B'

A''

A'''

B''

B'''

1

0

0

0

0

1

0

0

0

0

0

0

0

0

0

-17206.4161

-17.4666

9205.2331

0.9086

0.0036

0.0029

1.5377

0.0002

2

0

0

2

-2

2

0

0

0

0

0

0

0

0

0

-1317.0906

-0.1675

573.0336

-0.3015

-1.3696

0.0012

-0.4587

-0.0003

3

0

0

2

0

2

0

0

0

0

0

0

0

0

0

-227.6413

-0.0234

97.8459

-0.0485

0.2796

0.0002

0.1374

-0.0001

 

1번 행의 경우,

(Argument)1 = l × 0 + l' × 0 + F × 0 + D × 0 + Om × 1 + Lm × 0 + Lv × 0 + Le × 0 + Lma × 0 + Lj × 0 + Ls × 0 + Lu × 0 + Ln × 0 + Pa × 0

2번 행의 경우,

(Argument)2 = l × 0 + l' × 0 + F × 2 + D × -2 + Om × 2+ Lm × 0 + Lv × 0 + Le × 0 + Lma × 0 + Lj × 0 + Ls × 0 + Lu × 0 + Ln × 0 + Pa × 0

 

그리고 1번 행의 장동 운동량은 이렇게 계산합니다. 계산을 할 때에는 단위에 주의해야 합니다. (1)번 과정에서 계산한 값은 모두 초(arcsecond, as) 단위로 나타낸 식이지만 위의 표나 첨부한 수표는 모두 밀리초(milliarcsecond, mas) 단위로 값이 표시되어 있습니다.

Δφ1 = (-17206.4161 + -17.4666 × t) sin((Argument)1) + (0.0036 + 0.0039 × t) cos((Argument)1)

Δε1 = (9205.2331 + 0.9086 × t) cos((Argument)1) + (1.5377 + 0.0002 × t) sin((Argument)1)


 

Δφ와 Δε는 이들 각각의 계산 값을 모두 더한 값입니다.

Δφ = Δφ1 + Δφ2 + Δφ3 + Δφ4 + ...

Δε = Δε1 + Δε2 + Δε3 + Δε4 + ...


 

(3) IAU2006 세차 모형과 함께 쓰기

IAU2000 장동 모형은 IAU2006 세차 모형과 함께 쓰려면 간단한 변환을 통해 서로의 좌표를 맞추어 주어야 합니다.

변환식은 다음과 같습니다.

ΔφIAU2006 = ΔφIAU2000A + (0.4697×10-6 - 2.7774×10-6 t)ΔφIAU2000A

ΔεIAU2006 = ΔεIAU2000A - 2.7774×10-6 t ΔεIAU2000A


 

(4) 첨부 파일

첨부 파일은 IAU2000 장동 모형에 쓰이는 수표입니다. 세 가지 형태로 제공됩니다.

 

엑셀 문서(엑셀 97 이상에서 열 수 있습니다): nut.xls

ODF 문서(MS 오피스 및 오픈오피스 계열에서 열 수 있습니다): nut.ods

텍스트 문서(각각의 항목은 탭으로 구분되어 있습니다): nut.txt

 

 

 **이 계산식은 성도, 달력 프로그램에는 아직 적용되어 있지 않습니다.

?

  1. No Image

    [프로그램 소스]행성의 위치를 계산(DE404)하는 DLL 파일

    성도 프로그램에서 쓰고 있는 행성, 달 위치 계산 라이브러리입니다. 비주얼 C++ 2008에서 컴파일 할 수 있고, 프로젝트 파일을 새로 만들면 비주얼 C++ 6에서도 쓸 수 있습니다. (첨부 파일 안에 들어있는 DLL 파일은 비주얼 C++ 2008에서 컴파일 한 것으로 ...
    Date2009.08.12 Category천문 계산 By창환 Views8934
    Read More
  2. [책]역서

    (사진은 2008년판, 출처: 천문연구원) -제목: 역서 -편찬: 한국천문연구원 -발행사: 남산당 -언어: 한국어 한국천문연구원에서 펴낸 역서입니다. 현재 매년 새로운 역서를 발행하고 있습니다. 이 책에는 해당 해의 일력 자료, 해와 달, 행성, 혜성, 소행성의 ...
    Date2009.08.12 Category천문 계산 By창환 Views6511
    Read More
  3. [책]The Astronomical Almanac

    (사진은 2008년판, 출처는 구글 검색) -제목: Astronomical Almanac for the Year 2010 and Its Companion, The Astronomical Almanac Online: Data for Astronomy, Space Sciences, Geodesy, Surveying, Navigation, and other applications (제목 속의 연도...
    Date2009.08.12 Category천문 계산 By창환 Views4217
    Read More
  4. [책]Explanatory Supplement to the Astronomical Almanac

    (책 표지 사진은 구글에서 가져왔습니다.) -제목: Explanatory Supplement to the Astronomical Almanac -저자: P. Kenneth Seidelmann -출판사, 펴낸 날: University Science Books, Revised edition (2005년 8월 15일) -분량: 752쪽 -언어: 영어 -ISBN-10: ...
    Date2009.08.12 Category천문 계산 By창환 Views4955
    Read More
  5. [책]Astronomy with your personal computer

    (책 표지 사진은 구글에서 가져왔습니다.) -제목: Astronomy with your personal computer -저자: Peter Duffett-Smith -출판사, 펴낸 날: Cambridge University Press, 2nd edition (1990년 6월 29일) -분량: 272쪽 -언어: 영어 -ISBN-10: 052138995X -ISBN-...
    Date2009.08.12 Category천문 계산 By창환 Views4173
    Read More
  6. [책]Practical Astronomy with your Calculator

    (책 표지 사진은 구글에서 가져왔습니다.) -제목: Practical Astronomy with your Calculator -저자: Peter Duffett-Smith -출판사, 펴낸 날: Cambridge University Press, 3rd edition (1988년 3월 31일) -분량: 약 200쪽 -언어: 영어 -ISBN-10: 0521356997 ...
    Date2009.08.12 Category천문 계산 By창환 Views4060
    Read More
  7. [책]Astronomical Algorithms

    (책 표지 사진은 구글에서 가져왔습니다.) -제목: Astronomical Algorithms -저자: Jean Meeus -출판사, 펴낸 날: Willmann-Bell, 2nd edition (1998년 12월) -분량: 477쪽 -언어: 영어 -ISBN-10: 0943396611 -ISBN-13: 978-0943396613 이 책은 천문 계산에 ...
    Date2009.08.12 Category천문 계산 By창환 Views5976
    Read More
  8. [책]역법의 원리분석

    -제목: 역법의 원리분석 -저자: 이은성 -출판사, 펴낸 날: 정음사 (1985년 2월 1일 초판 발행) -분량: 444쪽 -언어: 한국어 -ISBN-13: 2005853002988 '역법의 원리분석'은 달력을 만드는 방법에 대하여 포괄적으로 설명하고 있는 책입니다. 달력 제작에 어떤 ...
    Date2009.08.12 Category천문 계산 By창환 Views10568
    Read More
  9. No Image

    IAU2000 장동 모델

    IAU 2000 장동 모델 IAU2000 장동 모델은 국제천문연맹에서 정한 장동 계산식으로 MHB2000 모델에 기반을 두고 있습니다. 수 세기의 범위에서 마이크로초(microarcsecond, μas) 단위의 정밀도를 가지는 IAU2000A 모델과 이보다 조금 낮은 정밀도를 가지지만 항...
    Date2009.08.12 Category천문 계산 By창환 Views5567
    Read More
  10. No Image

    IAU2006 세차 모델

    IAU2006 세차 모델 IAU에서는 2006년에 IAU2000에서 결정한 세차운동 식을 대체할 새로운 세차운동 식을 채택했습니다. 현재 IAU에서는 IAU2000A에서 결정한 세차운동식 대신 2006년에 채택한 새로운 식을 사용할 것을 권하고 있고, 미해군천문대에서는 2009년...
    Date2009.08.12 Category천문 계산 By창환 Views4375
    Read More
Board Pagination Prev 1 ... 12 13 14 15 16 17 18 19 20 21 Next
/ 21
Powered by XE