fastkernel {Rwave}

Kernel for Reconstruction from Wavelet Ridges

Description

Computes the cost from the sample of points on the estimated ridge and the matrix used in the reconstruction of the original signal, using simple trapezoidal rule for integrals.

Usage

fastkernel(node, phinode, nvoice, x.inc=1, x.min=node[1],
x.max=node[length(node)], w0=2 * pi, plot=FALSE)

Arguments

node

values of the variable b for the nodes of the ridge.

phinode

values of the scale variable a for the nodes of the ridge.

nvoice

number of scales within 1 octave.

x.inc

step unit for the computation of the kernel

x.min

minimal value of x for the computation of Q_2.

x.max

maximal value of x for the computation of Q_2.

w0

central frequency of the wavelet

plot

if set to TRUE, displays the modulus of the matrix of Q_2.

Details

Uses trapezoidal rule (instead of Romberg's method) to evaluate the kernel.

Value

matrix of the Q_2 kernel.

References

See discussions in the text of “Practical Time-Frequency Analysis”.

See Also

kernel, rkernel, gkernel, zerokernel.