%%%%以下提取小波包重構信號的故障診斷
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