[論文レビュー] The generalised Lomb-Scargle periodogram. A new formalism for the floating-mean and Keplerian periodograms

M. Zechmeister, M. Kürster|ArXiv.org|Jan 16, 2009

Stellar, planetary, and galactic studies参考文献 22被引用数 689

ひとこと要約

本論文は、不規則に分布したデータに対してオフセットおよび重みを含む完全な正弦波をフィッティングするための統一的で解析的なフレームワーク、一般化されたロンブ＝スカルグル周期ogram（GLS）を導入する。古典的なロンブ＝スカルグル法の欠点を克服し、周波数推定の正確性を向上させ、アリスティングを低減し、重み付き最小二乗法によるフィッティングによってスペクトル強度の決定を改善する。計算コストは元の手法と同等のままである。

ABSTRACT

The Lomb-Scargle periodogram is a common tool in the frequency analysis of unequally spaced data equivalent to least-squares fitting of sine waves. We give an analytic solution for the generalisation to a full sine wave fit, including an offset and weights ($χ^{2}$ fitting). Compared to the Lomb-Scargle periodogram, the generalisation is superior as it provides more accurate frequencies, is less susceptible to aliasing, and gives a much better determination of the spectral intensity. Only a few modifications are required for the computation and the computational effort is similar. Our approach brings together several related methods that can be found in the literature, viz. the date-compensated discrete Fourier transform, the floating-mean periodogram, and the "spectral significance" estimator used in the SigSpec program, for which we point out some equivalences. Furthermore, we present an algorithm that implements this generalisation for the evaluation of the Keplerian periodogram that searches for the period of the best-fitting Keplerian orbit to radial velocity data. The systematic and non-random algorithm is capable of detecting eccentric orbits, which is demonstrated by two examples and can be a useful tool in searches for the orbital periods of exoplanets.

研究の動機と目的

• 古典的なロンブ＝スカルグル周期ogramを、浮動平均と重み付きデータフィッティングを含む一般化した統一的で解析的な形式を構築すること。
• ゼロ平均を仮定し、測定誤差を無視する古典的手法の限界を克服すること。
• 日付補正離散フーリエ変換、浮動平均周期ogram、およびSigSpecにおけるスペクトル有意性推定器といった既存手法の等価性を統合・明確化すること。
• ラジアル速度データにおける軌道周期の検出を可能にするケプラー周期ogramへの形式の拡張。特に、離心軌道を含むもの。
• ラジアル速度時系列データから系外惑星の軌道周期を検出するための、計算的に効率的で体系的かつ非ランダムなアルゴリズムを提供すること。

提案手法

• オフセット $ c $ と重み付きデータを含む完全な正弦波モデル $ y(t) = a\cos\omega t + b\sin\omega t + c $ の最小二乗フィッティングに対する解析的解を導出する。
• 測定誤差を組み込むために、正規化された重み $ w_i = \frac{1}{\sigma_i^2} / \sum \frac{1}{\sigma_i^2} $ を用いた重み付き和を用いる。
• ラグランジュ乗数法を適用し、最小 $ \chi^2 $ 条件を導出し、パラメータ $ a $, $ b $, $ c $ に対する3つの一次方程式の連立式を得る。
• 連立方程式を解析的に解き、周期ogramのパワー $ p(\omega) $ を重み付き和 $ Y, C, S, CC, SS, CS $ および行列式 $ D = CC\cdot SS - CS^2 $ で表す。
• 周波数依存の位相シフト $ \tau $ を $ \tan 2\omega\tau = \frac{2CS}{CC - SS} $ で導入し、古典的手法の $ \hat{\tau} $ を一般化する。
• ケプラー軌道モデルをラジアル速度データにフィットすることで、ケプラー周期ogramへの形式の適用を行い、離心軌道の検出を可能にする。

実験結果

リサーチクエスチョン

• RQ1 ロンブ＝スカルグル周期ogramを、浮動平均と測定誤差の重みを含むように一般化しつつ、解析的効率を維持する方法は何か？
• RQ2 一般化されたロンブ＝スカルグル周期ogramと、日付補正離散フーリエ変換やSigSpecにおけるスペクトル有意性推定器といった既存手法との関係は何か？
• RQ3 一般化された形式は、不規則にサンプリングされた時系列における周波数検出の正確性を向上させ、アリスティングを低減できるか？
• RQ4 GLSは、特に離心軌道を持つ系外惑星の軌道周期検出にどのように応用できるか？
• RQ5 古典的手法および浮動平均周期ogramと比較して、GLSの計算コストと性能向上はどの程度か？

主な発見

• 一般化されたロンブ＝スカルグル周期ogram（GLS）は、オフセットと重み付きデータを含む完全な正弦波フィッティングを可能にすることで、古典的手法よりも優れた代替手段を提供し、より正確な周波数推定が可能になる。
• GLSの定式化により、アリスティングへの感受性が低減され、標準的手法よりもスペクトル強度の決定が改善される。
• GLSは、古典的手法のアルゴリズムに対するわずかな修正で実現され、元のロンブ＝スカルグル周期ogramと同等の計算効率を維持する。
• GLSの形式は、日付補正離散フーリエ変換、浮動平均周期ogram、およびSigSpecで用いられるスペクトル有意性推定器の間の等価性を統合・明確化する。
• GLSにより、離心軌道を含むラジアル速度データにおける周期の検出が可能な体系的かつ非ランダムなケプラー周期ogramの構築が可能となり、2つの実例でその有効性が示された。
• GLS周期ogramの解析的解は、重み付き最小二乗フィッティングを用いて導出され、周期ogramのパワーは $ p(\omega) = \frac{1}{YY \cdot D} \left[ SS \cdot YC^2 + CC \cdot YS^2 - 2CS \cdot YC \cdot YS \right] $ で表され、ここで $ D = CC\cdot SS - CS^2 $ である。