当前位置: 首页>>代码示例>>用法及示例精选>>正文


python scipy signal.freqs_zpk用法及代码示例

用法:

scipy.signal.freqs_zpk(z, p, k, worN=200)

计算模拟滤波器的频率响应。

给定滤波器的零点z,极点p和增益k,计算其频率响应:

(jw-z[0]) * (jw-z[1]) * ... * (jw-z[-1])
H(w) = k * ----------------------------------------
           (jw-p[0]) * (jw-p[1]) * ... * (jw-p[-1])

参数:

zarray_like

线性滤波器的零点

parray_like

线性滤波器的极点

kscalar

线性滤波器的增益

worN{None, int, array_like}, 可选参数

如果为“无”,则以响应曲线的有趣部分(由pole-zero位置确定)周围的200个频率进行计算。如果是单个整数,则以那么多的频率进行计算。否则,计算在worN中给出的角频率(例如rad /s)下的响应。

返回值:

wndarray

计算h的角频率。

hndarray

频率响应。

注意:

0.19.0版中的新函数。

例子:

>>> from scipy.signal import freqs_zpk, iirfilter
>>> z, p, k = iirfilter(4, [1, 10], 1, 60, analog=True, ftype='cheby1',
...                     output='zpk')
>>> w, h = freqs_zpk(z, p, k, worN=np.logspace(-1, 2, 1000))
>>> import matplotlib.pyplot as plt
>>> plt.semilogx(w, 20 * np.log10(abs(h)))
>>> plt.xlabel('Frequency')
>>> plt.ylabel('Amplitude response [dB]')
>>> plt.grid()
>>> plt.show()
../_images/scipy-signal-freqs_zpk-1.png

源码:

scipy.signal.freqs_zpk的API实现见:[源代码]

相关用法

注:本文由纯净天空筛选整理自 scipy.signal.freqs_zpk。非经特殊声明,原始代码版权归原作者所有,本译文的传播和使用请遵循“署名-相同方式共享 4.0 国际 (CC BY-SA 4.0)”协议。