케플러 방정식 풀이 (평균 → 이심 → 진근점이각)

케플러 방정식 풀이 쉽게 이해하기 — 개념과 원리

케플러 방정식은 타원 궤도에서 시간과 위치를 이어 주는 식입니다. 시간에 비례해 일정하게 증가하는 평균근점이각 M과 이심률 e를 넣으면, 이심근점이각 E, 근점에서 실제로 돌아간 각도인 진근점이각 ν, 그리고 긴반지름 a에 대한 거리의 비 r/a를 구해 줍니다. 행성, 위성, 혜성, 인공위성이 특정 시각에 어디에 있는지 예측할 때 씁니다.

계산 방법

• 케플러 방정식: M = E − e sin E (각도는 라디안). E를 식 하나로 풀어 낼 수 없으므로 E₀ = M에서 시작해 뉴턴 방법 Eₙ₊₁ = Eₙ − (Eₙ − e sin Eₙ − M) / (1 − e cos Eₙ)을 반복합니다.
• 진근점이각: ν = 2 atan2(√(1+e) sin(E/2), √(1−e) cos(E/2)). 반각 관계식 tan(ν/2) = √((1+e)/(1−e)) tan(E/2)를 사분면이 틀리지 않도록 쓴 형태입니다.
• 거리: r/a = 1 − e cos E.

평균근점이각은 도(°)로 입력하고, 이심률은 0(원) 이상 1 미만이어야 합니다.

원리

케플러 제2법칙에 따르면 행성은 같은 시간 동안 같은 넓이를 쓸고 지나갑니다. 평균근점이각 M = n(t − T)는 평균 운동 n = 2π/P로 일정하게 도는 가상의 천체가 가진 각도입니다. 이심근점이각은 타원을 둘러싼 보조원 위에서 잰 기하학적 보조 각도이며, 근점부터 쓸고 지나간 넓이를 계산하면 정확히 M = E − e sin E가 나옵니다. E를 알면 타원의 기하로 ν와 r을 바로 구할 수 있습니다. e < 1이면 도함수 1 − e cos E가 0이 되지 않아 뉴턴 방법이 매우 빨리 수렴합니다.

예제 풀이

1. 입력값: M = 45° = 0.785398 rad, e = 0.2입니다.
2. E₀ = M에서 시작하면 1회 후 약 54.4375°, 2회 후 약 54.30658°가 되고 3회째에 값이 더 바뀌지 않습니다.
3. 수렴값은 E = 54.30655693°(0.947828 rad)입니다. 검산: 0.947828 − 0.2 × 0.812150 = 0.785398 rad.
4. 진근점이각은 ν = 64.27172656°입니다.
5. 거리는 1 − 0.2 × cos 54.30656° = 0.8833103455로, 긴반지름의 88.3% 위치에 있습니다.

결과 해석

원 궤도(e = 0)에서는 세 각이 모두 같습니다. 타원에서는 근점 부근에서 가장 빠르게 움직이므로 근점에서 원점까지 가는 동안 ν가 M보다 앞서고, E는 그 사이에 놓입니다(45° → 54.3° → 64.3°). r/a는 근점에서 1 − e, 원점에서 1 + e입니다. 지구(e ≈ 0.0167)는 차이가 작지만 e가 0.9를 넘는 혜성은 차이가 매우 큽니다.

흔한 실수

• M = E − e sin E 안에서 도와 라디안을 섞는 경우. 이 식은 라디안에서만 성립합니다.
• atan2 대신 atan을 써서 ν의 사분면을 잘못 구하는 경우.
• 포물선·쌍곡선 궤도(e ≥ 1)에 그대로 적용하는 경우.

자주 묻는 질문

왜 수치 계산으로 풀어야 하나요?

E가 그대로 한 번, 사인 함수 안에 한 번 들어 있는 초월방정식이어서 대수적으로 분리할 수 없기 때문입니다.

시간으로 평균근점이각을 구하려면 어떻게 하나요?

M = 360° × (t − T)/P를 씁니다. T는 근점 통과 시각, P는 공전 주기입니다.

E와 ν는 무엇이 다른가요?

ν는 초점(태양)에서 본 실제 각도이고, E는 보조원의 중심에서 잰 계산용 중간 각도입니다.

출처: D. A. Vallado, Fundamentals of Astrodynamics and Applications, Microcosm Press · C. D. Murray and S. F. Dermott, Solar System Dynamics, Cambridge University Press