Trapezoid-equivalent rectangular integration of PSD over a band.
Usage
morie_dsp_band_power(psd, freqs, f_low, f_high)
Arguments
- psd
PSD vector.
- freqs
Matching frequency vector.
- f_low
Lower edge (Hz).
- f_high
Upper edge (Hz).
Value
Band power (units of PSD x Hz).
References
Rangayyan & Krishnan (2015), Ch. 6.
Examples
freqs <- seq.int(0, 100) / 100
psd <- rep(1, 101)
morie_dsp_band_power(psd, freqs, 0.2, 0.5) # ~0.3
#> [1] 0.31
