多尺度的连续小波变换分解后时间尺度图
墨西哥Marr小波函数:
%delta 小波变换尺度
%N 小波函数的长度
%s 原始信号
%g 原始信号某个尺度下的小波变化系数
function g=Singularity_Detection(delta,N,s);
%保持信号长度
n=length(s);
% 构造墨西哥小帽小波函数
for index_x=1:N;
x=index_x-(N+1)/2;
phi_x(index_x)=((pi^(-1/4))*(2/sqrt(3)))*(1-x.*x/(delta^2))*exp(-(x.*x)/(2*delta^2));
end;
phi_x=phi_x/norm(phi_x); %能量归一化
g=conv(s,phi_x); %卷积
g=wkeep(g,n); %保持信号长度
绘图:
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