【故障诊断分析】小波变换外圈轴承故障诊断【含Matlab源码 1678期】

💥💥💥💥💥💥💥💥💞💞💞💞💞💞💞💞💞Matlab武动乾坤博客之家💞💞💞💞💞💞💞💞💞💥💥💥💥💥💥💥💥
🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚤🚀🚀🚀🚀🚀🚀🚀🚀🚀🚀

🔊博主简介：985研究生，Matlab领域科研开发者；

🚅座右铭：行百里者，半于九十。

🏆代码获取方式：
CSDN Matlab武动乾坤—代码获取方式

更多Matlab信号处理仿真内容点击👇
①Matlab信号处理（进阶版）

⛳️关注CSDN Matlab武动乾坤，更多资源等你来！！

⛄一、简介

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

1 小波包理论



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

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

⛄二、部分源代码

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

%消噪，将高频成分中小于阈值的数置零，规则自行改变

figure(2);
subplot(5,1,1);plot(ca1);ylabel(‘V10’);title(‘第一层近似系数’)
subplot(5,1,2);plot(ca2);ylabel(‘V9’);title(‘第二层近似系数’)
subplot(5,1,3);plot(ca3);ylabel(‘V8’);title(‘第三层近似系数’)
subplot(5,1,4);plot(ca4);ylabel(‘V7’);title(‘第四层近似系数’)
subplot(5,1,5);plot(ca5);ylabel(‘V6’);title(‘第五层近似系数’)

figure(3);
subplot(6,1,1);plot(ca6);ylabel(‘V5’);title(‘第六层近似系数’)
subplot(6,1,2);plot(ca7);ylabel(‘V4’);title(‘第七层近似系数’)
subplot(6,1,3);plot(ca8);ylabel(‘V3’);title(‘第八层近似系数’)
subplot(6,1,4);plot(ca9);ylabel(‘V2’);title(‘第九层近似系数’)
subplot(6,1,5);plot(ca10);ylabel(‘V1’);title(‘第十层近似系数’)
subplot(6,1,6);plot(ca11);ylabel(‘V0’);title(‘第十一层近似系数’)

figure(4);
subplot(6,1,1);plot(cd1);ylabel(‘W10’);title(‘第一层小波细节系数’)%细节系数
subplot(6,1,2);plot(cd2);ylabel(‘W9’);title(‘第二层小波细节系数’)
subplot(6,1,3);plot(cd3);ylabel(‘W8’);title(‘第三层小波细节系数’)
subplot(6,1,4);plot(cd4);ylabel(‘W7’);title(‘第四层小波细节系数’)
subplot(6,1,5);plot(cd5);ylabel(‘W6’);title(‘第五层小波细节系数’)
subplot(6,1,6);plot(cd6);ylabel(‘W5’);title(‘第六层小波细节系数’)
figure(5);
subplot(5,1,1);plot(cd7);ylabel(‘W4’);title(‘第七层小波细节系数’)
subplot(5,1,2);plot(cd8);ylabel(‘W3’);title(‘第八层小波细节系数’)
subplot(5,1,3);plot(cd9);ylabel(‘W2’);title(‘第九层小波细节系数’)
subplot(5,1,4);plot(cd10);ylabel(‘W1’);title(‘第十层小波细节系数’)
subplot(5,1,5);plot(cd11);ylabel(‘W0’);title(‘第十一层小波细节系数’)

figure(6)
subplot(2,1,1);plot(s_rec);xlim([0 N]);title(‘进行小波分解后重组信号时域图’)
axis([0,15000,0,0.22]);
subplot(2,1,2);plot(0:2fs/N:2(N/2-1)*fs/N,f_rec(1:N/2));title(‘外圈故障小波分解重组信号FFT频谱图’)
axis([0,1000,0,0.005]);

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.
[4]张国新,刘祚时.基于小波包分析的滚动轴承故障诊断[J].江西理工大学学报. 2007,(01)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合