Python 中的 spec.pgram 等效项

Question

我正在使用 R 编写一些遗留代码。在代码中，周期图的频率和强度被提取并用于一些下游分析：

spec.pgram(tdata,pad=0.3,tap=0.3,span=2,plot=F,detrend=TRUE,demean=TRUE)

tdata 是一维时间序列数据。

问题是我想将此部分移植到 Python，但我无法重现此命令的输出。

到目前为止，我尝试使用 scipy 实现：

scipy.signal.periodogram(tdata)

以及pyspectrum包的实现：

p = spectrum.Periodogram(tdataA1n)
freq = p.frequencies()
intensities = p.psd

不幸的是，这两种方法都产生了相似的“模式”，但值的幅度明显较低（例如，R 中的值范围为 0-300，而 R 中的值为 0-600）。我假设这是基于 R 中的一些参数，但我不知道如何将它们转换为 python。

我们将不胜感激。

Answer 1

如果有人在寻找解决方案，我找到了一个 GitHub repo，它带有一个 python 版本的 spec.pgram，这对我来说就像一个魅力。不幸的是它不包含明确的许可证，所以我猜原始实现的 GPL 许可证是适用的。