성도 관련 문서

계산 방법에 관한 정보

by 창환 posted Aug 11, 2009
?

단축키

Prev이전 문서

Next다음 문서

ESC닫기

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

 

<계산 방법에 관한 정보>


이 프로그램에서 각종 천문 계산에 사용한 방법은 다음과 같습니다.

  • 각종 좌표의 변환: 주로 Meeus의 방법을 썼지만 Peter Duffett-Smith의 방법을 쓴 것도 있고 직접 개발하여 쓴 것도 있습니다. 사람에 따라 각각의 계산 방법은 약간의 차이가 있습니다만 결과는 모두 같습니다.
  • 황도경사각의 계산: Meeus의 방법 가운데 정밀한 방법(식 22.3)을 썼습니다.
  • 세차 운동: Meeus의 방법 가운데 정밀한 방법(식 21.3)을 썼습니다.
  • 장동: 장동 운동은 Meeus의 책 144쪽에 소개되어 있는 간단한 방법을 쓰고 있습니다. 계산 정밀도는 황경 방향으로 0.5초각, 황위 방향으로 0.1초각 이하입니다.
  • DeltaT: NASA의 일식 페이지(Polynomial Expressions for Delta T)에 있는 식을 쓰고 있습니다. 조석으로 인한 가속도는 -26.0 arcsec/cy^2로 계산한 값입니다. 이에 대한 설명은 NASA의 일식 페이지에 상세히 나와 있습니다(http://eclipse.gsfc.nasa.gov/SEhelp/deltatpoly2004.html).
  • 천체의 출몰 시각: 출몰 시각 계산은 직접 개발한 방법을 쓰고 있습니다. 정밀도는 10초 이내입니다.
  • 항성시의 계산: Meeus가 소개한 방법을 쓰고 있습니다.
  • 24기 입기 시각 및 정삭망 시각, 일월식 예측: 전부 자체 개발한 방법을 쓰고 있습니다. 계산 정밀도는 5분 이내입니다(광행차 및 장동 보정이 반영되어 있으며 DeltaT의 불확실성에 따라 오차가 조금 더 커질 수 있습니다).
  • 음력 계산: 음력 계산 역시 자체 개발한 방법을 쓰고 있습니다. 유효 기간은 기원전 1350년에서 기원후 3000년입니다. 이 범위 밖에서도 계산은 가능하지만 오차가 커질 수 있습니다. 양력 날짜의 경우 1582년 10월 4일까지는 율리우스력으로, 그 이후에는 그레고리력으로 표시합니다. 음력 계산에서 주의할 점은 과거의 역 계산법을 전혀 반영하지 않았다는 점입니다. 예를 들어 과거에 쓰였던 역법 중 하나인 원가력은 평기법과 평삭법을 쓰고 있지만 이 프로그램에서는 정기법과 정삭법만 씁니다. 따라서 음력의 날짜에 하루나 이틀 정도의 차이가 생길 수 있으며 윤달 계산에서도 다른 결과가 나올 수 있습니다. 이 때에는 일진을 비교하여 날짜를 결정하는 것이 좋습니다.
  • 이슬람력 계산: 이슬람력의 계산은 Meeus가 소개한 방법을 쓰고 있습니다.
  • 혜성, 소행성의 위치 계산: 계산법은 Meeus의 책 33~35장에 쓰인 방법을 바탕으로 제작하였습니다. 궤도 요소의 분점 이동 역시 계산에 포함되어 있습니다.
  • 지구의 그림자 위치: 직접 만든 방법을 쓰고 있습니다. DE404에 기반을 두고 있습니다.
  • 달, 행성의 위치 계산: NASA/JPL에서 만든 천체력인 DE404에 기반을 두고 있습니다(http://www.moshier.net). 계산 정밀도는 내행성은 기원전 1350년에서 기원후 3000년까지 0.5초각 이내이며, 외행성은 기원전 3000년에서 기원후 3000년 사이에서 1초각 이내입니다(일심 황도좌표 기준, 성도 프로그램에서 실제로 출력하는 결과는 세차 및 장동 보정의 오차, DeltaT의 불확실성, 광행차 보정 생략으로 인한 오차가 포함되어 있으므로 오차 범위가 조금 더 큽니다.). 달은 기원전 3000년에서 기원후 3000년 사이에서 황위 방향으로 5초각, 황경 방향으로 7초각 이내의 정밀도로 계산합니다. 지구-달의 질량 중심(Earth-Moon Barycenter) 보정 계산에 포함되어 있습니다. 빛의 속도 한계로 인한 시간차는 1차 보정까지만 계산하고 있습니다. 좀 더 정밀한 보정은 추후 반영토록 하겠습니다.
  • 균시차: VSOP87(태양의 평균 황경)과 DE404(태양의 진황경)를 결합한 방법을 쓰고 있으며 오차는 1분 미만으로 추정됩니다.