scipy.signal.

bilinear_zpk

scipy.signal.bilinear_zpk(z, p, k, fs)

使用双线性变换从模拟类型的滤波器返回一个数字 IIR 滤波器。

使用 Tustin 方法将一组极点和零点从模拟 s 平面变换到数字 z 平面，该方法用 2*fs*(z-1) / (z+1) 代替 s，以保持频率响应的形状。

参数:

zarray_like

模拟滤波器传递函数的零点。

parray_like

模拟滤波器传递函数的极点。

kfloat

模拟滤波器传递函数的系统增益。

fsfloat

作为普通频率（如赫兹）的采样率。此函数中不进行预翘曲。

返回:

zndarray

转换后的数字滤波器传递函数的零点。

pndarray

转换后的数字滤波器传递函数的极点。

kfloat

转换后的数字滤波器的系统增益。

另请参见

lp2lp_zpk、lp2hp_zpk、lp2bp_zpk、lp2bs_zpk

双线性

备注

版本 1.1.0 中新增。

示例

在您的浏览器中试用！

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

>>> fs = 100
>>> bf = 2 * np.pi * np.array([7, 13])
>>> filts = signal.lti(*signal.butter(4, bf, btype='bandpass', analog=True,
...                                   output='zpk'))
>>> filtz = signal.lti(*signal.bilinear_zpk(filts.zeros, filts.poles,
...                                         filts.gain, fs))
>>> wz, hz = signal.freqz_zpk(filtz.zeros, filtz.poles, filtz.gain)
>>> ws, hs = signal.freqs_zpk(filts.zeros, filts.poles, filts.gain,
...                           worN=fs*wz)
>>> plt.semilogx(wz*fs/(2*np.pi), 20*np.log10(np.abs(hz).clip(1e-15)),
...              label=r'$|H_z(e^{j \omega})|$')
>>> plt.semilogx(wz*fs/(2*np.pi), 20*np.log10(np.abs(hs).clip(1e-15)),
...              label=r'$|H(j \omega)|$')
>>> plt.legend()
>>> plt.xlabel('Frequency [Hz]')
>>> plt.ylabel('Magnitude [dB]')
>>> plt.grid(True)

返回在新选项卡中打开