singularity_detection.m

采用墨西哥帽子小波函数

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);