边际谱相关讨论的阶段性总结

zhlong

呵呵，这里讨论的是HHT，不是FFT呀！

williamzhaos

那位高手能不能详细的说一下求出时频图以后求边际谱的流程，我看过一些论文，好像说的都不是很详细。另外，在时频图中，是将所有的IMF分量的时频图做在一张图上好呢还是分开好，请赐教！谢谢！

zhangnan3509

咱们版里有做边际谱的程序就在置顶贴里，至于时频图我想还是把所有的成分作在一张图上好一些

williamzhaos

又是老兄，呵呵，谢谢了！
不过我看到论坛里的帖子上基本都是关于法国那个源程序来的，我现在是自己用一些简单的方法实现了hht，通过上午老兄和其他高手的帮助，我现在也已完成了时频图部分。而关于边际谱部分还是有些迷茫，不知老兄可否点拨一二，不胜感激！

[ 本帖最后由 williamzhaos 于 2007-11-8 14:36 编辑 ]

zhangnan3509

简单？你仅仅想做时频图吗？rilling的程序里已经有了那部分，好好看一下就能找到，你最后做出来时频图的程序和rilling的如出一辙。边际谱这个帖子里面已经有了程序，我也不太清楚该怎么补充

williamzhaos

我只需要做出时频图和边际谱就行了，主要是要讲程序加到我的仿真系统中去，实现谐波检测的作用！

Snikas

破凰的原程序关于频率点数设置处有这么一段：
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs*1/tspan;  
end
f=(0:N-3)/N*(fs/2);
通过将E矩阵对时间积分后得到bjp是一个行向量，其中每一列对应的是此频率在整个时间段对时间积分的结果，在画HHT边际谱时用到plot(f,bjp);语句，f应该和bjp是相同长度大小的向量吧，我认为边际谱的横坐标既然是频率那肯定f矩阵的长度应该对应的是E矩阵的行向量的数量（正好是HH谱中频率轴上的点数），不知道为什么楼主将信号在时域内的采样点数对应为边际谱的频率点数？我用楼主的频率设置方式画边际谱时总提示Vectors的数量不匹配，后来我将画边际谱的程序改成了：
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs*1/tspan;  
end
H=size(E,1);
f=(0:H-1)/H*(fs/2);
figure(3);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
就可以画出图来了。

lingyunzhi

学习一下，很有启发，谢谢楼主

nkdtxf

学习了，初学，还要一段时间才能体会啊

zt6130

非常之感谢啊。。。。。。。。。。

zhouli19860517

还不错啊

lixiaolong08

回复 7 # 破凰 的帖子

您真是一语点破我了思路，做边际谱时，幅值一直不对。原来我的采样时间只有0.16秒，所以边际谱幅值比实际幅值低了很多。

bzp小鹏

好东西啊，谢谢啊

飞驰吧少年

赞一个！很有用

奋斗奋斗

很实用。