y = fft(ephat, nfft, axis=0)