是否有用於在 Python 中計算長期平均頻譜 (LTAS) 的軟件包？

Question

我在 Python 3 中搜索類似於https://de.mathworks.com/matlabcentral/fileexchange/55212-long-term-average-spectrum的包。我找到了https://homepage.univie.ac.at/christian。 Herbst/python/index.html#praatUtilDemo ，但是它似乎不適用於 Python 3，我無法讓它運行。 另外，我嘗試自己實現它，但頻譜看起來很奇怪。 我的步驟是：

1. 從我的波形文件中獲取 SPL。
2. 將音頻序列分成 20 ms 間隔（非重疊）
3. 使用 numpy 計算每個時間幀的傅立葉變換timeframe = np.fft.rfft(timeframe-np.mean(timeframe)) 。
4. 使用np.abs(timeframe)處理結果的虛部。
5. 計算所有區間的平均頻譜。

Answer 1

Welch 方法可用於計算時間序列的平均頻譜。 它執行您在步驟 2-5 中打算執行的操作。

查看 SciPy 信號包和 welch 方法： https ://docs.scipy.org/doc/scipy-0.14.0/reference/generated/scipy.signal.welch.html