关于瞬时频率图形的解释
执行下面的代码得到附图x=;
y=sin(2*x*pi);
hy=hilbert(y);
=instfreq(hy‘);
plot(t,f);
请教以下问题:
1、这个图该如何解释?就是由这个图可以看出原信号的何种性质
2、这个图为什么是对称的,而且,变化幅度这么大?
3、这种方法求瞬时频率是否正确,所求得的结果和HHT求得的瞬时频率是否相同?
[ 本帖最后由 buzhiyao 于 2008-3-13 20:51 编辑 ] 可以看出,该信号的是各频谱分量很平滑,因为在你的图中纵坐标不是从0开始的,
建议你将第1点变为0,再看plot
当然低频成分比高频还是多一些。
对称是由于在FFT的性质,通过公式可以推导出(一般信号处理书上都有解释)
瞬时频率这就要和你要求计算的精度有关了, 本帖最后由 wdhd 于 2016-9-9 13:26 编辑
原帖由 yaoguo2 于 2008-3-15 11:04 发表
可以看出,该信号的是各频谱分量很平滑,因为在你的图中纵坐标不是从0开始的,
建议你将第1点变为0,再看plot
当然低频成分比高频还是多一些。
对称是由于在FFT的性质,通过公式可以推导出(一般信号处理书上都 ...
\
先谢过:handshake ,
我其实是对瞬时频率不太懂,instfreq是和fft一回事吗?
另外要是和fft差不多的话,这个信号就一个频率,为什么出来 两个差不多高的峰值(FFT中的泄漏也不是这样的)?
还有就是横、纵坐标的单位是什么?你说的各频谱分量指的是什么,能不能帮我在图中标出来,同时给予解释?:loveliness:
[ 本帖最后由 eight 于 2008-3-17 13:10 编辑 ] 关于楼上的问题我也一直很迷茫,希望能帮忙解释一下 很 Sorry 啊,对 HHT 我不太了解,结果却回答了这个问题。。。。。
刚才也查找了下资料,希望能有所帮助吧。
你可以看一下以前讨论的一个帖子
http://forum.vibunion.com/thread-41401-1-1.html
上面有解释和介绍为什么会出现波动的问题
主要的问题是信号离散化后非整数采样是会对HHT结果产生影响。
将你的程序改成下面这样
x=;
y=sin(2*x*pi);
hy=hilbert(y);
=instfreq(hy');
plot(t,f);
就可以了
把"1:0.1:10"改为"1:0.1:10-0.01",0.01为你设定的一个采样点长度。因为你的步长,△t 为0.1,△t = 1/fs,假案你的fs = 10(为N=10的倍数),一个采样点长度就为0.01, fs=100,△t =0.001。
所以你改成"1:0.1:10-0.001"也是可以的 先感谢yaogou2的指点:handshake ,和棒棒糖662的顶贴捧场:lol
yaogou2给的链接已看过,很有帮助
知道了几个问题:victory: :
1、该图的横坐标是时间,纵坐标是频率
2、hilbert变换是调用fft的,所以泄漏状况相似
3、周期信号的瞬时频率谱是条线
[ 本帖最后由 buzhiyao 于 2008-3-18 20:21 编辑 ]
页:
[1]