关于瞬时频率图形的解释

buzhiyao

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

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

yaoguo2

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

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

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

buzhiyao

原帖由 yaoguo2 于 2008-3-15 11:04 发表
可以看出，该信号的是各频谱分量很平滑，因为在你的图中纵坐标不是从0开始的，
建议你将第1点变为0，再看plot
当然低频成分比高频还是多一些。

对称是由于在FFT的性质，通过公式可以推导出（一般信号处理书上都 ...

\


先谢过:handshake ，
我其实是对瞬时频率不太懂，instfreq是和fft一回事吗？
另外要是和fft差不多的话，这个信号就一个频率，为什么出来 两个差不多高的峰值（FFT中的泄漏也不是这样的）？
还有就是横、纵坐标的单位是什么？你说的各频谱分量指的是什么，能不能帮我在图中标出来，同时给予解释？:loveliness:

[ 本帖最后由 eight 于 2008-3-17 13:10 编辑 ]

棒棒糖662

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

yaoguo2

很 Sorry 啊，对 HHT 我不太了解，结果却回答了这个问题。。。。。
刚才也查找了下资料，希望能有所帮助吧。
你可以看一下以前讨论的一个帖子
http://forum.vibunion.com/thread-41401-1-1.html
上面有解释和介绍为什么会出现波动的问题

主要的问题是信号离散化后非整数采样是会对HHT结果产生影响。

将你的程序改成下面这样
    x=[1:0.1:10-0.01];
    y=sin(2*x*pi);
    hy=hilbert(y);
    [f,t]=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

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

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