xiaobofenjie.m

对轴承故障振动信号的matlab小波分析程序

s=load('内圈损伤.TXT');
N=length(s);n=N-1;
fs=15360;f=(0:n)*fs/N;
t=0:(1/fs):(n/fs);
s=(s-mean(s))/std(s,1);
figure(1);
plot(t,s);

[c,l]=wavedec(s,4,'db10');
d4=wrcoef('d',c,l,'db10',4);
d3=wrcoef('d',c,l,'db10',3);
d2=wrcoef('d',c,l,'db10',2);
d1=wrcoef('d',c,l,'db10',1);
figure(2)
subplot(4,1,1);
plot(d4,'LineWidth',2);
ylabel('d4');
subplot(4,1,2);
plot(d3,'LineWidth',2);
ylabel('d3');
subplot(4,1,3);
plot(d2,'LineWidth',2);
ylabel('d2');
subplot(4,1,4);
plot(d1,'LineWidth',2);
ylabel('d1');
xlabel('样本序号 n/个');

y1=hilbert(d3);
my1=abs(y1);
py1=fft(my1-mean(my1));
mpy1=abs(py1);
figure;
subplot(211);plot(t,my1);title('滤波后的包络信号');
subplot(212);plot(f(1:(N/2)),mpy1(1:(N/2)));title('滤波信号的包络谱');