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

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

用法:

scipy.signal.dfreqresp(system, w=None, n=10000, whole=False)#

计算discrete-time 系统的频率响应。

参数 ：：

system： dlti 类的实例或说明系统的元组。

下面给出元组中元素的数量和解释：

• 1 (instance of dlti)

• 2 (numerator, denominator, dt)

• 3 (zeros, poles, gain, dt)

• 4 (A, B, C, D, dt)

w： 数组，可选

频率数组(以弧度/样本为单位)。为该数组中的每个值计算幅度和相位数据。如果没有给出合理的设置，将被计算。

n： 整数，可选

如果未给出 w，则要计算的频点数。 n 个频率以对数方式隔开，间隔选择包括系统极点和零点的影响。

whole： 布尔型，可选

通常，如果没有给出‘w’，则从 0 到奈奎斯特频率计算频率，pi 弧度/样本(upper-half of unit-circle)。如果整体为真，则计算从 0 到 2*pi 弧度/样本的频率。

返回 ：：

w： 一维ndarray

频率阵列 [弧度/样本]

H： 一维ndarray

复数幅度值数组

注意：

如果 (num, den) 为 system 传入，则分子和分母的系数应按 index 降序指定(例如 z^2 + 3z + 5 将表示为 [1, 3, 5] )。

例子：

生成传递函数的奈奎斯特图

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

构造采样时间为 0.05 秒的传递函数\(H(z) = \frac{1}{z^2 + 2z + 3}\)：

>>> sys = signal.TransferFunction([1], [1, 2, 3], dt=0.05)

>>> w, H = signal.dfreqresp(sys)

>>> plt.figure()
>>> plt.plot(H.real, H.imag, "b")
>>> plt.plot(H.real, -H.imag, "r")
>>> plt.show()