24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tony1230

银虫 (小有名气)

应助: 5 (幼儿园)
金币: 511.9
帖子: 72
在线: 73.3小时
虫号: 1793983

[交流] 求高手解释HHT变换求频率等问题

利用HHT变换分解信号以及求瞬时频率等，但是为何不同IMF分量求出来的瞬时频率是一样的呢，搞不清楚。下面是matlab代码求频率的，高手看看是不是存在什么问题。
hx = hilbert(C);
xr=real(hx);xi=imag(hx);
sz=sqrt(xr.^2 + xi.^2);
sx=atan(xi./xr);
dt = diff(t);
dx = diff(sx);
sp = dx./dt;
figure
plot(t(1:N-1),sp)
title(['InstantFreq imf ' num2str(i)])
end

%计算HHT时频谱和边际谱
[A,infeq,tt] = hhspectrum(C,t,1,1);
[im,tt1,Cenf] = toimage(A,infeq,tt,length(tt));
disp_hhs(im,tt1) %二维图显示HHT时频谱，E是求得的HHT谱

这里由两种瞬时频率，一种是根据HHT定义的根据相位时间导数求的（sp），还有下面的根据matlab内部函数instfreq求的瞬时频率（infeq）。sp有正有负，infeq严格在0-0.5之间。不知道哪种方法是合理的，本人不是学这专业的，最近用到这个工具，但是有些搞不懂，请高手指点。

另外，上面还有一个求解中心频率的（Cenf），但是出来的结果是一个从0-0.5线性序列。而且出来只有一个序列，看EMD分解可看到有5个IMF分量。看到有文章说，中心频率就是每个IMF瞬时频率的平均值，如果是这样算的话倒是可以理解了。

下面三幅图分别是EMD分解的结果，HHT频谱图，还有个边际谱图。虽然做了但是看不懂太深的东西。请大家帮忙分析指点。

marginal_sp.png

EMD.png

HHT_spectrum.png

[ Last edited by tony1230 on 2013-5-6 at 13:42 ]

回复此楼