morlet.m

基于morlet小波的最小二乘匹配追踪算法

function [psi,x] = morlet(lb,ub,nb,center_time,delt,amp,f) 
x = linspace(lb,ub,nb);        % wavelet support.
%psi = amp*exp(-((x-center_time).^2)/delt) .* cos(2*pi*(x-center_time));
%%%correct

psi = amp*exp(-((x-center_time).^2)*(f^2)*log(2)/delt) .* cos(2*pi*f*(x-center_time));
%psi = amp*exp(-((x-center_time).^2)*f/delt);