dyadicwaveletdecomp.m

使用程序detsigularity.m调用多孔算法dyadicwaveletdecomp

function [a,d]=dyadicWaveletDecomp(I,wname,leveln)
[Lo_D,Hi_D] = wfilters(wname,'d');
lengthOfSignal=length(I);
sigLow=I;
for m=1:leveln
    filterLow=upsample(Lo_D,m);
    filterHigh=upsample(Hi_D,m);
    temp=wkeep(conv(sigLow,filterLow),lengthOfSignal,'c');
    sigHigh((m-1)*lengthOfSignal+1:m*lengthOfSignal)=wkeep(conv(sigLow,filterHigh),lengthOfSignal,'c');;
    sigLow=temp;
    if m==leveln
        a=sigLow;
    end
end
    d=sigHigh;