# Here are four simplified power spectra, each with 9 frequency bins:
s = list(
c(rep(0, 4), 1, rep(0, 4)), # a single peak in spectrum
c(0, 0, 1, 0, 0, .75, 0, 0, .5), # perfectly periodic, with 3 harmonics
rep(0, 9), # a silent frame
rep(1, 9) # white noise
)
# Weiner entropy is ~0 for periodic, NA for silent, 1 for white noise
sapply(s, function(x) round(getEntropy(x), 2))
# Shannon entropy is ~0 for periodic with a single harmonic, moderate for
# periodic with multiple harmonics, NA for silent, highest for white noise
sapply(s, function(x) round(getEntropy(x, type = 'shannon'), 2))
# Normalized Shannon entropy - same but forced to be 0 to 1
sapply(s, function(x) round(getEntropy(x,
type = 'shannon', normalize = TRUE), 2))
Run the code above in your browser using DataLab