凯斯西储大学数据处理边际谱不对

hht123

我用凯斯西储大学的实际轴承故障数据进行处理，选用外环故障144.mat,但是边际谱图出来却是在3000hz左右，附件里是边际谱图，这是怎么回事，程序如下：
Nstd =0.05;
NE=100; % # of ensemble
numImf =10; % # of imf数量
runCEEMD =1;
maxSift =20;
typeSpline =2;
toModifyBC =2;
randType =2;
seedNo =0;
checksignal=1;
Imf= rcada_eemdgai(Sig,Nstd,NE,numImf); % run EEMD [3]


%% ********************* Marginal spectrum *******************论坛里的边际谱程序
N=length(t);
[A,fa,tt]=hhspectrum(Imf);
[E,tt1]=toimage(A,fa,tt,length(tt));
E=flipud(E);
for k=1:size(E,1)
%  bjp(k)=sum(E(k,:))*1/fs;
bjp(k)=sum(E(k,:));
end
f=(0:N-3)/N*(fs/2);
figure;
plot(f,bjp);
xlabel('Freq/ Hz');
ylabel('Amplitude');
Imf=Imf';
%give samplerate
samplerate=fs;
%give frequency-axis resolution for hilbert-spectrum
freqsol=500;
%give time-axis resolution for hilbert-spectrum
timesol=1000;
%give frequency-axis maximun value
hsp_fre1=fs/2;
au=size(Imf);
nIMF=au(2);
nPT=au(1);
totalt=nPT/samplerate;
Xlow=1/totalt;
[nta,taa,faa]=nnspa(Imf(1:nPT,1:nIMF), 0, totalt, freqsol, timesol, 0, hsp_fre1,0,totalt);
q=fspecial('gaussian', 7, 0.6);     
nsa=filter2(q, nta);
nsua=filter2(q, nsa);
%Plot the HESP and MS  
figure(3);
imagesc(taa,faa,nsua.^.5);
axis xy;set(gca,'FontSize',6);
title('Hilbert Amplitude Spectrum- for EEMD','FontSize',9);
colorbar;
xlabel('Time(sec)','FontSize',7,'VerticalAlignment','middle');
ylabel('Freq(hz)','FontSize',7);
ylim([Xlow hsp_fre1]);
clear  samplerate freqsol timesol hsp_fre1 au totalt Xlow Xhig nta taa faa nte tae fae

hht123

用仿真信号出来的结果很好啊，实际数据怎么出现这么大的问题？我已经去噪了啊。

hht123

顺便问一下，论文里单个imf分量的hilbert谱怎么画，论坛里画边际谱的程序输入必须是两行。

jianjianliu076

不懂帮顶

hht123

jianjianliu076 发表于 2015-4-6 16:35
不懂帮顶

谢谢！！！
我用仿真信号：
Sig1=(t>=0&t<=0.0853).*(1+sin(2*pi*2000*t)).*cos(2*pi*100*t+sin(2*pi*2000*t));
Sig2=(t>=0.0853&t<=0.1706).*(1+sin(2*pi*2500*t)).*cos(2*pi*200*t+sin(2*pi*2500*t));
x=Sig1+Sig2;
载波频率100hz和200hz都能在边际谱里显示出来啊，为什么实际数据就不行了？

hht123

答：需要解调

狠青春2012

不懂，学习学习

woshiqiao

不懂，前来学习。。。

sqanna

楼主 我也在做相关研究 能发我一份数据吗 15682070660@163.com