4、利用instfreq函数求取瞬时频率时出现的问题.docVIP

4、利用instfreq函数求取瞬时频率时出现的问题 最近利用HHT做突变检测，遇到了如下问题：首先根据定义求取的瞬时频率：s1=hilbert(imf(1,:));s2(:,1)=s1(1,:);instphase=angle(s2);unwrapinstphase=unwrap(instphase);instanglefrequency=diff(unwrapinstphase);realistfre=instanglefrequency/(2*pi);利用函数直接求取的瞬时频率int=instfreq(s2);figure(1);subplot(2,1,1);plot(realistfre);title(自己定义的瞬时频率);subplot(2,1,2);plot(inp);title(老外定义的瞬时频率);结果如下： 为什么会出现负频率现象，根据HHT理论，IMF满足窄带信号且是单分量信号，利用解析信号法求取瞬时频率应该不会有负频率，我查看了那个instfreq函数，是将相位的差分取了绝对值，当然不会出现负值，可是我不知道这是什么含义，而且它难道不会影响原信号的频谱特征吗？希望大家能够帮忙！ 答：看了几天的文献，终于搞明白了情况，下面和大家说一下，基本上做完EMD分解后，剩下的就是如何得到瞬时频率了，HUANG最一开始是由解析信