www.gusucode.com > wavelet工具箱matlab源码程序 > wavelet/wavelet/+wavelet/+internal/admConstant.m
function cpsi = admConstant(wavname,morseparams) % Admissibility constants are not magnitude-squared here because % we are using a delta distribution as the reconstruction wavelet % Morlet's single integral formula. % Lilly,J.M. & Ohlede, S. (2009) Higher-order properties of analytic % wavelets, IEEE Transactions on Signal Processing, 57(1),146-160. if strcmpi(wavname,'bump') mu = 5; sigma = 0.6; bumpwav = @(w)abs(2*exp(1)*exp(-1./(1-((w-mu)/sigma).^2)))./w; cpsi = integral(bumpwav,mu-sigma,mu+sigma); elseif strcmpi(wavname,'morse') ga = morseparams(1); be = morseparams(2); anorm = wavelet.internal.morsenormconstant(ga,be); cpsi = anorm/ga*gamma(be/ga); % LilLy and Ohlede, 2009 % Equivalent to the following but has better convergence % morsewav = @(w)abs(anorm*w.^be.*exp(-w.^ga))./w; % cpsi = integral(morsewav,0,Inf); elseif strcmpi(wavname,'amor') cf = 6; morlwav = @(w)abs(2*exp(-(w-cf).^2/2))./w; cpsi = integral(morlwav,0,Inf); end