buzhiyao 发表于 2008-3-13 20:47

关于瞬时频率图形的解释

执行下面的代码得到附图
x=;
y=sin(2*x*pi);
hy=hilbert(y);
=instfreq(hy‘);
plot(t,f);
请教以下问题:
1、这个图该如何解释?就是由这个图可以看出原信号的何种性质
2、这个图为什么是对称的,而且,变化幅度这么大?
3、这种方法求瞬时频率是否正确,所求得的结果和HHT求得的瞬时频率是否相同?

[ 本帖最后由 buzhiyao 于 2008-3-13 20:51 编辑 ]

yaoguo2 发表于 2008-3-15 11:04

可以看出,该信号的是各频谱分量很平滑,因为在你的图中纵坐标不是从0开始的,
建议你将第1点变为0,再看plot
当然低频成分比高频还是多一些。

对称是由于在FFT的性质,通过公式可以推导出(一般信号处理书上都有解释)

瞬时频率这就要和你要求计算的精度有关了,

buzhiyao 发表于 2008-3-16 18:50

本帖最后由 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 编辑 ]

棒棒糖662 发表于 2008-3-18 15:14

关于楼上的问题我也一直很迷茫,希望能帮忙解释一下

yaoguo2 发表于 2008-3-18 16:41

很 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"也是可以的

buzhiyao 发表于 2008-3-18 20:13

先感谢yaogou2的指点:handshake ,和棒棒糖662的顶贴捧场:lol
yaogou2给的链接已看过,很有帮助
知道了几个问题:victory: :
1、该图的横坐标是时间,纵坐标是频率
2、hilbert变换是调用fft的,所以泄漏状况相似
3、周期信号的瞬时频率谱是条线

[ 本帖最后由 buzhiyao 于 2008-3-18 20:21 编辑 ]
页: [1]
查看完整版本: 关于瞬时频率图形的解释