[천문학] 균시차 곡선 뿌시기

균시차 곡선은 아래와 같다.
균시차는 [시태양시-평균태양시]이며, [평균태양적경-시태양적경]이다.
 
예를 들어 균시차가 +10분이라는 것은
평균 태양시가 정오에 남중할 때 12시이므로 , 시태양시가 12시 10분이라는 것을 의미한다.
시태양이 먼저 남중했다고 해석이 가능하며 시태양이 먼저 남중하고 10분 후에 평균태양이 남중하게 된다.
 

 

타원 공전 궤도에 의한 균시차

 
오른쪽 그림에서 볼 수 있듯이 근일점(perihelion)에서 원일점(aphelion)으로 반시계 방향으로 태양이 공전하는 경우에는 시태양(real sun)이 평균 태양(mean sun)보다 빠르다. (=평균 태양이 시태양보다 느리다)
근일점에서부터 시태양이 빨랐다가 원일점에서 평균 태양에게 따라 잡히는데, 근일점과 원일점 사이에서 두 태양의 차이가 최대로 나타난다. 시태양이 빠른 경우, 균시차가 (-)값을 가지며 왼쪽 그래프에서 보듯이 아랫쪽으로 차이가 점점 커져 중간에서 최대로 나타나고 이후 점점 차이가 줄어들어 원일점에서 균시차가 0이 되는 것이다.
 

황도 경사에 의한 균시차

 
춘분점에서 하지점 사이의 경우, 평균 태양보다 시태양이 느리다. 시태양이 느린 경우 균시차가 (+)값을 가지며 이 효과는 춘분점과 하지점 사이에서 최대로 나타난다. 이후 시태양이 빨라지며 평균 태양을 하지점에서 따라 잡아 하지점에서 균시차는 0이 된다.
 

-----------------------------------------------------------------------------------------------------------------------------------------------------------------

 
균시차 곡선을 직접 그려보고 싶어 이것저것 검색하다가
github에서 균시차와 관련한 자료를 접하게 되었고, csv에 있는 균시차 자료를 발견하게 되었다.
날짜별 균시차 곡선을 구현해내었지만, 황도 경사 효과와 타원 궤도 효과를 따로 구현해내지는 못한다.

 
csv를 다운로드한 github에 균시차 코드 설명을 보니, 특정 경도에서 평균 태양시와 실제 태양시를 구할 수 있다고 한다.
그러면 csv에서 B열의 값이 시태양시-평균태양시(=균시차)이기 때문에 값들을 분해해 볼 수 있을 것 같다.
GitHub - behrouzz/eot: 시간 방정식

GitHub - behrouzz/eot: Equation of time

 
그런데 황도 경사 효과와 타원 궤도 효과는 어떻게 구현해 낼 수 있을까?