用凯斯西储大学轴承故障数据画边际谱不对怎么解

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

                               
登录/注册后可看大图

zhangnan3509

您可以考虑用别的处理方式观察是不是可以看出故障频率。如果可以那就是处理方式上有问题 要修正。

hht123

zhangnan3509 发表于 2015-4-7 09:11
您可以考虑用别的处理方式观察是不是可以看出故障频率。如果可以那就是处理方式上有问题 要修正。

我用仿真信号：
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

zhangnan3509 发表于 2015-4-7 09:11
您可以考虑用别的处理方式观察是不是可以看出故障频率。如果可以那就是处理方式上有问题 要修正。

实际数据的故障特征频率是载波频率还是调制频率，难道是调制频率？
我的仿真信号出来的是载波频率，那个3000是载波频率还是调制频率？

hht123

hht123 发表于 2015-4-7 12:29
实际数据的故障特征频率是载波频率还是调制频率，难道是调制频率？
我的仿真信号出来的是载波频率，那个 ...

自问自答一下：
论文"滚动轴承故障特征提取与应用研究"中提到：
轴承故障信号多为幅值调制信号，即故障特征频率调制轴承系统某阶固有频率。
包络解调是指对轴承故障振动信号取其幅值包络，去除掉信号中的高频载波频率，然后做幅值谱来查看解调出来的故障频率。
好像是：3000hz是系统固有的载波频率，故障特征频率是调制频率。