Median frequency from PSD

Frequency at which the cumulative spectrum reaches half the total power. Robust to high-frequency outliers vs. the mean frequency.

Usage

morie_dsp_median_frequency(psd, freqs)

Arguments

psd

PSD vector.

freqs

Matching frequency vector.

Value

Scalar (Hz).

References

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

Examples

freqs <- seq.int(0, 100) / 100
psd <- rep(1, 101)
morie_dsp_median_frequency(psd, freqs)
#> [1] 0.5