新建 文本文档 (2).txt

采用墨西哥帽子小波函数

>> load vonkoch
vonkoch=vonkoch(1:400);
S_Min=1;S_Max=32;
index=0;
for scale=S_Max:-1:S_Min;
index=index+1;
cwt_coef(index,:)=Singularity_Detection(scale,32*(scale),vonkoch);
end
cwtcoef_abs=abs(cwt_coef);
for index=S_Min:S_Max
       max_coef=max(cwtcoef_abs(index,:));
       min_coef=min(cwtcoef_abs(index,:));
       ext=max_coef-min_coef;
       cwtcoefabs(index,:)=64*(cwtcoef_abs(index,:)-min_coef)/ext;
      end
figure(1)
subplot(2,1,1);
plot(vonkoch);
xlabel('shijian')
ylabel('fudu')
title('fenxingxinhao')
axis([1 400 0 0.002])
subplot(2,1,2)
colormap(pink(64));
image(cwtcoef_abs)
set(gca,'YTick',2:3:32)
set(gca,'YTickLabel',32:-3:2)
title('lianxuxiaobobianhuanshijantu')
xlabel('shijian')
vonkoch=vonkoch(1:400);
S_Min=1;S_Max=32;
index=0;
for scale=S_Max:-1:S_Min;
index=index+1;
cwt_coef(index,:)=Singularity_Detection(scale,32*(scale),vonkoch);
end
cwtcoef_abs=abs(cwt_coef);
for index=S_Min:S_Max
       max_coef=max(cwtcoef_abs(index,:));
       min_coef=min(cwtcoef_abs(index,:));
       ext=max_coef-min_coef;
       cwtcoefabs(index,:)=64*(cwtcoef_abs(index,:)-min_coef)/ext;
      end
figure(1)
subplot(2,1,1);
plot(vonkoch);
xlabel('shijian')
ylabel('fudu')
title('fenxingxinhao')
subplot(2,1,2)
colormap(pink(64));
image(cwtcoef_abs)
set(gca,'YTick',2:3:32)
set(gca,'YTickLabel',32:-3:2)
title('lianxuxiaobobianhuanshijantu')
xlabel('shijian')
ylabel('chidu')
>> vonkoch=vonkoch(1:400);
S_Min=1;S_Max=32;
index=0;
for scale=S_Max:-1:S_Min;
index=index+1;
cwt_coef(index,:)=Singularity_Detection(scale,32*(scale),vonkoch);
end
cwtcoef_abs=abs(cwt_coef);
for index=S_Min:S_Max
       max_coef=max(cwtcoef_abs(index,:));
       min_coef=min(cwtcoef_abs(index,:));
       ext=max_coef-min_coef;
       cwtcoefabs(index,:)=64*(cwtcoef_abs(index,:)-min_coef)/ext;
      end
figure(1)
subplot(2,1,1);
plot(vonkoch);
xlabel('shijian')
ylabel('fudu')
title('fenxingxinhao')
axis([1 400 0 0.002])
subplot(2,1,2)
colormap(pink(64));
image(cwtcoef_abs)
set(gca,'YTick',2:3:32)
set(gca,'YTickLabel',32:-3:2)
title('lianxuxiaobobianhuanshijantu')
xlabel('shijian')
vonkoch=vonkoch(1:400);
S_Min=1;S_Max=32;
index=0;
for scale=S_Max:-1:S_Min;
index=index+1;
cwt_coef(index,:)=Singularity_Detection(scale,32*(scale),vonkoch);
end
cwtcoef_abs=abs(cwt_coef);
for index=S_Min:S_Max
       max_coef=max(cwtcoef_abs(index,:));
       min_coef=min(cwtcoef_abs(index,:));
       ext=max_coef-min_coef;
       cwtcoefabs(index,:)=64*(cwtcoef_abs(index,:)-min_coef)/ext;
      end
figure(1)
subplot(2,1,1);
plot(vonkoch);
xlabel('shijian')
ylabel('fudu')
title('fenxingxinhao')

subplot(2,1,2)
colormap(pink(64));
image(cwtcoef_abs)
set(gca,'YTick',2:3:32)
set(gca,'YTickLabel',32:-3:2)
title('lianxuxiaobobianhuanshijantu')
xlabel('shijian')
ylabel('chidu')