[논문 리뷰] Solving Kepler's equation CORDIC-like

M. Zechmeister|arXiv (Cornell University)|2018. 08. 21.

Scientific Research and Discoveries참고 문헌 12인용 수 5

한 줄 요약

이 논문은 타원 궤도와 쌍곡선 궤도의 케플러 방정식을 해결하기 위해 런타임 초월함수 호출 없이 CORDIC 유사 알고리즘을 제안한다. 사전에 계산된 각도 단계를 사용하여 이동적으로 근위도를 조정하고 부호 기반 방향 제어를 통해, 단일 정밀도로는 29회 반복, 이중 정밀도로는 55회 반복 내에 정확도를 확보하며, 뉴턴-랩슨 방법보다 하드웨어에 유리하고 수치적으로 안정적인 대안을 제공한다.

ABSTRACT

Context. Many algorithms to solve Kepler's equations require the evaluation of trigonometric or root functions. Aims. We present an algorithm to compute the eccentric anomaly and even its cosine and sine terms without usage of other transcendental functions at run-time. With slight modifications it is applicable for the hyperbolic case, too. Methods. Based on the idea of CORDIC, it requires only additions and multiplications and a short table. The table is independent of eccentricity and can be hardcoded. Its length depends on the desired precision. Results. The code is short. The convergence is linear for all mean anomalies and eccentricities e (including e = 1). As a stand-alone algorithm, single and double precision is obtained with 29 and 55 iterations, respectively. One half or two third of the iterations can be saved in combination with Newton's or Halley's method at the cost of one division.

연구 동기 및 목표

• 실행 시 sine, cosine 또는 기타 초월함수 평가를 피하는 케플러 방정식 해결 방법을 개발하기 위해.
• CORDIC 알고리즘의 반복적 회전 원리를, 특히 $ E - e\sin E = M $ 형태의 비선형적인 케플러 방정식에 적응시키기 위해.
• 모든 이심률 범위에서 수치적으로 안정되고 고정밀도인 해를 제공하여 $ e = 1 $ (직선 궤도) 포함하여 유효하게 하기 위해.
• 이심률에 독립적인 고정된 작은 룩업 테이블과 반복당 최소한의 연산을 통해 효율적인 하드웨어 구현을 가능하게 하기 위해.
• 뉴턴-랩슨 또는 할리 방법과 같은 고차수 반복 해법의 강력한 초기 추정치를 제공하여 총 반복 횟수를 줄이기 위해.

제안 방법

• 알고리즘은 사전에 계산된 기저 각도 $ \alpha_n = \pi / 2^n $ 를 사용하여 반복적 각도 덧셈을 수행하며, 초기값은 $ E_0 = 0 $ 으로 설정한다.
• 각 단계에서 현재 $ M_n = E_n - e\sin E_n $ 를 목표 평균 이면 $ M $ 과 비교하여 반복의 방향($ \sigma_{n+1} = \pm 1 $)을 결정한다.
• 직접 sine/cosine 평가를 피하기 위해, 삼각함수 덧셈 공식에서 소규모 각도 근사치를 사용하여 $ \sin E $ 와 $ \cos E $ 를 동시에 계산한다.
• 회전 논리를 수정하고 하이퍼볼릭 함수를 유사하게 적용하여, 이론적으로 쌍곡선 케이스로 확장한다.
• 효율성 향상을 위해, CORDIC 결과를 초기 추정치로 사용하여 뉴턴 또는 할리 방법과 조합할 수 있으며, 총 반복 횟수를 최대 2/3까지 줄일 수 있다.
• 각도와 삼각함수 값의 테이블은 사전에 계산되어 하드코딩되어 있어 이심률과 무관하며 임베디드 시스템에 적합하다.

실험 결과

연구 질문

• RQ1 실행 시 sine 또는 cosine 함수를 호출하지 않고도, 초월항 $ e\sin E $ 를 포함하는 케플러 방정식을 CORDIC 알고리즘으로 적응시킬 수 있는가?
• RQ2 모든 이심률 범위, 특히 $ e = 1 $ 에서도 수렴 특성과 정확도를 유지할 수 있는가?
• RQ3 소규모 각도에 대한 단순화된 테일러 전개를 사용하여 $ \sin E $ 와 $ \cos E $ 를 동시에 계산할 수 있는 정도는 어느 정도이며, 정확도를 유지할 수 있는가?
• RQ4 이 방법이 뉴턴-랩슨 또는 할리 방법과 같은 고차수 반복 해법의 효과적인 초기 추정치로 기능할 수 있는가? 전체 계산 비용을 줄일 수 있는가?
• RQ5 이 방법은 타원 궤도와 쌍곡선 궤도가 혼합된 상황에서의 통합 케플러 방정식으로 확장 가능한가?

주요 결과

• 알고리즘은 단일 정밀도로는 29회 반복, 이중 정밀도로는 55회 반복 내에 정확도를 확보하며, 모든 평균 이면과 이심률 범위에서 선형 수렴 특성을 보인다.
• 특히 $ e = 1 $ 이고 $ M = 0 $ 인 경우, 이중 정밀도 정확도는 약 $ 10^{-15} \sqrt{\frac{2}{1-e}} $ 수준으로 제한되며, 이는 직선 궤도 한계에서의 특이성을 반영한다.
• 이 알고리즘은 실행 시 초월함수 호출 없이도 $ E $, $ \sin E $, $ \cos E $ 를 동시에 계산할 수 있으며, 덧셈, 곱셈 및 사전에 계산된 테이블 외에는 추가 연산이 필요로 하지 않는다.
• CORDIC 결과를 초기 추정치로 사용할 경우, 뉴턴 또는 할리 방법을 적용하여 반복 횟수를 최대 2/3까지 줄일 수 있으며, 추가로 하나의 나눗셈만 필요하다.
• 알고리즘은 타원 궤도와 쌍곡선 궤도 모두에 적용 가능하며, $ e = 1 $ 의 경우 극한 경우로 간주하여 수치적 안정성을 유지한다.
• 반복당 계산 복잡도가 낮고 예측 가능한 구조를 지니므로, 현장 프로그래머블 게이트 어레이(FPGA)와 같은 하드웨어 구현에 매우 적합하다.