小波变换matlab源代码,连续小波变换和逆小波变换

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')