【数字信号】基于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
标签:load,abs,时频,源码,matlab,K4,xp,txt,xrms 来源: https://blog.csdn.net/TIQCmatlab/article/details/117191872