Complement to Continuous Wavelet Transform
scale = CWTscale(sig,cwt);
1-d signal
continous wavelet transform of sig, Output from CWT(sig)
complement of continuous wavelet Transform of sig
scale = fft(sig) - sum(columns of cwt)
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)); | ![]() | ![]() |
Mallat, "Wavelet Signal Processing"; 4.?.?