A. T. Walden, D. B. Percival and E. J. McCoy (1998), `Spectrum Estimation by Wavelet Thresholding of Multitaper Estimators,' IEEE Transactions on Signal Processing, 46, no. 12, pp. 3153-65.
Current methods for power spectrum estimation by wavelet thresholding use the empirical wavelet coefficients derived from the log periodogram. Unfortunately, the periodogram is a very poor estimate when the true spectrum has a high dynamic range and/or is rapidly varying. Also, because the distribution of the log periodogram is markedly non-Gaussian, complicated wavelet-dependent thresholding schemes are needed. These difficulties can be bypassed by starting with a multitaper spectrum estimator. The logarithm of this estimator is close to Gaussian distributed provided a moderate number (5 or more) of tapers are used. In contrast to the log periodogram, log multitaper estimates are not approximately pairwise uncorrelated at the Fourier frequencies, but the form of the correlation can be accurately and simply approximated. For scale-independent thresholding the correlation acts in accordance with the wavelet shrinkage paradigm to strongly suppress `noise spikes' while leaving informative coarse scale coefficients relatively unattenuated. This simple approach to spectrum estimation is demonstrated to work very well in practice. Additionally, the progression of the variance of wavelet coefficients with scale can be accurately calculated, allowing the use of scale-dependent thresholds. This more complex approach also works well in practice, but is not uniformly preferable to the scale-independent approach.
Go to next summary or home page for Don Percival