Matlab_Marr小波尺度变换图

多尺度的连续小波变换分解后时间尺度图

墨西哥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