Python SciPy signal.lp2bp用法及代码示例

本文简要介绍 python 语言中 scipy.signal.lp2bp 的用法。

用法:

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

将低通滤波器原型转换为带通滤波器。

返回具有中心频率的模拟 band-pass 滤波器我和带宽体重来自具有统一截止频率的模拟 low-pass 滤波器原型，以传递函数 (‘ba’) 表示。

参数 ：：

b： array_like

分子多项式系数。

a： array_like

分母多项式系数。

wo： 浮点数

所需的通带中心，作为角频率(例如，rad/s)。默认为无变化。

bw： 浮点数

所需的通带宽度，作为角频率(例如，rad/s)。默认为 1。

返回 ：：

b： array_like

变换后的band-pass 滤波器的分子多项式系数。

a： array_like

变换后的band-pass 滤波器的分母多项式系数。

注意：

这源自 s-plane 替换

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

这是“wideband” 变换，产生一个关于 wo 几何(对数频率)对称的通带。

例子：

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

>>> lp = signal.lti([1.0], [1.0, 1.0])
>>> bp = signal.lti(*signal.lp2bp(lp.num, lp.den))
>>> w, mag_lp, p_lp = lp.bode()
>>> w, mag_bp, p_bp = bp.bode(w)

>>> plt.plot(w, mag_lp, label='Lowpass')
>>> plt.plot(w, mag_bp, label='Bandpass')
>>> plt.semilogx()
>>> plt.grid(True)
>>> plt.xlabel('Frequency [rad/s]')
>>> plt.ylabel('Magnitude [dB]')
>>> plt.legend()