Python 多窗口法mtm功率谱实现

callmezoro

pmtm是matlab里面多窗口法的函数，但是python里使用就不像welch那么方便，我找到了一个包 spectrum, url: https://pyspectrum.readthedocs.io/en/latest/ 它里面可以提供pmtm的方法，但是问题是他返回的是复数Sk，权重和特征值三个变量，这三个变量前两个是二维的，特征值是一维的，无法直接进行画图处理，然而这个函数自带一个功能可以直接出图，却不给你画图的数据，但是所以我就写了个小函数，只要把pmtm返回的结果放到这个函数里就可以得到画图的数据，和自带出图功能比较是完全一致的也可证明其准确性

ps:spectrum.mtm.pmtm这个函数默认的NFFT和通常的默认方式不同，所以要想和原数据一样，最好自己输入NFFT

2020.1201更新：计算的时候写完就传上来了，结果结果不太对，检查的时候发现是没有进行fftshift，然后顺便又加上了计算频率/波数的部分，现在函数返回的结果可以直接用来画图了

xiaoyun11070

搞一个效果图出来撒

雨里1998

请问大佬们有没有python的多窗谱相干方法的代码？