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