c09caf 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.
c09caf computes the one-dimensional DWT of a given input data array,
, for
,
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
or
is returned by the initialization routine
c09aaf.
If on entry
or
, 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.
None.