beshe3.m

以上是我大学毕业设计

%本程序经过对信号的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