基于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