scipy.signal.

freqresp#

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

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

参数:
systemlti 类的实例或描述系统的元组。

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

  • 1 ( lti 的实例)

  • 2 (num, den)

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

  • 4 (A, B, C, D)

warray_like, 可选

频率数组(单位为 rad/s)。计算此数组中每个值的大小和相位数据。如果未给出,将计算一组合理的值。

nint, 可选

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

返回:
w一维 ndarray

频率数组 [rad/s]

H一维 ndarray

复数幅度值的数组

注释

如果将 (num, den) 传入 system,则分子和分母的系数都应按降幂顺序指定(例如,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