Code covered by the BSD License  

Highlights from
AR filter + Minimum Entropy Deconvolution for Bearing Fault Diagnosis

4.0

4.0 | 1 rating Rate this file 40 Downloads (last 30 days) File Size: 652 KB File ID: #41614
image thumbnail

AR filter + Minimum Entropy Deconvolution for Bearing Fault Diagnosis

by

 

AR Filter by YuleWalker method combined with Minimum Entropy Deconvolution for bearing fault diagnos

| Watch this File

File Information
Description

This function AR_MED_FILTER takes input SIGNAL with Sampling Frequency, Fs, and applies the Yule Walker method based AR filter. The order if the filter is found by Maximum kurtosis. After the application od AR filter, the signal is passed through Minimum Entropy Deconvolution. This combined AR+MED method brings out the Bearing faults hidden in Noise.

The function plots two figures for AR alone and another for AR+MED
Example:
   load('s4.mat');
   signal=s4;
   Fs=12000;
   ar_med_filter(signal,Fs);
The File 's4.mat' is the vibration signal recorded from a OR faulty bearing with a sampling frequency of 12000Hz. The Fault frequency is 161 Hz and is brought out.

This program isa based on the paper:
 Sawalhi N, Randall RB and Endo H (2007) The enhancement of fault detection and diagnosis in rolling element bearings using minimum entropy deconvolution combined with spectral kurtosis. Mechanical Systems and Signal Processing. 21:2616-2633

This function is basically written for Bearing fault diagnosis from Vibration signal.

Dont forget to rate or comment on the matlab central site
http://www.mathworks.in/matlabcentral/fileexchange/authors/258518

Author:Santhana Raj.A
https://sites.google.com/site/santhanarajarunachalam/

Acknowledgements

Minimum Entropy Deconvolution (Med 1 D And 2 D) inspired this file.

Required Products Filter Design Toolbox
Signal Processing Toolbox
MATLAB
MATLAB release MATLAB 7.14 (R2012a)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
26 Jul 2013 BK

Thank you for this interesting code. On line 53 ([~,index]=sort(kurt,1,'descend')), the variable kurt is row vector. Thus, the sorting processing should be applied on dimension 2 to have an effect.

14 May 2013 tk

thank you for your wonderful works, which is what I need.

Contact us