对凯斯西储大学的轴承故障信号进行小波变换特征提取

最新推荐文章于 2024-07-25 10:45:22 发布

麻吉霍霍

最新推荐文章于 2024-07-25 10:45:22 发布

阅读量2.1w

点赞数 17

分类专栏： matlab 文章标签：故障重构 matlab 小波变换

本文链接：https://blog.csdn.net/sayaitachi/article/details/24877261

版权

本文探讨了对凯斯西储大学轴承故障信号的分析，通过小波变换提取特征。正常信号、外圈故障（164Hz及倍频）、内圈故障（105.5Hz及倍频）的分析结果显示，滚珠体故障频率理论值141Hz存在疑问，同时35.16Hz的出现令人困惑，未发现预期的转频29.9Hz。作者疑惑是否代码有误或小波变换不足，期待专家指导。

摘要由CSDN通过智能技术生成

最近本菜鸟对凯斯西储大学的轴承信号做了时域参数和小波变换，希望能分辨出其轴承出现故障与否和故障类型—外圈故障，内圈故障，滚珠体故障。以下是用matlab2011b处理信号后的结果及遇到的问题。先贴代码

close all;
clear all;
clc;
%采样频率   
fs=12000;
%采样间隔
Ts=1/fs;
%采样点数
N=10000;n=1/N;
%轴承信号
 load 97.mat;%故障
 xdata=X097_DE_time(1:10000);
% xdata=BF(1:10000);
xdata=(xdata-mean(xdata))/std(xdata,1);
%信号时域波形
figure(1);
plot(1:N,xdata);
xlabel('时间 t/n');
ylabel('电压 V/v');
%使用db10小波进行尺度为5时的分解
[c,l]=wavedec(xdata,5,'db10');
%从小波分解结构[c,l]重构信号xdata
a0=waverec(c,l,'db10');

figure(2);
subplot(3,1,1);
plot(xdata);
title('原始信号')
subplot(3,1,2);
plot(a0);
title('重构信号')
subplot(3,1,3);
plot(xdata-a0);
title('误差信号')
err=max(abs(xdata-a0))

%重构第1~5层高频细节信号
d5=wrcoef('d',c,l,'db10',5);
d4=wrcoef('d',c,l,'db10',4);
d3=wrcoef('d',c,l,'db10',3);
d2=wrcoef(