I am reading this paper by Han et al. (2014). In this article, the authors extract detailed information from geomagnetic data sampled at 1Hz using a daubechies 5 (db5) wavelet: they reconstruct the signal using the coefficients from the 6th level only. The authors explain that the central frequency of the resulting reconstructed signal is 0.01Hz.
I have reprocessed similar data following the same steps and have convinced myself that the central frequency of the signal reconstructed using only coefficients from the 6th level is indeed around 0.01Hz. Below are the PSD plots for both the original and reconstructed signals, generated using Python:
But is there a way to analytically determine the central frequency of the reconstructed signal?
- On one hand, I have found explanations stating that, considering a SR of 1Hz, the central frequency of the signals reconstructed using the 6th level should lie within [0.5, 1]/2**6, i.e. [0.0078125, 0.015625]Hz, which is the case. Is this the correct answer? If this was the case, it would mean that the central frequency of the analyzed signal would be "wavelet-independent"...
- On the other hand, I know about the formula from Abry, 1997, stating that the pseudo-frequency of a wavelet at a given scale Fa is the central frequency of the wavelet Fc divided by the scale a: Fa=Fc/a. But my undestanding is that this formula yields the pseudo frequency of the wavelet, not the central frequency of the analyzed signal!
Any clarification would be appreciated. Thanks!