python - 如何在obspy中设置零相位带通滤波器的阶数

Question

我有一个关于在 obspy 中设置零相位带通滤波器的顺序的问题。

我想应用一个二阶零相位滤波器，转角频率为 0.05Hz 和 0.1Hz。

我认为输入'st.filter（'bandpass'，freqmin = 0.05，freqmax = 0.1，corners = 2，zerophase = True）'是正确的。

但是，我有点困惑，因为根据下面的链接，设置“zerophase = True”会导致过滤器顺序的两倍。

https://docs.obspy.org/_modules/obspy/signal/filter.html

所以我想问一下在我的情况下设置'corners = 2'或'corners = 1'等是否正确。

如果有人给出答案，这将对我有很大帮助。

谢谢你。

Answer 1

您应该指定corners = 1.

实现零相位滤波器的方式是使数据通过滤波器两次，一次向前，一次向后。这导致幅度响应加倍，但抵消了相位响应。因此，如果您期望阻带中的衰减为 12 dB/倍频程，就像在正常的二阶滤波器中一样，则指定corners = 1. 这将导致前向通过 6 dB/八度音程衰减，然后在反向通过时另一个 6 dB/八度音程衰减。