computes the one-dimensional discrete wavelet transform (DWT) at a single level. The initialization routine c09aaf
must be called first to set up the DWT options.
computes the one-dimensional DWT of a given input data array,
at a single level. For a chosen wavelet filter pair, the output coefficients are obtained by applying convolution and downsampling by two to the input,
. The approximation (or smooth) coefficients,
, are produced by the low pass filter and the detail coefficients,
, by the high pass filter. To reduce distortion effects at the ends of the data array, several end extension methods are commonly used. Those provided are: periodic or circular convolution end extension, half-point symmetric end extension, whole-point symmetric end extension or zero end extension. The number
, of coefficients
is returned by the initialization routine c09aaf
If on entry
, explanatory error messages are output on the current error message unit (as defined by x04aaf
The accuracy of the wavelet transform depends only on the floating-point operations used in the convolution and downsampling and should thus be close to machine precision.