scipy周期图和自我实现的功率谱密度之间的差异

Question

我正在尝试评估几个信号的频域。为此，我使用了this answer中给出的PSD实现。作为比较，我使用了scipy中提供的signal.periodogram函数：

from scipy.signal import tukey
import scipy as sp
f, Pxx_den = sp.signal.periodogram(a_gtrend_orig,12,window=tukey( len(a_gtrend_orig) ))

但是，当我将其绘制在自实现的PSD旁边时，它们看起来却大不相同：

由于使用了相同的窗口函数，并且周期图函数还应使用FFT，这种差异来自何处？

Answer 1

您要与之比较的示例是绘制每个频点处的振幅，即abs（fft（））

周期图产生一个功率频谱密度，这意味着它是每个频点处的振幅的平方。

标签“ windowed psd”来自早期的编辑，后来被更正。