IMF的周期

yingziyingzi

各位大侠，我用HHT分解得到了一系列IMF，但不知道如何求其周期。
如以下IMF1和IMF2，时间t=1:5:3655，我用如下程序：
Y=fft(x);
N = length(Y);
Y(1) = [];
power = (sqrt(Y.*conj(Y))/(N/2)).^2;
nyquist = (N-1)/N;
freq = (1:N-1)./(N-1)*nyquist./dt;

period = 1./freq;
[mp,index] = max(power);
T=period(index);

得到两个序列的周期为178.83和5.1422，第二个周期明显不对，请问，应该怎么作呢。
谢谢！！