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()