Papers › A Fast Multitaper Power Spectrum Estimation in Nonuniformly Sampled Time Series
A Fast Multitaper Power Spectrum Estimation in Nonuniformly Sampled Time Series
Jie Cui, Benjamin H. Brinkmann, Gregory A. Worrell
Nonuniformly sampled signals are prevalent in real-world applications, but their power spectra estimation, usually from a finite number of samples of a single realization, presents a significant challenge. The optimal solution, which uses Bronez Generalized Prolate Spheroidal Sequence (GPSS), is computationally demanding and often impractical for large datasets. This paper describes a fast nonparametric method, the MultiTaper NonUniform Fast Fourier Transform (MTNUFFT), capable of estimating power spectra with lower computational burden. The method first derives a set of optimal tapers through cubic spline interpolation on a nominal analysis band. These tapers are subsequently shifted to other analysis bands using the NonUniform FFT (NUFFT). The estimated spectral power within the band is the average power at the outputs of the taper set. This algorithm eliminates the needs for time-consuming computation to solve the Generalized Eigenvalue Problem (GEP), thus reducing the computational load from O(N⁴) to O(N logN + N log(1/ϵ)), which is comparable with the NUFFT. The statistical properties of the estimator are assessed using Bronez GPSS theory, revealing that the bias of estimates and variance bound of the MTNUFFT estimator are identical to those of the optimal estimator. Furthermore, the degradation of bias bound may serve as a measure of the deviation from optimality. The performance of the estimator is evaluated using both simulation and real-world data, demonstrating its practical applicability. The code of the proposed fast algorithm is available on GitHub (https://github.com/jiecui/mtnufft).
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Tasks
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Methods
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections