A Fast Template Periodogram for Detecting Non-sinusoidal Fixed-shape Signals in Irregularly Sampled Time Series

01/29/2021
by   John Hoffman, et al.
0

Astrophysical time series often contain periodic signals. The large and growing volume of time series data from photometric surveys demands computationally efficient methods for detecting and characterizing such signals. The most efficient algorithms available for this purpose are those that exploit the 𝒪(Nlog N) scaling of the Fast Fourier Transform (FFT). However, these methods are not optimal for non-sinusoidal signal shapes. Template fits (or periodic matched filters) optimize sensitivity for a priori known signal shapes but at a significant computational cost. Current implementations of template periodograms scale as 𝒪(N_f N_obs), where N_f is the number of trial frequencies and N_obs is the number of lightcurve observations, and due to non-convexity, they do not guarantee the best fit at each trial frequency, which can lead to spurious results. In this work, we present a non-linear extension of the Lomb-Scargle periodogram to obtain a template-fitting algorithm that is both accurate (globally optimal solutions are obtained except in pathological cases) and computationally efficient (scaling as 𝒪(N_flog N_f) for a given template). The non-linear optimization of the template fit at each frequency is recast as a polynomial zero-finding problem, where the coefficients of the polynomial can be computed efficiently with the non-equispaced fast Fourier transform. We show that our method, which uses truncated Fourier series to approximate templates, is an order of magnitude faster than existing algorithms for small problems (N≲ 10 observations) and 2 orders of magnitude faster for long base-line time series with N_obs≳ 10^4 observations. An open-source implementation of the fast template periodogram is available at https://www.github.com/PrincetonUniversity/FastTemplatePeriodogram.

READ FULL TEXT

page 1

page 2

page 3

page 4

research
10/01/2021

pyFFS: A Python Library for Fast Fourier Series Computation and Interpolation with GPU Acceleration

Fourier transforms are an often necessary component in many computationa...
research
09/29/2022

Improving Power Spectral Estimation using Multitapering: Precise asteroseismic modeling of stars, exoplanets, and beyond

Asteroseismic time-series data have imprints of stellar oscillation mode...
research
09/18/2023

Scaling the time and Fourier domains to align periodically and their convolution

This note shows how to align a periodic signal with its the Fourier tran...
research
04/01/2020

From Fourier to Koopman: Spectral Methods for Long-term Time Series Prediction

We propose spectral methods for long-term forecasting of temporal signal...
research
02/28/2018

Multichannel Interpolation for Periodic Signals via FFT, Error Analysis and Image Scaling

This paper describes a new method for the multichannel interpolation of ...
research
12/08/2020

Fourier-Informed Knot Placement Schemes for B-Spline Approximation

Fitting B-splines to discrete data is especially challenging when the gi...
research
12/02/2021

Time-Series Estimation from Randomly Time-Warped Observations

We consider the problem of estimating a signal from its warped observati...

Please sign up or login with your details

Forgot password? Click here to reset