continuous wavelet transform and inverse
UPDATED FEB 25 2011
This collection of files is an extension of the wavelet software package originally written by Torrence and Compo (http://paos.colorado.edu/research/wavelets/software.html)
The main functions are:
1. contwt.m: (continuous wavelet transform). This is essentially Torrence and Compo's wavelet.m with a few modifications (more inputs and outputs for easier access)
2. invcwt.m: inverse continuous wavelet transform.
Please see the help in each function for details and usage.
Example usage:
%make a test signal
dt = 0.1
t = 0:dt:10;
x = cos(2*t)
%compute the CWT
[wave, period, scale, coi, dj,paramout, k] = contwt(x,dt);
%reconstruct the original signal
Xrec = invcwt(wave, 'MORLET', scale, paramout, k);
%Plot results
figure;
plot(t,x); %original signal
hold on;
plot(t, Xrec, 'r--')
legend('Original signal', 'Reconstructed Signal')
xlabel('Time');
ylabel('Signal (arbitrary units)')
%plot wavelet coeffs
figure;imagesc(abs(wave))
xlabel('Time (integer index)')
ylabel('Scale')