.BG
.FN cwtp
.TL
Continuous Wavelet Transform with Phase Derivative
.DN
Computes the continuous wavelet transform with (complex-valued) 
Morlet wavelet and its phase derivative.
.CS
cwtp(input, noctave, nvoice=1, w0=2 * pi, twoD=TRUE, plot=TRUE)
.RA
.AG input
input signal (possibly complex-valued)
.AG noctave
number of powers of 2 for the scale variable
.OA
.AG nvoice
number of scales in each octave
(i.e. between two consecutive powers of 2).
.AG w0
central frequency of the wavelet.
.AG twoD
logical variable set to T to organize the output as a 2D array 
(signal\_size x nb\_scales), otherwise, the output is a 3D array 
(signal\_size x noctave x nvoice).
.AG plot
if set to T, display the modulus of the
continuous wavelet transform on the graphic device.
.RT
list containing the continuous (complex) wavelet transform
and the phase derivative

wt:array of complex numbers for the values of the continuous
wavelet transform

f:array of the same dimensions containing the values of the
derivative of the phase of the continuous wavelet transform
.SE
.DT
.SH REFERENCE
See discussions in the text of "Practical Time-Frequency Analysis".
.SA
"cgt", "cwt", "cwtTh", "DOG" for wavelet transform, and "gabor"
for continuous Gabor transform.
.EX
