Python SciPy signal.lp2bs用法及代碼示例

本文簡要介紹 python 語言中 scipy.signal.lp2bs 的用法。

用法:

scipy.signal.lp2bs(b, a, wo=1.0, bw=1.0)#

將低通濾波器原型轉換為帶阻濾波器。

返回具有中心頻率的模擬 band-stop 濾波器我和帶寬體重來自具有統一截止頻率的模擬 low-pass 濾波器原型，以傳遞函數 (‘ba’) 表示。

參數 ：：

b： array_like

分子多項式係數。

a： array_like

分母多項式係數。

wo： 浮點數

所需的阻帶中心，作為角頻率(例如，rad/s)。默認為無變化。

bw： 浮點數

所需的阻帶寬度，作為角頻率(例如，rad/s)。默認為 1。

返回 ：：

b： array_like

變換後的band-stop 濾波器的分子多項式係數。

a： array_like

變換後的band-stop 濾波器的分母多項式係數。

注意：

這源自 s-plane 替換

\[s \rightarrow \frac{s \cdot \mathrm{BW}}{s^2 + {\omega_0}^2}\]

這是“wideband” 變換，產生一個關於 wo 幾何(對數頻率)對稱的阻帶。

例子：

>>> from scipy import signal
>>> import matplotlib.pyplot as plt

>>> lp = signal.lti([1.0], [1.0, 1.5])
>>> bs = signal.lti(*signal.lp2bs(lp.num, lp.den))
>>> w, mag_lp, p_lp = lp.bode()
>>> w, mag_bs, p_bs = bs.bode(w)
>>> plt.plot(w, mag_lp, label='Lowpass')
>>> plt.plot(w, mag_bs, label='Bandstop')
>>> plt.semilogx()
>>> plt.grid(True)
>>> plt.xlabel('Frequency [rad/s]')
>>> plt.ylabel('Magnitude [dB]')
>>> plt.legend()