Reconstructs signal from a ``regularly sampled'' ridge, in the wavelet case, from a precomputed kernel.
regrec2(siginput, cwtinput, phi, nbnodes, noct, nvoice, Q2,
epsilon=0.5, w0=2 * pi, plot=FALSE)input signal.
wavelet transform, output of cwt.
unsampled ridge.
number of samples on the ridge
number of octaves (powers of 2)
number of different scales per octave
second term of the reconstruction kernel
coefficient of the \(Q_2\) term in reconstruction kernel
central frequency of Morlet wavelet
if set to TRUE, displays original and reconstructed signals
Returns a list containing:
reconstruction from a ridge.
<wavelets,dualwavelets> matrix.
coefficients of dual wavelets in reconstructed signal.
array containing the dual wavelets.
array containing the wavelets on sampled ridge.
wavelet transform of sol, restricted to the ridge.
wavelet transform of signal, restricted to the ridge.
number of nodes used for the reconstruction.
The computation of the kernel may be time consuming. This function avoids recomputing it if it was computed already.
See discussions in the text of ``Practical Time-Frequency Analysis''.