Blackman-Tukey Spectrum
The Blackman-Tukey spectrum (or correlogram) is defined as follows
\[S_{BT}(\omega)=\sum_{k=-M+1}^{M-1} w_k \hat r_k e^{-i \omega k}\]where $w_k$ are window weights and $\hat r_k$ are standard sample covariance (or correlation) estimates
Implementation
The Blackman-Tukey periodogram is implemented in the demetra.data.SmoothedPeriodogram
class.
Details
The estimation of the spectrum is defined as follows:
- The series is corrected for mean effect ($y_t^c=y_t-\overline y$)
- Optionally, tapering is applied on the series
- The $M-1$ first sample correlations are computed
- The spectrum is computed, using the provided smoothing window (Tukey window by default)
- it is evaluated at frequencies defined by $2 \pi k / \left(r\left(2M-1\right)\right)$, where $r$ is the relative (in $]0,1]$) resolution defined by the user (.5 by default)
Example
SmoothedPeriodogram periodogram = SmoothedPeriodogram.builder()
.data(DoubleSequence.of(data))
.taper(new TukeyHanningTaper(.1))
.windowLength(85)
.relativeResolution(1)
.windowFunction(DiscreteWindowFunction.Tukey)
.build();