An algorithm for Morlet wavelet transform based on generalized discrete Fourier transform

摘要

Continuous wavelet transform (CWT) is a linear convolution of signal and wavelet function for a fixed scale. This paper studies the algorithm of CWT with Morlet wavelet as mother wavelet by using nonzero-padded linear convolution. The time domain filter, which is a non-causal filter, is the sample of wavelet function. By making generalized discrete Fourier transform (GDFT) and inverse transform for this filter, we can get a geometrically weighted periodic extension of the filter when evaluated outside its original support. From this extension of the time domain filter, we can get a causal filter. In this paper, GDFT-based algorithm for CWT, which has a more concise form than that of linear convolution proposed by Jorge Martinez, is constructed by using this causal filter. The analytic expression of the GDFT of this filter, which is essential for GDFT-based algorithm for CWT, is deduced in this paper. The numerical experiments show that the calculation results of GDFT-based algorithm are stable and reliable; the running speed of GDFT-based algorithm is faster than that of the other two algorithms studied in our previous work.