%外圈故障特征頻率?理論值?107.36Hz????
clc
clear?all
load?OR007@12_0_drive;???%outer-race?fault（drive）
data1=X156_DE_time;??
load?normal_0_1797;?????%Normal?baseline?Data（drive）
data2=X097_DE_time;
N=4096;????????????????%采樣點數
fs=12000;?
t=1:N;
t=t/fs;
xdata1=data1（1:4096）;xdata_outer=（xdata1-mean（xdata1）/std（xdata11））;
xdata2=data2（1:4096）;xdata_normal=（xdata2-mean（xdata2）/std（xdata21））;%零均值處理
?%%?====================故障信號小波包=====================%%
wpt1=wpdec（xdata_outer5‘db1‘‘shannon‘）;
%%?=================故障信號峭度值========
for?i=1:32
????c1=wpcoef（wpt1[5i-1]）;
????K1（i）=kurtosis（c1）;
end
figure;
subplot（311）;
bar（K1）;
title（‘外圈故障信號小波包分解各節點峭度值‘）;
%%?==============正常信號小波包===============%%
wpt2=wpdec（xdata_normal5‘db1‘‘shannon‘）;
%%?===============正常信號峭度值?=============

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1540??2015-04-23?16:35??outer_normal_kurtosis.m
?????文件?????2935152??2015-01-02?15:28??OR007@12_0_drive.mat