Spectral crest for audio signals and auditory spectrograms

`crest = spectralCrest(x,f)`

`crest = spectralCrest(x,f,Name,Value)`

`[crest,spectralPeak,spectralMean] = spectralCrest(___)`

specifies options using one or more `crest`

= spectralCrest(`x`

,`f`

,`Name,Value`

)`Name,Value`

pair arguments.

`[`

returns the spectral peak and spectral mean.`crest`

,`spectralPeak`

,`spectralMean`

] = spectralCrest(___)

The spectral crest is calculated as described in [1]:

$$\text{crest}=\frac{\mathrm{max}\left({s}_{k\in [{b}_{1},{b}_{2}]}\right)}{\frac{1}{{b}_{2}-{b}_{1}}{\displaystyle \sum _{k={b}_{1}}^{{b}_{2}}{s}_{k}}}$$

where

*s*is the spectral value at bin_{k}*k*.*b*_{1}and*b*_{2}are the band edges, in bins, over which to calculate the spectral crest.

[1] Peeters, G. "A Large Set of Audio Features for Sound Description (Similarity and Classification) in the CUIDADO Project." Technical Report; IRCAM: Paris, France, 2004.