resample_mvfft: Resample an fft using varying numbers of sine tapers

In abarbour/psd: Adaptive, Sine-Multitaper Power Spectral Density and Cross Spectrum Estimation

Resample an fft using varying numbers of sine tapers

Description

Produce an un-normalized psd based on an fft and a vector of optimal sine tapers.

Usage

resample_mvfft(fftz, tapers, verbose = TRUE, dbl = TRUE, tapcap = 10000L)

Arguments

fftz	complex; a matrix representing the dual-length fft; see also the dbl argument
tapers	integer; a vector of tapers
verbose	logical; should messages be given?
dbl	logical; should the code assume fftz is dual-length or single-length?
tapcap	integer; the maximum number of tapers which can be applied; note that the length is automatically limited by the length of the series.

Details

To produce a psd estimate with our adaptive spectrum estimation method, we need only make one fft calculation initially and then apply the weighting factors given by parabolic_weights, which this function does.

Value

list that includes the auto and cross-spectral density, and the number of tapers

See Also

riedsid

Examples

fftz <- complex(real=1:8, imaginary = 1:8)
taps <- 1:4
try(resample_mvfft(fftz, taps))

abarbour/psd documentation built on Aug. 15, 2023, 8:56 a.m.