scipy.signal.

freqresp#

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

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

参数:

system一个 lti 类实例或描述该系统的元组。

下面给出元组中元素的数量及其解释

1 ( lti 的实例)

2 (num, den)

3 (zeros, poles, gain)

4 (A, B, C, D)

w类数组, 可选

频率数组 (单位: 弧度/秒)。将为该数组中的每个值计算幅度和相位数据。如果未给出，将计算一组合理的值。

n整数, 可选

要计算的频率点数，如果 w 未给出。这 n 个频率在选定区间内按对数间隔，该区间包括系统零点和极点的影响。

返回:

备注

如果为 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()