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（scale32*（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（211）;
plot（vonkoch）;
xlabel（‘Time‘）;
ylabel（‘Amplitude‘）;
title（‘Sigal‘）;
axis（[1?510?0?0.02]）

subplot（212）
colormap（pink（64））;
image（cwtcoef_abs）
set（gca‘YTick‘2:3:32）
set（gca‘YTickLabel‘32:-3:2）
title（‘連續(xù)小波變換時間尺度圖‘）
xlabel（‘時間‘）
ylabel（‘尺度‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????747??2015-10-10?12:41??Marr小波\CWT.m
?????文件??????129147??2015-10-10?12:43??Marr小波\MarrCWT.fig
?????文件???????42059??2015-10-10?12:46??Marr小波\MarrCWT.jpg
?????文件?????????437??2015-10-10?12:28??Marr小波\Singularity_Detection.m