基于小波分析的轴承故障频率检测matlab程序

daihong

1. %外圈故障特征频率：107.67Hz
   
2. clear
   
3. clc
   
4. load 156.mat
   
5. sig=X156_DE_time;fs=12000;
   
6. N=12000;Ts=1/fs;
   
7. sig=sig(1:N);
   
8. t=0:Ts:(N-1)*Ts;
   
9. %归一化
   
10. %sig=(sig-mean(sig))/std(sig,1);
    
11. plot(t,sig);title('外圈故障轴承的归一化时域波形');grid on;
    
12. xlabel('时间 t/s');
    
13. ylabel('幅值 A/V');
    
14. %%%%%%%%%%%%%
    
15. %功率谱
    
16. nfft=16384;
    
17. y=abs(psd(sig,nfft));
    
18. figure(2)
    
19. plot((0:nfft/2-1)/nfft*fs,y(1:nfft/2));
    
20. xlabel('频率 f/Hz');
    
21. ylabel('归一化功率谱 S/V');
    
22. %db10小波进行4层分解
    
23. %一维小波分解
    
24. [c,l] = wavedec(sig,4,'db10');
    
25. %重构第1-4层细节信号
    
26. d4 = wrcoef('d',c,l,'db10',4);
    
27. d3 = wrcoef('d',c,l,'db10',3);
    
28. d2 = wrcoef('d',c,l,'db10',2);
    
29. d1 = wrcoef('d',c,l,'db10',1);
    
30. %显示细节信号
    
31. figure;
    
32. subplot(4,1,1);
    
33. plot(d4,'LineWidth',2);
    
34. ylabel('d4');
    
35. subplot(4,1,2);
    
36. plot(d3,'LineWidth',2);
    
37. ylabel('d3');
    
38. subplot(4,1,3);
    
39. plot(d2,'LineWidth',2);
    
40. ylabel('d2');
    
41. subplot(4,1,4);
    
42. plot(t,d1,'LineWidth',2);
    
43. ylabel('d1');
    
44. xlabel('样本序号 n/个');
    
45. %第1层细节信号的包络谱
    
46. y=hilbert(d1);
    
47. ydata=abs(y);
    
48. ydata=ydata-mean(ydata);
    
49. nfft=16384;
    
50. p=abs(fft(ydata,nfft));
    
51. figure;
    
52. plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
    
53. xlabel('频率 f/Hz');
    
54. ylabel('功率谱 P/W');
    
55. %%%%%%故障检测%%%%%%
    
56. P=p(1:nfft/2);
    
57. [M,f1]=max(P)
    
58. disp('故障特征频率为:');
    
59. f1=f1*fs/nfft-1
    
60. 
    
61. 
    

复制代码