资源简介

小波包能量谱在滚动轴承故障诊断中的应用,文件包括轴承内圈、外圈、滚珠故障振动数据原始和正常工作振动数据等,编写小波包能量谱的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


评论

共有 条评论