Inverse discrete wavelet transform
"Reconstruct" reconstructs a signal from its wavelet decomposition
(discrete wavelet transform).
The only parameter is the
If we have the approximation and detail
signal at octave i (see the manual section of Decompose) we define new signals
- number n of octaves of the discrete wavelet transform.
Thus the new signals are upsampled versions of the original signals.
The approximation at octave i-1 is
where mathtextmt2conv denotes a convolution operation.
Starting at the last octave n, this iteration is repeated until we finally
get the reconstructed signal .
First, decompose a signal by applying Decompose , confirm the number of
octaves suggested; then pass the result to "Reconstruct" , leaving the number of
The result from reconstruction should be equal to the original signal.
(Don't mind the padding that will occur if the signal length was not a
power of 2.)
y = WaveletReconst(x,n);
Note that x must be a discrete wavelet transform (DWT) signal.
Define wavelet, Load wavelet,
Save wavelet, Decompose,
Daubechies , Rioul/Vetterli