Matlab Marr小波 墨西哥草帽小波变换尺度图程序

clc;clear load vonkoch vonkoch=vonkoch(1:510); 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; cwtcoef_abs(index,:)=64*(cwtcoef_abs(index,:)-min_coef)/ext; end figure(1) subplot(2,1,1); plot(vonkoch); xlabel('Time'); ylabel('Amplitude'); title('Sigal'); axis([1 510 0 0.02]) subplot(2,1,2) colormap(pink(64)); image(cwtcoef_abs) set(gca,'YTick',2:3:32) set(gca,'YTickLabel',32:-3:2) title('连续小波变换时间尺度图') xlabel('时间') ylabel('尺度')