我正在尝试评估几个信号的频域。为此,我使用了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,这种差异从何而来?在
将其与之比较的示例是绘制每个频点的振幅,即abs(fft())
周期图产生了一个功率谱密度,这意味着它是每个频点处振幅的平方。在
“窗口psd”标签来自早期编辑,后来被更正。在
相关问题 更多 >
编程相关推荐