tjucruiser 发表于 2008-9-25 09:48

求边际谱的问题，频率不对应

fs=1000;    %fs为采样频率;
N=1000;   %采样点数
t=1/fs:1/fs:1;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=; %IMF集
%%%%%%%%%%%%%求边际谱
=hhspectrum(y);
=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
    bjp(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');
求出的边际谱如下图，为什么得到的频率不是30和45呢？

吃书的老虎 发表于 2008-9-25 11:24

HHT中频率是指瞬时频率，我想问的是f=(0:N-3)/N*(fs/2)是指瞬时频率吗？好像不对吧，如果这样画边际谱肯定不对的。

news 发表于 2008-9-26 15:22

回复 楼主 tjucruiser 的帖子

E=flipud(E);
把这个语句去掉，新版的程序不需要这个！

tjucruiser 发表于 2008-9-26 15:27

多谢您了，我刚试了，果然是这样。
页: [1]
查看完整版本: 求边际谱的问题，频率不对应