I am in the process of implementing a sampled continuous wavelet transform.
I have coded two wavelets (Morlet and a Hanning based). I require the phase
information to be indicated on the output as Vetterli and Herley showed in
the article IEEE Trans on Sig.Proc Vol 40 No 9. They showed how the magnitude
and phase info for a burst sine can be shown on the one diagram.

I would like a formula for a "real" wavelet (something as a function of t).
Vetterli's derivation is a discrete one so it does not fit in well with the
code I have put together.

B.T.W How have other people who require sub octave resolution implemented
a wavelet transform. I am using a fast convolution of filters approach which
is reasonably quick but are there better ways. (The Fast DWT works on an
octave basis. Is there a way to do the DWT on a sub octave approach).