轴承故障诊断小波分析.7z

clear all
clc

load 10k_32768_5.dat
s=X10k_32768_5‘;

fs=10000;                                                 %采样频率
N=32768;                                                  %采样点数


[cl]=wavedec（s11‘haar‘）;  %11层分解
ca1=appcoef（cl‘haar‘1）;   %V10
ca2=appcoef（cl‘haar‘2）;   %V9
ca3=appcoef（cl‘haar‘3）;   %V8
ca4=appcoef（cl‘haar‘4）;   %V7
ca5=appcoef（cl‘haar‘5）;   %V6
ca6=appcoef（cl‘haar‘6）;   %V5
ca7=appcoef（cl‘haar‘7）;   %V4
ca8=appcoef（cl‘haar‘8）;   %V3
ca9=appcoef（cl‘haar‘9）;   %V2
ca10=appcoef（cl‘haar‘10）;   %V1
ca11=appcoef（cl‘haar‘11）;   %V0

cd1=detcoef（cl1）;   %W10
cd2=detcoef（cl2）;   %W9
cd3=detcoef（cl3）;   %W8
cd4=detcoef（cl4）;   %W7
cd5=detcoef（cl5）;   %W6
cd6=detcoef（cl6）;   %W5
cd7=detcoef（cl7）;   %W4
cd8=detcoef（cl8）;   %W3
cd9=detcoef（cl9）;   %W2
cd10=detcoef（cl10）;   %W1
cd11=detcoef（cl11）;   %W0

%消噪，将高频成分中小于阈值的数置零，规则自行改变
k=100;                             %设置一个阈值自行改变这个阈值，以达到最佳效果

d11=cd11;

d10=cd10;

d9=cd9;

d8=cd8;

d7=cd7;

d6=cd6;
d6（find（d6
d5=cd5;
d5（find（d5
d4=cd4;
d4（find（d4
d3=cd3;
d3（find（d3
d2=cd2;
d2（find（d2
d1=cd1;
d1（find（d1
d1=zeros（1length（d1））;
c=[ca11 d11 d10 d9 d8 d7 d6 d5 d4 d3 d2 d1];
s_rec=waverec（cl‘haar‘）;
f_rec=abs（fft（hilbert（s_rec）））*2/N;   %希尔伯特变换

figure（1）
plot（s）;xlim（[0 N]）;title（‘外圈故障原始信号‘）
axis（[01500000.5]）;

figure（2）;
subplot（511）;plot（ca1）;ylabel（‘V10‘）;title（‘第一层近似系数‘）
subplot（512）;plot（ca2）;ylabel（‘V9‘）;title（‘第二层近似系数‘）
subplot（513）;plot（ca3）;ylabel（‘V8‘）;title（‘第三层近似系数‘）
subplot（514）;plot（ca4）;ylabel（‘V7‘）;title（‘第四层近似系数‘）
subplot（515）;plot（ca5）;ylabel（‘V6‘）;title（‘第五层近似系数‘）

figure（3）;
subplot（611）;plot（ca6）;ylabel（‘V5‘）;title（‘第六层近似系数‘）
subplot（612）;plot（ca7）;ylabel（‘V4‘）;title（‘第七层近似系数‘）
subplot（613）;plot（ca8）;ylabel（‘V3‘）;title（‘第八层近似系数‘）
subplot（614）;plot（ca9）;ylabel（‘V2‘）;title（‘第九层近似系数‘）
subplot（615）;plot（ca10）;ylabel（‘V1‘）;title（‘第十层近似系数‘）
subplot（616）;plot（ca11）;ylabel（‘V0‘）;title（‘第十一层近似系数‘）

figure（4）;
subplot（611）;plot（cd1）;ylabel（‘W10‘）;title（‘第一层小波细节系数‘）%细节系数
subplot（612）;plot（cd2）;ylabel（‘W9‘）;title（‘第二层小波细节系数‘）
subplot（613）;plot（cd3）;ylabel（‘W8‘）;title（‘第三层小波细节系数‘）
subplot（614）;plot（cd4）;ylabel（‘W7‘）;title（‘第四层小波细节系数‘）
subplot（615）;plot（cd5）;ylabel（‘W6‘）;title（‘第五层小波细节系数‘）
subplot（616）;plot（cd6）;ylabel（‘W5‘）;title（‘第六层小波细节系数‘）
figure（5）;
subplot（511）;plot（cd7）;ylabel（‘W4‘）;title（‘第七层小波细节系数‘）
subplot（512）;plot（cd8）;ylabel（‘W3‘）;title（‘第八层小波细节系数‘）
subplot（513）;plot（cd9）;ylabel（‘W2‘）;title（‘第九层小波细节系数‘）
subplot（514）;plot（cd10）;ylabel（‘W1‘）;title（‘第十层小波细节系数‘）
subplot（515）;plot（cd11）;ylabel（‘W0‘）;title（‘第十一层小波细节系数‘）

figure（6）
subplot（211）;plot（s_rec）;xlim（[0 N]）;title（‘进行小波分解后重组信号时域图‘）
axis（[01500000.22]）;
subplot（212）;plot（0:2*fs/N:2*（N/2-1）*fs/Nf_rec（1:N/2））;titl