【数字信号】基于matlab轴承故障分析（时频图+功率谱）【含Matlab源码 922期】

一、简介

基于matlab轴承故障分析（时频图+功率谱）

二、源代码

clc
clear all
x1=load('sig1.txt');
x2=load('sig2.txt');
x3=load('sig3.txt');
x4=load('sig4.txt');

x=x1;%sig1.txt
fs=15360;
N=length(x);
t=(0:N-1)/fs;
figure();
subplot(211);
plot(t,x);
title('sig1波形图');
xlabel('时间 ');ylabel('振幅 ');
axis([0,max(t),min(x)*1.1,max(x)*1.1]);
f=(0:N/2-1)*fs/N;
y=abs(fft(x-mean(x),N));
ffty=y(1:N/2)*2/N;
subplot(212);
plot(f,ffty);
title('sig1频谱图');
xlabel('频率');ylabel('幅值');
[xp,xrms,k,C,I,K4]=canshu(x,N);
shu1=[xp,xrms,k,C,I,K4];

%%平方解调&hilbert包络解调
z=hilbert(x);
z1=abs(z);
z2=z1-mean(z1);
z3=abs(fft(z2,N));
Hilbertz=z3(1:N/2)*2/N;
xx=x.^2;
function [xp,xrms,k,C,I,K4]=canshu(x,N)%峰值、有效值、峭度、峰值指标、脉冲指标、峭度指标
xp=max(abs(x));%峰值
 xrms=sqrt(1/N*sum(x.^2));%有效值（均方根值）
 k=sum(x.^4)/N;%峭度
 K4=kurtosis(x);%峭度指标
 C=xp./xrms;%峰值指标：峰值/有效值
 u=1/N*sum(x);%均值

三、运行结果

四、备注

版本：2014a