请选择 进入手机版 | 继续访问电脑版

Hello Mat

 找回密码
 立即注册
查看: 3947|回复: 0

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

[复制链接]

1064

主题

1239

帖子

13

金钱

管理员

Rank: 9Rank: 9Rank: 9

积分
21393
发表于 2018-7-26 23:11:26 | 显示全部楼层 |阅读模式
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. % peaks = findpeaks(bjp);
  38. % plot(peaks);

  39. %% 移除路径
  40. rmpath(genpath('./data'))
  41. rmpath(genpath('./Hilbert'))
  42. rmpath(genpath('./libsvm-3.22'))
复制代码
GUI如下:


本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
算法QQ  3283892722
群智能算法链接http://halcom.cn/forum.php?mod=forumdisplay&fid=73
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Python|Opencv|MATLAB|Halcom.cn

GMT+8, 2022-7-3 14:38 , Processed in 0.184980 second(s), 22 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表