|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
关于这段程序,有点儿不明白啊,f=(0:N-3)/N*(fs/2);
这里fs/2是什么意思啊?
还有做出来的图怎么显示频率是比450还大
可是这两个正弦信号的频率分别是30和45啊
HHT.fig
(3.92 KB, 下载次数: 40)
%破凰版边界谱分析
clear;
fs=1000; %fs为采样频率;
N=1000; %采样点数
t=1/fs:1/fs:1;
y1=2*sin(60*pi*t);
y2=5*sin(90*pi*t);
y=[y1;y2;zeros(size(y1))]; %IMF集
%%%%%%%%%%%%%求边际谱
[A,fa,tt]=hhspectrum(y);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
b(k)=sum(E(k,:))*1/fs;
end
f=(0:N-3)/N*(fs/2);
plot(f,b);
xlabel('频率 / Hz');
ylabel('幅值');
[ 本帖最后由 serendipity_w 于 2007-7-30 16:40 编辑 ] |
|