4-基于小波变换的故障诊断（MATLAB视频）

Halcom

4-基于小波变换的故障诊断（MATLAB视频）百度网盘链接：
链接：https://pan.baidu.com/s/1COCeo6kL5fDORHIyW_csoA
具体链接在halcom.cn论坛，联系人QQ：3283892722
该论坛是一个学习交流平台，我会逐一的和大家分享学习。
欢迎大家录制视频，并提交给我，我来设置视频，你可在论坛进行打赏分享。
视频专用播放器：http://halcom.cn/forum.php?mod=viewthread&tid=258&extra=page%3D1

主程序如下：

1. addpath(genpath('./data'))
   
2. addpath(genpath('./Hilbert'))
   
3. addpath(genpath('./libsvm-3.22'))
   
4. %% 随机分析一个信号
   
5. load('data.mat')  % 加载数据
   
6. x = RC304(:, 2);  % 信号
   
7. Ts = 7.8125e-5;   % 采样时间
   
8. Fs = 1/Ts;        % 采样频率
   
9. %% 小波滤波
   
10. [thr,sorh,keepapp,crit] = ddencmp('cmp','wp',x);
    
11. [xc,wpt,perf0,perfl2] = wpdencmp(x,sorh,3,'db2',crit,thr,keepapp);
    
12. [c,l] = wavedec(xc,3,'db1');   % 多尺度小波分解
    
13. %% FFT分析
    
14. L = length(c);
    
15. NFFT = 2^nextpow2(L); % Next power of 2 from length of y
    
16. c_fft = fft(c, NFFT)/L;
    
17. figure(1)
    
18. plot( Fs/2*linspace(0, 1, NFFT/2+1), abs( c_fft(1:NFFT/2+1) ) )
    
19. %% Hilbert谱
    
20. t0 = 0;    % 起始时刻
    
21. t1 = 2.5;  % 结束时间
    
22. fres = 100;  % 频率分辨率
    
23. tres = 100;  % 时间分辨率
    
24. fw0 = 0;     % 样本最小频率
    
25. fw1 = 200;   % 样本最大频率
    
26. tw0 = t0;    % 最小时间缩放
    
27. tw1 = t1;    % 最大时间缩放
    
28. [ns, ts, fs] = nnspe( c, t0, t1, fres, tres, fw0,fw1,tw0,tw1,'hilbert','pchip',3 );
    
29. filter_1D = fspecial('gaussian',9, 0.9);
    
30. nss = filter2(filter_1D, ns);
    
31. % 边际谱
    
32. bjp = sum( nss, 2 )./ Fs;
    
33. figure(2)
    
34. plot(fs, bjp);
    
35. xlabel('频率'); ylabel('幅值');
    
36. title('Hilbert边际谱')
    
37. 
    
38. % peaks = findpeaks(bjp);
    
39. % plot(peaks);
    
40. 
    
41. %% 移除路径
    
42. rmpath(genpath('./data'))
    
43. rmpath(genpath('./Hilbert'))
    
44. rmpath(genpath('./libsvm-3.22'))

复制代码

GUI如下：