%本程序經過對信號的cwt變換，從變換系數中提取高頻和低頻部分的信號；高頻系數可以確定缺陷的位置，
%低頻系數可以顯示底拼波形
load c.txt%裝載待分解的信號
s=b;
n1=1;
n2=11;
a=(n1:n2);%a為尺度因子，與頻率大致成反比，所以a低對應高頻，a高對應低頻
%wavelets='cmor1-0.5';
w=cwt(s,a,'cmor1-0.5');%1為頻帶寬；0.5為中心頻率。兩個都可以調節
w=w';
w=angle(w);
%plot(w)；%顯示系數整體
w=w';
for i=1:(n2-n1);%通過循環顯示相應的頻段
   si=w(i,:);
  subplot((n2-n1)/2,2,i);
  i=i+1;
  plot(si);
 end