杭州锐达数字技术有限公司
查看: 847|回复: 5

[小波] 多小波变换用于复合故障特征的分类提取

[复制链接]
发表于 2016-1-18 15:46 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有帐号?我要加入

x
有没有人研究过多小波变换,进来讨论讨论!!!
回复
分享到:

使用道具 举报

发表于 2016-1-20 10:24 | 显示全部楼层
怎样讨论,哈哈
 楼主| 发表于 2016-1-20 13:20 | 显示全部楼层

你接触过多小波变换吗?我多小波分解出两个通道的信号 然后用包络解调分析 得出的频率特征和之前计算出来的对应不上 你了解吗?
发表于 2017-3-7 10:14 | 显示全部楼层
可以交流下复合故障诊断的相关知识吗
发表于 2017-3-8 08:39 | 显示全部楼层
发表于 2017-3-10 09:35 | 显示全部楼层

  1. %t=0.001:0.001:1;
  2. t=1:1000;
  3. s1=sin(2*pi*50*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t));

  4. for t=1:500;
  5. s2(t)=sin(2*pi*50*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t));
  6. end

  7. for t=501:1000;
  8. s2(t)=sin(2*pi*200*t*0.001)+sin(2*pi*120*t*0.001)+rand(1,length(t));
  9. end


  10. subplot(9,2,1)
  11. plot(s1)
  12. title('原始信号')
  13. ylabel('S1')

  14. subplot(9,2,2)
  15. plot(s2)
  16. title('故障信号')
  17. ylabel('S2')

  18. wpt=wpdec(s1,3,'db1','shannon');

  19. %plot(wpt);
  20. s130=wprcoef(wpt,[3,0]);
  21. s131=wprcoef(wpt,[3,1]);
  22. s132=wprcoef(wpt,[3,2]);
  23. s133=wprcoef(wpt,[3,3]);
  24. s134=wprcoef(wpt,[3,4]);
  25. s135=wprcoef(wpt,[3,5]);
  26. s136=wprcoef(wpt,[3,6]);
  27. s137=wprcoef(wpt,[3,7]);
  28. s10=norm(s130);
  29. s11=norm(s131);
  30. s12=norm(s132);
  31. s13=norm(s133);
  32. s14=norm(s134);
  33. s15=norm(s135);
  34. s16=norm(s136);
  35. s17=norm(s137);

  36. st10=std(s130);
  37. st11=std(s131);
  38. st12=std(s132);
  39. st13=std(s133);
  40. st14=std(s134);
  41. st15=std(s135);
  42. st16=std(s136);
  43. st17=std(s137);

  44. disp('正常信号的特征向量');
  45. snorm1=[s10,s11,s12,s13,s14,s15,s16,s17]
  46. std1=[st10,st11,st12,st13,st14,st15,st16,st17]

  47. subplot(9,2,3);plot(s130);
  48. ylabel('S130');
  49. subplot(9,2,5);plot(s131);
  50. ylabel('S131');
  51. subplot(9,2,7);plot(s132);
  52. ylabel('S132');
  53. subplot(9,2,9);plot(s133);
  54. ylabel('S133');
  55. subplot(9,2,11);plot(s134);
  56. ylabel('S134');
  57. subplot(9,2,13);plot(s135);
  58. ylabel('S135');
  59. subplot(9,2,15);plot(s136);
  60. ylabel('S136');
  61. subplot(9,2,17);plot(s137);
  62. ylabel('S137');


  63. wpt=wpdec(s2,3,'db1','shannon');

  64. %plot(wpt);
  65. s230=wprcoef(wpt,[3,0]);
  66. s231=wprcoef(wpt,[3,1]);
  67. s232=wprcoef(wpt,[3,2]);
  68. s233=wprcoef(wpt,[3,3]);
  69. s234=wprcoef(wpt,[3,4]);
  70. s235=wprcoef(wpt,[3,5]);
  71. s236=wprcoef(wpt,[3,6]);
  72. s237=wprcoef(wpt,[3,7]);

  73. s20=norm(s230);
  74. s21=norm(s231);
  75. s22=norm(s232);
  76. s23=norm(s233);
  77. s24=norm(s234);
  78. s25=norm(s235);
  79. s26=norm(s236);
  80. s27=norm(s237);

  81. st20=std(s230);
  82. st21=std(s231);
  83. st22=std(s232);
  84. st23=std(s233);
  85. st24=std(s234);
  86. st25=std(s235);
  87. st26=std(s236);
  88. st27=std(s237);

  89. disp('故障信号的特征向量');
  90. snorm2=[s20,s21,s22,s23,s24,s25,s26,s27]
  91. std2=[st20,st21,st22,st23,st24,st25,st26,st27]


  92. subplot(9,2,4);plot(s230);
  93. ylabel('S230');
  94. subplot(9,2,6);plot(s231);
  95. ylabel('S231');
  96. subplot(9,2,8);plot(s232);
  97. ylabel('S232');
  98. subplot(9,2,10);plot(s233);
  99. ylabel('S233');
  100. subplot(9,2,12);plot(s234);
  101. ylabel('S234');
  102. subplot(9,2,14);plot(s235);
  103. ylabel('S235');
  104. subplot(9,2,16);plot(s236);
  105. ylabel('S236');
  106. subplot(9,2,18);plot(s237);
  107. ylabel('S237');

  108. %fft

  109. figure
  110. y1=fft(s1,1024);
  111. py1=y1.*conj(y1)/1024;
  112. y2=fft(s2,1024);
  113. py2=y2.*conj(y2)/1024;

  114. y130=fft(s130,1024);
  115. py130=y130.*conj(y130)/1024;
  116. y131=fft(s131,1024);
  117. py131=y131.*conj(y131)/1024;
  118. y132=fft(s132,1024);
  119. py132=y132.*conj(y132)/1024;
  120. y133=fft(s133,1024);
  121. py133=y133.*conj(y133)/1024;
  122. y134=fft(s134,1024);
  123. py134=y134.*conj(y134)/1024;
  124. y135=fft(s135,1024);
  125. py135=y135.*conj(y135)/1024;
  126. y136=fft(s136,1024);
  127. py136=y136.*conj(y136)/1024;
  128. y137=fft(s137,1024);
  129. py137=y137.*conj(y137)/1024;

  130. y230=fft(s230,1024);
  131. py230=y230.*conj(y230)/1024;
  132. y231=fft(s231,1024);
  133. py231=y231.*conj(y231)/1024;
  134. y232=fft(s232,1024);
  135. py232=y232.*conj(y232)/1024;
  136. y233=fft(s233,1024);
  137. py233=y233.*conj(y233)/1024;
  138. y234=fft(s234,1024);
  139. py234=y234.*conj(y234)/1024;
  140. y235=fft(s235,1024);
  141. py235=y235.*conj(y235)/1024;
  142. y236=fft(s236,1024);
  143. py236=y236.*conj(y236)/1024;
  144. y237=fft(s237,1024);
  145. py237=y237.*conj(y237)/1024;

  146. f=1000*(0:511)/1024;
  147. subplot(1,2,1);
  148. plot(f,py1(1:512));
  149. ylabel('P1');
  150. title('原始信号的功率谱')
  151. subplot(1,2,2);
  152. plot(f,py2(1:512));
  153. ylabel('P2');
  154. title('故障信号的功率谱')

  155. figure

  156. subplot(4,2,1);
  157. plot(f,py130(1:512));
  158. ylabel('P130');
  159. title('S130的功率谱')
  160. subplot(4,2,2);
  161. plot(f,py131(1:512));
  162. ylabel('P131');
  163. title('S131的功率谱')
  164. subplot(4,2,3);
  165. plot(f,py132(1:512));
  166. ylabel('P132');
  167. subplot(4,2,4);
  168. plot(f,py133(1:512));
  169. ylabel('P133');
  170. subplot(4,2,5);
  171. plot(f,py134(1:512));
  172. ylabel('P134');
  173. subplot(4,2,6);
  174. plot(f,py135(1:512));
  175. ylabel('P135');
  176. subplot(4,2,7);
  177. plot(f,py136(1:512));
  178. ylabel('P136');
  179. subplot(4,2,8);
  180. plot(f,py137(1:512));
  181. ylabel('P137');
  182. figure

  183. subplot(4,2,1);
  184. plot(f,py230(1:512));
  185. ylabel('P230');
  186. title('S230的功率谱')
  187. subplot(4,2,2);
  188. plot(f,py231(1:512));
  189. ylabel('P231');
  190. title('S231的功率谱')
  191. subplot(4,2,3);
  192. plot(f,py232(1:512));
  193. ylabel('P232');
  194. subplot(4,2,4);
  195. plot(f,py233(1:512));
  196. ylabel('P233');
  197. subplot(4,2,5);
  198. plot(f,py234(1:512));
  199. ylabel('P234');
  200. subplot(4,2,6);
  201. plot(f,py235(1:512));
  202. ylabel('P235');
  203. subplot(4,2,7);
  204. plot(f,py236(1:512));
  205. ylabel('P236');
  206. subplot(4,2,8);
  207. plot(f,py237(1:512));
  208. ylabel('P237');
  209. figure
  210. plottree(wpt)

复制代码


您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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