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

本文簡要介紹 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])

參數 ：：

z： array_like

線性濾波器的零點

p： array_like

線性濾波器的極點

k： 標量

線性濾波器的增益

worN： {無，int，數組}，可選

如果沒有，則計算響應曲線感興趣部分周圍的 200 個頻率(由 pole-zero 位置確定)。如果是單個整數，則以那麽多頻率進行計算。否則，計算工作中給出的角頻率(例如，rad/s)的響應。

返回 ：：

w： ndarray

計算 h 的角頻率。

h： ndarray

頻率響應。

注意：

例子：

>>> import numpy as np
>>> 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(True)
>>> plt.show()