hht123 发表于 2015-4-5 20:23

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

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

hht123 发表于 2015-4-5 21:44

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

hht123 发表于 2015-4-6 15:34

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

jianjianliu076 发表于 2015-4-6 16:35

不懂帮顶{:{20}:}

hht123 发表于 2015-4-6 18:54

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 发表于 2015-4-8 14:30

答:需要解调

狠青春2012 发表于 2015-4-13 13:51

不懂,学习学习

woshiqiao 发表于 2015-5-22 10:42

不懂,前来学习。。。{:{39}:}

sqanna 发表于 2017-4-25 16:15

楼主 我也在做相关研究 能发我一份数据吗 15682070660@163.com
页: [1]
查看完整版本: 凯斯西储大学数据处理边际谱不对