We first generated simulated time-series data of
length from an auto-regressive process
of order with a peaked spectrum.
We then computed the DFT, calculated the magnitude-squared
of the bins, keeping the unique bins as our
"raw data" .
We then created the
matrix that computed the
first three auto-correlation lags (0 through 2),
then computed the feature (5.1).