Rufter2007 发表于 2010-1-21 12:15

对HHT边际谱中幅值的一些疑问

看了前边有关hht边际谱的讨论，收获很大，但仍有几个问题想请教下高手。
（1）边际谱线的高度表征该频率点在所有时刻对应的幅值总和，我们设置的频率点数NN，实际上是将整个频域划分为NN+1个区间，将落入某个区间的所有瞬时频率点对应的幅值求和，即得该区域对应的边际谱幅值。那么，当NN取值越大，则所划区域范围越大，结果应是求和的幅值越大。那么，这个NN的值的设置，有什么依据么？如原论坛贴中，有取信号长度的。
（2）既然是对各区域的瞬时频率点对应幅值求和，那程序中，为什么还要对幅值求和结果进行“/fs”处理呢？
以上问题对应的程序如下（之前论坛的程序）：
clear;
fs=1000;    %fs为采样频率;

t=1/fs:1/fs:1;N=length(t);   %采样点数y1=5*sin(2*pi*241*t);
y2=3*sin(2*pi*73*t);
y=; %IMF集
%%%%%%%%%%%%%求边际谱
=hhspectrum(y);NN=500;%设置频率点数，有什么依据么？
=toimage(A,fa,tt,NN);E=flipud(E);%
for k=1:NN
    bjp(k)=sum(E(k,:))*1/fs;为什么要将幅值求和结果都除fs？
end
f=(0:NN-1)/NN*(fs/2);
plot(f,bjp);
xlabel('频率 / Hz');
ylabel('幅值');

h526306 发表于 2010-5-1 20:59

NN=500,是因为fs=1000，要满足采样定理
至于除以fs，感觉没有什么物理意义吧，但是既然都除以fs，对结果是没有任何影响的

dzkt 发表于 2010-5-2 14:24

你把fs去掉，再看看图，纵轴能一样么?

cboboc 发表于 2010-5-2 16:40

回复 沙发 h526306 的帖子

我有一个问题想请教您，你的意思是满足采样定理才将nn设置为500，但是这是采样点数，而不是数据中所含成分的最大频率吧？

dzkt 发表于 2010-5-2 23:13

NN就是采样点数，与采样频率无关，楼主理解有误，信号中的频率必须满足采样定理。

pandacute 发表于 2010-8-24 21:26

NN=500好像是人为想要显示的频率范围吧，跟什么采样频率无关的，至于乘以1/fs是因为积分嘛，离散不就是一个一个小矩形的面积嘛，这样理解应该是对的！

欢迎大家讨论{:{39}:}
页: [1]
查看完整版本: 对HHT边际谱中幅值的一些疑问