Code covered by the BSD License  

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

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

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

AR filter + Minimum Entropy Deconvolution for Bearing Fault Diagnosis


Santhana Raj (view profile)


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

| Watch this File

File Information

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
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

Author:Santhana Raj.A


Med2d(X,Filter Size,Term Iter,Term Delta,Overlap Mode,Plot Mode) inspired this file.

Required Products Filter Design Toolbox
Signal Processing Toolbox
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

BK (view profile)

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.

Comment only
14 May 2013 tk

tk (view profile)

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

Contact us