基于小波包能量分析的轴承故障诊断

小波包能量谱在滚动轴承故障诊断中的应用，文件包括轴承内圈、外圈、滚珠故障振动数据原始和正常工作振动数据等，编写小波包能量谱的Matlab的源程序，经测试仿真可用，供从事这方面研究的学者参考。

%%%%以下提取小波包重构信号的故障诊断
clc
clear
fs=12e3;
load zhengchang.mat;%加载正常信号
load 98.mat;
s1=X098_DE_time（19001:1:200001）;% 正常情况下驱动端数据
% s1=X098_FE_time（1:1:100001）;% 正常情况下驱动端数据
load in.mat;%加载内圈故障信号
% load fan_in.mat;%加载风扇端内圈故障信号
load 106.mat;
s2=X106_DE_time（19001:1:200001）;% 内圈故障信号情况下驱动端数据
% s2=X279_FE_time（1:1:100001）;% 内圈故障信号情况下驱动端数据
% load fan_ball.mat;%加载滚动体故障信号
load ball.mat;%加载滚动体故障信号
load 119.mat;
s3=X119_DE_time（19001:1:200001）;% 滚动体故障情况下驱动端数据
% s3=X283_FE_time（1:1:100001）;% 滚动体故障情况下驱动端数据
% load fan_out_12.mat;%加载外圈12点位置故障信号
load out_12.mat;%加载外圈12点位置故障信号
load 158.mat;
s4=X158_DE_time（14001:1:150001）;% 外圈12点位置情况下驱动端数据
% s4=X305_FE_time（1:1:100001）;% 外圈12点位置情况下驱动端数据
n=3;
ji=‘db5‘;
wpt=wpdec（s1nji）;%使用db5小波包分解正常信号到底三层，使用shannon熵
% plot（wpt）;
rcfs10=wprcoef（wpt[n0]）;%利用相应的频段进行重构
rcfs11=wprcoef（wpt[n1]）;
rcfs12=wprcoef（wpt[n2]）;
rcfs13=wprcoef（wpt[n3]）;
rcfs14=wprcoef（wpt[n4]）;
rcfs15=wprcoef（wpt[n5]）;
rcfs16=wprcoef（wpt[n6]）;
rcfs17=wprcoef（wpt[n7]）;
cfs10=wpcoef（wpt[n0]）;%提取小波包分解系数
cfs11=wpcoef（wpt[n1]）;
cfs12=wpcoef（wpt[n2]）;
cfs13=wpcoef（wpt[n3]）;
cfs14=wpcoef（wpt[n4]）;
cfs15=wpcoef（wpt[n5]）;
cfs16=wpcoef（wpt[n6]）;
cfs17=wpcoef（wpt[n7]）;
wpt1=wpdec（s2nji）;
rcfs20=wprcoef（wpt1[n0]）;
rcfs21=wprcoef（wpt1[n1]）;
rcfs22=wprcoef（wpt1[n2]）;
rcfs23=wprcoef（wpt1[n3]）;
rcfs24=wprcoef（wpt1[n4]）;
rcfs25=wprcoef（wpt1[n5]）;
rcfs26=wprcoef（wpt1[n6]）;
rcfs27=wprcoef（wpt1[n7]）;
cfs20=wpcoef（wpt1[n0]）;
cfs21=wpcoef（wpt1[n1]）;
cfs22=wpcoef（wpt1[n2]）;
cfs23=wpcoef（wpt1[n3]）;
cfs24=wpcoef（wpt1[n4]）;
cfs25=wpcoef（wpt1[n5]）;
cfs26=wpcoef（wpt1[n6]）;
cfs27=wpcoef（wpt1[n7]）;
wpt2=wpdec（s3nji）;
rcfs30=wprcoef（wpt2[n0]）;
rcfs31=wprcoef（wpt2[n1]）;
rcfs32=wprcoef（wpt2[n2]）;
rcfs33=wprcoef（wpt2[n3]）;
rcfs34=wprcoef（wpt2[n4]）;
rcfs35=wprcoef（wpt2[n5]）;
rcfs36=wprcoef（wpt2[n6]）;
rcfs37=wprcoef（wpt2[n7]）;
cfs30=wpcoef（wpt2[n0]）;
cfs31=wpcoef（wpt2[n1]）;
cfs32=wpcoef（wpt2[n2]）;
cfs33=wpcoef（wpt2[n3]）;
cfs34=wpcoef（wpt2[n4]）;
cfs35=wpcoef（wpt2[n5]）;
cfs36=wpcoef（wpt2[n6]）;
cfs37=wpcoef（wpt2[n7]）;
wpt3=wpdec（s4nji）;
rcfs40=wprcoef（wpt3[n0]）;
rcfs41=wprcoef（wpt3[n1]）;
rcfs42=wprcoef（wpt3[n2]）;
rcfs43=wprcoef（wpt3[n3]）;
rcfs44=wprcoef（wpt3[n4]）;
rcfs45=wprcoef（wpt3[n5]）;
rcfs46=wprcoef（wpt3[n6]）;
rcfs47=wprcoef（wpt3[n7]）;
cfs40=wpcoef（wpt3[n0]）;
cfs41=wpcoef（wpt3[n1]）;
cfs42=wpcoef（wpt3[n2]）;
cfs43=wpcoef（wpt3[n3]）;
cfs44=wpcoef（wpt3[n4]）;
cfs45=wpcoef（wpt3[n5]）;
cfs46=wpcoef（wpt3[n6]）;
cfs47=wpcoef（wpt3[n7]）;
for i=1:2^n
E（i）=（norm（wpcoef（wpt[ni-1]）2））^2;%求第i个节点的范数平方，其实也就是平方和，即能量
end
E_total=sum（E）; %求总能量
for i=1:2^n
pfir（i）= E（i）/E_total;%若求每个节点的概率，即所占的能量比，则为E（i）/E_total
end

for i=1:2^n
E1（i）=（norm（wpcoef（wpt1[ni-1]）2））^2;%求第i个节点的范数平方，其实也就是平方和，即能量
end
E1_total=sum（E1）; %求总能量
for i=1:2^n
pfir1（i）= E1（i）/E1_total;%若求每个节点的概率，即所占的能量比

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件    2928192  2015-09-23 22:26  00004xiaobobaonengliangfenxi\106.mat

     文件    2914248  2015-09-23 22:29  00004xiaobobaonengliangfenxi\119.mat

     文件    2928192  2015-09-23 22:30  00004xiaobobaonengliangfenxi\158.mat

     文件    7742720  2015-09-23 22:03  00004xiaobobaonengliangfenxi\98.mat

     文件       5362  2016-05-26 16:30  00004xiaobobaonengliangfenxi\ball.mat

     文件       5354  2016-05-26 16:29  00004xiaobobaonengliangfenxi\in.mat

     文件       5376  2016-05-26 16:30  00004xiaobobaonengliangfenxi\out_12.mat

     文件       6699  2016-05-26 22:14  00004xiaobobaonengliangfenxi\xiaobobaonengliangfenxi.m

     文件       5376  2016-05-26 16:26  00004xiaobobaonengliangfenxi\zhengchang.mat

     目录          0  2016-05-26 22:15  00004xiaobobaonengliangfenxi

----------- ---------  ---------- -----  ----

             16541519                    10