rwtscale.sci

小波分解源代码

function scale = RWTscale(sig,cwt)
// RWTscale is equal to CWTscale
// CWTscale -- Complement to Continuous Wavelet Transform
//  Usage
//    scale = CWTscale(sig,cwt);
//  Inputs
//    sig      1-d signal
//    cwt      continous wavelet transform of sig, Output from CWT(sig)
//  Outputs
//    scale    complement of continuous wavelet Transform of sig
//  Side Effects   
//  Description
//     scale = fft(sig) - sum(columns of cwt)
//  Algorithm
//     
//  Examples
//     N = 1024;     // signal length; 
//     nvoice = 12;  
//     // create a Brownian
//     CantorMeasure = MakeFractal(N,3,'Deterministic',[.5 0 .5]);
//     Devil  = cumsum(CantorMeasure); t = (.5:(N-.5))./N;
//     figure; plot(t, Devil); title(sprintf(' Devil Staircase Signal'));
//     // make CWT
//     Devil_cwt = CWT(Devil,nvoice,'Sombrero');
//     // complement to CWT
//     scale = CWTscale(Devil,Devil_cwt);
//     plot(real(scale));
//  See Also
//    CWT   iCWT
//  References
//    Mallat, "Wavelet Signal Processing"; 
//    
//  Copyright Aldo I Maalouf

 sz = size(cwt);
 
 fftsig = mtlb_fft(sig(:));
 for i = 1:sz(2),
	fftsig = fftsig - mtlb_fft(cwt(:,i));
 end;
 scale = (mtlb_ifft(fftsig));
    
endfunction