有关实测齿轮箱信号的边际谱

王子期待

最近在用HHT做齿轮箱故障诊断的内容，实测信号采集正常和具有故障（局部损伤或磨损状态）的齿轮振动信号加速度信号。齿轮箱试验台振动数据单位为mm/s2，采样点数为4096点，采样频率为3838.77 Hz，齿轮啮合频率为307Hz，转轴I的频率为10Hz，转轴II的频率为7Hz。其中一组局部损伤状态的数据输入后程序如下：
clc;
N=4096;   %采样点数
n=1:N;
fs=3838.77;   %采样频率
t=n/fs;
imf = emd(x2);  %x2为实测的振动信号

[m,z]=size(imf);
figure(1);
subplot(m+1,1,1);
plot(t,x2)
ylabel('signal');   

for i=1:m-1
subplot(m+1,1,i+1);
plot(t,imf(i,:),'k')
ylabel(['imf',int2str(i)])
end

subplot(m+1,1,m+1);
plot(t,imf(m,:),'k')
ylabel(['r',int2str(m-1)])
xlabel('Time')



[A, fa, tt] = hhspectrum(imf);
[E, tt1] = toimage(A,fa,tt,length(tt));
disp_hhs(E,t,-20);                               %求 Hilbert谱


for k=1:size(E,1)
bjp(k) = sum(E(k,:))*1/fs;
end
f = (0:N-3)/N*(fs/2);
figure(2);
plot(f,bjp);                              % 求边际谱
xlabel('频率 / Hz');
ylabel('幅值');
pause;
plot(0:fs/N:fs/2-fs/N, s_fft(1:end/2))       % FFT变换
得出来的图形
C:\Users\Administrator.VOOOWPOF137KJAV\Desktop\IMF分解图.png
C:\Users\Administrator.VOOOWPOF137KJAV\Desktop\边际谱图.png
首先得出来的分解图怎么会这样呢，乱七八糟的，哪几个才是有效的信号啊，哪些是虚假的或者是噪声的啊。
还有边际谱怎么跟啮合频率或者他的倍频一点联系没有呢。。
求高人指点...

王子期待

图好像没贴好，我重新插下附件。

王子期待

大婶们快来帮忙解答哇。。小弟不胜感激！！

shuihai707

首先你得对你的数据很了解，其次对故障情况比较了解，最后直接将数据用EMD处理，有时效果并不好，所以你从图上看不出什么，建议相应的预处理下，比如去噪什么的。具体数据还要具体分析。

王子期待

shuihai707 发表于 2013-11-17 17:19

                               
登录/注册后可看大图

首先你得对你的数据很了解，其次对故障情况比较了解，最后直接将数据用EMD处理，有时效果并不好，所以你从图 ...

数据就是采集到的齿轮具有局部损伤的振动信号啊，上面的程序没有问题把。预处理的话用小波可以么，要不要消趋势项。

shuihai707

王子期待 发表于 2013-11-17 18:34

                               
登录/注册后可看大图

数据就是采集到的齿轮具有局部损伤的振动信号啊，上面的程序没有问题把。预处理的话用小波可以么，要不要 ...

有时候工程实践和理论分析对不上，就要多尝试了，尝试多了，发现效果好了，就知道问题出在哪里。可以去除趋势项试试效果。

mayaview

最近EMD很热啊，总感觉EMD不大靠谱，而且分解速度那叫一个慢啊！
LZ的问题感觉能解释的希望很小，分解的信号很难讲有什么含义，某些书上马后炮的解释一通，感觉就是胡扯。
PS： Matlab里现有那几个EMD的函数包感觉边界处理都不咋地，很头疼啊！

狠青春2012

本人小白，想问下边际谱即为包络谱吗