|
楼主 |
发表于 2008-10-30 20:34
|
显示全部楼层
该信号的HHT谱图和边际谱如附件所示,信号的取样率为50hz,N=36001,为什么在20-25hz这一部分,HHT谱图和边际谱图不对应呢,出现了没有值,这是为什么呀?sig为信号。
t=sig(:,1);
pm=sig(:,2);
N=length(t);
fs=50;
T=N*1/fs;
% imf = emd(pm);
[imf,ort,nbits] = emd(pm);
emd_visu(pm,t,imf);
[A,f,tt] = hhspectrum(imf(1:end-1,:));
[im,tt,ff] = toimage(A,f,tt);
figure;
for k=1:size(im,2);
bjp1(k)=sum(im(:,k))*1/fs*1/T;
end
NN=size(im,2);
f=(0:N-3)/N*(fs/2);
plot(f,bjp1);
title('边际谱');
xlabel('频率HZ');
ylabel('幅值'); |
-
HHT谱图
-
边际谱图
|