hht123 发表于 2015-4-6 22:17

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

我用凯斯西储大学的实际轴承故障数据进行处理,选用外环故障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


%% ********************* Marginal spectrum *******************论坛里的边际谱程序
N=length(t);
=hhspectrum(Imf);
=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;
=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();
clearsamplerate freqsol timesol hsp_fre1 au totalt Xlow Xhig nta taa faa nte tae fae


http://forum.vibunion.com/forum.php?mod=attachment&aid=Njg0OTN8MDI5ZGUxMjd8MTQyODMyOTgwNXwyMjMyNDZ8MTM2MDYz

zhangnan3509 发表于 2015-4-7 09:11

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

hht123 发表于 2015-4-7 10:09

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 发表于 2015-4-7 12:29

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

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

hht123 发表于 2015-4-7 17:06

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

自问自答一下:
论文"滚动轴承故障特征提取与应用研究"中提到:
轴承故障信号多为幅值调制信号,即故障特征频率调制轴承系统某阶固有频率。
包络解调是指对轴承故障振动信号取其幅值包络,去除掉信号中的高频载波频率,然后做幅值谱来查看解调出来的故障频率。
好像是:3000hz是系统固有的载波频率,故障特征频率是调制频率。
页: [1]
查看完整版本: 用凯斯西储大学轴承故障数据画边际谱不对怎么解