scipy.signal.

freqresp#

scipy.signal.freqresp(system, w=None, n=10000)[源代码][源代码]#

计算连续时间系统的频率响应。

参数:
系统 : lti 类的一个实例或描述系统的元组。一个实例

以下给出了元组中元素的数量及其解释:

  • 1 (lti 的实例)

  • 2 (分子, 分母)

  • 3 (零点, 极点, 增益)

  • 4 (A, B, C, D)

w类似数组, 可选

频率数组(以 rad/s 为单位)。幅度和相位数据是为该数组中的每个值计算的。如果未提供,将计算一个合理的集合。

nint, 可选

如果未给出 w,则计算的频率点数。n 个频率在包含系统极点和零点影响的区间内对数间隔分布。

返回:
w1D ndarray

频率数组 [rad/s]

H1D ndarray

复数幅值数组

注释

如果为 system 传递了 (num, den),则分子和分母的系数应以降幂顺序指定(例如,s^2 + 3s + 5 应表示为 [1, 3, 5])。

示例

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

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

构建传递函数 \(H(s) = \frac{5}{(s-1)^3}\)

>>> s1 = signal.ZerosPolesGain([], [1, 1, 1], [5])
>>> w, H = signal.freqresp(s1)
>>> plt.figure()
>>> plt.plot(H.real, H.imag, "b")
>>> plt.plot(H.real, -H.imag, "r")
>>> plt.show()
../../_images/scipy-signal-freqresp-1.png