k-th spectral moment — morie_dsp_spectral

m_k = sum(freqs^k * psd) * df. Used to derive mean, median, edge, and higher-order frequency descriptors.

Usage

morie_dsp_spectral_moment(psd, freqs, order = 0L)

Arguments

psd: PSD vector.
freqs: Matching frequency vector.
order: Moment order. Default 0.

Value

Scalar moment.

References

Rangayyan & Krishnan (2015), Ch. 6, sec. 6.6.

Examples

fs <- 1024; t <- seq.int(0, 1023) / fs
out <- morie_dsp_psd_periodogram(cos(2 * pi * 50 * t), fs = fs)
morie_dsp_spectral_moment(out$psd, out$freqs, order = 0L)
#> [1] 0.5