利用小波变换检测轴承滚动体

200810100659

轴承故障为滚动体局部剥落。轴承故障特征频率为外环通过频率及其倍频，理论值为219HZ

1. %采样频率
   
2. fs=10000;
   
3. 
   
4. %轴承滚动体故障信号
   
5. fid=fopen('bearingroll.dat','r');%故障
   
6. N=1024;
   
7. xdata=fread(fid,N,'int16');
   
8. fclose(fid);
   
9. xdata=(xdata-mean(xdata))/std(xdata,1);
   
10. 
    
11. %时域波形
    
12. figure(1);
    
13. plot(1:N,xdata);
    
14. xlabel('样本序号 n/个');
    
15. ylabel('电压 V/v');
    
16. 
    
17. %db10小波进行4层分解
    
18. %一维小波分解
    
19. [c,l] = wavedec(xdata,4,'db10');
    
20. 
    
21. %重构第1-4层细节系数
    
22. d4 = wrcoef('d',c,l,'db10',4);
    
23. d3 = wrcoef('d',c,l,'db10',3);
    
24. d2 = wrcoef('d',c,l,'db10',2);
    
25. d1 = wrcoef('d',c,l,'db10',1);
    
26. 
    
27. %显示细节信号
    
28. figure(2)
    
29. subplot(4,1,1);
    
30. plot(d4,'LineWidth',2);
    
31. ylabel('d4');
    
32. subplot(4,1,2);
    
33. plot(d3,'LineWidth',2);
    
34. ylabel('d3');
    
35. subplot(4,1,3);
    
36. plot(d2,'LineWidth',2);
    
37. ylabel('d2');
    
38. subplot(4,1,4);
    
39. plot(d1,'LineWidth',2);
    
40. ylabel('d1');
    
41. xlabel('样本序号 n/个');
    
42. 
    
43. %第1层细节信号的包络谱
    
44. y=hilbert(d1);
    
45. ydata=abs(y);
    
46. y=y-mean(y);
    
47. nfft=1024;
    
48. p=abs(fft(ydata,nfft));
    
49. figure(3);
    
50. plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
    
51. xlabel('频率 f/Hz');
    
52. ylabel('功率谱 P/W');
    

复制代码

200810100659

如果想将小波和双谱分析结合起来，该怎么做呢

hanyou

'bearingroll.dat'
请楼主把这个文件传一传吧。