【故障诊断分析】基于小波变换实现外圈轴承故障诊断含Matlab源码-CSDN博客

本文链接：https://blog.csdn.net/m0_57702748/article/details/140031384

1 简介

在滚动轴承的故障诊断时，传统的频谱分析法通常采用共振解调技术，具有良好的效果，但当内圈、滚动体或多点故障时，解调谱线却很难分辨故障类型小波包是小波理论在信号处理应用领域的又一重大发展，它在小波理论的基础上引人了最优基准准则，目前已在通信、图像处理、分形、信号检测等方面获得成功应用小波包分解是多分辨率分析的推广，是一种比小波分析更为精细的信号时频分析方法，它将频带进行多层次划分，对小波分析没有细分的高频部分也进行了分解，并能根据被分析信号的特征，自适应地选择频带，使之与信号频谱相匹配，从而提高时频分辨率，因此它是一种非常有效的信号处理方法。本文将用小波包分析方法构造滚动轴承状态检测信号的能量特征向量，再用神经网络方法对滚动轴承的状态进行识别。

1 小波包理论

【故障诊断分析】基于小波变换实现外圈轴承故障诊断含Matlab源码_小波分析

【故障诊断分析】基于小波变换实现外圈轴承故障诊断含Matlab源码_小波分析_02

【故障诊断分析】基于小波变换实现外圈轴承故障诊断含Matlab源码_信号处理_03

2 特征频率估算 滚动轴承各部件发生单点故障时，其特征信号频率为

【故障诊断分析】基于小波变换实现外圈轴承故障诊断含Matlab源码_小波分析_04

3 构造能量特征向量 当轴承在负载旋转过程中出现疲劳、压痕或局部腐蚀等故障时，将产生冲击振动，而这些振动反复互激又会引起更大的振颤，用加速度传感器获取运动中轴承的振动加速度信号，由于轴承的故障使得该振动信号在某些频带内的信号能量减小，而使另外一些频带内的信号能量增强，因此，在各频率成分信号的能量中，包含了丰富的故障信息，某种或某几种频率成分能量的改变代表一种故障，基于这点，由轴承振动信号的各频段能量相对值构造轴承故障特征向量。

【故障诊断分析】基于小波变换实现外圈轴承故障诊断含Matlab源码_小波分析_05

2 部分代码


          
          
           
           clear all
           
           
clc
           
           
load 10k_32768_5.dat
           
           
s=X10k_32768_5';
           
           
fs=10000;                                                 %采样频率
           
           
N=32768;                                                  %采样点数
           
           
[c,l]=wavedec(s,11,'haar');  %11层分解
           
           
ca1=appcoef(c,l,'haar',1);   %V10
           
           
ca2=appcoef(c,l,'haar',2);   %V9
           
           
ca3=appcoef(c,l,'haar',3);   %V8
           
           
ca4=appcoef(c,l,'haar',4);   %V7
           
           
ca5=appcoef(c,l,'haar',5);   %V6
           
           
ca6=appcoef(c,l,'haar',6);   %V5
           
           
ca7=appcoef(c,l,'haar',7);   %V4
           
           
ca8=appcoef(c,l,'haar',8);   %V3
           
           
ca9=appcoef(c,l,'haar',9);   %V2
           
           
ca10=appcoef(c,l,'haar',10);   %V1
           
           
ca11=appcoef(c,l,'haar',11);   %V0
           
           
cd1=detcoef(c,l,1);   %W10
           
           
cd2=detcoef(c,l,2);   %W9
           
           
cd3=detcoef(c,l,3);   %W8
           
           
cd4=detcoef(c,l,4);   %W7
           
           
cd5=detcoef(c,l,5);   %W6
           
           
cd6=detcoef(c,l,6);   %W5
           
           
cd7=detcoef(c,l,7);   %W4
           
           
cd8=detcoef(c,l,8);   %W3
           
           
cd9=detcoef(c,l,9);   %W2
           
           
cd10=detcoef(c,l,10);   %W1
           
           
cd11=detcoef(c,l,11);   %W0
           
           
%消噪，将高频成分中小于阈值的数置零，规则自行改变
           
           
k=100;                             %设置一个阈值,自行改变这个阈值，以达到最佳效果
           
           
d11=cd11;
           
           
d10=cd10;
           
           
d9=cd9;
           
           
d8=cd8;
           
           
d7=cd7;
           
           
d6=cd6;
           
           
d6(find(d6<max(d6)/k))=0;
           
           
d5=cd5;