tsa
power spectral density using fftshift
currently returns two-sided according to fft frequencies, use first half
ArmaFft.spdshift