Welch PSD estimate (delegated to signal::pwelch / specgram)

Thin wrapper that prefers the signal package's pwelch-style routine for the production estimator; falls back to a pure-R Hamming-windowed averaged periodogram when signal is unavailable.

Usage

morie_dsp_psd_welch(x, fs = 1, nperseg = 256L, noverlap = NULL)

Arguments

x

Numeric vector.

fs

Sampling frequency (Hz). Default 1.

nperseg

Segment length. Default 256.

noverlap

Overlap in samples. Default nperseg %/% 2.

Value

List with freqs and psd.

References

Rangayyan & Krishnan (2015), Ch. 6, sec. 6.4; Welch (1967).