0
$\begingroup$

I'm trying to find the instantaneous frequencies for one of the Intrinsic Mode Functions (IMFs) I extracted from the original signal using Hilbert-Huang Transform. So, I got the following spectrum for that IMF:

enter image description here

However, when I tried to get the instantaneous frequencies for this IMF, I got the following (I'm using this library):

[x,_]=inst_freq(imfs[0])
x*fs
Out[17]: 
array([100., 200., 100., 100., 200., 200., 100., 100., 100., 100., 100.,
       100., 100., 100., 100., 100., 100., 100., 200., 100., 100., 100.,
       100., 100., 100., 100., 100., 200., 100., 100., 200., 100., 100.,
       100., 200., 100., 200., 100., 100., 100., 100., 200., 200., 200.,
       100., 200., 100., 100., 100., 100., 100., 100., 100., 100., 100.,
       100., 100., 100., 100., 200., 100., 100., 200., 100., 100., 100.,
       200., 100., 200., 200., 100., 100., 200., 100., 100., 100., 100.,
       100., 100., 200., 100., 100., 100., 100., 100., 100., 200., 200.,
       200., 100., 100., 200., 100., 100., 100., 100., 200., 100., 100.,
       100., 100., 100., 100., 200., 100., 100., 200., 200., 100., 100.,
       100., 100., 200., 100., 100., 100., 100., 100., 100., 200., 200.,
       100., 100., 100., 200., 100., 100., 100., 200., 200., 200., 100.,
       100., 200., 200., 100., 100., 100., 100., 200., 100., 100., 100.,
       100., 100., 100., 200., 200., 200., 100., 100., 200., 200., 200.,
       100., 100., 100., 100.])

From the instantaneous frequencies, one could imagine that the spectrum contains only the 100 and 200 Hz frequencies. Am I missing something here?

$\endgroup$
  • $\begingroup$ Do you think you can turn your comment into an answer and "self-accept" that answer? This will close this question gracefully and it will stop it from circulating the board as unanswered. All the best. $\endgroup$ – A_A Jun 20 '18 at 10:13
0
$\begingroup$

This is a major bug in this library. see http://github.com/jaidevd/pyhht/issues/43#issuecomment-398388106 I think it's better to use more robust libraries that are well-tested.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.