B-Spline Wavelet Transform

Computes the stationary wavelet transform of a signal using B-spline filters. This specific implementation is optimized for 1D signals like otolith radial distances.

Usage

fwaveletspl_3(ra0, sjm, det = FALSE)

Arguments

ra0

A numeric vector representing the input signal (must be a power of 2 length for some applications, though here it uses the signal length).

sjm

An integer specifying the number of decomposition scales.

det

A logical value. If FALSE (default), it returns the approximation coefficients (low-frequency). If TRUE, it returns the detail coefficients (high-frequency).

Value

A matrix of dimensions (sjm x length(ra0)) containing the coefficients for each scale.

Examples

# Decompose a signal into 5 scales
sig <- rnorm(512)
sjm <- 9 # Number of scales
wavelets <- fwaveletspl_3(sig, sjm)
# plot(wavelets[1,], type="l") # Plot first scale