声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1887|回复: 1

[综合讨论] 实际轴承故障信号倒谱分析

[复制链接]
发表于 2012-6-9 22:43 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
有人用实际故障信号做过倒谱分析吗,为什么我用倒谱分析,结果出来不理想。我看很多文献都说用倒谱分析能够提取出故障特征频率,但是我试了好久了 还是搞不出来。。。
%倒谱:
clc
clear all
fs=12000;
N=10000;
nfft=8192;
t=0:1/fs:(N-1)/fs;
nn=1:N/4;
load('F:\美国西储大学轴承数据\12k Drive End Bearing Fault Data\IR007_0.mat');
d=X105_DE_time(1:N,1);
figure(1)
subplot(211)
f=0:fs/nfft:fs/2-fs/nfft;
y=abs(fft(d,nfft));          %FFT
plot(f(nn),y(nn))
xlabel('Frequency');
ylabel('|FFT)|')
subplot(212)
z=rceps(d);
plot(t(nn),abs(z(nn)));
title('倒频谱')
xlabel('时间(s)');
ylabel('幅值(dB)');
set(gca, 'YGrid', 'on');
xlim([0 0.02])
imf=emd(d);%EMD分解
%对每个imf分量进行倒谱分析
figure(2)
z=rceps(imf(1,:));
plot(t(nn),abs(z(nn)));
title('imf1倒频谱')
xlabel('时间(s)');
ylabel('幅值(dB)');
set(gca, 'YGrid', 'on');
%对imf2分量进行倒谱分析
figure(3)
z=rceps(imf(2,:));
plot(t(nn),abs(z(nn)));
title('imf2倒频谱')
xlabel('时间(s)');
ylabel('幅值(dB)');
set(gca, 'YGrid', 'on');
%对imf3分量进行倒谱分析
figure(4)
z=rceps(imf(3,:));
plot(t(nn),abs(z(nn)));
title('imf3倒频谱')
xlabel('时间(s)');
ylabel('幅值(dB)');
set(gca, 'YGrid', 'on');
%对imf4分量进行倒谱分析
figure(5)
z=rceps(imf(4,:));
plot(t(nn),abs(z(nn)));
title('imf4倒频谱')
xlabel('时间(s)');
ylabel('幅值(dB)');
set(gca, 'YGrid', 'on');
%对imf5分量进行倒谱分析
figure(6)
z=rceps(imf(5,:));
plot(t(nn),abs(z(nn)));
title('imf5倒频谱')
xlabel('时间(s)');
ylabel('幅值(dB)');
set(gca, 'YGrid', 'on');






11.png
22.png
33.png
44.png
55.png
66.png
回复
分享到:

使用道具 举报

发表于 2013-5-16 21:45 | 显示全部楼层
做了FFT后可以做倒谱吗?和直接做倒谱有什么区别啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 17:49 , Processed in 0.102478 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表