불규칙 샘플 데이터의 빠른 카이제곱 주기 탐색

초록

본 논문은 불규칙하게 샘플링되고 오차가 비균일한 시계열 데이터에서, 원하는 고조파 수만큼의 조화 성분을 포함한 주기 신호를 효율적으로 찾기 위한 Fast χ² 알고리즘을 제안한다. FFT를 이용해 χ² 최소값을 주파수별로 계산함으로써 O(N log N)의 시간 복잡도를 달성하고, 기존 Lomb‑Scargle 방식보다 통계적 효율성과 계산 속도에서 우수함을 보인다. 구현 코드는 논문에 함께 제공된다.

상세 분석

Fast χ² 알고리즘은 전통적인 주기 탐색 방법이 갖는 두 가지 근본적인 한계를 동시에 해결한다. 첫째, 불규칙 샘플링과 비균일 오차를 가진 데이터에 대해 가중 최소제곱(Weighted Least Squares) 형태의 χ² 함수를 정의하고, 이를 주파수와 고조파 차수(K)별로 최소화한다. 기존 Lomb‑Scargle 방법은 기본적으로 단일 사인·코사인 모델에 국한되며, 가중치 적용이 복잡하거나 고조파를 직접 포함하기 어렵다. 반면 Fast χ²는 모델을
( m(t)=\sum_{k=1}^{K}