Improved Methods for Power Spectrum Modelling of Red Noise
Abstract
We discuss techniques of modelling the power spectrum of variability in cases where the spectral power is continuous but diverges strongly to low frequencies (socalled `red noise'), as is seen, for example, in the erratic variability of active galaxies and Xray binaries. First we review the sampling properties of the periodogram and traditional smoothed periodogram estimates of the power spectral density function. Such estimates are biased, are of unknown variance, and have a strongly nonGaussian distribution, and so are inappropriate for a leastsquares `goodnessoffit' test. We suggest a new method based on grouping estimates of the logarithm of spectral power density. We show that these estimates are of known variance and require 23 times less smoothing to have distributions close to Gaussian. For spectral density functions that diverge exactly as a power law, these estimates are unbiased, and for any strongly diverging power spectrum will be less biased than traditional estimates. These estimates are therefore much superior for the purposes of fitting analytical models to power spectra and of assessing these models.
 Publication:

Monthly Notices of the Royal Astronomical Society
 Pub Date:
 April 1993
 DOI:
 10.1093/mnras/261.3.612
 Bibcode:
 1993MNRAS.261..612P
 Keywords:

 methods: data analysis