关于破凰网友的求边际谱的程序
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=; %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('幅值');
[ 本帖最后由 eight 于 2008-4-8 16:51 编辑 ] 你的效果我猜不好，因为边际效应比较严重，还会有模态混叠
页:
[1]