Python中的多锥函数？

Posted 2023-02-25

【中文标题】Python中的多锥函数？

【英文标题】：Multitaper function in Python?

【发布时间】：2015-10-27 09:03:54

【问题描述】：

我正在 Python 中搜索一个函数，它返回多锥功率谱密度，例如 Matlab 的 pmtm (http://fr.mathworks.com/help/signal/ref/pmtm.html)

有人知道吗？我尝试在 Python 3 上安装 mtspec 但没有成功。还有其他选择吗？

谢谢

【参考方案1】：

您可以使用来自spectrum 的pmtm，它返回多锥度谱估计

使用 pip 安装：

pip install spectrum

并像文档示例一样使用：

from spectrum import *

data = data_cosine(N=2048, A=0.1, sampling=1024, freq=200)
# If you already have the DPSS windows
[tapers, eigen] = dpss(2048, 2.5, 4)
res = pmtm(data, e=tapers, v=eigen, show=False)
# You do not need to compute the DPSS before end
res = pmtm(data, NW=2.5, show=False)
res = pmtm(data, NW=2.5, k=4, show=True)