Spectral kurtosis from PSD

Standardised fourth central moment of frequency under the PSD treated as a probability mass.

Usage

morie_dsp_spectral_kurtosis(psd, freqs)

Arguments

psd

PSD vector.

freqs

Matching frequency vector.

Value

Scalar.

References

Rangayyan & Krishnan (2015), Ch. 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_kurtosis(out$psd, out$freqs)
#> [1] 2.184481e+28