Optimization of Weighting Factors for Multiple Window Spectrogram of Event-Related Potentials

This paper concerns the mean square error optimal weighting factors for multiple window spectrogram of different stationary and nonstationary processes. It is well known that the choice of multiple windows is important, but here we show that the weighting of the different multiple window spectrograms in the final average is as important to consider and that the equally averaged spectrogram is not mean square error optimal for non-stationary processes. The cost function for optimization is the normalized mean square error where the normalization factor is the multiple window spectrogram. This means that the unknown weighting factors will be present in the numerator as well as in the denominator. A quasi-Newton algorithm is used for the optimization. The optimization is compared for a number of well-known sets of multiple windows and common weighting factors and the results show that the number and the shape of the windows are important for a small mean square error. Multiple window spectrograms using these optimal weighting factors, from ElectroEncephaloGram data including steady-state visual evoked potentials, are shown as examples.

[1]  Richard G. Baraniuk,et al.  Multiple window time-frequency analysis , 1996, Proceedings of Third International Symposium on Time-Frequency and Time-Scale Analysis (TFTS-96).

[2]  D. Thomson,et al.  Spectrum estimation and harmonic analysis , 1982, Proceedings of the IEEE.

[3]  Maria Hansson,et al.  A multiple window method for estimation of peaked spectra , 1997, IEEE Trans. Signal Process..

[4]  R. Fletcher Practical Methods of Optimization , 1988 .

[5]  William J. Williams,et al.  Spectrogram decompositions of time-frequency distributions , 2001, Proceedings of the Sixth International Symposium on Signal Processing and its Applications (Cat.No.01EX467).

[6]  P. Welch The use of fast Fourier transform for the estimation of power spectra: A method based on time averaging over short, modified periodograms , 1967 .

[7]  Patrick J. Loughlin,et al.  Multiple window time-varying spectral analysis , 2001, IEEE Trans. Signal Process..

[8]  Maria Hansson,et al.  Kernels and Multiple Windows for Estimation of the Wigner-Ville Spectrum of Gaussian Locally Stationary Processes , 2007, IEEE Transactions on Signal Processing.

[9]  Donald B. Percival,et al.  The variance of multitaper spectrum estimates for real Gaussian processes , 1994, IEEE Trans. Signal Process..

[10]  Behrouz Farhang-Boroujeny Prolate Filters for Nonadaptive Multitaper Spectral Estimators With High Spectral Dynamic Range , 2008, IEEE Signal Processing Letters.

[11]  Magnus Lindgren,et al.  Multiple-window spectrogram of peaks due to transients in the electroencephalogram , 2001, IEEE Transactions on Biomedical Engineering.

[12]  Kurt S. Riedel,et al.  Minimum bias multiple taper spectral estimation , 2018, IEEE Trans. Signal Process..

[13]  Thomas P. Bronez,et al.  On the performance advantage of multitaper spectral analysis , 1992, IEEE Trans. Signal Process..

[14]  Jun Xiao,et al.  Multitaper Time-Frequency Reassignment for Nonstationary Spectrum Estimation and Chirp Enhancement , 2007, IEEE Transactions on Signal Processing.

[15]  Richard A. Silverman,et al.  Locally stationary random processes , 2018, IRE Trans. Inf. Theory.

[16]  Louis L. Scharf,et al.  Toeplitz and Hankel kernels for estimating time-varying spectra of discrete-time random processes , 2001, IEEE Trans. Signal Process..

[17]  Maria Hansson,et al.  Optimization of weighting factors for multiple window time-frequency analysis , 2009, 2009 17th European Signal Processing Conference.

[18]  M. Hansson Optimized weighted averaging of peak matched multiple window spectrum estimators , 1999, IEEE Trans. Signal Process..