pub fn psd(data: &[f64]) -> Vec<f64>
Compute the power spectral density: |FFT(x)|^2 for each frequency bin.
|FFT(x)|^2