有关实测齿轮箱信号的边际谱
最近在用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为实测的振动信号
=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')
= hhspectrum(imf);
= 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
首先得出来的分解图怎么会这样呢，乱七八糟的，哪几个才是有效的信号啊，哪些是虚假的或者是噪声的啊。
还有边际谱怎么跟啮合频率或者他的倍频一点联系没有呢。。
求高人指点... 图好像没贴好，我重新插下附件。 大婶们快来帮忙解答哇。。小弟不胜感激！！ 首先你得对你的数据很了解，其次对故障情况比较了解，最后直接将数据用EMD处理，有时效果并不好，所以你从图上看不出什么，建议相应的预处理下，比如去噪什么的。具体数据还要具体分析。 shuihai707 发表于 2013-11-17 17:19 static/image/common/back.gif
首先你得对你的数据很了解，其次对故障情况比较了解，最后直接将数据用EMD处理，有时效果并不好，所以你从图 ...
数据就是采集到的齿轮具有局部损伤的振动信号啊，上面的程序没有问题把。预处理的话用小波可以么，要不要消趋势项。 王子期待 发表于 2013-11-17 18:34 static/image/common/back.gif
数据就是采集到的齿轮具有局部损伤的振动信号啊，上面的程序没有问题把。预处理的话用小波可以么，要不要 ...
有时候工程实践和理论分析对不上，就要多尝试了，尝试多了，发现效果好了，就知道问题出在哪里。可以去除趋势项试试效果。 最近EMD很热啊，总感觉EMD不大靠谱，而且分解速度那叫一个慢啊！
LZ的问题感觉能解释的希望很小，分解的信号很难讲有什么含义，某些书上马后炮的解释一通，感觉就是胡扯。
PS： Matlab里现有那几个EMD的函数包感觉边界处理都不咋地，很头疼啊！ 本人小白，想问下边际谱即为包络谱吗
页:
[1]