声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1798|回复: 5

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

[复制链接]
发表于 2010-1-21 12:15 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
看了前边有关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=[y1;y2;zeros(size(y1))]; %IMF

%%%%%%%%%%%%%
求边际谱
[A,fa,tt]=hhspectrum(y);
NN=500;%设置频率点数,有什么依据么?
[E,tt1]=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('
幅值');

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2010-5-1 20:59 | 显示全部楼层
NN=500,是因为fs=1000,要满足采样定理
至于除以fs,感觉没有什么物理意义吧,但是既然都除以fs,对结果是没有任何影响的
发表于 2010-5-2 14:24 | 显示全部楼层
你把fs去掉,再看看图,纵轴能一样么?
发表于 2010-5-2 16:40 | 显示全部楼层

回复 沙发 h526306 的帖子

我有一个问题想请教您,你的意思是满足采样定理才将nn设置为500,但是这是采样点数,而不是数据中所含成分的最大频率吧?
发表于 2010-5-2 23:13 | 显示全部楼层
NN就是采样点数,与采样频率无关,楼主理解有误,信号中的频率必须满足采样定理。
发表于 2010-8-24 21:26 | 显示全部楼层
NN=500好像是人为想要显示的频率范围吧,跟什么采样频率无关的,至于乘以1/fs是因为积分嘛,离散不就是一个一个小矩形的面积嘛,这样理解应该是对的!

欢迎大家讨论
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-5-4 08:13 , Processed in 0.062379 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表